Sequence data are available at the NCBI Sequence Read Archive (SRA) under the identifier SRA: SRP069841.

Introduction {#sec001}
============

Organisms in earth's photic zone experience diurnal cycles that influence physiology, metabolism, and growth in different ways. Photosynthetic organisms are particularly affected by light cycles since they depend on energy acquired from light harvesting to drive the biosynthesis of macromolecules essential for growth and metabolic demands in the dark. For these organisms, survival depends not only on sufficient illumination, but also on adequate supply of carbon dioxide and both macro and micronutrients. In marine ecosystems, availability of the micronutrient iron (Fe) plays a major role in controlling primary productivity, as Fe limitation is involved in regulating phytoplankton ecophysiology (at least seasonally) in 30% to 60% of sunlit ocean waters \[[@pgen.1006490.ref001]--[@pgen.1006490.ref004]\]. Large-scale Fe enrichment experiments typically result in diatom dominance of the phytoplankton biomass within days \[[@pgen.1006490.ref005]\], implying that low density diatom stocks are maintained at low Fe concentrations \[[@pgen.1006490.ref006]\], and indicating that certain diatoms have unique adaptive mechanisms to cope with low Fe poising them for rapid growth and ecological dominance upon Fe addition.

Fe is essential for the metabolism and growth of all organisms and diatoms respond in a variety of ways to Fe-limitation. For example, open-ocean species display lower concentrations of photosystem I, cytochrome b~6~f \[[@pgen.1006490.ref007]\] and cellular pigments at the cost of light capture efficiency \[[@pgen.1006490.ref008]\] and/or substitute Fe-containing enzymes with Fe-free equivalents, such as ferredoxin with flavodoxin or cytochrome b6 with plastocyanin \[[@pgen.1006490.ref009]--[@pgen.1006490.ref013]\]. Additionally, diatoms appear to induce reductive \[[@pgen.1006490.ref014]--[@pgen.1006490.ref016]\] and potentially non-reductive \[[@pgen.1006490.ref017]\] high-affinity transport systems to acquire Fe. Some Fe-acquisition and homeostasis mechanisms have been shown to be dependent on the light:dark cycle in phytoplankton \[[@pgen.1006490.ref018]\]. Fe stress also seems to impose a strain on overall proteome health that cells counteract by induction of proteolysis, chaperone expression and various reactive oxygen species (ROS) defense mechanisms \[[@pgen.1006490.ref019]--[@pgen.1006490.ref023]\]. The unique redox capabilities of Fe-containing metalloproteins are particularly essential in both photosynthetic and respiratory electron transport chains \[[@pgen.1006490.ref024]--[@pgen.1006490.ref026]\]. Therefore, decreases in Fe would be expected to impair cellular energetic processes during both light phases (photosynthesis) and dark phases (catabolism, respiration). Previous studies on the low Fe physiology of diatoms have been largely conducted under continuous illumination and/or sampled single time points and it is currently unknown how Fe-limitation affects cellular processes including cell division and metabolic shifts over diel cycles.

High resolution sampling has revealed the extent to which physiological states vary over diel cycles and how these state changes are accompanied and regulated by transcriptome shifts in phytoplankton \[[@pgen.1006490.ref027]--[@pgen.1006490.ref030]\]. Here, we provide the first global cellular insights into metabolic shifts between diel cycles and at varying Fe concentrations in the model pennate diatom *Phaeodactylum tricornutum* using a time course sampling the physiology, metabolome, and transcriptome at high resolution. We investigate the effect of steady-state Fe-limitation on the temporal expression of genes involved in diverse cellular processes including Fe acquisition, endocytosis, cell cycle progression and division, ribosome biogenesis, light harvesting complex assembly, and primary carbon and nitrogen metabolism. We greatly expand the catalog of highly Fe sensitive genes putatively involved in novel Fe acquisition mechanisms \[[@pgen.1006490.ref019]\], and identify specific transcription factors and gene regulatory mechanisms utilized by the model diatom *P*. *tricornutum* to adapt to Fe-limitation. Though *P*. *tricornutum* is not ecologically dominant, it is a useful model for Fe studies in diatoms because it has a low susceptibility to Fe limitation, grows well at ecologically relevant Fe concentrations, and has the highly iron responsive gene set possessed by oceanic but not coastal diatoms \[[@pgen.1006490.ref014],[@pgen.1006490.ref020]\]. Taken together, these data provide an in-depth view of the global cellular processes operating over the diel cycle, leading to physiological acclimation to low Fe conditions.

Results and Discussion {#sec002}
======================

Fe-limited response of *P*. *tricornutum* {#sec003}
-----------------------------------------

Steady-state cultures of *P*. *tricornutum* were grown at three different iron (Fe) concentrations, two of which were limiting at 20 pM and 40 pM Fe′ (low Fe) in biological duplicate, and non-limiting Fe conditions at 400 pM Fe′ (high Fe) in biological triplicate over a 12:12 light:dark cycle. Cultures were sampled over a complete diel cycle, every 4 hr for 28 hr (7 time points) resulting in a total of 49 time points sampled ([S1 Table](#pgen.1006490.s023){ref-type="supplementary-material"}). Relative to high Fe cultures, Fe limitation reduced growth rates by 25 and 28% at 40 pM Fe′ and 20 pM Fe′ respectively ([Table 1](#pgen.1006490.t001){ref-type="table"}, [S1](#pgen.1006490.s006){ref-type="supplementary-material"}--[S7](#pgen.1006490.s012){ref-type="supplementary-material"} Figs). At all Fe conditions, cell division peaked during the dark phase (2--6 AM; [Fig 1A](#pgen.1006490.g001){ref-type="fig"}). Generally, cell division in phytoplankton is entrained to light:dark cycles, however several diatom species including *P*. *tricornutum*, have been difficult to synchronize based on light:dark cycles alone \[[@pgen.1006490.ref030]--[@pgen.1006490.ref036]\]. These observations are consistent with the relatively small degree of growth synchronization observed in our cultures where percentage of dividing cells reached a maximum of c. 20% ([Fig 1A](#pgen.1006490.g001){ref-type="fig"}).

![Physiological response of *P*. *tricornutum* to low Fe across diel cycles.\
Diel patterns in cell division and fluctuations of F~v~/F~m~ in *P*. *tricornutum* grown at 20 pM, 40 pM, and 400 pM Fe′ over 12:12 L:D cycles, lights on = 9am, lights off = 9pm. (**A**) Bar plot of the time point averaged percentage of dividing cells shown with standard deviation across biological replicates. (**B**) Time point averaged F~v~/F~m~. Legend is shared for both plots and shaded background indicates dark time points.](pgen.1006490.g001){#pgen.1006490.g001}

10.1371/journal.pgen.1006490.t001

###### Physiological response of *P*. *tricornutum* to varying Fe and illumination.

Average growth rates (biological duplicate for 20 pM and 40 pM Fe′, biological triplicate for 400 pM Fe′ with standard deviation), average divisions per day, and F~v~/F~m~.

![](pgen.1006490.t001){#pgen.1006490.t001g}

                                                 Low Fe            High Fe           
  ---------------------------------------------- ----------------- ----------------- -----------------
  **Specific growth rate (μ)**                   0.57 (± 0.14)     0.81 (±0.02)      1.07 (±0.08)
  **Divisions per day**                          0.82              1.16              1.54
  **Average F**~**v**~**/F**~**m**~ **(**± SD)   0.412 (± 0.070)   0.512 (± 0.075)   0.642 (± 0.011)

In addition to reduction in growth rates, Fe availability impacted F~v~/F~m~, (the maximum potential quantum efficiency of PSII). Average F~v~/F~m~ of *P*. *tricornutum* grown at 400 pM Fe′ was near the theoretical maximum (F~v~/F~m~ = 0.642), indicating high photosynthetic efficiencies, which decreased on average by 20% and 36% (in 20 and 40 pM Fe′ respectively; [Table 1](#pgen.1006490.t001){ref-type="table"}).

Here, reductions in F~v~/F~m~ result from both the reduced capacity to synthesize Fe rich components of the photosynthetic electron transport chain (PSII, cytochrome *b*~*6*~*f*, cytochrome oxidase, PSI and ferredoxin), and the production of photoprotective chlorophyll-protein complexes that are energetically disconnected from the photosynthetic electron transport chain particularly under non-limiting macronutrient concentrations \[[@pgen.1006490.ref004], [@pgen.1006490.ref037]--[@pgen.1006490.ref039]\]. Diel fluctuations in F~v~/F~m~ were consistent in all Fe conditions with a mid-day minimum (likely due to photoinhibition), maxima at 10PM (2 h after transition to darkness), followed by a significant decline by 2AM ([Fig 1B](#pgen.1006490.g001){ref-type="fig"}). Nocturnal decreases in F~v~/F~m~ are caused by chlororespiration mediated reduction in the plastoquinone (PQ) pool and resulting back-transfer of electrons which decreases the overall efficiency of electron transport at PSII \[[@pgen.1006490.ref038]\], corroborated by significant rates of chlororespiration previously measured in *P*. *tricornutum* \[[@pgen.1006490.ref040]--[@pgen.1006490.ref043]\]. There was a modest recovery of F~v~/F~m~ before dawn across Fe concentrations, which occurs as cells deplete the pool of respiratory substrates \[[@pgen.1006490.ref038]\]. Together, growth rate measurements and physiological fluorescence confirm that *P*. *tricornutum* grown at 20--40 pM Fe′ experience Fe limitation, and are energetically impaired relative to the Fe replete (400 pM Fe′) cultures. These physiological patterns are further explored in the context of diel patterns in metabolite and transcript profiles.

Global transcriptome dynamics and diel expression of highly Fe responsive genes {#sec004}
-------------------------------------------------------------------------------

The transcriptome was analyzed using two principal methods, expression clustering and differential expression (DE) binning (Methods). First, we explored genome wide co-expression by using Weighted Gene Co-Expression Network Analysis of RPKMs (WGCNA; \[[@pgen.1006490.ref044]--[@pgen.1006490.ref046]\]). WGCNA yielded 27 co-expression clusters or "modules" from a filtered set of the active transcriptome (n = 8499) and containing only nuclear-encoded genes, with modules ranging in size from 81 to 937 genes ([Fig 2](#pgen.1006490.g002){ref-type="fig"}, [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}). For DE binning, genes were determined to be expressed as a function of Fe′ or illumination status by assignment to 8 different response types, characterized by elevated transcript abundance under those conditions. These response types were "low Fe" (20 pM and 40 pM Fe′), "high Fe" (400 pM Fe′), "light" (10 AM, 2 PM, 6 PM), and "dark\" (10 PM, 2 AM, 6 AM) with the four remaining bins being genes assigned to both an Fe or light condition ("dark, low Fe", "dark, high Fe", "light, low Fe", "dark, low Fe"). The two DE approaches utilized (non-parametric Skillings-Mack (SM) test on reads per kilobase of transcript per million (RPKM) values and the LR test in EdgeR, see [methods](#sec011){ref-type="sec"}) returned different total numbers of genes, but several genes (n = 2,389) were identically binned by both tests yielding response type calls we considered to be high confidence ([Table 2](#pgen.1006490.t002){ref-type="table"}, [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}). The main difference in the gene sets returned between the different tests were that the SM test returned more significantly DE genes in the dark and at low Fe, and the LR test returned more significantly DE genes in the light and at high Fe, differences which are likely attributed to variation in the normalization procedures (methods). The percentage of each module assigned to response types by the Skillings-Mack test are shown ([Fig 2A and 2B](#pgen.1006490.g002){ref-type="fig"}).

![Genome wide transcriptome dynamics and transcript level expression of Fe status biomarker genes.\
The transcriptome of semi-continuous cultures of *P*. *tricornutum* grown over diel cycles (12:12) at Fe limited (20 and 40 pM Fe′) and Fe replete (400 pM Fe′) conditions is shown. (**A**) Heat map shows standardized transcript abundance (standardized RPKM across Fe conditions) for each gene. Genes organized into co-expression modules (n = 8418 genes, labeled with module color names) and ordered by hierarchical clustering results. (**B**) Percentage of genes from each expression module belonging to response type groups (light, dark, low Fe, or high Fe determined with Skillings-Mack) indicated. Black bars indicate genes previously identified as Fe-responsive in \[[@pgen.1006490.ref019]\]. Identities of genes and module and response type assignments can be found in [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}. (**C**) Time point averaged transcript abundance (RPKM) is shown, plotted on a log scale, error bars omitted for clarity (**D**) Standardized transcript abundance (calculated independently for each Fe condition) is also plotted to illustrate dark phase expression maxima. Legend includes gene names and Phatr3 PIDs. Background shading indicates dark time points. \*Binned low Fe by Skillings-Mack, †Binned low Fe by EdgeR](pgen.1006490.g002){#pgen.1006490.g002}

10.1371/journal.pgen.1006490.t002

###### Differentially expressed gene totals.

Total number of genes determined to be significantly differentially expressed (p \< 0.05) by the Skillings-Mack test and the LR test in EdgeR (EdgeR) and by both tests. Total number of genes determined to be differentially expressed only by one method or the other (Total Unique) are also shown. The three most abundant unique response types for each test are shown in bold. Ratios of Low Fe and High Fe genes expressed in the dark relative to light are calculated for each test.

![](pgen.1006490.t002){#pgen.1006490.t002g}

                            Total Significant   Total Unique                       
  ------------------------- ------------------- -------------- -------- ---------- ---------
  **Dark only**             1155                2616           1390     **1461**   235
  **Light only**            722                 897            1549     175        **827**
  **Low Fe only**           250                 1240           322      **990**    72
  **High Fe only**          100                 147            440      47         **340**
  **Dark, Low Fe**          122                 732            135      **610**    13
  **Dark, High Fe**         3                   27             31       24         28
  **Light, Low Fe**         18                  99             43       81         25
  **Light, High Fe**        19                  24             299      5          **280**
  ***Total*:**              *2389*              *5782*         *4209*   *3061*     *1447*
  **(\# Dark, Low Fe)**     6.7x                7.4x           3.1x                
  **(\# Light, Low Fe)**                                                           
  **(\# Dark, High Fe)**    0.2x                1.1x           0.1x                
  **(\# Light, High Fe)**                                                          

Though there were differences in total number and response type distribution of DE genes returned with both statistical approaches, there were consistently more "dark, low Fe" genes relative to "light, low Fe" genes: a pattern that did not hold at high Fe ([Table 2](#pgen.1006490.t002){ref-type="table"}). This points to the dark as a particularly important time for transcriptional output in cells grown at low Fe. Genes which had elevated dark phase gene expression at low Fe are likely to be involved in nocturnal metabolic and cellular processes required for growth in these conditions. Known Fe responsive genes such as Fe starvation inducible proteins (ISIPs), flavodoxin, ferric reductases (FRE) were expressed strongly at low Fe and also during the dark phase at all Fe conditions suggesting the promoters of these genes are sensitive to both Fe status and illumination phase ([Fig 2C and 2D](#pgen.1006490.g002){ref-type="fig"}). There were Fe-specific differences in the timing of nocturnal induction. At low and med Fe, elevated transcript levels were detected at 2AM, preceding maximum induction of the same genes at 6AM in Fe-replete cultures, possibly as a response to reduced cellular Fe quotas at Fe-limitation ([Fig 2D](#pgen.1006490.g002){ref-type="fig"}).

The significance of nighttime induction of Fe-responsive genes is unknown. Previously it has been shown that ISIP2a has a role in Fe uptake at the cell surface in *P*. *tricornutum*, and ISIP2a protein levels and Fe-uptake rates are elevated at night \[[@pgen.1006490.ref017]\]. This, combined with diel regulation of several Fe-responsive genes points to the dark phase as an important time for Fe metabolism. Diel patterns in Fe metabolism have been observed in other marine phototrophs. For example, the marine diazotroph *Crocosphaera watsonii* recycles intracellular Fe reserves for synthesis of day or night specific metalloenzymes ("hotbunking") as a strategy to reduce its overall Fe requirement, \[[@pgen.1006490.ref047]\]. Also, *Ostreococcus tauri* expresses ferritin nocturnally in order to store intracellular Fe recovered from damaged/oxidized iron-binding proteins \[[@pgen.1006490.ref018]\]. Ferritin is used for iron storage in some bloom-forming diatoms \[[@pgen.1006490.ref048]\] and may be similarly nocturnally expressed in these species, though no diel pattern was detected for ferritin in our experiments with *P*. *tricornutum*. However, nocturnal expression of other Fe-sensitive genes (like ISIPs) in *P*. *tricornutum* may occur in order to scavenge Fe from damaged proteins (like *O*. *tauri)*, and/or to reuse cellular Fe during diel protein turnover (as in *C*. *watsonii*), even under Fe-replete conditions.

In order to identify other candidate genes involved in specific mechanisms to acquire and maintain cellular Fe and gain additional insight into the significance of nocturnal expression of Fe-sensitive genes, the complete low Fe response type gene set (significant by either SM test or EdgeR, n = 2090) was extracted, clustered hierarchically, and surveyed for gene functions ([S2 Dataset](#pgen.1006490.s002){ref-type="supplementary-material"}, [S8 Fig](#pgen.1006490.s013){ref-type="supplementary-material"}). Improved detection of low Fe sensitive transcripts in this experiment over previous efforts can be attributed to both improvements in transcriptome sequencing coverage, and increased resolution of sampling frequency \[[@pgen.1006490.ref019]\]. This gene set contained several annotated genes not previously known to be Fe sensitive, with putative roles in Fe sequestration and transport, including a highly expressed ZIP Zinc transporter (ZIP1; Phatr3_J52343), and two genes that encode heavy metal-transporting P-type ATPases (Phatr3_EG02342, Phatr3_Jdraft391). One of the most highly induced low Fe genes encoded a protein with a hemerythrin-like domain and a CHY zinc finger domain at the C-terminus (Phatr3_J12097). Hemerythrin domains are found in Fe-binding regulator proteins in both animals and plants, raising the intriguing possibility that this gene has a similar role in *P*. *tricornutum* \[[@pgen.1006490.ref049]--[@pgen.1006490.ref050]\]. Phytochelatin synthase (PCS; Phatr3_J48175) was among the genes highly co-expressed with the ISIPs. Phytochelatins are peptides that play an important role in metal detoxification in plants and fungi and are synthesized by PCS from glutathione and glycine \[[@pgen.1006490.ref051]\]. Nighttime induction of this enzyme in Fe-limited *P*. *tricornutum* could aid in contributing to defense against oxidative stress; directly or through sequestration, storage, and intracellular buffering of Fe, reducing free radical generation from temporal fluctuations in excess free Fe \[[@pgen.1006490.ref052]\]. This would be particularly important for cells adapted to Fe-limitation, which are already heavily impacted by oxidative stress from impaired photosynthetic and mitochondrial electron transport \[[@pgen.1006490.ref004]\].

Other genes found within the highly Fe-responsive set and strongly induced at the transcript level at night and at low Fe were two cathepsins (Phatr3_J25433, Phatr3_J4936) and a v-SNARE protein (Phatr3_J40521/Phatr3_EG02577). Cathepsins are proteases that are active in the low pH of lysosomes, while v-SNARE proteins are involved in intracellular vesicle transport. Vesicle trafficking is a major cellular activity that occurs during normal cellular function to transport proteins between organelles, but it is particularly important during endocytosis, exocytosis, and autophagy, when molecules are transported to the lysosome for degradation. Genes with specific roles in these processes are well-studied in other organisms and homologs of many of these genes were identified in the *P*. *tricornutum* genome. Homologs of genes involved in endocytosis, exocytosis, and autophagy including genes of the exocyst complex, cytosokeleton components, and additional SNARE proteins, were upregulated at low Fe during the dark phase \[[@pgen.1006490.ref053]\] ([S3 Dataset](#pgen.1006490.s003){ref-type="supplementary-material"}). Other transcripts for genes with endocytosis and vesicular sorting functions, such as endosomal coat protein subunits, adaptins, and putative GTPases from the dynamin and Rab family, and members of the ESCRT complex were upregulated in a similar fashion, as were transcripts with established homology to genes with known roles in autophagy. Finally, genes with roles in glycosylation, oligosaccharide processing, and phosphoinositide metabolism, all important in membrane trafficking and sorting, were also upregulated at low Fe ([S3 Dataset](#pgen.1006490.s003){ref-type="supplementary-material"}). Overall, this illustrates that the processes of endocytosis, exocytosis, and autophagy and associated intracellular transport of myriad cellular components occurs mostly during the dark phase in *P*. *tricornutum*, reflecting at least in part a higher degree of lysosomal protein degradation. ƒ trafficking and adaptive significance of these processes for diatoms grown under these conditions. Lommer et al. \[[@pgen.1006490.ref020]\] hypothesized that endocytosis has a role in Fe acquisition, based on the identification of a conserved endocytosis motif in ISIP1 from *Thalassiosia oceanica*. Our results provide additional support for this hypothesis and point to this cellular process as an evolutionarily conserved strategy among diatoms to regulate cellular Fe. Further, our data identify candidate genes specifically involved in that process, such as the cathepsins and v-SNAREs.

There were several transcription factors (TF) identified among the low Fe responsive gene set \[[@pgen.1006490.ref054]\]([S2 Dataset](#pgen.1006490.s002){ref-type="supplementary-material"}). Of the three most highly induced TFs, one was from the C2H2 Zn finger family (CCHH11; Phatr3_J38018), one was from the Myb family (Myb1R_SHAQKYF3; Phatr3_J44256), and another was a sigma factor (sigma70.1a; Phatr3_J14599) predicted to be targeted to the chloroplast, presumably to regulate the chloroplast genome ([S2 Fig](#pgen.1006490.s007){ref-type="supplementary-material"}, Dataset). Several other TFs were also induced at low Fe but to a lesser extent ([S9 Fig](#pgen.1006490.s014){ref-type="supplementary-material"}). We hypothesize at least the most highly upregulated TFs could be involved in regulating downstream cellular responses to low Fe. CCHH11 was the most strongly upregulated gene at low Fe with an average of 453 reads per kilobase of transcript per million mapped (RPKM), 171 RPKM, and 10 RPKM detected at 20 pM, 40 pM, and 400 pM Fe′ respectively ([S2 Fig](#pgen.1006490.s007){ref-type="supplementary-material"}), but lacked clear diel expression patterns. A transmembrane domain was detected at the C-terminus of CCHH11 that did not overlap with the predicted DNA binding domain ([S2 Fig](#pgen.1006490.s007){ref-type="supplementary-material"}). Often, TFs which possess a transmembrane domain are cleaved proteolytically to become liberated and active in regulating gene transcription in a process known as regulated intramembrane proteolysis (Rip; \[[@pgen.1006490.ref055]\]). Dramatic upregulation of transmembrane-containing CCHH11 at low Fe, combined with our observed elevated expression of endocytotic and lysosomal genes at low Fe, could implicate Rip as a previously unknown mechanism for *P*. *tricornutum* to respond and cope with low Fe.

This dataset expands the catalog of putative Fe acquisition and Fe regulatory gene candidates, but the cellular response to Fe-limited growth is not limited to Fe acquisition and homeostasis. The more global cellular response to low Fe, evidenced by the large number and diverse patterns of genes assigned to the low Fe response type, is also important to investigate low Fe adaptation at the genome scale ([S2 Dataset](#pgen.1006490.s002){ref-type="supplementary-material"}).

Transcriptomic hallmarks of reduced growth at low Fe and successional expression of cell cycle regulators {#sec005}
---------------------------------------------------------------------------------------------------------

To characterize Fe and diel-modulated differences in expression, biological roles of genes within co-expression modules and response types were surveyed by calculating over and under-representation of Gene Ontology (GO) annotations ([S2 Table](#pgen.1006490.s024){ref-type="supplementary-material"}). Several modules were disproportionately enriched in genes attributed to various cellular processes, such as growth and division ([S2 Table](#pgen.1006490.s024){ref-type="supplementary-material"}). Core histone transcripts (H2A/B, H2, H4), main nucleosome components, were strikingly abundant in the tan module ([S2 Table](#pgen.1006490.s024){ref-type="supplementary-material"}, [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}). In *Drosophila*, expression of histone mRNAs is tightly constrained to a brief pulse during S-phase \[[@pgen.1006490.ref056]\], and a similar correlation has been observed in the diatom *Thalassiosira pseudonana* \[[@pgen.1006490.ref057]\]. Manual inspection of annotated genes in the tan module revealed that several genes were functionally associated with S-phase, indicating that peak timing of expression in this module may correspond to a time in the diel cycle when the majority of the *P*. *tricornutum* population is progressing through S-phase ([S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}). Key transcripts of mitosis-related genes were found in the darkgreen module ([S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}), including the cell division cycle protein 20 (CDC20, Phatr3_J12783): a highly conserved key activator of the anaphase promoting complex specifically expressed during the G2-to-M-phase in *P*. *tricornutum* \[[@pgen.1006490.ref058]\]. Peak expression of the darkgreen module lagged shortly behind S-phase expression (selected tan module genes, [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}), as would be expected as the population progresses through the cell cycle ([Fig 3](#pgen.1006490.g003){ref-type="fig"}). Predominantly, these S-phase and M-phase genes were expressed during the light to dark transition correlating with the timeframe when cells were dividing ([Fig 1B](#pgen.1006490.g001){ref-type="fig"} and [Fig 3](#pgen.1006490.g003){ref-type="fig"}). Fe-replete cultures expressed cell cycle phase-specific genes earlier and over a shorter period of time (2PM, 12hr) than what was detected in Fe-limited cultures (6PM, 16hr) ([Fig 3A](#pgen.1006490.g003){ref-type="fig"}). This indicates that at high Fe, the population is both able to begin dividing sooner and is more highly synchronized than cells grown at low Fe. We hypothesize that this delay in timing of cell cycle progression in low Fe cultures, inferred from mRNA levels of phase marker genes, is a consequence of reduced energetic status in Fe-impoverished and photosynthetically compromised cells.

![Expression of cell cycle progression-related genes, ribosome biogenesis, and cyclins and cyclin-dependent kinases (CDKs).\
(**A**) Data show average standardized transcript abundance (across Fe conditions) in WGCNA co-expression modules enriched for gene functions associated with cell cycle progression (tan--selected genes see [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}, darkgreen) and ribosome biogenesis (darkgrey) genes at 20, 40, and 400 pM Fe′ with shaded area indicating one standard deviation ([S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"} for gene IDs). (**B**) Cyclin and CDK transcript abundance. Waves of co-expressed cyclins and CDKs were identified with hierarchical clustering ([S10 Fig](#pgen.1006490.s015){ref-type="supplementary-material"}). Standardized cluster averages (within each Fe condition, positive values only) are shown. (**C**) Transcript abundance (RPKM) of wave 1 cyclins and dsCYC2-regulating transcription factors (bZIP10 and Aureo1a) are shown. Dashed lines indicate data are plotted on the secondary y-axis.](pgen.1006490.g003){#pgen.1006490.g003}

Cell growth (leading to division) depends on protein synthesis by ribosomes, and reductions in ribosome assembly and translation are often among the first and most drastic responses when cells experience environmental stress \[[@pgen.1006490.ref059],[@pgen.1006490.ref060]\]. In our experiments, two modules (turquoise, darkgrey) with similar overall patterns of expression were enriched in genes required for ribosome biogenesis ([S2 Table](#pgen.1006490.s024){ref-type="supplementary-material"}). The darkgrey module was nearly completely comprised of transcripts for almost all of the ribosomal proteins (RPs) encoded in the *P*. *tricornutum* genome ([S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}). Peak expression of RPs occurred largely during the illuminated phase of the diel cycle, and in an inverse pattern to expression of cell cycle progression genes, illustrating that the cell segregates expression of genes essential for growth and division temporally ([Fig 3](#pgen.1006490.g003){ref-type="fig"}). RPs are some of the most highly expressed mRNAs in the *P*. *tricornutum* genome, and though they are highly expressed overall, they were less highly expressed at low Fe ([S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}). Reduction of ribosome biogenesis, and subsequent decreased protein biosynthesis could explain reduced growth rates at low Fe.

Progression through the cell cycle is regulated by cyclins and cyclin-dependent kinases (CDKs). Cyclins and CDKs are evolutionarily conserved in eukaryotes, and cyclins have expanded in diatoms, presumably to regulate cell cycle progression in response to varied environmental signals \[[@pgen.1006490.ref061]--[@pgen.1006490.ref062]\]. Over diel cycles, four distinct waves in expression of cyclins and CDKs were identified with expression of one wave of these cell cycle regulators limited to the day (wave 1) while most other cyclins and CDKs (waves 2, 3, 4) were more highly expressed during the light-dark transition period. Cyclin/CDK wave 1 contained the G1 marker cyclin CYCP6 and the diatom specific cyclin dsCYC2 known to be involved in regulating the light-dependent onset of G1 phase in *P*. *tricornutum*. This finding, combined with the observation of decreased proportion of dividing cells detected during the light indicate that the majority of the population is in the G1 phase at this time. Dark phase expression of the other waves of cyclins/CDKs is consistent with observations of increased cell division and the conclusion that this time period is particularly active for cell cycle progression ([Fig 3B](#pgen.1006490.g003){ref-type="fig"}, [S10 Fig](#pgen.1006490.s015){ref-type="supplementary-material"}). Like the S- and M-phase modules, the window of peak cyclin/CDK expression (waves 2--4) was shorter in Fe-replete cultures than in Fe-limited cultures, corroborating a higher degree of synchrony in Fe-replete populations. The timing of cyclin/CDK waves 2 and 3 correlated with the S-phase and M-phase enriched modules respectively, leading us to investigate whether these cyclin/CDK waves were cell cycle stage or transition specific. Indeed, cyclin/CDK wave 2 contained several regulators previously shown to be expressed at the G1-S boundary in *P*. *tricornutum* such as CYCL, dsCYC6, hCDK2, and hCDK5 and cyclin/CDK wave 3 contained marker genes for the G2-M transition such as CDKA2 and CYCB1 \[[@pgen.1006490.ref062]--[@pgen.1006490.ref063]\]. However, each wave also contained genes previously assigned to the opposite phase transition ([S10 Fig](#pgen.1006490.s015){ref-type="supplementary-material"}, \[[@pgen.1006490.ref062]\]). In the previous study, Huysman et al. \[[@pgen.1006490.ref062]\] obtained a much higher degree of synchrony with prolonged dark exposure, so the observed discrepancy is most likely attributed to heterogeneity in our mildly synchronized populations. However, the temporal segregation in G1-S and G2-M cyclin expression we observed suggests their transcript level expression dynamics are more nuanced than has been demonstrated previously, illustrating our limited understanding of cell cycle regulation, and underscoring the need for an increased number of observations and experiments to continue to elucidate cyclin/CDK specificity and roles in diatom cell cycle.

Diatom cyclin and CDK expression is affected by nutrient supply, with previous work indicating that dsCYC10 is sensitive to phosphate availability \[[@pgen.1006490.ref062]\]. In this study, dsCYC10 (Phatr3_50251) was also up-regulated at low Fe (with high confidence), supporting the idea that it may sense or relay the status of multiple nutrients during the cell cycle \[[@pgen.1006490.ref062]\]. Other cell cycle regulators were also sensitive to low Fe, however few dramatic Fe-specific modulations in transcript abundance were detected ([S10 Fig](#pgen.1006490.s015){ref-type="supplementary-material"}). At low Fe, these genes may be modestly upregulated in order to drive cell cycle progression forward in particularly energetically impaired cells.

Transitions in illumination state are strong signals that are associated with transcriptome state shifts in diatoms \[[@pgen.1006490.ref028]\]. In *P*. *tricornutum*, light is a driver of cell cycle progression through the expression of light sensitive cyclins such as dsCYC2 (cyclin/CDK wave 1), regulated transcriptionally by AUREOCHROME1a (AUREO1a, Phatr3_J8113) and bZIP10 (Phatr3_J43744) transcription factors \[[@pgen.1006490.ref064]\]. Strong temporal dynamics in the transcript abundance of AUREO1a and bZIP10 were observed under diel cycles ([Fig 3B](#pgen.1006490.g003){ref-type="fig"}). Transcripts of AUREO1a steadily accumulate throughout the diel cycle, reaching maximum abundance prior to the onset of illumination indicating that the dark-light transition represses transcription of AUREO1a. In contrast, bZIP10 is at maximum transcript abundance directly following the dark-light transition, and is minimally expressed during the dark. Opposing diel expression patterns by the two TFs that regulate dsCYC2 production, could indicate a functionally distinct role to coordinate cell division with both the light-dark and dark-light transitions to ensure cell cycle progression is synchronized fully to light cycles, and not just the onset of illumination. Despite the detectible diel cycling of dsCYC2, and the previously described role of this gene in cell cycling, we noted that overall mRNA levels of this gene were extremely low (\<5 RPKM) relative to other similarly expressed cyclins ([Fig 3B](#pgen.1006490.g003){ref-type="fig"}). This suggests that cyclins other than dsCYC2, such as CYCP5 and CYCP6, which are drastically modulated in our experiments, also have important roles in linking cell cycle progression to the dark-light transition. Though *P*. *tricornutum* does not synchronize tightly to diel cycles, there is a subtle natural synchrony evidenced by the dark phase increase in proportion of dividing cells, and peaks in expression of the majority of cell cycle regulators (including cyclins and CDKs). The degree of synchronization is dependent of the Fe status of the population, with rapidly growing (Fe-replete) populations not only dividing earlier, but over a shorter period of time in the diel cycle than Fe-limited populations. Overall, these results provide additional details and mechanistic insight into how nutrients and light drive populations to become synchronous with the light cycle \[[@pgen.1006490.ref065]--[@pgen.1006490.ref066]\].

It is possible that there are biogeochemical consequences to this slowed growth at low Fe. Generally, slower growth results in a reduced rate of progression through the cell cycle, and as a consequence cells spend a longer period of time in M-phase. This is biogeochemically and ecologically significant, since M phase is the primary period of silica deposition \[[@pgen.1006490.ref067]\]. Some diatoms grown in Fe-limiting conditions with sufficient silicic acid availability (e.g., HNLC regions) have more heavily silicified cell walls and are therefore more efficient at long term carbon sequestration from the surface ocean \[[@pgen.1006490.ref012], [@pgen.1006490.ref068]--[@pgen.1006490.ref072]\]. Although *P*. *tricornutum* is a lightly silicified diatom, the Fe-specific and diel regulation of these (likely cell cycle driven) processes are likely to be similar within more heavily silicified diatoms.

Regulation of the light harvesting complex and partial downregulation at low Fe {#sec006}
-------------------------------------------------------------------------------

Up to 80% of the cellular Fe requirement is associated with photosynthesis in photoautotrophic cells, where Fe is bound to various cytochromes, and Fe-S cluster proteins of the photosynthetic electron transport chain \[[@pgen.1006490.ref026]\]. Fe depletion in phototrophs impairs photosynthetic efficiencies, leading to increased production of damaging reactive oxygen species (ROS) when photosynthetic light harvesting exceeds the capacity of photosynthetic electron transport \[[@pgen.1006490.ref073]\]. Reorganization of the light harvesting machinery is expected be an important aspect of acclimation to shifting Fe and light conditions. To better define how this is accomplished, the temporal expression dynamics of key light harvesting genes were investigated.

Transcripts encoding for both pigment biosynthesis and light harvesting complex (LHC) proteins were found nearly exclusively in the similarly expressed cyan (n = 199) and darkorange (n = 81) modules, which peak in the late afternoon (2PM--6PM; [Fig 2](#pgen.1006490.g002){ref-type="fig"}). The cyan module contained the complete pathway for porphyrin and chlorophyll metabolism, and several steps in terpenoid backbone biosynthesis, including genes such as polyprenyl synthetase (Phatr3_J19000), geranylgeranyl reductase (Phatr3_J31683), and chlorophyll synthase (Phatr3_J12807). Together, the cyan and darkorange modules contained nearly all annotated antenna proteins of the LHC superfamily which bind pigment and transfer absorbed light energy to the reaction center complexes. These modules also contained transcripts for nuclear-encoded components of photosystem II (PSII) such as PsbM (Phatr3_J55057), PsbO (Phatr3_J20331), and PetJ (Phatr3_J44056). A large proportion of the darkorange module (LHC transcripts), was downregulated at low Fe, consistent with findings from previous transcriptome studies of Fe-limited phytoplankton \[[@pgen.1006490.ref019]--[@pgen.1006490.ref020]\]. In contrast, several genes contained within the cyan module (pigment biosynthesis genes) were upregulated in Fe limiting conditions. The overall downregulation of the LHCs and upregulation of pigment (and specifically carotenoid biosynthesis genes) is consistent with a transcriptional effort to remodel the LHC, towards increased photoprotection, for example by production of DLHCs; expected in Fe-limited, macronutrient replete conditions \[[@pgen.1006490.ref004]\].

The expression patterns of the full suite of LHC proteins identified in our dataset (n = 48) were analyzed to obtain a complete view of how these genes are modulated over diel cycles and in response to Fe. Phylogenetically, these proteins belong to 4 distinct groups: the chromophyte-specific fucoxantin Chl *a/c* LHCFs, the red algal-like LHCRs, LHCZ proteins and LHCX proteins \[[@pgen.1006490.ref074]\]. Hierarchical clustering of the LHC transcript abundance patterns revealed that the majority are co-expressed (mirroring results obtained for the WGCNA module analysis) and belong to the LHCF/LHCR subtype ([Fig 4](#pgen.1006490.g004){ref-type="fig"}, [S11 Fig](#pgen.1006490.s016){ref-type="supplementary-material"}). Strong temporal dynamics were conserved across Fe conditions in this group. LHC transcript abundance was minimal directly following illumination (10AM) in all Fe conditions with peak mRNA detected in the late afternoon ([Fig 4](#pgen.1006490.g004){ref-type="fig"}). LHC transcripts were also downregulated in *P*. *tricornutum* when illuminated following a 48hr dark period \[[@pgen.1006490.ref075]\]. This suggests that cellular signals generated upon illumination are involved in repressing LHC transcription. It has been shown that this could be effected by PtCPF1, a cryptochrome/photolyase \[[@pgen.1006490.ref076]\] which is maximally expressed in our experiments at 10AM. LHC expression is induced in the late afternoon, directly preceding the detectible onset of transcripts enriched for cell cycle progression in all Fe conditions ([Fig 4](#pgen.1006490.g004){ref-type="fig"}). Expression of light harvesting genes has been coupled to chloroplast ontogeny at the G2/M transition in the diatom *Seminavis robusta*, and at S-phase in both *Thalassiosira pseudonana* and *P*. *tricornutum* \[[@pgen.1006490.ref057],[@pgen.1006490.ref062],[@pgen.1006490.ref077]\]. Taken together our data further substantiates that transcriptional upregulation of the LHCs occurs at least in part as preparation for chloroplast replication during cell cycle division.

![Light harvesting antenna protein expression.\
(**A**) Heat map shows min/max transcript abundance (RPKM across Fe conditions, plotted on a log scale). Dendogram shows hierarchical clustering (Pearson correlation) of RPKM values with corresponding module assignment from the WGCNA clustering analysis. Expression clades are labeled to indicate Fe-response types (high Fe, low Fe). Dominant protein family composition of expression clades (Lhcf/Lhcr and Lhcx/Lhcz) is also indicated. Gene identifiers (Phatr3 protein IDs) can be found in [S3 Fig](#pgen.1006490.s008){ref-type="supplementary-material"}. (**B**) Selected expression patterns of Fe-sensitive light harvesting complex antenna proteins from the Lhcx, ac2, and Lhcz families. Lhcz (Phatr3_J16481) is from the red module and is indicated on the heat map with a red arrow.](pgen.1006490.g004){#pgen.1006490.g004}

Genes of the photoprotective LHCX and LHCZ families had a different pattern of expression, with day-time maxima; the majority of which do not have any Fe-modulated expression ([Fig 4](#pgen.1006490.g004){ref-type="fig"}). However, three LHC genes, from the LHCF, LHCZ, and LHCX families were significantly upregulated under low Fe ([Fig 4B](#pgen.1006490.g004){ref-type="fig"}). Transcripts encoding two of these proteins (Lhcx2; Phatr2_54065, and Lhcz; Phatr3_16481) were previously shown to be Fe responsive \[[@pgen.1006490.ref019]\]. LHCs, and especially members from the stress-response LHCX family, are thought to play an essential role in compensating for decreased levels of PSI units in Fe-deplete cells and to enhance non-photochemical quenching (NPQ) both under light and Fe stress \[[@pgen.1006490.ref078]--[@pgen.1006490.ref079]\]. A photoprotective function for LHCX1 (Phatr3_J27278) has been confirmed in *P*. *tricornutum* \[[@pgen.1006490.ref079]--[@pgen.1006490.ref080]\].

In addition to LHC remodeling, there are other changes to the light harvesting antenna that may have photoadaptive significance in Fe limiting conditions. For example, one of the most highly day-time induced transcripts at low Fe encoded a chloroplast-targeted heme oxygenase (Phatr3_J5902). Heme oxygenase catalyzes the rate limiting step in degradation of heme to bilirubin, producing ferrous Fe. In *Arabidopsis thaliana*, heme oxygenase is required to synthesize phytochrome photoreceptors \[[@pgen.1006490.ref081]\]. It is possible that *P*. *tricornutum* heme oxygenase is involved in plastid Fe homeostasis, and may be upregulated to increase supply of ferrous Fe under these conditions, perhaps through recovery of Fe from damaged porphyrins. Further, membrane lipid desaturation in cyanobacteria lead to reductions in photoinhibition of PSII \[[@pgen.1006490.ref082]--[@pgen.1006490.ref084]\] and three genes encoding fatty acid desaturases showed a strong light-sensitive induction at low Fe in *P*. *tricornutum*, potentially as a strategy to protect the photosystem in a similar manner (Phatr3_J41570; Phatr3_J55137; Phatr3_EG02619).

Diel oscillations in transcripts of primary metabolism {#sec007}
------------------------------------------------------

When illuminated, phototrophic organisms acquire light energy that is then stored within macromolecular chemical bonds. In the dark, stored chemical energy is then catabolized to meet cellular demands. To obtain a comprehensive view of how *P*. *tricornutum* regulates metabolic shifts between photoautotrophic (anabolic) and autoheterotrophic (catabolic) states, the diel expression of genes involved in several major primary metabolic pathways for acquiring, partitioning, storing, and catabolizing chemical energy was evaluated. A filtered set of 1,187 genes with either hypothesized or known direct or intersecting roles in central carbon metabolism and nitrogen assimilation were identified ([S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"}) and clustered hierarchically (methods). Strong, coordinated diel regulation of transcript abundance was observed for several essential genes in carbon metabolism. Broadly, pathways that were upregulated during the day included the Calvin Benson (CB) cycle, gluconeogenesis, carbohydrate biosynthesis, amino acid biosynthesis, mitochondrial glycolysis, fatty acid biosynthesis, and nitrogen assimilation. Nighttime expressed pathways included cytosolic and chloroplast glycolysis, carbohydrate degradation, the pentose phosphate pathway, fatty acid degradation (β-oxidation), and the TCA cycle.

Genes in the CB cycle were among the most highly expressed metabolism-related transcripts in *P*. *tricornutum*, consistent with the central importance of carbon fixation for a photoautotroph ([S12 Fig](#pgen.1006490.s017){ref-type="supplementary-material"}). Most CB transcripts reached peak expression at 10AM across Fe conditions, consistent with previous observations in *P*. *tricornutum* ([Fig 5A](#pgen.1006490.g005){ref-type="fig"}) \[[@pgen.1006490.ref085]\]. However, not all chloroplast-localized isoforms of putative CB enzymes peaked in expression during the light phase. Fructose bisphosphatase (FBP; Phatr3_J2793), phosphoribulokinase (PRK; Phatr3_EG02325), and a fructose/seduheptulose bisphosphatase (FBP/SBP; Phatr3_EG02409) were maximally expressed during the dark phase. Nocturnal expression indicates that these isozymes catalyze anapleurotic reactions, and are required during conditions in which the chloroplast is being resupplied with carbon skeletons in preparation for daytime carbon fixation. Co-expression of these genes and several putative chrysolaminarin genes with two triosephosphate transporters (TPT2; Phatr3_J46451, TPT4b; Phatr3_J24610), suggests that these transporters could mediate plastid carbon import under these conditions \[[@pgen.1006490.ref086]\] ([S13 Fig](#pgen.1006490.s018){ref-type="supplementary-material"}).

![Coordinated expression of central carbon metabolic pathways.\
(**A**) Biological replicate averaged transcript abundance (for each time point) shown for genes from the Calvin Benson cycle, mitochondrial glycolysis pathway, cytosolic glycolysis pathway, and TCA cycle for each Fe condition (error bars omitted for clarity). Series labeled with gene names (see [S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"} for abbreviations) and Phatr3 PID. (**B**) Schematic representation of reactions catalyzed by GAPDH and PGK across the cytosol, mitochondria, and chloroplast (shown with 4 membranes). NADPH and ATP produced from light harvesting depicted in green, and NAD(P)H and ATP regenerated in the cytosol or mitochondria depicted in magenta. Yellow/blue highlighted arrows indicate day/night upregulation respectively. Solid arrows indicate likely transport processes, dashed arrow indicates hypothesized carbon skeleton transport. GAP: glyceraldehyde 3-phosphate, 3PG: 3-phosphoglycerate, FAB: fatty acid biosynthesis. (**C**) Transcript abundance for chloroplast, mitochondrial, and cytosolic targeted isoforms of glyceraldehyde phosphate dehydrogenase (GAPDH) and phosphoglycerate kinase (PGK) are plotted (TPI-GAPDH: triosephosphate isomerase-GAPDH fusion protein).](pgen.1006490.g005){#pgen.1006490.g005}

Carbohydrate biosynthesis genes (UGP-PGM; Phatr3_EG02613, SKN1; Phatr3_J50238, UGP; Phatr3_J23639, BGS1; Phatr3_J9617) were upregulated during the light when synthesis of the diatom storage carbohydrate chrysolaminarin occurs ([S13B Fig](#pgen.1006490.s018){ref-type="supplementary-material"}) \[[@pgen.1006490.ref087]\]. Conversely, putative chrysolaminarin and TCA cycle genes were upregulated at night reflecting elevated carbon catabolism ([Fig 5A](#pgen.1006490.g005){ref-type="fig"}, [S13B Fig](#pgen.1006490.s018){ref-type="supplementary-material"}). Photosynthetically fixed carbon is partitioned toward carbohydrates or the TCA cycle by glycolysis and gluconeogenesis, which are the central energetic and carbon metabolic pathways. Since many of the glycolysis/gluconeogenesis reactions are reversible, they are incompatible metabolic processes and in order for them to occur simultaneously they must be segregated either temporally (as in *Chlamydomonas)* or by intracellular compartmentation \[[@pgen.1006490.ref030],[@pgen.1006490.ref088]\]. Diatoms possess several glycolysis/gluconeogenesis pathways localized to the chloroplast, cytosol, and notably (for stramenopiles) to the mitochondria \[[@pgen.1006490.ref089]--[@pgen.1006490.ref091]\]. The unique compartmentation of glycolysis in the mitochondria of diatoms (and other stramenopiles) has been identified previously, but the purpose of this organization is unknown; it does not exist in land plants or other eukaryotes \[[@pgen.1006490.ref090]--[@pgen.1006490.ref093]\]. Since this novel organization is characteristic of diatoms, diel expression of transcripts encoding these glycolytic/gluconeogenic enzymes was analyzed in order gain new insight into its significance.

Diatoms possess several isozymes for many steps in the cytosolic glycolysis/gluconeogenesis pathway and it remains unclear whether they are simply functionally redundant (i.e. extra copies) or whether they might be differentiated \[[@pgen.1006490.ref090]--[@pgen.1006490.ref091]\]. Some isozymes were more highly expressed during the day, while others were expressed at night indicating there could be more specialization than previously known. Specifically, two isozymes of GPI and FBA (GPI, Phatr3_EG02208 and FBA3 II, Phatr3_J29014) are co-expressed with the committed gluconeogenesis enzyme FBP (Phatr3_J23247). Selected carbon transporters (TPT1; Phatr3_J50742, TPT4a; Phatr3_J54017) were also diurnally co-expressed. TPT4a has been experimentally localized to the chloroplast inner membrane where it is hypothesized to play an important role in export or import of carbon to the cytosol \[[@pgen.1006490.ref086]\]. Diurnal expression of these isozymes and their co-expression with carbohydrate biosynthesis is evidence that they may be specialized for plastid carbon export and gluconeogenesis. Other cytosolic glycolysis isozymes were most highly expressed in the dark phase and co-expressed with the TCA cycle, which is consistent with a role for this pathway in supplying pyruvate from catabolized carbohydrate for nocturnal respiration ([Fig 5A](#pgen.1006490.g005){ref-type="fig"}). In contrast to cytosolic glycolysis, mitochondrial glycolysis isozymes were expressed maximally during the light phase and were coordinated with the CB cycle indicating that operation of this pathway may be elevated in the day, perhaps to directly process triose phosphates exported from the chloroplast ([Fig 5](#pgen.1006490.g005){ref-type="fig"}). Though glycolysis/gluconeogenesis pathways were already known to be segregated spatially (compartmentalized) in diatoms, temporal segregation of the expression of these different glycolytic pathways into phases of the diel cycle is strong evidence that they are more dramatically functionally differentiated than previously understood.

An advantage of localizing glycolysis to the mitochondria and operating this pathway in the day is that diatoms can simultaneously synthesize carbohydrates (via cytosolic gluconeogenesis) and direct carbon towards the TCA cycle as pyruvate (mitochondrial glycolysis), processes which are otherwise incompatible if localized to the same compartment. While the actual fate of the pyruvate produced in mitochondrial glycolysis is unclear, it could either be exported from the mitochondria for biosynthetic processes in other compartments, or could enter the TCA cycle. The TCA cycle is both a degradative and biosynthetic, since its intermediates are drawn away as precursors for a wide variety of molecules such as amino acids. Based on the daytime expression of mitochondrial glycolysis and co-expression with other anabolic pathways, we hypothesize that this pathway may be biosynthetic and utilized to sustain daytime synthesis of amino acids and other important cellular molecules. A key consequence of mitochondrial glycolysis is that ATP and NAD(P)H fixed into GAP during the CB cycle are regenerated directly in the mitochondria, effectively shuttling ATP and NAD(P)H directly from the chloroplast to supply daytime mitochondrial energy demands ([Fig 5B](#pgen.1006490.g005){ref-type="fig"}). Energetic coupling between diatom mitochondria and chloroplasts has been demonstrated, and *in silico* metabolic construction in *P*. *tricornutum* predicted that the nature of this connection is the mitochondrial glycolysis pathway \[[@pgen.1006490.ref093], [@pgen.1006490.ref094]\]. Our data validate a strong functional connection between CB and mitochondrial glycolysis, thus clarifying the metabolic significance of the latter pathway.

Few genes in these primary metabolic pathways were strongly Fe sensitive as defined by assignment to either the low Fe or high Fe response type (methods). Only one CB gene, fructose bisphosphate aldolase (FBAC5; Phatr3_J41423), previously identified to be highly sensitive to low Fe was upregulated in Fe limited cultures with high confidence \[[@pgen.1006490.ref019],[@pgen.1006490.ref095]\]. However, some primary metabolism genes were more highly expressed at high Fe than in the more energy impoverished low Fe acclimated cultures. The CB cycle Ribose 5-P isomerase (RPI; Phatr3_J13382) and mitochondrial enolase (ENO; Phatr3_Jdraft1192) were both upregulated in high Fe-adapted *P*. *tricornutum* with high confidence. Other CB genes had elevated maximum transcript abundance in high Fe relative to low Fe at specific time points. For example, CB cycle transcripts of isoforms of GAPDH (Phatr3_J22122) and PGK (Phatr3_J29157) were significantly elevated at 2AM and 6AM in high Fe cultures ([Fig 5C](#pgen.1006490.g005){ref-type="fig"}). Upregulation under high Fe reflects energetically replete carbon fixation conditions in *P*. *tricornutum* with unimpaired photosynthetic electron transport. Increased mRNA for the mitochondrial glycolysis genes TPI-GAPDH and PGK was observed in high Fe conditions at specific time points, similar to observations for the CB cycle ([Fig 5C](#pgen.1006490.g005){ref-type="fig"}). At low Fe, *P*. *tricornutum* decreases expression of these genes, presumably to reduce the concentration of enzymes that partition photosynthetic energy to this pathway thus matching the reduced supply from impaired photosynthetic electron transport. Overall, this interpretation is consistent with physiological data showing that in diatoms with reduced growth rates, the proportion of the transient carbon pool belonging to cellular reductants is also reduced \[[@pgen.1006490.ref096]\]. Cytosolic GAPDH and PGK transcripts were either elevated at or unaffected by low Fe ([Fig 5C](#pgen.1006490.g005){ref-type="fig"}). We interpret the modest upregulation of cytosolic glycolysis in Fe limited cells as an attempt to increase cellular energy levels by drawing more carbon through this pathway to sustain production of ATP and NAD(P)H in the cytosol, and supply the TCA cycle and oxidative phosphorylation at night.

This pathway-level analysis of carbon metabolism across Fe conditions clarifies the functional differentiation between glycolytic pathways targeted to different compartments in *P*. *tricornutum*, and reveals that transcript abundance is modulated in a manner reflecting energy status at the different Fe conditions. At high Fe, when cells are energetically replete, there is an increased expression of genes involved in generating chemical energy, through carbon acquisition and storage in the light phase. These are some of the only genes in the *P*. *tricornutum* genome that are upregulated at high Fe along with genes of the light harvesting antenna and ribosomal proteins. At low Fe, *P*. *tricornutum* increases the expression of genes involved in acquiring energy from intracellular carbon (cytosolic glycolysis during the dark phase).

Dark phase proteolysis and catabolism in *P*. *tricornutum* {#sec008}
-----------------------------------------------------------

The night appears to be a particularly active phase for cellular processes and metabolism, especially in low Fe cultures ([Table 2](#pgen.1006490.t002){ref-type="table"}). As discussed previously, dark expression of Fe acquisition genes, vesicle transport and autophagy genes, and partial synchronization of division with the dark phase means that a significant proportion of the differentially expressed transcriptome is involved in these processes. Additionally, the dark phase is particularly important for autotrophs since energy demands must be met by catabolizing macromolecules stored during the daytime. Reduced cellular nutritional quotas for *P*. *tricornutum* cultivated at low Fe would likely affect the magnitude of nighttime catabolic processes. To explore this, we analyzed diel fluctuations in the metabolite content and composition of *P*. *tricornutum* along with the functional roles of metabolic genes active in the nighttime transcriptome, and particularly those that were upregulated at the transcript level under low Fe.

Diel fluctuations were documented for several metabolites (n = 33 total detected) with most reaching their maxima during the daytime ([Fig 6](#pgen.1006490.g006){ref-type="fig"}). Generally, metabolites were more abundant in cells grown at 400 pM Fe′ than in cells grown at lower Fe concentrations, with only threonate elevated significantly at low Fe ([Fig 6](#pgen.1006490.g006){ref-type="fig"}). Threonate is a product of ascorbate degradation, indicating that Fe limited *P*. *tricornutum* has higher ascorbate levels and may rely on its anti-oxidative properties to cope with elevated reactive oxygen stress \[[@pgen.1006490.ref097]--[@pgen.1006490.ref099]\]. Ascorbic acid likely functions in preventing photo-oxidative damage in Fe-deplete cells. On the other hand, if transported to the outer membrane, nighttime elevation of ascorbic acid can reflect an important role in the reduction of ferric to ferrous Fe during Fe-uptake; which has shown to be primarily segregated to the night phase in other phytoplankton and microbes \[[@pgen.1006490.ref100]--[@pgen.1006490.ref102]\].

![Fluctuations in metabolite composition of *P*. *tricornutum*.\
Abundant and/or condition-modulated metabolites are shown. Heat map colors scale to metabolite relative minimum and maximum across Fe conditions. Data are clustered hierarchically (Pearson correlation). Abundant metabolites (upper 50th percentile of average relative abundance) are indicated with black bars. Metabolites significantly more abundant at high Fe/low Fe are indicated (High Fe: 400pM Fe′, Low Fe: 20 & 40 pM Fe′ combined). Metabolites significantly more abundant during the light/dark phases are labeled (Light or Dark).](pgen.1006490.g006){#pgen.1006490.g006}

Daytime metabolite composition was dominated by amino acids, but also included glycerate (and glycerate 3-P) formed in the first step of carbon fixation. Detection of these metabolites during the day is in agreement with general upregulation of essential transcripts for amino acid biosynthesis ([S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"}). At night, the metabolite composition of *P*. *tricornutum* shifts. The relative abundance of sugars (maltose, galactose) increases and is correlated with increases in mRNA for chrysolaminarin degradation genes, suggesting the observed accumulation is a result of increased catabolism of storage carbon. Urea and spermidine are also abundant at night, as are lysine and the branched chain amino acid isoleucine.

The global transcriptome was examined to identify metabolic drivers of these diel and Fe-specific metabolite shifts in *P*. *tricornutum*. Co-expression modules (lightcyan, royal blue) were identified which were strongly enriched with genes involved in proteolysis ([S2 Table](#pgen.1006490.s024){ref-type="supplementary-material"}). Proteolysis may remediate the damage caused by elevated levels of ROS generated by an Fe-impoverished electron-transport chain \[[@pgen.1006490.ref019]\]. Additionally, proteolysis could be upregulated in order to recycle fixed carbon and nitrogen to decrease dependence on extracellular input and the subsequent demand for additional reductant \[[@pgen.1006490.ref019]--[@pgen.1006490.ref020],[@pgen.1006490.ref103]\]. However, induction of genes in the lightcyan and royal blue modules was dark phase specific and conserved across Fe conditions, suggesting an Fe-independent function of this gene repertoire in the remodeling of a day-time to a night-time proteome, which could further serve as a dark phase source of amino acids.

In addition to proteome remodeling, dark phase regulation of metabolic pathways likely affected metabolite composition. Several catabolic pathways are upregulated at night in *P*. *tricornutum* with varying degrees of Fe sensitivity. Fatty acid breakdown occurs through β-oxidation in the peroxisome and mitochondria to produce acetyl-CoA and propionyl-CoA that fuels nighttime respiration. Peroxisomal production of acetyl-coA is processed through the glyoxylate cycle prior to import into to the mitochondria as malate. Both β-oxidation and the glyoxylate cycle were strongly induced at night, and did not display sensitivity to Fe status ([Fig 7](#pgen.1006490.g007){ref-type="fig"}). However, some transcripts encoding dark catabolism genes were highly upregulated in low Fe conditions, such as both isozymes of methylcrotonyl-CoA carboxylase (MCC1; Phatr3_J843, MCC2; Phatr3_J19329), and branched chain alpha-keto acid dehydrogenase E1, beta subunit (BCKDC, Phatr3_J11021), which are important for degradation of branched chain amino acids (BCAA; [Fig 7](#pgen.1006490.g007){ref-type="fig"}). Recently, MCC1 has been shown to have a role in lipid metabolism \[[@pgen.1006490.ref104]\]. Seemingly, *P*. *tricornutum* increases the capacity to catabolize a different pool of metabolites (BCAA) under low Fe compared to Fe replete conditions. Upregulation of BCAA catabolism in low Fe-adapted *P*. *tricornutum* could be in response to elevated influx of BCAA from the enhanced proteolysis occurring under these conditions. Beyond this, it is possible that the upregulation of BCAA catabolism is a response to having reduced carbohydrate and lipid stores under Fe limitation, forcing cells to catabolize less labile sources of carbon skeletons and chemical energy. Interestingly, acetolactate synthase (ALS; Phatr3_J37341), the first step in BCAA synthesis, was upregulated under low Fe, though it was maximally expressed during the daytime/anabolic phase. Apparently, under low Fe, *P*. *tricornutum* not only increases degradation of BCAA, but also increases the synthesis of these molecules, indicating there is a regulated shift in amino acid composition with some adaptive significance under low Fe conditions. The advantage of shifting towards enhanced BCAA metabolism in response low Fe remains unclear.

![Dark phase expression of catabolism and pyruvate hub enzymes.\
(**A**) Biological replicate-averaged transcript abundance (for each time point) shown for genes from the β-oxidation (mitochondrial and peroxisomal combined), glyoxylate cycle, branched chain amino acid (BCAA) degradation pathways and mitochondrial pyruvate hub (error bars omitted for clarity). Series labeled with gene names (see [S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"} for abbreviations) and Phatr3 PID. (**B**) Pathway diagram showing reactions and co-reactants (H^+^ omitted) of the pyruvate hub. Each arrow represents a unique isozyme. PEP = phosphoenolpyruvate, OAA = oxaloacetate, PK = pyruvate kinase, ME = malic enzyme, MDH = malate dehydrogenase, PEPCK = phosphoenolpyruvate carboxykinase, PEPC = phosphoenolpyruvate carboxylase.](pgen.1006490.g007){#pgen.1006490.g007}

Carbon skeletons catabolically generated in the dark phase enter mitochondrial metabolism as a variety of 3C/4C molecules that intersect with the TCA cycle at different points ([Fig 7B](#pgen.1006490.g007){ref-type="fig"}). The fate of 3C/4C skeletons may be determined by the activity of TCA enzymes as well as several other enzymes that catalyze intersecting reactions. Such enzymes include malic enzyme (ME), pyruvate carboxylase (PC), phospho*enol*pyruvate carboxykinase (PEPCK), and phosphoenolpyruvate carboxylase (PEPC), together referred to as the pyruvate hub \[[@pgen.1006490.ref091],[@pgen.1006490.ref105]\]. Though these enzymes likely have important roles in cellular metabolism, though their specific functions remain largely unknown \[[@pgen.1006490.ref090]--[@pgen.1006490.ref091],[@pgen.1006490.ref105]--[@pgen.1006490.ref108]\]. Most of the genes in the mitochondrial pyruvate hub are strongly upregulated in the dark, consistent with enhanced processing of 3C/4C skeletons during this phase ([Fig 7](#pgen.1006490.g007){ref-type="fig"}). Perhaps this is to regulate metabolite abundance and flux into other intersecting pathways as a function of the availability of cofactors required for these reactions, such as NAD(P)H, ATP, and GTP ([Fig 7](#pgen.1006490.g007){ref-type="fig"}). Notably, one of the mitochondrial isoforms of malic enzyme (ME; Phatr3_J54082) is strongly upregulated in the day, in stark contrast to the other pyruvate hub enzymes, suggesting it may have a daytime role in processing malate uncoupled from dark production in the glyoxylate cycle ([Fig 7](#pgen.1006490.g007){ref-type="fig"}).

Due to the central importance of the TCA cycle in cellular energy generation, it is generally subject to post-transcriptional (i.e. allosteric, substrate supply, redox, post-translational) rather than transcriptional (coarse) regulation. Despite this, there were marked day: night fluctuations that were consistent within the greater catabolic context occurring in *P*. *tricornutum* in the dark phase, and indicate there is a clear component of transcript level regulation in facilitating metabolic flux shifts between illuminated and non-illuminated states.

Diel regulation of nitrogen assimilation and the ornithine-urea cycle {#sec009}
---------------------------------------------------------------------

The acquisition of nitrogen is an essential process for maintenance and growth in phytoplankton, and nitrogen assimilation genes are some of the most highly expressed metabolism-related genes in *P*. *tricornutum* ([S12 Fig](#pgen.1006490.s017){ref-type="supplementary-material"}). Nitrogen assimilation, particularly from nitrate, is energetically expensive requiring ATP and NAD(P)H \[[@pgen.1006490.ref109]\]. Additionally, several steps in the pathway have a direct requirement for Fe. None of the genes involved in primary nitrogen assimilation were found to be sensitive to Fe at the transcript level statistically, but strong diel regulation of several key genes involved in primary nitrogen assimilation was observed ([Fig 8](#pgen.1006490.g008){ref-type="fig"}). The nitrate transporter (NRT; Phatr3_J26029), nitrate reductase (NR; Phatr3_54983), formate/nitrite transporter (FNT; Phatr3_J13076), and chloroplast-localized Fd-NIR (Phatr3_J12092), were all coordinately expressed with maximum transcript abundance detected at 10AM (turquoise module, [S14 Fig](#pgen.1006490.s019){ref-type="supplementary-material"}). Additionally, chloroplast-localized GS/GOGAT (GSII; Phatr3_J, GOGATchl; Phatr3_J24739) and mitochondrial GOGAT (GOGAT; Phatr3_J20342) were upregulated in the same way as NRT, NR, FNT, and Fd-NIR, though notably, mitochondrial GS (GSIII; Phatr3_J22357) was not coordinately expressed with these genes pointing to a different role in N metabolism for this enzyme ([S14 Fig](#pgen.1006490.s019){ref-type="supplementary-material"}). The lack of coordinated expression of GSIII with other nitrogen assimilation genes has been previously observed in *T*. *pseudonana*, where GSIII was also found to be upregulated at night \[[@pgen.1006490.ref110]\]. Also in low Fe cultivated *T*. *oceanica*, GS was found to be coordinately upregulated with isocitrate lyase (from the glyoxylate cycle), which is consistent with our observations \[[@pgen.1006490.ref020]\]. With the exception of GSIII, there was coordinate, diurnal expression of the complete pathway for assimilation of nitrate to glutamine in *P*. *tricornutum*, expressed independently of Fe availability. We did not observe the downregulation in nitrogen assimilation transcripts at low Fe observed by Allen et al. \[[@pgen.1006490.ref019]\] in Fe starved *P*. *tricornutum*. We attribute this difference to the fact that in the diel experiments, *P*. *tricornutum* growth was in steady state, indicating that expression of N assimilation transcripts is more affected by growth state than by Fe status alone.

![Diel expression of the N-assimilation pathway.\
Diagram shows the N-assimilation pathway and enzyme subcellular localizations with the general pattern of expression of genes indicated. Compartments are labeled, and the chloroplast is shown with four membranes. NRT = nitrate transporter, NR = nitrate reductase, NOD = nitric oxide dioxygenase, FNT = formate-nitrite transporter, Fd-NIR = ferredoxin-dependent nitrite reductase, NIR = NADH-dependent nitrite reductase, GS(II/III) = glutamine synthetase, GOGAT = glutamine oxoglutarate aminotransferase, NMO = nitronate monooxygenase. See [S6 Fig](#pgen.1006490.s011){ref-type="supplementary-material"} for PIDs and normalized transcript abundance.](pgen.1006490.g008){#pgen.1006490.g008}

Organic nitrogen is obtained through reduction of nitrate, but it can also be re-assimilated from nitric oxide. Nitric oxide is an important part of stress surveillance in diatoms, and a candidate gene involved in production of nitric oxide, nitric oxide synthase (PtNOA, Phatr3_J40200), has been identified in *P*. *tricornutum* \[[@pgen.1006490.ref111]--[@pgen.1006490.ref112]\]. In addition to production by PtNOA, nitric oxide can also be produced by abiotic photolysis of nitrite, or NR acting on nitrite as a substrate instead of nitrate \[[@pgen.1006490.ref113]\]. How and when *P*. *tricornutum* would produce nitric oxide via NR is unclear, but this activity of NR may be either a mechanism to produce this signaling molecule under certain cellular nitrogen and energetic conditions, a wasteful process, or both. Nitric oxide can be recovered to nitrite by nitric oxide dioxygenase (NOD; Phatr3_J45621), which we identify here in *P*. *tricornutum*. Over diel cycles, PtNOA, NR and NOD were all co-expressed with primary nitrogen assimilation genes, supporting active turnover of the nitric oxide pool in coordination with the nitrogen assimilation pathway ([S14 Fig](#pgen.1006490.s019){ref-type="supplementary-material"}). Specifically, this co-expression suggests that the relative activities of these enzymes (particularly NOD and NR) may function together in a system to regulate intracellular levels of nitrate, nitrite, and nitric oxide, and subsequently the magnitude of either N assimilation or NO production.

In marked contrast to most of the other nitrogen assimilation genes, chloroplast NADH-dependent nitrite reductase (NAD(P)H-NIR; Phatr3_EG02286) was upregulated during the dark ([Fig 8](#pgen.1006490.g008){ref-type="fig"}, [S14 Fig](#pgen.1006490.s019){ref-type="supplementary-material"}). NAD(P)H-NIR may sustain assimilation of nitrate during the dark phase, which has been observed to occur in diatoms \[[@pgen.1006490.ref114]--[@pgen.1006490.ref115]\], when the supply of reduced ferredoxin (from PET) is diminished in the dark. In addition to nitrite produced from nitrate (by NR) and from nitric oxide (by NOD), nitrite can be produced from the oxygenation of alkyl nitronates by nitronate monooxygenase (NMO). Two genes coding putative NMOs were identified in the *P*. *tricornutum* genome (Phatr3_J38755, Phatr3_J50353), and like NAD(P)H-NIR, these genes were expressed maximally at night ([Fig 8](#pgen.1006490.g008){ref-type="fig"}, [S14 Fig](#pgen.1006490.s019){ref-type="supplementary-material"}). Taken together, these data suggest that diatom nitrogen assimilation during the day and night phases may rely on temporally partitioned pathways and different nitrogen sources (day = nitrate/nitric oxide, night = alkyl nitronates), expanding on previous observations of diurnal and nitrate-limited expression of nitrogen assimilation genes in diatoms \[[@pgen.1006490.ref110],[@pgen.1006490.ref116]\]. Ammonia produced at night by NAD(P)H-NIR may be assimilated by existing GSII machinery in the chloroplast, or transported to the mitochondria and assimilated into glutamine by GSIII. Nocturnal upregulation of GSIII is consistent with a role in assimilating ammonia released from dark phase proteolysis and amino acid catabolism. *P*. *tricornutum* possesses several ammonia permeases \[[@pgen.1006490.ref117]\], and increases in nocturnal expression of these transporters (Phatr3_J10881, Phatr3_Jdraft972, Phatr3_J13418, Phatr3_J51516, Phatr3_J1862, Phatr3_J27877) indicates enhanced dark phase intracellular ammonia transport, further supporting this interpretation.

Diatoms possess an ornithine-urea cycle (OUC), considered to be an important metabolic hub at the intersection of cellular carbon and nitrogen, which is thought to have conferred a competitive advantage on diatoms \[[@pgen.1006490.ref118]--[@pgen.1006490.ref119]\]. However, little is known regarding the specific role of the OUC across diverse physiological states. Recent work by Levering et al. \[[@pgen.1006490.ref120]\] identified the existence of a chloroplast ornithine cycle, the function of which was hypothesized to exchange NAD(P)H and ATP between the chloroplast and mitochondria. To gain insight into the biological role of these pathways and potential interactions between them, the diel and Fe-specific expression of genes and key metabolites involved in the OUC cycle were analyzed.

Most of the core genes comprising the first four reactions of the urea cycle (unCPS; Phatr3_J24195, OTC; Phatr3_J30514, ASuS; Phatr3_J21116, ASL; Phatr3_J34526) were coordinately regulated with maximal expression detected during the light phase ([Fig 9](#pgen.1006490.g009){ref-type="fig"}, [S14 Fig](#pgen.1006490.s019){ref-type="supplementary-material"}). Arginase (ARG; Phatr3_J38509), however, which catalyzes the final step in the urea cycle was the exception and was most strongly expressed during the dark phase ([Fig 9](#pgen.1006490.g009){ref-type="fig"}). The first four reactions of the urea cycle work together to synthesize the amino acid arginine, which accumulated during the day, and to much higher levels in *P*. *tricornutum* grown at 400 pm Fe′ ([Fig 6](#pgen.1006490.g006){ref-type="fig"} and [Fig 9](#pgen.1006490.g009){ref-type="fig"}). Arginase then catalyzes the formation of ornithine and urea from arginine. Elevated arginase expression was indeed accompanied by a decrease in arginine and concomitant increase in both ornithine and urea, implicating this enzyme in modulating the relative abundance of these metabolites over diel cycles and under balanced growth ([Fig 9](#pgen.1006490.g009){ref-type="fig"}). It is not clear why the complete urea cycle is not coordinately expressed, but it is consistent with earlier observations of urea cycle expression in *P*. *tricornutum* grown on alternate N sources \[[@pgen.1006490.ref118]\], in which arginase transcript abundance did not correlate with expression of the other urea cycle enzymes. We hypothesize that the urea cycle functions primarily in the day to synthesize arginine, and that the fumarate produced in this process is returned to the TCA cycle as malate ([Fig 7](#pgen.1006490.g007){ref-type="fig"} and [Fig 9](#pgen.1006490.g009){ref-type="fig"}). To sustain arginine biosynthesis, an input of ornithine to the mitochondria in the absence of arginase activity is required. Considering the pattern of urea cycle expression and metabolic shifts observed over day: night cycles, we propose a conceptual model in which the chloroplast ornithine cycle and urea cycle are coordinated to regulate metabolite shifts based on cellular energy, carbon, and nitrogen status ([Fig 9](#pgen.1006490.g009){ref-type="fig"}).

![Urea cycle and chloroplast ornithine cycle expression and metabolites.\
(**A**) Pathway diagram shows the localization of reactions in the urea cycle, distributed between the mitochondria and the cytosol, and the chloroplast ornithine cycle \[[@pgen.1006490.ref116],[@pgen.1006490.ref119]\]. Compartments are labeled, and the chloroplast is shown with four membranes. Dashed arrow indicates hypothesized transport. unCPS = carbamoyl-phosphate synthase, OTC = ornithine transcarbamylase, ASuS = arginosuccinate synthase, ASL = arginosuccinate lyase, FUM = fumarase, ARG = arginase, ACOAT = acetylornithine aminotransferase, NAGS = N-acetyl glutamate synthase, GACT = glutamate acetyltranferase, AGK = acetylglutamate kinase, AGPR = acetyl-glutamyl phosphate reductase. (**B**) Metabolite abundance (replicate averaged) over diel cycles in iron-replete conditions, shown normalized to maximum value within each metabolite. Average transcript abundance (across Fe conditions) for arginase is also plotted (with standard deviation).](pgen.1006490.g009){#pgen.1006490.g009}

Unlike the urea cycle, the chloroplast ornithine cycle is capable of synthesizing ornithine *de novo* with N-acetylglutamate synthase (NAGS) catalyzing the first step in which glutamate and acetyl-CoA are utilized to synthesize N-acetylglutamate (NAG). Two putative NAGS were identified in the *P*. *tricornutum* genome, which lacked high confidence plastid-targeting predictions but still may have N-terminal targeting sequences (NAGS1; Phatr3_J10902, NAGS2; Phatr3_J44492, [S15 Fig](#pgen.1006490.s020){ref-type="supplementary-material"}). Under conditions in the chloroplast when glutamate, acetyl-CoA, NADPH, and ATP are elevated (i.e. illuminated, N-replete), ornithine production is favored \[[@pgen.1006490.ref120]\]. Under this scenario, ornithine is then exported from the chloroplast, supplying the urea cycle for arginine biosynthesis, and effectively signaling to the mitochondria that the chloroplast is metabolically "charged". In bacteria and unicellular eukaryotes, there is negative feedback of arginine levels on NAGS activity, thereby modulating arginine biosynthesis \[[@pgen.1006490.ref121]\]. The NAGS isozymes that *P*. *tricornutum* possesses are most similar to sequences from other diatoms, red algae, and oomycetes, suggesting that they were acquired during the secondary endosymbiosis and supporting the idea that they could be regulated by the arginine feedback inhibition characteristic of NAGS from these lineages \[[@pgen.1006490.ref121]\]. Simultaneously in the mitochondria, the first phase of the urea cycle is favored towards arginine biosynthesis when NH~4~^+^, HCO~3~^-^, aspartate, and ATP are abundant, and there is a reliable supply of ornithine imported from the chloroplast, which again occurs during energetically favorable and N-replete conditions. Abundant glutamate favors the production of NAG in the mitochondria (by ACCOAT). In humans, NAG activates the urea cycle by allosterically interacting with CPSI, in order to signal urea production and excretion at high cellular nitrogen levels \[[@pgen.1006490.ref121]\]. Phylogenetically, diatom mitochondrial CPS (unCPS) clades with human CPSI \[[@pgen.1006490.ref118]\]; therefore, it is likely that the form found in diatoms (unCPS) is also activated by NAG. However, the ultimate outcome in diatoms would be to favor arginine biosynthesis at high nitrogen, rather than urea production. In the dark, conditions for arginine biosynthesis are no longer favorable, and arginase is upregulated, converting arginine synthesized in the daytime to ornithine and urea. An increase in mitochondrial ornithine decreases NAG production, slowing the urea cycle and reducing the flux of precursors towards the arginine biosynthesis pathway. Urease (Phatr3_J29702) expression is strongly upregulated at night, correlating with a decrease in cellular urea levels, and increasing the nocturnal supply of nitrogen in the form of ammonia. Other genes involved in the OUC, such as agmatinase (Phatr3\_ J40880), ornithine cyclodeaminase (Phatr3_J35643), and ornithine aminotransferase (OAT; Phatr3_J27726), were also upregulated at night.

This model of nitrogen metabolism explains mechanistically why at low Fe there are statistically significant decreases in metabolites including arginine and urea ([Fig 7](#pgen.1006490.g007){ref-type="fig"}). *P*. *tricornutum* is not limited for N at any of the Fe′ concentrations examined, but the biosynthesis of arginine and urea rely on availability of carbon skeletons (from acetyl-CoA), and sufficient ATP and NAD(P)H, which have been seen to be limiting under low Fe acclimated *P*. *tricornutum*. Further, this model highlights the degree to which the chloroplast and mitochondria are metabolically and energetically connected, and able to communicate their metabolic status through the dynamic exchange of carbon skeletons and nitrogen groups. Despite the essential roles of transporters in regulating pathway connectivity and intracellular metabolite fluxes, these transporters remain conspicuously absent from this and many models of diatom metabolism \[[@pgen.1006490.ref090]--[@pgen.1006490.ref092],[@pgen.1006490.ref108],[@pgen.1006490.ref120]\]. In the *P*. *tricornutum* genome, there are over 350 putative transporters (intracellular and extracellular) that are dynamically expressed over day: night cycles ([S16 Fig](#pgen.1006490.s021){ref-type="supplementary-material"}). The majority of transporters are induced in the dark phase, perhaps to accommodate elevated shuttling of metabolites within cellular compartments during enhanced division, proteolysis, autophagy, and catabolism observed at this time. Effective interpretation of the significance of fluctuating transporter transcripts requires better biochemical characterization and localization of these transporters. Such knowledge will be essential for construction of accurate metabolic models, required to further improve our understanding of the significance of fluctuating transcripts over diel cycles.

Conclusions {#sec010}
===========

Oceanic phytoplankton must possess adaptations that confer survivability in adverse conditions, and strategies used to cope with varying levels of resource availability provide insight into forces that structure phytoplankton populations across different ocean regimes. Initial molecular insights into adaptations diatoms possess to cope with limiting Fe, such as Fe-acquisition machinery, have begun to explain why this group is so successful in Fe-limited regions of the ocean, but have generally not been studied in a dynamic physiological context. In addition to nutrient acquisition, primary producers must cope with vast fluctuations in energy inputs across day: night cycles. Mechanistically, how phytoplankton are able maintain homeostasis and coordinate growth under these dramatic metabolic and energetic shifts is just beginning to be understood in detail.

*P*. *tricornutum* exhibits massive and coordinated shifts in the global transcriptome, metabolome, and physiology over day: night cycles revealing the extent to which fluctuations in light drive major cellular changes. Through incorporating time course expression measurements of Fe-sensitive genes, we identify the night as a particularly important phase for Fe metabolism in both Fe-limited and Fe-replete cultures. At night, Fe-limited cultures also express genes associated with catabolism, proteolysis, autophagy, and endocytosis more highly than Fe-replete cultures, suggesting that efficient recycling of intracellular Fe and acquisition of Fe from the environment are essential for successful acclimation to low Fe. Functional investigation of this highly expanded Fe-sensitive gene catalog, including genes involved in endocytosis and gene regulation (i.e. transcription factors) will be important to further characterize the specific mechanisms diatoms possess to acquire, respond and acclimate to low Fe.

Diel expression of core metabolism reveals the strong influence of light:dark cycles in regulating primary metabolism at the transcriptional level. Expression of these metabolic pathways was largely unaffected by Fe-limitation. However subtle shifts in metabolism, such as the upregulation of cytosolic glycolysis and branched chain amino acid metabolism function to maintain the metabolic and energetic homeostasis required to support growth in low Fe cells. Downregulation of energy acquisition (light harvesting and carbon storage), reduction in growth rate, and upregulation of Fe acquisition appear to work together to "protect" core metabolic function across diel cycles, revealing the metabolic resilience of *P*. *tricornutum* to low Fe. These insights explain further how diatoms persist in Fe-limited areas of the global ocean, and highlight specific adaptations which may have contributed to their ecological dominance in these regions.

Methods {#sec011}
=======

Cultivation and physiological measurements {#sec012}
------------------------------------------

Steady-state cultures of axenic *P*. *tricornutum* (CCMP2561) were maintained under a diel cycle of 12h Light: 12h darkness (320 μmol photons·m^−2^·s^−1^; lights ON at 9 AM; lights OFF at 9PM) at 18°C and at three different concentrations of Fe′ (Fe′: sum of all Fe species not complexed to EDTA): a low concentration (20 pM Fe′), intermediate concentration (40 pM Fe′) and at Fe-replete levels (400 pM Fe′). Cultures were grown in a modification of Aquil medium \[[@pgen.1006490.ref122]\] containing 100, 10 and 100 **μ**M nitrate, phosphate and silicates, standard "f/2" vitamin concentrations, and 300 μM EDTA with free Zn, Mn, Co, Cu, MoO~4~ and SeO~3~ concentrations (expressed as negative logarithms) of 10.93, 8.03, 10.77, 12.63, 7.00 and 8.00. Total Fe concentrations added were 15, 30 and 600 nM Fe. Fe′ was calculated assuming a 2 nM background Fe in our Aquil preparation, and the equations in Sunda et al. \[[@pgen.1006490.ref122]\] considering temperature, pH and diel-integrated light levels. Cells were acclimated to each Fe condition in 250 mL polycarbonate bottles, and growth rates were calculated as the natural log of cell density over time, using a Beckman-Coulter counter. After growth rates changed by less than 10% in subsequent transfers (usually by 30 generations), 100--200 mL of each culture was inoculated into 20 L cultures, and growth was monitored as before. At each time point, 250 mL of cultures was harvested at an average density of 4.84x10^5^ cells mL^-1^ (sd = 1.15x10^5^ cells mL^-1^), and fresh (and equilibrated with respect to metal speciation) media was added to maintain modest cell densities throughout the diel sampling period. Samples were taken three times during the day (10AM, 2PM, 6PM) and during the night (10PM, 2AM, 6AM). Three distinct diel experiments were run, including one with high Fe conditions only and an additional two with all three experimental treatments, leading to independent duplicate diel data sets for the 20 and 40 pM Fe′ treatments and triplicate diel data sets for the 400 pM Fe′ treatment, totaling 49 samples ([S1 Table](#pgen.1006490.s023){ref-type="supplementary-material"}). Growth curves for each experiment are included ([S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"}--[S7 Fig](#pgen.1006490.s012){ref-type="supplementary-material"}). Due to a computational correction, the Fe′ values should be 27, 54, and 540 pM instead of 20, 40, and 400 pM.

At each time point, triplicate determinations of F~v~/F~m~ were made using the DCMU method \[[@pgen.1006490.ref123]\]. The significance of nocturnal drop in F~v~/F~m~ value was tested by a paired t-test of the diurnal maxima vs. the nocturnal minima in either low Fe conditions or Fe-replete conditions.

Percentage dividing cells was microscopically determined by observations of culture samples using a Zeiss Axioplan 200 equipped with differential interference contrast (DIC) module. Therefore, part of the sample at each time point was fixed with formaldehyde at 0.5% final concentration using freshly prepared paraformaldehyde in PBS. Discrimination of dividing vs. interphase cells was achieved at 600--1000× total magnification, based on the quantification of doublet cells. Doublet cells characterized by a daughter cell wall that is only present in mitotic cells and for each sample a minimum number of 200 cells were observed. This quantification resulted in a conservative estimation of M-phase cells, since only part of the M-phase (actual cytokinesis), was considered, thereby excluding nuclear division during karyokinesis.

RNA extraction and sequencing {#sec013}
-----------------------------

During sampling, an average of 1.21 × 10^8^ cells ml^-1^ was harvested (250 mL) on 0.45 **μ**m Sterivex filters and flash frozen in liquid nitrogen. Total RNA was extracted using Trizol reagent: (Thermo Fisher Scientific), enriched and amplified (15 cycles) in mRNA using MessageAmp II aRNA Amplification Kit (Ambion). 80 ng of amplified mRNA was used for SOLiD library construction using SOLiD Total RNA-Seq Kit (Applied Biosystems). SOLiD-based transcriptome sequencing of poly-A enriched mRNA was performed on all 49 experimental time points. Sequence data (n = 49 samples) were deposited and are available at the NCBI Sequence Read Archive (SRA) under the identifier SRA: SRP069841.

Read-mapping and RPKM calculations {#sec014}
----------------------------------

Paired-end SOLiD reads (50 bp forward and 35 bp reverse) were mapped to the *Phaeodactylum tricornutum* genome via Life Technologies LifeScope software (which handles reads in SOLiD color space format directly). This mapping yielded a BAM alignment, which, combined with a file detailing the genomic locations of each *P*. *tricornutum* exon, was converted to a set of read counts per exon using the samtools package. An in-house script was used to aggregate gene counts from exons and to compute RPKM values, where gene RPKM = gene_read_count / (total_mapped_reads (Mb) \* gene_transcript_length (kb))

RNA-Seq reads for each of the 49 libraries were mapped to both the *P*. *tricornutum* Phatr2 assembly (ASM15095v2), which includes 34 finished chromosomes and 55 unassembled scaffolds \[[@pgen.1006490.ref094]\] and the re-annotation of the *P*. *tricornutum* genome (Phatr3) (<http://protists.ensembl.org/Phaeodactylum_tricornutum/>). The re-annotation effort was augmented by 90 new RNAseq libraries, in addition to the 49 considered in this study, leading to the discovery of additional gene models and improvement of gene model boundaries. Quantitative temporal transcript levels were calculated for both versions of the *P*. *tricornutum* genome, yielding RPKM values for 10,577 gene models from Phatr2 and 12,294 gene models from Phatr3. These totals include transcript levels for genes on the chloroplast and mitochondrial genomes (n = 160). Though generally gene models were improved in Phatr3, reads were mapped back to both versions of the gene models to obtain transcript levels for gene models that were erroneously deleted or called incorrectly in Phatr3. A comparison between Phatr2 and the updated Phatr3 gene models, and RPKMs for version of *P*. *tricornutum* gene models are included in [S5 Dataset](#pgen.1006490.s005){ref-type="supplementary-material"}.

Subsequent genome-scale analyses utilized the transcriptome mapped to Phatr3 gene models, and all reported gene expression is from Phatr3 except when specifically denoted. Significance of differences of mean vs. median RPKM ratios were verified using a Wilcoxon rank-sum test performed in Matlab using group-based ratios for all genes in the active transcriptome. This led to the identification of 10,588 genes in the *P*. *tricornutum* genome that were transcriptionally active under the tested experimental conditions. The number of Phatr3 genes which had an average count value greater than 10 and minimum count greater than 1 were scored, resulting in the identification of 10,557 genes which were considered transcriptionally active. Of these, 158 were genes from the mitochondrial and chloroplast genomes, meaning 10,399 nuclear genes were actively transcribed. In the manuscript, both RPKM and standardized RPKM \[[@pgen.1006490.ref124]\] are reported.

Global Expression Analysis {#sec015}
--------------------------

### WGCNA clustering {#sec016}

Correlation network analysis of gene expression diel profiles was performed using the WGCNA package in R, following the principles laid out by Langfelder and Horvath. We specifically applied the comparative consensus network option to compare data sets across conditions \[[@pgen.1006490.ref045]\]. The consensus network was constructed using the 20 pM Fe′ and the 400 pM Fe′ conditions, but to increase performance only data from most highly expressed nuclear genes (8,499 genes; Median ≥ 10 RPKM) were used as input. Application of a *signed* network was crucial to allow clustering of oppositely correlated gene expression values in distinct modules. In addition, by raising the absolute value of the pairwise gene expression correlations to the soft-thresholding parameter *β* = 30, the scale-free topology criterion for WGCNA was achieved, thereby emphasizing high correlations at the expense of low correlations \[[@pgen.1006490.ref045]\]. Selection of *β* occurred as described by Langfelder and Horvath \[[@pgen.1006490.ref045]\] through testing its dependence (5\<*β* \<40) on the scale-free topology fit model for both Fe conditions ([S17 Fig](#pgen.1006490.s022){ref-type="supplementary-material"}). Thus, for construction of the consensus network, the function *blockwiseConsensusModules* was used using the following parameter settings: *β* = 30, *networkType* = \"signed\", *TOMType* = \"signed\". Subsequently, modules were detected on basis of the Topological Overlap Measure (TOM) using the following parameters: *minModuleSize* = 40, *deepSplit* = 2, *maxBlockSize* = 10000, *mergeCutHeight* = 0.20. These parameters resulted in detection of 27 color-coded modules that described distinct units of correlated gene expression in both Fe conditions. The color-code 'grey' ME is reserved for background genes that could not be clustered.

### Response type binning {#sec017}

Genes with significant responses to Fe and illumination state were identified by binning replicate time points per Fe and light condition and applying the non-parametric Skillings-Mack (SM) test for 2-way experimental design in the TM4 gene expression software suite \[[@pgen.1006490.ref125]--[@pgen.1006490.ref129]\]. With significance threshold set at *P* ≤ 0.01, 55% percent of the active transcriptome (5,803 genes; 45% of the *P*. *tricornutum* genome) was sorted into four main expression groups (*P* ≤ 0.01). Genes were categorized as possessing higher transcript abundance in the dark (Dark) or in the light (Light), and at low Fe (20--40 pM Fe′) conditions or at high Fe (400 pM Fe′) relative to the opposite condition. The direction of Fe or light-dependent modulation, either an increase or decrease of transcript abundance in each condition, was identified by calculation of the log~2~ of the average RPKM ratio within the Low Fe vs. High Fe and Light vs. Dark samples, and is reported in [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"}, [S2 Dataset](#pgen.1006490.s002){ref-type="supplementary-material"}, [S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"}, and [S5 Dataset](#pgen.1006490.s005){ref-type="supplementary-material"} as "Iron Depl. vs. Iron Repl" and "Light vs. Dark".

The differential expression of genes responsive to Fe and illumination state was validated using EdgeR with a generalized linear model using the likelihood ratio test (LR) in a two factor analysis \[[@pgen.1006490.ref130]\]. P-values were adjusted by the Benjamini-Hochberg procedure for multiple testing correction to control for false discovery rate at a level of adjusted p \< 0.05. Response type bins for each gene are shown in [S1 Dataset](#pgen.1006490.s001){ref-type="supplementary-material"} with genes found to be significant for the same response category by both methods called "high confidence". The larger set of genes found to be significant by the SM test are considered the primary response type bins unless otherwise specified, including those that are high confidence. This approach to binning was chosen to represent a more inclusive set and was, in this dataset, found to be more consistent with genes known to be responsive to Fe. Among the iron starvation inducible protein (ISIP) genes, all were called significant in low Fe conditions by SM, but ISIP2a was not called significant by EdgeR.

### GO functional enrichment calculations {#sec018}

Gene functional enrichments (FDR-corrected P ≤ 0.05) were identified in differentially modulated gene groups against the whole genome as a reference. Significance was assessed using a Binomial test and Benjamini & Hochberg\'s FDR multiple testing correction using the BiNGO Cytoscape plugin \[[@pgen.1006490.ref131]\].

Identification and curation of metabolic pathways {#sec019}
-------------------------------------------------

Genes were assigned to metabolic pathways and transport processes by manual curation in Phatr3, using a combination of approaches and criteria. First, genes were identified with keyword searches of both existing and *de novo* gene catalog annotations, and published pathway annotations specific to *P*. *tricornutum* \[[@pgen.1006490.ref085]--[@pgen.1006490.ref086],[@pgen.1006490.ref090]--[@pgen.1006490.ref092],[@pgen.1006490.ref132]\]. To identify all desired isozymes of individual steps in desired pathways, we used BLASTP against the *P*. *tricornutum* genome (draft version at NCBI, Phatr2 version at JGI, and/or Phatr3 version at ENSEMBL) utilizing query sequences from validated genes from model organisms, or *P*. *tricornutum* isozymes. Identified genes were then curated and functionally assigned based on a variety of criteria. These criteria were the presence of specific conserved domains (protein queries to the NCBI Conserved Domain Database), subcellular bioinformatics-based targeting predictions, and expert knowledge of biochemistry and biochemical pathways. Subcellular targeting predictions were improved in Phatr3 with extensions to gene models to include full N-terminal ends where targeting peptides are found. When high-confidence subcellular localization predictions were important for pathway curation, these predictions were further refined by using SignalP 3.0 \[[@pgen.1006490.ref133]\], SignalP 4.1 \[[@pgen.1006490.ref134]\], ASAFind v1.0.0 \[[@pgen.1006490.ref135]\], Mitoprot \[[@pgen.1006490.ref136]\], and TMHMM \[[@pgen.1006490.ref137]\].

Metabolite analysis {#sec020}
-------------------

Metabolites were extracted and derivatization and GC-MS analysis were carried out as described in Lisec et al. \[[@pgen.1006490.ref138]\], with minor modifications made for diatom analysis as defined in Allen et al. \[[@pgen.1006490.ref019]\]. Metabolites were identified through comparison with database entries of authentic MSRI libraries \[[@pgen.1006490.ref139]\]. The impact of illumination or Fe availability on the abundance of 35 metabolites was determined using T-test (*P*-value \< 0.05).

Supporting Information {#sec021}
======================

###### Active transcriptome and assignment of genes to WGCNA modules and response types.

(XLSX)

###### 

Click here for additional data file.

###### List of highly low Fe responsive genes.

(XLSX)

###### 

Click here for additional data file.

###### List of low Fe induced genes involved in endocytosis, exocytosis and autophagy.

(XLSX)

###### 

Click here for additional data file.

###### Metabolic pathway annotations.

(XLSX)

###### 

Click here for additional data file.

###### Phatr3 and Phatr2 comparison.

(XLSX)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment L1 at 20 pM Fe′.

The start of the diel is indicated by filled inverted triangles, times of carboy inoculation used for diel is indicated by hollow triangles, and time is given as ordinal day 2010. In (**A**) and (**C**), breaks in regression lines indicate culture transfer (prior to diel sampling) or dilution (during diel sampling). (**A**) Cell density versus time seven days prior to and during the diel sampling period. (**B**) Cumulative cell density computed from initial starting density (disregarding dilution events) in order to calculate and visualize growth history prior to and during diel experiments. Growth rates, as the slope of the natural log of cell density are shown for the period preceding and after carboy inoculation. (**C**) Cell densities determined by coulter counting during the diel. Vertical lines indicate transition between light and dark phases. Cell densities were maintained from 400,000--530,000 cells mL^-1^ during the diel.

(AI)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment L2 at 20 pM Fe′.

As in [S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"}, except in panel (**C**), cell densities were maintained from 400,000--526,000 cells mL^-1^ during the diel.

(AI)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment M1 at 40 pM Fe′

As in [S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"}, except in panel (**C**), cell densities were maintained from 415,000--840,000 cells mL^-1^ during the diel.

(AI)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment M2 at 40 pM Fe′.

As in [S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"} except in panel (**C**), cell densities were maintained from 375,000--517,000 cells mL^-1^ during the diel.

(AI)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment H1 at 400 pM Fe′.

As in [S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"}, except in panel (**C**), cell densities were maintained from 360,000--946,000 cells mL-1 during the diel.

(AI)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment H3 at 400 pM Fe′.

As in [S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"}, except in panel (**C**), cell densities were maintained from cells mL-1 during the diel.

(AI)

###### 

Click here for additional data file.

###### Natural log of cell density versus time preceding and during experiment H2 at 400 pM Fe′.

As in [S1 Fig](#pgen.1006490.s006){ref-type="supplementary-material"}, except in panel (**C**), cell densities were maintained from 340,000--540,000 cells mL^-1^ during the diel.

(AI)

###### 

Click here for additional data file.

###### Expansion of low Fe sensitive gene catalog in P. tricornutum.

Diel expression of low Fe response type genes (n = 2090, [S2 Dataset](#pgen.1006490.s002){ref-type="supplementary-material"}, methods). Transcript abundance patterns are hierarchically clustered (Pearson correlation). Heat map shows standardized transcript abundance (RPKM across Fe conditions) for each gene. Black bars indicate which genes were identified as low Fe responsive by each statistical approach (Skillings-Mack vs. EdgeR), as well as genes previously identified as low Fe responsive in \[[@pgen.1006490.ref019]\].

(AI)

###### 

Click here for additional data file.

###### Transcription factor expression pattern.

(**A**) Heat map shows standardized transcript abundance (across Fe conditions) for annotated transcription factors \[[@pgen.1006490.ref054]\] clustered hierarchically (Pearson correlation). Adjacent heat maps indicate the maximum transcript detected for each gene across experiments (RPKM), and the degree to which genes responded to Fe (see [methods](#sec011){ref-type="sec"}). Names and PIDs of the top ten most highly expressed genes and the top ten most highly iron-responsive genes are labeled, with overlap between these lists noted. Predicted organellar localization is indicated. (**B**) Transcript abundance of CCHH11 is indicated, and shown RPKM values from each experimental time point (not averages). (**C**) Output from TMHMM showing a C-terminal transmembrane domain in the CCHH11 gene model (Phatr3_J38018).

(AI)

###### 

Click here for additional data file.

###### Cyclin and CDK expression.

Heat map shows transcript abundance (RPKM) for *P*. *tricornutum* cyclins and CDKs \[[@pgen.1006490.ref062]\]. Genes binned as low Fe responsive are shown in bold (\*Binned by Skillings-Mack, †Binned by EdgeR). Dendogram shows hierarchical clustering (Pearson correlation) of RPKM values. Gene ID, gene names, and relative expression level (max RPKM for each gene normalized to max across entire gene set) included to illustrate which cyclins/CDKs transcripts are most abundant. Cyclin and CDK transcripts previously assigned to cell cycle phases \[[@pgen.1006490.ref062]\] are indicated.

(AI)

###### 

Click here for additional data file.

###### Light harvesting antenna protein expression.

Heat map shows transcript abundance (RPKM, plotted on a log scale). Dendogram shows hierarchical clustering (Pearson correlation) of RPKM values. Protein ID, gene names, WGCNA module and response type (SM: Skillings-Mack and EdgeR) are indicated for each gene.

(AI)

###### 

Click here for additional data file.

###### Most highly expressed primary metabolic pathways.

(**A**) Bar plot shows maximum and median RPKM detected for each gene across the entire experiment (all time points and Fe conditions) annotated with a primary metabolic function, organized by pathway ([S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"}). (**B**) List of the pathways included and the number and percentage of genes assigned to each pathway with a maximum RPKM greater than 500. (**C**) List of the genes within the most highly expressed pathways with a maximum RPKM greater than 500, with maximum and median values indicated. Genes binned as Fe-responsive (low Fe and high Fe) are indicated. \*Binned by Skillings-Mack, †Binned by EdgeR

(AI)

###### 

Click here for additional data file.

###### Carbohydrate biosynthesis and degradation.

(**A**) Pathway diagram shows the putative localization of glycolysis, gluconeogenesis, and carbohydrate biosynthesis and degradation enzymes. Compartments are labeled, and the chloroplast is shown with four membranes. Each line connecting metabolites indicates an isozyme. (**B**) Biological replicate averaged transcript abundance (for each time point) shown for genes from gluconeogenesis, carbohydrate biosynthesis, carbohydrate degradation, and dark phase Calvin Benson cycle. Error bars omitted for clarity, series labeled with gene names (see [S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"} for abbreviations) and Phatr3 PID.

(AI)

###### 

Click here for additional data file.

###### Expression patterns of nitrogen metabolism genes.

Average-normalized transcript abundance for genes involved in nitrogen assimilation, dark phase nitrogen assimilation, urea cycle, and chloroplast ornithine cycle is shown. Series labeled with gene names (see [S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"} for abbreviations) and Phatr3 PID.

(AI)

###### 

Click here for additional data file.

###### N-acetylglutamate synthase (NAGS) gene candidates.

(**A**) Table summarizing targeting predictions (see [methods](#sec011){ref-type="sec"}) for both NAGS1 and NAGS2 gene candidates. (**B**) Sequences used for gene targeting, with Phatr3 gene model ID and EST-based N-term extensions (NAGS1) indicated. (**C**) Domain architecture (from NCBI blastp) showing hits to PRK05279 (N-acetylglutamate synthase, validated).

(PDF)

###### 

Click here for additional data file.

###### Transporter expression patterns.

Heat map shows standardized transcript abundance (across Fe conditions) of all putative transporters clustered hierarchically (Pearson correlation, [S4 Dataset](#pgen.1006490.s004){ref-type="supplementary-material"}). Adjacent heat maps indicate the maximum transcript detected for each gene across experiments, and the degree to which genes responded to Fe (see [methods](#sec011){ref-type="sec"}).

(AI)

###### 

Click here for additional data file.

###### Scale-free topology model fit as function the soft thresholding power (x-axis).

Numbers in the plots indicate the corresponding soft thresholding powers. The plot indicates that the approximate scale-free topology is attained around the soft-thresholding power of 30 for both Fe conditions.

(AI)

###### 

Click here for additional data file.

###### Experimental design.

(PDF)

###### 

Click here for additional data file.

###### GO term enrichment in WGCNA modules.

(PDF)

###### 

Click here for additional data file.

The author would like to acknowledge Stephane Lefebvre, James McCarthy, and Jeff McQuaid for valuable discussions and feedback, which assisted with data interpretation and greatly enhanced the quality of this study.

[^1]: The authors have declared that no competing interests exist.

[^2]: **Conceptualization:** ABK AEA.**Data curation:** ABK JPM JHB TO ARF.**Formal analysis:** SRS JTFG JPM JHB HZ TO.**Funding acquisition:** ABK CLD AEA.**Investigation:** SRS JTFG ABK HZ TO ARF AMN.**Methodology:** SRS JTFG ABK JPM JHB ARF AEA.**Project administration:** ABK CLD AEA.**Resources:** ABK CLD ARF AEA.**Software:** JTFG JPM JHB.**Supervision:** ABK CLD AEA.**Validation:** SRS JTFG ABK AEA.**Visualization:** SRS JTFG JPM.**Writing -- original draft:** SRS JTFG ABK CLD ARF AEA.**Writing -- review & editing:** SRS JTFG ABK JPM JHB HZ CLD TO ARF AEA.

[^3]: Current address: Center for Cancer Research, Bethesda, Maryland, United States of America
