SSU-rRNA Gene Sequencing Survey of Benthic Microbial Eukaryotes from Guaymas Basin Hydrothermal Vent.

Microbial eukaryotes have important roles in marine food webs, but their diversity and activities in hydrothermal vent ecosystems are poorly characterized. In this study, we analyzed microbial eukaryotic communities associated with bacterial (Beggiatoa) mats in the 2,000 m deep-sea Guaymas Basin hydrothermal vent system using 18S rRNA gene high-throughput sequencing of the V4 region. We detected 6,954 distinct Operational Taxonomic Units (OTUs) across various mat systems. Of the sequences that aligned with known protistan phylotypes, most were affiliated with alveolates (especially dinoflagellates and ciliates) and cercozoans. OTU richness and community structure differed among sediment habitats (e.g. different mat types and cold sediments away from mats). Additionally, full-length 18S rRNA genes amplified and cloned from single cells revealed the identities of some of the most commonly encountered, active ciliates in this hydrothermal vent ecosystem. Observations and experiments were also conducted to demonstrate that ciliates were trophically active and ingesting fluorescent bacteria or Beggiatoa trichomes. Our work suggests that the active and diverse protistan community at the Guaymas Basin hydrothermal vent ecosystem likely consumes substantial amounts of bacterial biomass, and that the different habitats, often defined by distances of just a few 10s of cm, select for particular assemblages and levels of diversity.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone.

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

Effect of light and prey availability on gene expression of the mixotrophic chrysophyte, Ochromonas sp.

BACKGROUND: Ochromonas is a genus of mixotrophic chrysophytes that is found ubiquitously in many aquatic environments. Species in this genus can be important consumers of bacteria but vary in their ability to perform photosynthesis. We studied the effect of light and bacteria on growth and gene expression of a predominantly phagotrophic Ochromonas species. Axenic cultures of Ochromonas sp. were fed with heat-killed bacteria (HKB) and grown in constant light or darkness. RNA was extracted from cultures in the light or in the dark with HKB present (Light + HKB; Dark + HKB), and in the light after HKB were depleted (Light + depleted HKB). RESULTS: There were no significant differences in the growth or bacterial ingestion rates between algae grown in light or dark conditions. The availability of light led to a differential expression of only 8% of genes in the transcriptome. A number of genes associated with photosynthesis, phagotrophy, and tetrapyrrole synthesis was upregulated in the Light + HKB treatment compared to Dark + HKB. Conversely, the comparison between the Light + HKB and Light + depleted HKB treatments revealed that the presence of HKB led to differential expression of 59% of genes, including the majority of genes involved in major carbon and nitrogen metabolic pathways. Genes coding for unidirectional enzymes for the utilization of glucose were upregulated in the presence of HKB, implying increased glycolytic activities during phagotrophy. Algae without HKB upregulated their expression of genes coding for ammonium transporters, implying uptake of inorganic nitrogen from the culture medium when prey were unavailable. CONCLUSIONS: Transcriptomic results agreed with previous observations that light had minimal effect on the population growth of Ochromonas sp. However, light led to the upregulation of a number of phototrophy- and phagotrophy-related genes, while the availability of bacterial prey led to prominent changes in major carbon and nitrogen metabolic pathways. Our study demonstrated the potential of transcriptomic approaches to improve our understanding of the trophic physiologies of complex mixotrophs, and revealed responses in Ochromonas sp. not apparent from traditional culture studies.

De novo sequences of Haloquadratum walsbyi from Lake Tyrrell, Australia, reveal a variable genomic landscape.

Hypersaline systems near salt saturation levels represent an extreme environment, in which organisms grow and survive near the limits of life. One of the abundant members of the microbial communities in hypersaline systems is the square archaeon, Haloquadratum walsbyi. Utilizing a short-read metagenome from Lake Tyrrell, a hypersaline ecosystem in Victoria, Australia, we performed a comparative genomic analysis of H. walsbyi to better understand the extent of variation between strains/subspecies. Results revealed that previously isolated strains/subspecies do not fully describe the complete repertoire of the genomic landscape present in H. walsbyi. Rearrangements, insertions, and deletions were observed for the Lake Tyrrell derived Haloquadratum genomes and were supported by environmental de novo sequences, including shifts in the dominant genomic landscape of the two most abundant strains. Analysis pertaining to halomucins indicated that homologs for this large protein are not a feature common for all species of Haloquadratum. Further, we analyzed ATP-binding cassette transporters (ABC-type transporters) for evidence of niche partitioning between different strains/subspecies. We were able to identify unique and variable transporter subunits from all five genomes analyzed and the de novo environmental sequences, suggesting that differences in nutrient and carbon source acquisition may play a role in maintaining distinct strains/subspecies.

Virus-host and CRISPR dynamics in Archaea-dominated hypersaline Lake Tyrrell, Victoria, Australia.

The study of natural archaeal assemblages requires community context, namely, a concurrent assessment of the dynamics of archaeal, bacterial, and viral populations. Here, we use filter size-resolved metagenomic analyses to report the dynamics of 101 archaeal and bacterial OTUs and 140 viral populations across 17 samples collected over different timescales from 2007-2010 from Australian hypersaline Lake Tyrrell (LT). All samples were dominated by Archaea (75-95%). Archaeal, bacterial, and viral populations were found to be dynamic on timescales of months to years, and different viral assemblages were present in planktonic, relative to host-associated (active and provirus) size fractions. Analyses of clustered regularly interspaced short palindromic repeat (CRISPR) regions indicate that both rare and abundant viruses were targeted, primarily by lower abundance hosts. Although very few spacers had hits to the NCBI nr database or to the 140 LT viral populations, 21% had hits to unassembled LT viral concentrate reads. This suggests local adaptation to LT-specific viruses and/or undersampling of haloviral assemblages in public databases, along with successful CRISPR-mediated maintenance of viral populations at abundances low enough to preclude genomic assembly. This is the first metagenomic report evaluating widespread archaeal dynamics at the population level on short timescales in a hypersaline system.

New approaches indicate constant viral diversity despite shifts in assemblage structure in an Australian hypersaline lake.

It is widely stated that viruses represent the most significant source of biodiversity on Earth, yet characterizing the diversity of viral assemblages in natural systems remains difficult. Viral diversity studies are challenging because viruses lack universally present, phylogenetically informative genes. Here, we developed an approach to estimate viral diversity using a series of functional and novel conserved genes. This approach provides direct estimates of viral assemblage diversity while retaining resolution at the level of individual viral populations in a natural system. We characterized viral assemblages in eight samples from hypersaline Lake Tyrrell (LT), Victoria, Australia, using 39,636 viral contigs. We defined viral operational taxonomic units (OTUs) in two ways. First, we used genes with three different functional predictions that were abundantly represented in the data set. Second, we clustered proteins of unknown function based on sequence similarity, and we chose genes represented by three clusters with numerous members to define OTUs. In combination, diversity metrics indicated between 412 and 735 sampled populations, and the number of populations remained relatively constant across samples. We determined the relative representation of each viral OTU in each sample and found that viral assemblage structures correlate with salinity and solution chemistry. LT viral assemblages were near-replicates from the same site sampled a few days apart but differed significantly on other spatial and temporal scales. The OTU definition approach proposed here paves the way for metagenomics-based analyses of viral assemblages using ecological models previously applied to bacteria and archaea.

Temporal profiling resolves the drivers of microbial nitrogen cycling variability in coastal sediments.

Here we describe the potential for sediment microbial nitrogen-cycling gene (DNA) and activity (RNA) abundances to spatially resolve coastal areas impacted by seasonal variability in external nutrient inputs. Three sites were chosen within a nitrogen-limited embayment, Port Phillip Bay (PPB), Australia that reflect variability in both proximity to external nutrient inputs and the dominant form of available nitrogen. At three sediment depths (0-1; 1-5; 5-10 cm) across a 2 year study key genes involved in nitrification (archaeal amoA and bacterial ß-amoA), nitrite reduction (clade I nirS and cluster I nirK, archaeal nirK-a), anaerobic oxidation of ammonium (anammox 16S rRNA phylogenetic marker) and nitrogen fixation (nifH) were quantified. Sediments impacted by a dominance of organic nitrogen inputs were characterised at all time-points and to sediment depths of 10 cm by the highest transcript abundances of archaeal amoA and archaeal nirk-a. Proximity to a dominance of external nitrate inputs was associated with the highest transcript abundances of nirS which temporally co-varied with seasonal changes in sediment nitrate. Sediments isolated from external inputs displayed the greatest depth-specific decrease in quantifiable transcript abundances. In these isolated sediments bacterial ß-amoA transcripts were temporally associated with increased sediment ammonium levels. Across this nitrogen limited system variability in the abundance of bacterial ß-amoA, archaeal amoA, archaeal nirk-a or nirS transcripts from the sediment surface (0-1 and 5 cm) demonstrated a capacity to improve our ability to monitor coastal zones impacted by anthropogenic nitrogen inputs. Specifically, the spatial detection sensitivity of bacterial ß-amoA transcripts could be developed as a metric to determine spatiotemporal impacts of large external loading events. This temporal study demonstrates a capacity for microbial activity metrics to facilitate coastal management strategies through greater spatial resolution of areas impacted by external nutrient inputs.

Using metatranscriptomics to better understand the role of microbial nitrogen cycling in coastal sediment benthic flux denitrification efficiency.

Spatial and temporal variability in benthic flux denitrification efficiency occurs across Port Phillip Bay, Australia. Here, we assess the capacity for untargeted metatranscriptomics to resolve spatiotemporal differences in the microbial contribution to benthic nitrogen cycling. The most abundant sediment transcripts assembled were associated with the archaeal nitrifier Nitrosopumilus. In sediments close to external inputs of organic nitrogen, the dominant transcripts were associated with Nitrosopumilus nitric oxide nitrite reduction (nirK). The environmental conditions close to organic nitrogen inputs that select for increased transcription in Nitrosopumilus (amoCAB, nirK, nirS, nmo, hcp) additionally selected for increased transcription of bacterial nitrite reduction (nxrB) and transcripts associated with anammox (hzo) but not denitrification (bacterial nirS/nirk). In sediments that are more isolated from external inputs of organic nitrogen dominant transcripts were associated with nitrous oxide reduction (nosZ) and changes in nosZ transcript abundance were uncoupled from transcriptional profiles associated with archaeal nitrification. Coordinated transcription of coupled community-level nitrification-denitrification was not well supported by metatranscriptomics. In comparison, the abundance of archaeal nirK transcripts were site- and season-specific. This study indicates that the transcription of archaeal nirK in response to changing environmental conditions may be an important and overlooked feature of coastal sediment nitrogen cycling.

Comparative analysis of eukaryotic marine microbial assemblages from 18S rRNA gene and gene transcript clone libraries by using different methods of extraction.

Eukaryotic marine microbes play pivotal roles in biogeochemical nutrient cycling and ecosystem function, but studies that focus on the protistan biogeography and genetic diversity lag-behind studies of other microbes. 18S rRNA PCR amplification and clone library sequencing are commonly used to assess diversity that is culture independent. However, molecular methods are not without potential biases and artifacts. In this study, we compare the community composition of clone libraries generated from the same water sample collected at the San Pedro Ocean Time Series (SPOTs) station in the northwest Pacific Ocean. Community composition was assessed using different cell lysis methods (chemical and mechanical) and the extraction of different nucleic acids (DNA and RNA reverse transcribed to cDNA) to build Sanger ABI clone libraries. We describe specific biases for ecologically important phylogenetic groups resulting from differences in nucleic acid extraction methods that will inform future designs of eukaryotic diversity studies, regardless of the target sequencing platform planned.

Dynamic viral populations in hypersaline systems as revealed by metagenomic assembly.

Viruses of the Bacteria and Archaea play important roles in microbial evolution and ecology, and yet viral dynamics in natural systems remain poorly understood. Here, we created de novo assemblies from 6.4 Gbp of metagenomic sequence from eight community viral concentrate samples, collected from 12 h to 3 years apart from hypersaline Lake Tyrrell (LT), Victoria, Australia. Through extensive manual assembly curation, we reconstructed 7 complete and 28 partial novel genomes of viruses and virus-like entities (VLEs, which could be viruses or plasmids). We tracked these 35 populations across the eight samples and found that they are generally stable on the timescale of days and transient on the timescale of years, with some exceptions. Cross-detection of the 35 LT populations in three previously described haloviral metagenomes was limited to a few genes, and most previously sequenced haloviruses were not detected in our samples, though 3 were detected upon reducing our detection threshold from 90% to 75% nucleotide identity. Similar results were obtained when we applied our methods to haloviral metagenomic data previously reported from San Diego, CA: 10 contigs that we assembled from that system exhibited a variety of detection patterns on a timescale of weeks to 1 month but were generally not detected in LT. Our results suggest that most haloviral populations have a limited or, possibly, a temporally variable global distribution. This study provides high-resolution insight into viral biogeography and dynamics and it places "snapshot" viral metagenomes, collected at a single time and location, in context.

The Sorcerer II Global Ocean Sampling expedition: expanding the universe of protein families.

Metagenomics projects based on shotgun sequencing of populations of micro-organisms yield insight into protein families. We used sequence similarity clustering to explore proteins with a comprehensive dataset consisting of sequences from available databases together with 6.12 million proteins predicted from an assembly of 7.7 million Global Ocean Sampling (GOS) sequences. The GOS dataset covers nearly all known prokaryotic protein families. A total of 3,995 medium- and large-sized clusters consisting of only GOS sequences are identified, out of which 1,700 have no detectable homology to known families. The GOS-only clusters contain a higher than expected proportion of sequences of viral origin, thus reflecting a poor sampling of viral diversity until now. Protein domain distributions in the GOS dataset and current protein databases show distinct biases. Several protein domains that were previously categorized as kingdom specific are shown to have GOS examples in other kingdoms. About 6,000 sequences (ORFans) from the literature that heretofore lacked similarity to known proteins have matches in the GOS data. The GOS dataset is also used to improve remote homology detection. Overall, besides nearly doubling the number of current proteins, the predicted GOS proteins also add a great deal of diversity to known protein families and shed light on their evolution. These observations are illustrated using several protein families, including phosphatases, proteases, ultraviolet-irradiation DNA damage repair enzymes, glutamine synthetase, and RuBisCO. The diversity added by GOS data has implications for choosing targets for experimental structure characterization as part of structural genomics efforts. Our analysis indicates that new families are being discovered at a rate that is linear or almost linear with the addition of new sequences, implying that we are still far from discovering all protein families in nature.

The Sorcerer II Global Ocean Sampling expedition: northwest Atlantic through eastern tropical Pacific.

The world's oceans contain a complex mixture of micro-organisms that are for the most part, uncharacterized both genetically and biochemically. We report here a metagenomic study of the marine planktonic microbiota in which surface (mostly marine) water samples were analyzed as part of the Sorcerer II Global Ocean Sampling expedition. These samples, collected across a several-thousand km transect from the North Atlantic through the Panama Canal and ending in the South Pacific yielded an extensive dataset consisting of 7.7 million sequencing reads (6.3 billion bp). Though a few major microbial clades dominate the planktonic marine niche, the dataset contains great diversity with 85% of the assembled sequence and 57% of the unassembled data being unique at a 98% sequence identity cutoff. Using the metadata associated with each sample and sequencing library, we developed new comparative genomic and assembly methods. One comparative genomic method, termed "fragment recruitment," addressed questions of genome structure, evolution, and taxonomic or phylogenetic diversity, as well as the biochemical diversity of genes and gene families. A second method, termed "extreme assembly," made possible the assembly and reconstruction of large segments of abundant but clearly nonclonal organisms. Within all abundant populations analyzed, we found extensive intra-ribotype diversity in several forms: (1) extensive sequence variation within orthologous regions throughout a given genome; despite coverage of individual ribotypes approaching 500-fold, most individual sequencing reads are unique; (2) numerous changes in gene content some with direct adaptive implications; and (3) hypervariable genomic islands that are too variable to assemble. The intra-ribotype diversity is organized into genetically isolated populations that have overlapping but independent distributions, implying distinct environmental preference. We present novel methods for measuring the genomic similarity between metagenomic samples and show how they may be grouped into several community types. Specific functional adaptations can be identified both within individual ribotypes and across the entire community, including proteorhodopsin spectral tuning and the presence or absence of the phosphate-binding gene PstS.

A tale of two mixotrophic chrysophytes: Insights into the metabolisms of two Ochromonas species (Chrysophyceae) through a comparison of gene expression.

Ochromonas spp. strains CCMP1393 and BG-1 are phagotrophic phytoflagellates with different nutritional strategies. Strain CCMP1393 is an obligate phototroph while strain BG-1 readily grows in continuous darkness in the presence of bacterial prey. Growth and gene expression of strain CCMP1393 were investigated under conditions allowing phagotrophic, mixotrophic, or phototrophic nutrition. The availability of light and bacterial prey led to the differential expression of 42% or 45-59% of all genes, respectively. Data from strain CCMP1393 were compared to those from a study conducted previously on strain BG-1, and revealed notable differences in carbon and nitrogen metabolism between the 2 congeners under similar environmental conditions. Strain BG-1 utilized bacterial carbon and amino acids through glycolysis and the tricarboxylic acid cycle, while downregulating light harvesting and carbon fixation in the Calvin cycle when both light and bacteria were available. In contrast, the upregulation of genes related to photosynthesis, light harvesting, chlorophyll synthesis, and carbon fixation in the presence of light and prey for strain CCMP1393 implied that this species is more phototrophic than strain BG-1, and that phagotrophy may have enhanced phototrophy. Cellular chlorophyll a content was also significantly higher in strain CCMP1393 supplied with bacteria compared to those without prey. Our results thus point to very different physiological strategies for mixotrophic nutrition in these closely related chrysophyte species.

Autotrophic and heterotrophic acquisition of carbon and nitrogen by a mixotrophic chrysophyte established through stable isotope analysis.

Collectively, phagotrophic algae (mixotrophs) form a functional continuum of nutritional modes between autotrophy and heterotrophy, but the specific physiological benefits of mixotrophic nutrition differ among taxa. Ochromonas spp. are ubiquitous chrysophytes that exhibit high nutritional flexibility, although most species generally fall towards the heterotrophic end of the mixotrophy spectrum. We assessed the sources of carbon and nitrogen in Ochromonas sp. strain BG-1 growing mixotrophically via short-term stable isotope probing. An axenic culture was grown in the presence of either heat-killed bacteria enriched with 15N and 13C, or unlabeled heat-killed bacteria and labeled inorganic substrates (13C-bicarbonate and 15N-ammonium). The alga exhibited high growth rates (up to 2 divisions per day) only until heat-killed bacteria were depleted. NanoSIMS and bulk IRMS isotope analyses revealed that Ochromonas obtained 84-99% of its carbon and 88-95% of its nitrogen from consumed bacteria. The chrysophyte assimilated inorganic 13C-carbon and 15N-nitrogen when bacterial abundances were very low, but autotrophic (photosynthetic) activity was insufficient to support net population growth of the alga. Our use of nanoSIMS represents its first application towards the study of a mixotrophic alga, enabling a better understanding and quantitative assessment of carbon and nutrient acquisition by this species.

Single-cell transcriptomics of small microbial eukaryotes: limitations and potential.

Single-cell transcriptomics is an emerging research tool that has huge untapped potential in the study of microbial eukaryotes. Its application has been tested in microbial eukaryotes 50 µm or larger, and it generated transcriptomes similar to those obtained from culture-based RNA-seq. However, microbial eukaryotes have a wide range of sizes and can be as small as 1 µm. Single-cell RNA-seq was tested in two smaller protists (8 and 15 µm). Transcript recovery rate was much lower and randomness in observed gene expression levels was much higher in single-cell transcriptomes than those derived from bulk cultures of cells. We found that the reason of such observation is that the smaller organisms had much lower mRNA copy numbers. We discuss the application of single-cell RNA-seq in studying smaller microbial eukaryotes in the context of these limitations.

Probing the evolution, ecology and physiology of marine protists using transcriptomics.

Protists, which are single-celled eukaryotes, critically influence the ecology and chemistry of marine ecosystems, but genome-based studies of these organisms have lagged behind those of other microorganisms. However, recent transcriptomic studies of cultured species, complemented by meta-omics analyses of natural communities, have increased the amount of genetic information available for poorly represented branches on the tree of eukaryotic life. This information is providing insights into the adaptations and interactions between protists and other microorganisms and macroorganisms, but many of the genes sequenced show no similarity to sequences currently available in public databases. A better understanding of these newly discovered genes will lead to a deeper appreciation of the functional diversity and metabolic processes in the ocean. In this Review, we summarize recent developments in our understanding of the ecology, physiology and evolution of protists, derived from transcriptomic studies of cultured strains and natural communities, and discuss how these novel large-scale genetic datasets will be used in the future.

Microbiomes of Muricea californica and M. fruticosa: Comparative Analyses of Two Co-occurring Eastern Pacific Octocorals.

Octocorals are sources of novel but understudied microbial diversity. Conversely, scleractinian or reef-building coral microbiomes have been heavily examined in light of the threats of climate change. Muricea californica and Muricea fruticosa are two co-occurring species of gorgonian octocoral abundantly found in the kelp forests of southern California, and thus provide an excellent basis to determine if octocoral microbiomes are host specific. Using Illumina MiSeq amplicon sequencing and replicate samples, we evaluated the microbiomes collected from multiple colonies of both species of Muricea to measure both inter- and intra-colony microbiome variabilities. In addition, microbiomes from overlying sea water and nearby zoanthids (another benthic invertebrate) were also included in the analysis to evaluate whether bacterial taxa specifically associate with octocorals. This is also the first report of microbiomes from these species of Muricea. We show that microbiomes isolated from each sample type are distinct, and specifically, that octocoral species type had the greatest effect on predicting the composition of the Muricea microbiome. Bacterial taxa contributing to compositional differences include distinct strains of Mycoplasma associated with either M. californica or M. fruticosa, an abundance of Spirochaetes observed on M. californica, and a greater diversity of γ-Proteobacteria associated with M. fruticosa. Many of the bacterial taxa contributing to these differences are known for their presence in photosymbiont-containing invertebrate microbiomes.

Gene expression characterizes different nutritional strategies among three mixotrophic protists.

Mixotrophic protists, i.e. protists that can carry out both phototrophy and heterotrophy, are a group of organisms with a wide range of nutritional strategies. The ecological and biogeochemical importance of these species has recently been recognized. In this study, we investigated and compared the gene expression of three mixotrophic protists, Prymnesium parvum, Dinobyron sp. and Ochromonas sp. under light and dark conditions in the presence of prey using RNA-Seq. Gene expression of the obligately phototrophic P. parvum and Dinobryon sp. changed significantly between light and dark treatments, while that of primarily heterotrophic Ochromonas sp. was largely unchanged. Gene expression of P. parvum and Dinobryon sp. shared many similarities, especially in the expression patterns of genes related to reproduction. However, key genes involved in central carbon metabolism and phagotrophy had different expression patterns between these two species, suggesting differences in prey consumption and heterotrophic nutrition in the dark. Transcriptomic data also offered clues to other physiological traits of these organisms such as preference of nitrogen sources and photo-oxidative stress. These results provide potential target genes for further exploration of the mechanisms of mixotrophic physiology and demonstrate the potential usefulness of molecular approaches in characterizing the nutritional modes of mixotrophic protists.

Gene expression in the mixotrophic prymnesiophyte, Prymnesium parvum, responds to prey availability.

The mixotrophic prymnesiophyte, Prymnesium parvum, is a widely distributed alga with significant ecological importance. It produces toxins and can form ecosystem disruptive blooms that result in fish kills and changes in planktonic food web structure. However, the relationship between P. parvum and its prey on the molecular level is poorly understood. In this study, we used RNA-Seq technology to study changes in gene transcription of P. parvum in three treatments with different microbial populations available as potential prey: axenic P. parvum (no prey), bacterized P. paruvm, and axenic P. parvum with ciliates added as prey. Thousands of genes were differentially expressed among the three treatments. Most notably, transcriptome data indicated that P. parvum obtained organic carbon, including fatty acids, from both bacteria and ciliate prey for energy and cellular building blocks. The data also suggested that different prey provided P. parvum with macro- and micro-nutrients, namely organic nitrogen in the form of amino acids from ciliates, and iron from bacteria. However, both transcriptomic data and growth experiments indicated that P. parvum did not grow faster in the presence of prey despite the gains in nutrients, although algal abundances attained in culture were slightly greater in the presence of prey. The relationship between phototrophy, heterotrophy and growth of P. parvum is discussed.

Changes in gene expression of Prymnesium parvum induced by nitrogen and phosphorus limitation.

Prymnesium parvum is a globally distributed prymnesiophyte alga commonly found in brackish water marine ecosystems and lakes. It possesses a suite of toxins with ichthyotoxic, cytotoxic and hemolytic effects which, along with its mixotrophic nutritional capabilities, allows it to form massive Ecosystem Disruptive Algal Blooms (EDABs). While blooms of high abundance coincide with high levels of nitrogen (N) and phosphorus (P), reports of field and laboratory studies have noted that P. parvum toxicity appears to be augmented at high N:P ratios or P-limiting conditions. Here we present the results of a comparative analysis of P. parvum RNA-Seq transcriptomes under nutrient replete conditions, and N or P deficiency to understand how this organism responds at the transcriptional level to varying nutrient conditions. In nutrient limited conditions we found diverse transcriptional responses for genes involved in nutrient uptake, protein synthesis and degradation, photosynthesis, and toxin production. As anticipated, when either N or P was limiting, transcription levels of genes encoding transporters for the respective nutrient were higher than those under replete condition. Ribosomal and lysosomal protein genes were expressed at higher levels under either nutrient-limited condition compared to the replete condition. Photosynthesis genes and polyketide synthase genes were more highly expressed under P-limitation but not under N-limitation. These results highlight the ability of P. parvum to mount a coordinated and varied cellular and physiological response to nutrient limitation. Results also provide potential marker genes for further evaluating the physiological response and toxin production of P. parvum populations during bloom formation or to changing environmental conditions.