Trophic interactions with heterotrophic bacteria limit the range of Prochlorococcus.

Prochlorococcus is both the smallest and numerically most abundant photosynthesizing organism on the planet. While thriving in the warm oligotrophic gyres, Prochlorococcus concentrations drop rapidly in higher-latitude regions. Transect data from the North Pacific show the collapse occurring at a wide range of temperatures and latitudes (temperature is often hypothesized to cause this shift), suggesting an ecological mechanism may be at play. An often used size-based theory of phytoplankton community structure that has been incorporated into computational models correctly predicts the dominance of Prochlorococcus in the gyres, and the relative dominance of larger cells at high latitudes. However, both theory and computational models fail to explain the poleward collapse. When heterotrophic bacteria and predators that prey nonspecifically on both Prochlorococcus and bacteria are included in the theoretical framework, the collapse of Prochlorococcus occurs with increasing nutrient supplies. The poleward collapse of Prochlorococcus populations then naturally emerges when this mechanism of "shared predation" is implemented in a complex global ecosystem model. Additionally, the theory correctly predicts trends in both the abundance and mean size of the heterotrophic bacteria. These results suggest that ecological controls need to be considered to understand the biogeography of Prochlorococcus and predict its changes under future ocean conditions. Indirect interactions within a microbial network can be essential in setting community structure.

The dynamic trophic architecture of open-ocean protist communities revealed through machine-guided metatranscriptomics.

Intricate networks of single-celled eukaryotes (protists) dominate carbon flow in the ocean. Their growth, demise, and interactions with other microorganisms drive the fluxes of biogeochemical elements through marine ecosystems. Mixotrophic protists are capable of both photosynthesis and ingestion of prey and are dominant components of open-ocean planktonic communities. Yet the role of mixotrophs in elemental cycling is obscured by their capacity to act as primary producers or heterotrophic consumers depending on factors that remain largely uncharacterized. Here, we develop and apply a machine learning model that predicts the in situ trophic mode of aquatic protists based on their patterns of gene expression. This approach leverages a public collection of protist transcriptomes as a training set to identify a subset of gene families whose transcriptional profiles predict trophic mode. We applied our model to nearly 100 metatranscriptomes obtained during two oceanographic cruises in the North Pacific and found community-level and population-specific evidence that abundant open-ocean mixotrophic populations shift their predominant mode of nutrient and carbon acquisition in response to natural gradients in nutrient supply and sea surface temperature. Metatranscriptomic data from ship-board incubation experiments revealed that abundant mixotrophic prymnesiophytes from the oligotrophic North Pacific subtropical gyre rapidly remodeled their transcriptome to enhance photosynthesis when supplied with limiting nutrients. Coupling this approach with experiments designed to reveal the mechanisms driving mixotroph physiology provides an avenue toward understanding the ecology of mixotrophy in the natural environment.

Diel transcriptional oscillations of light-sensitive regulatory elements in open-ocean eukaryotic plankton communities.

The 24-h cycle of light and darkness governs daily rhythms of complex behaviors across all domains of life. Intracellular photoreceptors sense specific wavelengths of light that can reset the internal circadian clock and/or elicit distinct phenotypic responses. In the surface ocean, microbial communities additionally modulate nonrhythmic changes in light quality and quantity as they are mixed to different depths. Here, we show that eukaryotic plankton in the North Pacific Subtropical Gyre transcribe genes encoding light-sensitive proteins that may serve as light-activated transcription factors, elicit light-driven electrical/chemical cascades, or initiate secondary messenger-signaling cascades. Overall, the protistan community relies on blue light-sensitive photoreceptors of the cryptochrome/photolyase family, and proteins containing the Light-Oxygen-Voltage (LOV) domain. The greatest diversification occurred within Haptophyta and photosynthetic stramenopiles where the LOV domain was combined with different DNA-binding domains and secondary signal-transduction motifs. Flagellated protists utilize green-light sensory rhodopsins and blue-light helmchromes, potentially underlying phototactic/photophobic and other behaviors toward specific wavelengths of light. Photoreceptors such as phytochromes appear to play minor roles in the North Pacific Subtropical Gyre. Transcript abundance of environmental light-sensitive protein-encoding genes that display diel patterns are found to primarily peak at dawn. The exceptions are the LOV-domain transcription factors with peaks in transcript abundances at different times and putative phototaxis photoreceptors transcribed throughout the day. Together, these data illustrate the diversity of light-sensitive proteins that may allow disparate groups of protists to respond to light and potentially synchronize patterns of growth, division, and mortality within the dynamic ocean environment.

A Bayesian approach to modeling phytoplankton population dynamics from size distribution time series.

The rates of cell growth, division, and carbon loss of microbial populations are key parameters for understanding how organisms interact with their environment and how they contribute to the carbon cycle. However, the invasive nature of current analytical methods has hindered efforts to reliably quantify these parameters. In recent years, size-structured matrix population models (MPMs) have gained popularity for estimating division rates of microbial populations by mechanistically describing changes in microbial cell size distributions over time. Motivated by the mechanistic structure of these models, we employ a Bayesian approach to extend size-structured MPMs to capture additional biological processes describing the dynamics of a marine phytoplankton population over the day-night cycle. Our Bayesian framework is able to take prior scientific knowledge into account and generate biologically interpretable results. Using data from an exponentially growing laboratory culture of the cyanobacterium Prochlorococcus, we isolate respiratory and exudative carbon losses as critical parameters for the modeling of their population dynamics. The results suggest that this modeling framework can provide deeper insights into microbial population dynamics provided by size distribution time-series data.

Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus.

The Southern Ocean houses a diverse and productive community of organisms. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean.

Anthropogenic Asian aerosols provide Fe to the North Pacific Ocean.

Fossil-fuel emissions may impact phytoplankton primary productivity and carbon cycling by supplying bioavailable Fe to remote areas of the ocean via atmospheric aerosols. However, this pathway has not been confirmed by field observations of anthropogenic Fe in seawater. Here we present high-resolution trace-metal concentrations across the North Pacific Ocean (158°W from 25°to 42°N). A dissolved Fe maximum was observed around 35°N, coincident with high dissolved Pb and Pb isotope ratios matching Asian industrial sources and confirming recent aerosol deposition. Iron-stable isotopes reveal in situ evidence of anthropogenic Fe in seawater, with low δ56Fe (-0.23‰ > δ56Fe > -0.65‰) observed in the region that is most influenced by aerosol deposition. An isotope mass balance suggests that anthropogenic Fe contributes 21-59% of dissolved Fe measured between 35° and 40°N. Thus, anthropogenic aerosol Fe is likely to be an important Fe source to the North Pacific Ocean.

Glycine betaine uptake and metabolism in marine microbial communities.

Glycine betaine (GBT) is a compatible solute in high concentrations in marine microorganisms. As a component of labile organic matter, GBT has complex biochemical potential as a substrate for microbial use that is unconstrained in the environment. Here we determine the uptake kinetics and metabolic fate of GBT in two natural microbial communities in the North Pacific characterized by different nitrate concentrations. Dissolved GBT had maximum uptake rates of 0.36 and 0.56 nM h-1 with half-saturation constants of 79 and 11 nM in the high nitrate and low nitrate stations respectively. During multiday incubations, most GBT taken into cells was retained as a compatible solute. Stable isotopes derived from the added GBT were also observed in other metabolites, including choline, carnitine and sarcosine, suggesting that GBT was used for biosynthesis and for catabolism to pyruvate and ammonium. Where nitrate was scarce, GBT was primarily metabolized via demethylation to glycine. Gene transcript data were consistent with SAR11 using GBT as a source of methyl groups to fuel the methionine cycle. Where nitrate concentrations were higher, more GBT was partitioned for lipid biosynthesis by both bacteria and eukaryotic phytoplankton. Our data highlight unexpected metabolic pathways and potential routes of microbial metabolite exchange.

Two distinct pools of B12 analogs reveal community interdependencies in the ocean.

Organisms within all domains of life require the cofactor cobalamin (vitamin B12), which is produced only by a subset of bacteria and archaea. On the basis of genomic analyses, cobalamin biosynthesis in marine systems has been inferred in three main groups: select heterotrophic Proteobacteria, chemoautotrophic Thaumarchaeota, and photoautotrophic Cyanobacteria. Culture work demonstrates that many Cyanobacteria do not synthesize cobalamin but rather produce pseudocobalamin, challenging the connection between the occurrence of cobalamin biosynthesis genes and production of the compound in marine ecosystems. Here we show that cobalamin and pseudocobalamin coexist in the surface ocean, have distinct microbial sources, and support different enzymatic demands. Even in the presence of cobalamin, Cyanobacteria synthesize pseudocobalamin-likely reflecting their retention of an oxygen-independent pathway to produce pseudocobalamin, which is used as a cofactor in their specialized methionine synthase (MetH). This contrasts a model diatom, Thalassiosira pseudonana, which transported pseudocobalamin into the cell but was unable to use pseudocobalamin in its homolog of MetH. Our genomic and culture analyses showed that marine Thaumarchaeota and select heterotrophic bacteria produce cobalamin. This indicates that cobalamin in the surface ocean is a result of de novo synthesis by heterotrophic bacteria or via modification of closely related compounds like cyanobacterially produced pseudocobalamin. Deeper in the water column, our study implicates Thaumarchaeota as major producers of cobalamin based on genomic potential, cobalamin cell quotas, and abundance. Together, these findings establish the distinctive roles played by abundant prokaryotes in cobalamin-based microbial interdependencies that sustain community structure and function in the ocean.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

Confounding effects of oxygen and temperature on the TEX86 signature of marine Thaumarchaeota.

Marine ammonia-oxidizing archaea (AOA) are among the most abundant of marine microorganisms, spanning nearly the entire water column of diverse oceanic provinces. Historical patterns of abundance are preserved in sediments in the form of their distinctive glycerol dibiphytanyl glycerol tetraether (GDGT) membrane lipids. The correlation between the composition of GDGTs in surface sediment and the overlying annual average sea surface temperature forms the basis for a paleotemperature proxy (TEX86) that is used to reconstruct surface ocean temperature as far back as the Middle Jurassic. However, mounting evidence suggests that factors other than temperature could also play an important role in determining GDGT distributions. We here use a study set of four marine AOA isolates to demonstrate that these closely related strains generate different TEX86-temperature relationships and that oxygen (O2) concentration is at least as important as temperature in controlling TEX86 values in culture. All of the four strains characterized showed a unique membrane compositional response to temperature, with TEX86-inferred temperatures varying as much as 12 °C from the incubation temperatures. In addition, both linear and nonlinear TEX86-temperature relationships were characteristic of individual strains. Increasing relative abundance of GDGT-2 and GDGT-3 with increasing O2 limitation, at the expense of GDGT-1, led to significant elevations in TEX86-derived temperature. Although the adaptive significance of GDGT compositional changes in response to both temperature and O2 is unclear, this observation necessitates a reassessment of archaeal lipid-based paleotemperature proxies, particularly in records that span low-oxygen events or underlie oxygen minimum zones.

Light-driven synchrony of Prochlorococcus growth and mortality in the subtropical Pacific gyre.

Theoretical studies predict that competition for limited resources reduces biodiversity to the point of ecological instability, whereas strong predator/prey interactions enhance the number of coexisting species and limit fluctuations in abundances. In open ocean ecosystems, competition for low availability of essential nutrients results in relatively few abundant microbial species. The remarkable stability in overall cell abundance of the dominant photosynthetic cyanobacterium Prochlorococcus is assumed to reflect a simple food web structure strongly controlled by grazers and/or viruses. This hypothesized link between stability and ecological interactions, however, has been difficult to test with open ocean microbes because sampling methods commonly have poor temporal and spatial resolution. Here we use continuous techniques on two different winter-time cruises to show that Prochlorococcus cell production and mortality rates are tightly synchronized to the day/night cycle across the subtropical Pacific Ocean. In warmer waters, we observed harmonic oscillations in cell production and mortality rates, with a peak in mortality rate consistently occurring ∼6 h after the peak in cell production. Essentially no cell mortality was observed during daylight. Our results are best explained as a synchronized two-component trophic interaction with the per-capita rates of Prochlorococcus consumption driven either directly by the day/night cycle or indirectly by Prochlorococcus cell production. Light-driven synchrony of food web dynamics in which most of the newly produced Prochlorococcus cells are consumed each night likely enforces ecosystem stability across vast expanses of the open ocean.

Cryptic carbon and sulfur cycling between surface ocean plankton.

About half the carbon fixed by phytoplankton in the ocean is taken up and metabolized by marine bacteria, a transfer that is mediated through the seawater dissolved organic carbon (DOC) pool. The chemical complexity of marine DOC, along with a poor understanding of which compounds form the basis of trophic interactions between bacteria and phytoplankton, have impeded efforts to identify key currencies of this carbon cycle link. Here, we used transcriptional patterns in a bacterial-diatom model system based on vitamin B12 auxotrophy as a sensitive assay for metabolite exchange between marine plankton. The most highly up-regulated genes (up to 374-fold) by a marine Roseobacter clade bacterium when cocultured with the diatom Thalassiosira pseudonana were those encoding the transport and catabolism of 2,3-dihydroxypropane-1-sulfonate (DHPS). This compound has no currently recognized role in the marine microbial food web. As the genes for DHPS catabolism have limited distribution among bacterial taxa, T. pseudonana may use this sulfonate for targeted feeding of beneficial associates. Indeed, DHPS was both a major component of the T. pseudonana cytosol and an abundant microbial metabolite in a diatom bloom in the eastern North Pacific Ocean. Moreover, transcript analysis of the North Pacific samples provided evidence of DHPS catabolism by Roseobacter populations. Other such biogeochemically important metabolites may be common in the ocean but difficult to discriminate against the complex chemical background of seawater. Bacterial transformation of this diatom-derived sulfonate represents a previously unidentified and likely sizeable link in both the marine carbon and sulfur cycles.

Recognition cascade and metabolite transfer in a marine bacteria-phytoplankton model system.

The trophic linkage between marine bacteria and phytoplankton in the surface ocean is a key step in the global carbon cycle, with almost half of marine primary production transformed by heterotrophic bacterioplankton within hours to weeks of fixation. Early studies conceptualized this link as the passive addition and removal of organic compounds from a shared seawater reservoir. Here, we analysed transcript and intracellular metabolite patterns in a two-member model system and found that the presence of a heterotrophic bacterium induced a potential recognition cascade in a marine phytoplankton species that parallels better-understood vascular plant response systems. Bacterium Ruegeria pomeroyi DSS-3 triggered differential expression of >80 genes in diatom Thalassiosira pseudonana CCMP1335 that are homologs to those used by plants to recognize external stimuli, including proteins putatively involved in leucine-rich repeat recognition activity, second messenger production and protein kinase cascades. Co-cultured diatoms also downregulated lipid biosynthesis genes and upregulated chitin metabolism genes. From differential expression of bacterial transporter systems, we hypothesize that nine diatom metabolites supported the majority of bacterial growth, among them sulfonates, sugar derivatives and organic nitrogen compounds. Similar recognition responses and metabolic linkages as observed in this model system may influence carbon transformations by ocean plankton.

Scalable clustering algorithms for continuous environmental flow cytometry.

MOTIVATION: Recent technological innovations in flow cytometry now allow oceanographers to collect high-frequency flow cytometry data from particles in aquatic environments on a scale far surpassing conventional flow cytometers. The SeaFlow cytometer continuously profiles microbial phytoplankton populations across thousands of kilometers of the surface ocean. The data streams produced by instruments such as SeaFlow challenge the traditional sample-by-sample approach in cytometric analysis and highlight the need for scalable clustering algorithms to extract population information from these large-scale, high-frequency flow cytometers. RESULTS: We explore how available algorithms commonly used for medical applications perform at classification of such a large-scale, environmental flow cytometry data. We apply large-scale Gaussian mixture models to massive datasets using Hadoop. This approach outperforms current state-of-the-art cytometry classification algorithms in accuracy and can be coupled with manual or automatic partitioning of data into homogeneous sections for further classification gains. We propose the Gaussian mixture model with partitioning approach for classification of large-scale, high-frequency flow cytometry data. AVAILABILITY AND IMPLEMENTATION: Source code available for download at https://github.com/jhyrkas/seaflow_cluster, implemented in Java for use with Hadoop. CONTACT: hyrkas@cs.washington.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nitrosopumilus maritimus gen. nov., sp. nov., Nitrosopumilus cobalaminigenes sp. nov., Nitrosopumilus oxyclinae sp. nov., and Nitrosopumilus ureiphilus sp. nov., four marine ammonia-oxidizing archaea of the phylum Thaumarchaeota.

Four mesophilic, neutrophilic, and aerobic marine ammonia-oxidizing archaea, designated strains SCM1T, HCA1T, HCE1T and PS0T, were isolated from a tropical marine fish tank, dimly lit deep coastal waters, the lower euphotic zone of coastal waters, and near-surface sediment in the Puget Sound estuary, respectively. Cells are straight or slightly curved small rods, 0.15-0.26 µm in diameter and 0.50-1.59 µm in length. Motility was not observed, although strain PS0T possesses genes associated with archaeal flagella and chemotaxis, suggesting it may be motile under some conditions. Cell membranes consist of glycerol dibiphytanyl glycerol tetraether (GDGT) lipids, with crenarchaeol as the major component. Strain SCM1T displays a single surface layer (S-layer) with p6 symmetry, distinct from the p3-S-layer reported for the soil ammonia-oxidizing archaeon Nitrososphaera viennensis EN76T. Respiratory quinones consist of fully saturated and monounsaturated menaquinones with 6 isoprenoid units in the side chain. Cells obtain energy from ammonia oxidation and use carbon dioxide as carbon source; addition of an α-keto acid (α-ketoglutaric acid) was necessary to sustain growth of strains HCA1T, HCE1T, and PS0T. Strain PS0T uses urea as a source of ammonia for energy production and growth. All strains synthesize vitamin B1 (thiamine), B2 (riboflavin), B6 (pyridoxine), and B12 (cobalamin). Optimal growth occurs between 25 and 32 °C, between pH 6.8 and 7.3, and between 25 and 37 ‰ salinity. All strains have a low mol% G+C content of 33.0-34.2. Strains are related by 98 % or greater 16S rRNA gene sequence identity, sharing ~85 % 16S rRNA gene sequence identity with Nitrososphaera viennensis EN76T. All four isolates are well separated by phenotypic and genotypic characteristics and are here assigned to distinct species within the genus Nitrosopumilus gen. nov. Isolates SCM1T (=ATCC TSD-97T =NCIMB 15022T), HCA1T (=ATCC TSD-96T), HCE1T (=ATCC TSD-98T), and PS0T (=ATCC TSD-99T) are type strains of the species Nitrosopumilusmaritimus sp. nov., Nitrosopumilus cobalaminigenes sp. nov., Nitrosopumilus oxyclinae sp. nov., and Nitrosopumilus ureiphilus sp. nov., respectively. In addition, we propose the family Nitrosopumilaceae fam. nov. and the order Nitrosopumilales ord. nov. within the class Nitrososphaeria.

Development of a molecular-based index for assessing iron status in bloom-forming pennate diatoms.

Iron availability limits primary productivity in large areas of the world's oceans. Ascertaining the iron status of phytoplankton is essential for understanding the factors regulating their growth and ecology. We developed an incubation-independent, molecular-based approach to assess the iron nutritional status of specific members of the diatom community, initially focusing on the ecologically important pennate diatom Pseudo-nitzschia. Through a comparative transcriptomic approach, we identified two genes that track the iron status of Pseudo-nitzschia with high fidelity. The first gene, ferritin (FTN), encodes for the highly specialized iron storage protein induced under iron-replete conditions. The second gene, ISIP2a, encodes an iron-concentrating protein induced under iron-limiting conditions. In the oceanic diatom Pseudo-nitzschia granii (Hasle) Hasle, transcript abundance of these genes directly relates to changes in iron availability, with increased FTN transcript abundance under iron-replete conditions and increased ISIP2a transcript abundance under iron-limiting conditions. The resulting ISIP2a:FTN transcript ratio reflects the iron status of cells, where a high ratio indicates iron limitation. Field samples collected from iron grow-out microcosm experiments conducted in low iron waters of the Gulf of Alaska and variable iron waters in the California upwelling zone verify the validity of our proposed Pseudo-nitzschia Iron Limitation Index, which can be used to ascertain in situ iron status and further developed for other ecologically important diatoms.

Marine ammonia-oxidizing archaeal isolates display obligate mixotrophy and wide ecotypic variation.

Ammonia-oxidizing archaea (AOA) are now implicated in exerting significant control over the form and availability of reactive nitrogen species in marine environments. Detailed studies of specific metabolic traits and physicochemical factors controlling their activities and distribution have not been well constrained in part due to the scarcity of isolated AOA strains. Here, we report the isolation of two new coastal marine AOA, strains PS0 and HCA1. Comparison of the new strains to Nitrosopumilus maritimus strain SCM1, the only marine AOA in pure culture thus far, demonstrated distinct adaptations to pH, salinity, organic carbon, temperature, and light. Strain PS0 sustained nearly 80% of ammonia oxidation activity at a pH as low as 5.9, indicating that coastal strains may be less sensitive to the ongoing reduction in ocean pH. Notably, the two novel isolates are obligate mixotrophs that rely on uptake and assimilation of organic carbon compounds, suggesting a direct coupling between chemolithotrophy and organic matter assimilation in marine food webs. All three isolates showed only minor photoinhibition at 15 µE ⋅ m(-2) ⋅ s(-1) and rapid recovery of ammonia oxidation in the dark, consistent with an AOA contribution to the primary nitrite maximum and the plausibility of a diurnal cycle of archaeal ammonia oxidation activity in the euphotic zone. Together, these findings highlight an unexpected adaptive capacity within closely related marine group I Archaea and provide new understanding of the physiological basis of the remarkable ecological success reflected by their generally high abundance in marine environments.

Plastid proteome prediction for diatoms and other algae with secondary plastids of the red lineage.

The plastids of ecologically and economically important algae from phyla such as stramenopiles, dinoflagellates and cryptophytes were acquired via a secondary endosymbiosis and are surrounded by three or four membranes. Nuclear-encoded plastid-localized proteins contain N-terminal bipartite targeting peptides with the conserved amino acid sequence motif 'ASAFAP'. Here we identify the plastid proteomes of two diatoms, Thalassiosira pseudonana and Phaeodactylum tricornutum, using a customized prediction tool (ASAFind) that identifies nuclear-encoded plastid proteins in algae with secondary plastids of the red lineage based on the output of SignalP and the identification of conserved 'ASAFAP' motifs and transit peptides. We tested ASAFind against a large reference dataset of diatom proteins with experimentally confirmed subcellular localization and found that the tool accurately identified plastid-localized proteins with both high sensitivity and high specificity. To identify nucleus-encoded plastid proteins of T. pseudonana and P. tricornutum we generated optimized sets of gene models for both whole genomes, to increase the percentage of full-length proteins compared with previous assembly model sets. ASAFind applied to these optimized sets revealed that about 8% of the proteins encoded in their nuclear genomes were predicted to be plastid localized and therefore represent the putative plastid proteomes of these algae.

The evolution of silicon transporters in diatoms.

Diatoms are highly productive single-celled algae that form an intricately patterned silica cell wall after every cell division. They take up and utilize silicic acid from seawater via silicon transporter (SIT) proteins. This study examined the evolution of the SIT gene family to identify potential genetic adaptations that enable diatoms to thrive in the modern ocean. By searching for sequence homologs in available databases, the diversity of organisms found to encode SITs increased substantially and included all major diatom lineages and other algal protists. A bacterial-encoded gene with homology to SIT sequences was also identified, suggesting that a lateral gene transfer event occurred between bacterial and protist lineages. In diatoms, the SIT genes diverged and diversified to produce five distinct clades. The most basal SIT clades were widely distributed across diatom lineages, while the more derived clades were lineage-specific, which together produced a distinct repertoire of SIT types among major diatom lineages. Differences in the predicted protein functional domains encoded among SIT clades suggest that the divergence of clades resulted in functional diversification among SITs. Both laboratory cultures and natural communities changed transcription of each SIT clade in response to experimental or environmental growth conditions, with distinct transcriptional patterns observed among clades. Together, these data suggest that the diversification of SITs within diatoms led to specialized adaptations among diatoms lineages, and perhaps their dominant ability to take up silicic acid from seawater in diverse environmental conditions.