Subcellular proteomics for determining iron-limited remodeling of plastids in the model diatom Thalassiosira pseudonana (Bacillariophyta).

Diatoms are important primary producers in the world's oceans, yet their growth is constrained in large regions by low bioavailable iron (Fe). Low-Fe stress-induced limitation of primary production is due to requirements for Fe in components of essential metabolic pathways including photosynthesis and other chloroplast plastid functions. Studies have shown that under low-Fe stress, diatoms alter plastid-specific processes, including components of electron transport. These physiological changes suggest changes of protein content and in protein abundances within the diatom plastid. While in silico predictions provide putative information on plastid-localized proteins, knowledge of diatom plastid proteins remains limited in comparison to well-studied model photosynthetic organisms. To address this, we employed shotgun proteomics to investigate the proteome of subcellular plastid-enriched fractions from Thalassiosira pseudonana to gain a better understanding of how the plastid proteome is remodeled in response to Fe limitation. Using mass spectrometry-based peptide identification and quantification, we analyzed T. pseudonana grown under Fe-replete and -limiting conditions. Through these analyses, we inferred the relative quantities of each protein, revealing that Fe limitation regulates major metabolic pathways in the plastid, including the Calvin cycle. Additionally, we observed changes in the expression of light-harvesting proteins. In silico localization predictions of proteins identified in this plastid-enriched proteome allowed for an in-depth comparison of theoretical versus observed plastid-localization, providing evidence for the potential of additional protein import pathways into the diatom plastid.

Oligotrophic waters of the Northwest Atlantic support taxonomically diverse diatom communities that are distinct from coastal waters.

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems.

Emerging harmful algal blooms caused by distinct seasonal assemblages of a toxic diatom.

Diatoms in the Pseudo-nitzschia genus produce the neurotoxin domoic acid. Domoic acid bioaccumulates in shellfish, causing illness in humans and marine animals upon ingestion. In 2017, high domoic acid levels in shellfish meat closed shellfish harvest in Narragansett Bay, Rhode Island for the first and only time in history, although abundant Pseudo-nitzschia have been observed for over 60 years. To investigate whether an environmental factor altered endemic Pseudo-nitzschia physiology or new domoic acid-producing strain(s) were introduced to Narragansett Bay, we conducted weekly sampling from 2017 to 2019 and compared closure samples. Plankton-associated domoic acid was quantified by LC-MS/MS and Pseudo-nitzschia spp. were identified using a taxonomically improved high-throughput rDNA sequencing approach. Comparison with environmental data revealed a detailed understanding of domoic acid dynamics and seasonal multi-species assemblages. Plankton-associated domoic acid was low throughout 2017-2019, but recurred in fall and early summer maxima. Fall domoic acid maxima contained known toxic species as well as a novel Pseudo-nitzschia genotype. Summer domoic acid maxima included fewer species but also known toxin producers. Most 2017 closure samples contained the particularly concerning toxic species, P. australis, which also appeared infrequently during 2017-2019. Recurring Pseudo-nitzschia assemblages were driven by seasonal temperature changes, and plankton-associated domoic acid correlated with low dissolved inorganic nitrogen. Thus, the Narragansett Bay closures were likely caused by both resident assemblages that become toxic depending on nutrient status as well as the episodic introductions of toxic species from oceanographic and climatic shifts.

Metatranscriptome analyses indicate resource partitioning between diatoms in the field.

Diverse communities of marine phytoplankton carry out half of global primary production. The vast diversity of the phytoplankton has long perplexed ecologists because these organisms coexist in an isotropic environment while competing for the same basic resources (e.g., inorganic nutrients). Differential niche partitioning of resources is one hypothesis to explain this "paradox of the plankton," but it is difficult to quantify and track variation in phytoplankton metabolism in situ. Here, we use quantitative metatranscriptome analyses to examine pathways of nitrogen (N) and phosphorus (P) metabolism in diatoms that cooccur regularly in an estuary on the east coast of the United States (Narragansett Bay). Expression of known N and P metabolic pathways varied between diatoms, indicating apparent differences in resource utilization capacity that may prevent direct competition. Nutrient amendment incubations skewed N/P ratios, elucidating nutrient-responsive patterns of expression and facilitating a quantitative comparison between diatoms. The resource-responsive (RR) gene sets deviated in composition from the metabolic profile of the organism, being enriched in genes associated with N and P metabolism. Expression of the RR gene set varied over time and differed significantly between diatoms, resulting in opposite transcriptional responses to the same environment. Apparent differences in metabolic capacity and the expression of that capacity in the environment suggest that diatom-specific resource partitioning was occurring in Narragansett Bay. This high-resolution approach highlights the molecular underpinnings of diatom resource utilization and how cooccurring diatoms adjust their cellular physiology to partition their niche space.

The Phaeodactylum genome reveals the evolutionary history of diatom genomes.

Diatoms are photosynthetic secondary endosymbionts found throughout marine and freshwater environments, and are believed to be responsible for around one-fifth of the primary productivity on Earth. The genome sequence of the marine centric diatom Thalassiosira pseudonana was recently reported, revealing a wealth of information about diatom biology. Here we report the complete genome sequence of the pennate diatom Phaeodactylum tricornutum and compare it with that of T. pseudonana to clarify evolutionary origins, functional significance and ubiquity of these features throughout diatoms. In spite of the fact that the pennate and centric lineages have only been diverging for 90 million years, their genome structures are dramatically different and a substantial fraction of genes ( approximately 40%) are not shared by these representatives of the two lineages. Analysis of molecular divergence compared with yeasts and metazoans reveals rapid rates of gene diversification in diatoms. Contributing factors include selective gene family expansions, differential losses and gains of genes and introns, and differential mobilization of transposable elements. Most significantly, we document the presence of hundreds of genes from bacteria. More than 300 of these gene transfers are found in both diatoms, attesting to their ancient origins, and many are likely to provide novel possibilities for metabolite management and for perception of environmental signals. These findings go a long way towards explaining the incredible diversity and success of the diatoms in contemporary oceans.

Temporal Dynamics of Cyanobacterial Bloom Community Composition and Toxin Production from Urban Lakes.

With a long evolutionary history and a need to adapt to a changing environment, cyanobacteria in freshwater systems use specialized metabolites for communication, defense, and physiological processes. However, the role that these metabolites play in differentiating species, maintaining microbial communities, and generating niche persistence and expansion is poorly understood. Furthermore, many cyanobacterial specialized metabolites and toxins present significant human health concerns due to their liver toxicity and their potential impact to drinking water. Gaps in knowledge exist with respect to changes in species diversity and toxin production during a cyanobacterial bloom (cyanoHAB) event; addressing these gaps will improve understanding of impacts to public and ecological health. In the current project, we utilized a multiomics strategy (DNA metabarcoding and metabolomics) to determine the cyanobacterial community composition, toxin profile, and the specialized metabolite pool at three freshwater lakes in Providence, RI during summer-fall cyanoHABs. Species diversity decreased at all study sites over the course of the bloom event, and toxin production reached a maximum at the midpoint of the event. Additionally, LC-MS/MS-based molecular networking identified new toxin congeners. This work provokes intriguing questions with respect to the use of allelopathy by organisms in these systems and the presence of emerging toxic compounds that can impact public health.

Simulated upwelling and marine heatwave events promote similar growth rates but differential domoic acid toxicity in Pseudo-nitzschia australis.

Along the west coast of the United States, highly toxic Pseudo-nitzschia blooms have been associated with two contrasting regional phenomena: seasonal upwelling and marine heatwaves. While upwelling delivers cool water rich in pCO2 and an abundance of macronutrients to the upper water column, marine heatwaves instead lead to warmer surface waters, low pCO2, and reduced nutrient availability. Understanding Pseudo-nitzschia dynamics under these two conditions is important for bloom forecasting and coastal management, yet the mechanisms driving toxic bloom formation during contrasting upwelling vs. heatwave conditions remain poorly understood. To gain a better understanding of what drives Pseudo-nitzschia australis growth and toxicity during these events, multiple-driver scenario or 'cluster' experiments were conducted using temperature, pCO2, and nutrient levels reflecting conditions during upwelling (13 °C, 900 ppm pCO2, replete nutrients) and two intensities of marine heatwaves (19 °C or 20.5 °C, 250 ppm pCO2, reduced macronutrients). While P. australis grew equally well under both heatwave and upwelling conditions, similar to what has been observed in the natural environment, cells were only toxic in the upwelling treatment. We also conducted single-driver experiments to gain a mechanistic understanding of which drivers most impact P. australis growth and toxicity. These experiments indicated that nitrogen concentration and N:P ratio were likely the drivers that most influenced domoic acid production, while the impacts of temperature or pCO2 concentration were less pronounced. Together, these experiments may help to provide both mechanistic and holistic perspectives on toxic P. australis blooms in the dynamic and changing coastal ocean, where cells interact simultaneously with multiple altered environmental variables.

Draft Genome Sequence of Salegentibacter sp. Strain BDJ18, a Plankton-Associated Bacterium in the Northeast Atlantic Ocean.

Salegentibacter sp. strain BDJ18 was isolated from a plankton-associated seawater sample from the northeast Atlantic Ocean. We report its draft genome assembly, which includes genes potentially important for microbial interactions in the marine environment.

Empirical bayes analysis of sequencing-based transcriptional profiling without replicates.

BACKGROUND: Recent technological advancements have made high throughput sequencing an increasingly popular approach for transcriptome analysis. Advantages of sequencing-based transcriptional profiling over microarrays have been reported, including lower technical variability. However, advances in technology do not remove biological variation between replicates and this variation is often neglected in many analyses. RESULTS: We propose an empirical Bayes method, titled Analysis of Sequence Counts (ASC), to detect differential expression based on sequencing technology. ASC borrows information across sequences to establish prior distribution of sample variation, so that biological variation can be accounted for even when replicates are not available. Compared to current approaches that simply tests for equality of proportions in two samples, ASC is less biased towards highly expressed sequences and can identify more genes with a greater log fold change at lower overall abundance. CONCLUSIONS: ASC unifies the biological and statistical significance of differential expression by estimating the posterior mean of log fold change and estimating false discovery rates based on the posterior mean. The implementation in R is available at http://www.stat.brown.edu/Zwu/research.aspx.

Phylogenetic diversity of cyanobacterial narB genes from various marine habitats.

Nitrate, the most abundant combined, dissolved form of inorganic nitrogen in global oceans, is a common source of nitrogen (N) for phytoplankton including cyanobacteria. Using a nested polymerase chain reaction (PCR) method, the diversity of the cyanobacterial nitrate reductase gene, narB, was examined in plankton samples from a variety of marine habitats. A total of 480 narB gene fragment sequences were obtained from a coastal coral reef (Heron Island, Australia), open-ocean tropical and subtropical oceanic waters (Atlantic and Pacific Oceans) and a temperate N. Pacific Ocean site (34 degrees N, 129 degrees W). Phylogenetic analyses distinguished eight picocyanobacterial narB clades comprised of DNA sequences derived from the nutrient-replete coastal, nutrient-deplete pelagic and tidally influenced coral reef habitats. The phylogeny of recovered narB gene sequences was consistent with 16S rRNA and ITS sequence phylogenies, suggesting minimal horizontal gene transfer of the narB gene. Depending on sampled habitat, environmental narB sequence types segregated into three divisions: non-picocyanobacterial, coastal picocyanobacterial and open-ocean picocyanobacterial sequences. Using a reverse transcription PCR method, narB mRNA sequences were amplified from Heron Island samples, indicating that narB expression can be detected in environmental samples.

Characterization of cyanobacterial glnA gene diversity and gene expression in marine environments.

PCR primers were designed and used to amplify glnA, the gene that encodes glutamine synthetase, from pure cultures of cyanobacteria and four samples from different marine environments. The glnA phylogeny was similar to that of the 16S rRNA gene, indicating that glnA gene sequences can be used to identify cyanobacteria expressing the glnA gene. Diverse unicellular cyanobacteria glnA genes were recovered from the North Pacific Subtropical Gyre, Monterey Bay, Chesapeake Bay and waters off the New Jersey coast. The majority of sequences were closely related to sequences from Synechococcus strains (78-88% identical DNA sequences). A few sequences that clustered with Prochlorococcus glnA genes were recovered from Monterey Bay and the North Pacific Subtropical Gyre. The expression of glnA was assayed by reverse transcriptase PCR to determine if there was a daily pattern in gene expression of samples collected from New Jersey's Longterm Environmental Observatory site (LEO-15). glnA expression varied over the day, with different glnA sequence types exhibiting different daily cycles. Results showed that the glnA gene can be used to characterize the diversity of natural populations of cyanobacteria, and to characterize gene expression patterns of individual species or strains.

Preferential depletion of zinc within Costa Rica upwelling dome creates conditions for zinc co-limitation of primary production.

The Costa Rica Dome (CRD) is a wind-driven feature characterized by high primary production and an unusual cyanobacterial bloom in surface waters. It is not clear whether this bloom arises from top-down or bottom-up processes. Several studies have argued that trace metal geochemistry within the CRD contributes to the composition of the phytoplankton assemblages, since cyanobacteria and eukaryotic phytoplankton have different transition metal requirements. Here, we report that total dissolved zinc (Zn) is significantly depleted relative to phosphate (P) and silicate (Si) within the upper water column of the CRD compared with other oceanic systems, and this may create conditions favorable for cyanobacteria, which have lower Zn requirements than their eukaryotic competitors. Shipboard grow-out experiments revealed that while Si was a limiting factor under our experimental conditions, additions of Si and either iron (Fe) or Zn led to higher biomass than Si additions alone. The addition of Fe and Zn alone did not lead to significant enhancements. Our results suggest that the depletion of Zn relative to P in upwelled waters may create conditions in the near-surface waters that favor phytoplankton with low Zn requirements, including cyanobacteria.

Genetic indicators of iron limitation in wild populations of Thalassiosira oceanica from the northeast Pacific Ocean.

Assessing the iron (Fe) nutritional status of natural diatom populations has proven challenging as physiological and molecular responses can differ in diatoms of the same genus. We evaluated expression of genes encoding flavodoxin (FLDA1) and an Fe-starvation induced protein (ISIP3) as indicators of Fe limitation in the marine diatom Thalassiosira oceanica. The specificity of the response to Fe limitation was tested in cultures grown under Fe- and macronutrient-deficient conditions, as well as throughout the diurnal light cycle. Both genes showed a robust and specific response to Fe limitation in laboratory cultures and were detected in small volume samples collected from the northeast Pacific, demonstrating the sensitivity of this method. Overall, FLDA1 and ISIP3 expression was inversely related to Fe concentrations and offered insight into the Fe nutritional health of T. oceanica in the field. As T. oceanica is a species tolerant to low Fe, indications of Fe limitation in T. oceanica populations may serve as a proxy for severe Fe stress in the overall diatom community. At two shallow coastal locations, FLD1A and ISIP3 expression revealed Fe stress in areas where dissolved Fe concentrations were high, demonstrating that this approach may be powerful for identifying regions where Fe supply may not be biologically available.

Nitrogenase genes in PCR and RT-PCR reagents: implications for studies of diversity of functional genes.

Studies of the diversity of microorganisms in the environment have been facilitated by use of PCR and reverse transcription PCR (RT-PCR). Inhibition of the PCR by complex sample matrices and low abundance of some target microorganisms require the use of high-sensitivity amplification procedures, involving a large number of cycles or nested PCR methods. Using these methods, we frequently observed contamination of the amplification reagents, including polymerases, by genomic DNA containing nitrogenase (nifH) and rRNA genes. Contaminating genes were sequenced and found to belong to a variety of rRNA clades, but only three major nifH clades. These sequence types included a few nifH sequences reported in previous studies of the environment. Contamination could be reduced by restriction digestion and ultrafiltration of PCR reagents, but efficiency of amplification was also reduced. Our results suggest that studies relying on large numbers of PCR amplification cycles to assess environmental gene diversity should take precautions to assure that clone libraries generated from amplified PCR products are not the result of contaminated PCR reagents.

Thalassiosira spp. community composition shifts in response to chemical and physical forcing in the northeast Pacific Ocean.

Diatoms are genetically diverse unicellular photosynthetic eukaryotes that are key primary producers in the ocean. Many of the over 100 extant diatom species in the cosmopolitan genus Thalassiosira are difficult to distinguish in mixed populations using light microscopy. Here, we examine shifts in Thalassiosira spp. composition along a coastal to open ocean transect that encountered a 3-month-old Haida eddy in the northeast Pacific Ocean. To quantify shifts in Thalassiosira species composition, we developed a targeted automated ribosomal intergenic spacer analysis (ARISA) method to identify Thalassiosira spp. in environmental samples. As many specific fragment lengths are indicative of individual Thalassiosira spp., the ARISA method is a useful screening tool to identify changes in the relative abundance and distribution of specific species. The method also enabled us to assess changes in Thalassiosira community composition in response to chemical and physical forcing. Thalassiosira spp. community composition in the core of a 3-month-old Haida eddy remained largely (>80%) similar over a 2-week period, despite moving 24 km southwestward. Shifts in Thalassiosira species correlated with changes in dissolved iron (Fe) and temperature throughout the sampling period. Simultaneously tracking community composition and relative abundance of Thalassiosira species within the physical and chemical context they occurred allowed us to identify quantitative linkages between environmental conditions and community response.

Identifying reference genes with stable expression from high throughput sequence data.

Genes that are constitutively expressed across multiple environmental stimuli are crucial to quantifying differentially expressed genes, particularly when employing quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) assays. However, the identification of these potential reference genes in non-model organisms is challenging and is often guided by expression patterns in distantly related organisms. Here, transcriptome datasets from the diatom Thalassiosira pseudonana grown under replete, phosphorus-limited, iron-limited, and phosphorus and iron co-limited nutrient regimes were analyzed through literature-based searches for homologous reference genes, k-means clustering, and analysis of sequence counts (ASC) to identify putative reference genes. A total of 9759 genes were identified and screened for stable expression. Literature-based searches surveyed 18 generally accepted reference genes, revealing 101 homologs in T. pseudonana with variable expression and a wide range of mean tags per million. k-means analysis parsed the whole transcriptome into 15 clusters. The two most stable clusters contained 709 genes, but still had distinct patterns in expression. ASC analyses identified 179 genes that were stably expressed (posterior probability < 0.1 for 1.25 fold change). Genes known to have a stable expression pattern across the test treatments, like actin, were identified in this pool of 179 candidate genes. ASC can be employed on data without biological replicates and was more robust than the k-means approach in isolating genes with stable expression. The intersection of the genes identified through ASC with commonly used reference genes from the literature suggests that actin and ubiquitin ligase may be useful reference genes for T. pseudonana and potentially other diatoms. With the wealth of transcriptome sequence data becoming available, ASC can be easily applied to transcriptome datasets from other phytoplankton to identify reference genes.

The transcriptome and proteome of the diatom Thalassiosira pseudonana reveal a diverse phosphorus stress response.

Phosphorus (P) is a critical driver of phytoplankton growth and ecosystem function in the ocean. Diatoms are an abundant class of marine phytoplankton that are responsible for significant amounts of primary production. With the control they exert on the oceanic carbon cycle, there have been a number of studies focused on how diatoms respond to limiting macro and micronutrients such as iron and nitrogen. However, diatom physiological responses to P deficiency are poorly understood. Here, we couple deep sequencing of transcript tags and quantitative proteomics to analyze the diatom Thalassiosira pseudonana grown under P-replete and P-deficient conditions. A total of 318 transcripts were differentially regulated with a false discovery rate of <0.05, and a total of 136 proteins were differentially abundant (p<0.05). Significant changes in the abundance of transcripts and proteins were observed and coordinated for multiple biochemical pathways, including glycolysis and translation. Patterns in transcript and protein abundance were also linked to physiological changes in cellular P distributions, and enzyme activities. These data demonstrate that diatom P deficiency results in changes in cellular P allocation through polyphosphate production, increased P transport, a switch to utilization of dissolved organic P through increased production of metalloenzymes, and a remodeling of the cell surface through production of sulfolipids. Together, these findings reveal that T. pseudonana has evolved a sophisticated response to P deficiency involving multiple biochemical strategies that are likely critical to its ability to respond to variations in environmental P availability.

Characterization of putative iron responsive genes as species-specific indicators of iron stress in thalassiosiroid diatoms.

Iron (Fe) availability restricts diatom growth and primary production in large areas of the oceans. It is a challenge to assess the bulk Fe nutritional health of natural diatom populations, since species can differ in their physiological and molecular responses to Fe limitation. We assayed expression of selected genes in diatoms from the Thalassiosira genus to assess their potential utility as species-specific molecular markers to indicate Fe status in natural diatom assemblages. In this study, we compared the expression of the photosynthetic genes encoding ferredoxin (a Fe-requiring protein) and flavodoxin (a Fe-free protein) in culture experiments with Fe replete and Fe stressed Thalassiosira pseudonana (CCMP 1335) isolated from coastal waters and Thalassiosira weissflogii (CCMP 1010) isolated from the open ocean. In T. pseudonana, expression of flavodoxin and ferredoxin genes were not sensitive to Fe status but were found to display diel periodicities. In T. weissflogii, expression of flavodoxin was highly responsive to iron levels and was only detectable when cultures were Fe limited. Flavodoxin genes have been duplicated in most diatoms with available genome data and we show that T. pseudonana has lost its copy related to the Fe-responsive copy in T. weissflogii. We also examined the expression of genes for a putative high affinity, copper (Cu)-dependent Fe uptake system in T. pseudonana. Our results indicate that genes encoding putative Cu transporters, a multi-Cu oxidase, and a Fe reductase are not linked to Fe status. The expression of a second putative Fe reductase increased in Fe limited cultures, but this gene was also highly expressed in Fe replete cultures, indicating it may not be a useful marker in the field. Our findings highlight that Fe metabolism may differ among diatoms even within a genus and show a need to validate responses in different species as part of the development pipeline for genetic markers of Fe status in field populations.

Spatial-temporal variability in diazotroph assemblages in Chesapeake Bay using an oligonucleotide nifH microarray.

The distribution of nitrogen-fixing microorganisms in the Chesapeake Bay was investigated using fingerprints from a nifH microarray comprised of 706 60-mer oligonucleotide nifH probes representing cultivated organisms and environmental clones from different nifH clusters. Diverse nifH targets, amplified from samples using degenerate nifH primers, were detected in water column and sediment samples collected in April and October, 2001-2002. Total nifH richness and diversity (Simpson's and Shannon indices) were highest at the most riverine, oligohaline North Bay station. In most samples, the highest diversity was in nifH Cluster 3, which includes many anaerobes, while Cluster 1 (alpha-, beta- gamma- Proteobacteria, Cyanobacteria) targets had the greatest microarray signal intensities. In a multidimensional scaling analysis, deep water communities from April and October were similar within each of the sampling sites, while the surface communities had more variability. Diazotroph communities in the water column in the North Bay were distinct from the Mid- and South Bay communities, and there was a gradual change in sediment diazotroph assemblages from the North to the South Bay. Diazotrophic assemblages from the majority of the water column samples from the Mid- and South Bay clustered with the sediment assemblage in Mid-Bay. Dissolved inorganic nitrogen, salinity, dissolved organic carbon and dissolved organic phosphorus had a significant relationship with the diazotrophic bacterioplankton community. Higher diversity in the freshwater end of the system may reflect variability in disturbance rates and environmental conditions such as forms and concentrations of organic matter, nutrients and oxygen.

Cyanobacterial assimilatory nitrate reductase gene diversity in coastal and oligotrophic marine environments.

Cyanobacteria are important primary producers in many marine ecosystems and their abundances and growth rates depend on their ability to assimilate various nitrogen sources. To examine the diversity of nitrate-utilizing marine cyanobacteria, we developed PCR primers specific for cyanobacterial assimilatory nitrate reductase (narB) genes. We obtained amplification products from diverse strains of cultivated cyanobacteria and from several marine environments. Phylogenetic trees constructed with the narB gene are congruent with those based on ribosomal RNA genes and RNA polymerase genes. Analysis of sequence library data from coastal and oligotrophic marine environments shows distinct groups of Synechococcus sp. in each environment; some of which are represented by sequences from cultivated organisms and others that are unrelated to known sequences and likely represent novel phylogenetic groups. We observed spatial differences in the distribution of sequences between two sites in Monterey Bay and differences in the vertical distribution of sequence types at the Hawai'i Ocean Time-series Station ALOHA, suggesting that nitrogen assimilation in Synechococcus living in different ecological niches can be followed with the nitrate reductase gene.