Proximity proteomics in a marine diatom reveals a putative cell surface-to-chloroplast iron trafficking pathway
Abstract
Iron is a biochemically critical metal cofactor in enzymes involved in photosynthesis, cellular respiration, nitrate assimilation, nitrogen fixation, and reactive oxygen species defense. Marine microeukaryotes have evolved a phytotransferrin-based iron uptake system to cope with iron scarcity, a major factor limiting primary productivity in the global ocean. Diatom phytotransferrin is endocytosed; however, proteins downstream of this environmentally ubiquitous iron receptor are unknown. We applied engineered ascorbate peroxidase APEX2-based subcellular proteomics to catalog proximal proteins of phytotransferrin in the model marine diatom Phaeodactylum tricornutum. Proteins encoded by poorly characterized iron-sensitive genes were identified including three that are expressed from a chromosomal gene cluster. Two of them showed unambiguous colocalization with phytotransferrin adjacent to the chloroplast. Further phylogenetic, domain, and biochemical analyses suggest their involvement in intracellular iron processing. Proximity proteomics holds enormous potential to glean new insights into iron acquisition pathways and beyond in these evolutionarily, ecologically, and biotechnologically important microalgae.
Introduction
Iron (Fe) likely played an important role in the origin of life (Bonfio et al., 2017; Jin et al., 2018) and is fundamental to almost all extant metabolisms (Aguirre et al., 2013) serving as a cofactor in enzymes involved in DNA replication, photosynthesis, cellular respiration, nitrate assimilation, nitrogen fixation, and reactive oxygen species defense (Crichton, 2016). Early anoxic oceans were rich in readily bioavailable ferrous, Fe2+, iron, but with the rise of oxygenic photosynthesis and the Great Oxygenation Event (GOE) in the Paleoproterozoic ~2.3 billion years ago followed by the Neoproterozoic Oxygenation Event (NOE) ~1.5 billion years later, most of the ferrous iron oxidized into insoluble ferric, Fe3+, minerals which are not bioavailable (Camacho et al., 2017; Knoll et al., 2016; Och and Shields-Zhou, 2012). Conceivably, these large global shifts in ocean chemistry could have had a major role in driving the evolution of novel iron uptake mechanisms.
Iron availability limits primary productivity in high-nutrient, low-chlorophyll (HNLC) regions which cover ~25% of the modern ocean habitat as demonstrated by numerous large-scale Fe fertilization experiments and natural Fe upwelling events invariably resulting in diatom-dominated phytoplankton blooms (Ardyna et al., 2019; de Baar, 2005; Martin et al., 1994). Although dissolved Fe in marine environments is primarily found complexed to organic ligands (Hutchins and Boyd, 2016; Tagliabue et al., 2017) such as bacterially produced siderophores and hemes (Boiteau et al., 2016; Hogle et al., 2014), low (pM) amounts of unchelated labile Fe, Fe’, are crucial for eukaryotic phytoplankton, particularly diatoms (Morel et al., 2008), widespread single-celled microalgae critical to oceanic primary production (Benoiston et al., 2017). Diatoms have convergently evolved phytotransferrin (pTF), which serves as the basis of a high-affinity ferric, Fe3+, iron binding, and acquisition pathway (McQuaid et al., 2018; Morrissey et al., 2015). Phytotransferrin sequences have a broad taxonomic distribution and are abundant in marine transcriptomic datasets (Bertrand et al., 2015; Marchetti et al., 2012).
Phytotransferrins are estimated to have emerged concurrently with the NOE exemplifying the link between large environmental changes and molecular innovation (McQuaid et al., 2018). They are phylogenetically related to Fe-assimilation (FEA) domain-containing proteins found in freshwater and marine single-celled algae, and represent functional analogues of iron delivery proteins transferrins (Behnke and LaRoche, 2020; Cheng et al., 2004; McQuaid et al., 2018; Scheiber et al., 2019). pTF (Ensembl ID: Phatr3_J54465, UniProt ID: B7FYL2) from the model marine diatom Phaeodactylum tricornutum (Bowler et al., 2008) was first identified as an iron deficiency-sensitive gene and named ISIP2a (iron starvation induced protein 2a) (Allen et al., 2008), although its expression is relatively high in iron-replete conditions as well (Smith et al., 2016). pTF localizes to the cell surface and intracellular puncta, presumably endosomal vesicles (McQuaid et al., 2018), which is further supported by the presence of an endocytosis motif in the C-terminus of the protein (Lommer et al., 2012). In diatom pTF, dissolved labile ferric iron is coordinated synergistically with carbonate anion (CO32-), a common feature of all transferrins (Cheng et al., 2004; McQuaid et al., 2018; Zak et al., 2002). Thus, decline in seawater carbonate concentrations due to ocean acidification caused by elevated atmospheric CO2 may negatively impact this prevalent iron uptake system in diatoms and other marine phytoplankton (McQuaid et al., 2018).
While accessing and binding dilute dissolved labile ferric iron on the cell surface is important, its subsequent internalization and delivery to target sites within complex cellular milieu are also critical. Fe3+ is highly insoluble and conversion between Fe3+ and Fe2+ can lead to toxic reactive oxygen species causing damage to proteins, lipids, and nucleic acids (Cheng et al., 2004). Precise control of iron internalization and intracellular trafficking in either of its redox states is thus crucial for maintenance of cellular homeostasis (Philpott and Jadhav, 2019; Wang and Pantopoulos, 2011). Human ferric iron-laden transferrin (Tf) bound to transferrin receptor is internalized via clathrin-mediated endocytosis. Endosome acidification leads to carbonate protonation and Tf conformational change resulting in iron release. Fe3+ is then reduced by six-transmembrane epithelial antigen of prostate 3 (STEAP3), exported to the cytoplasm through divalent metal transporter 1 (DMT1; also NRAMP2) or Zrt- and Irt-like protein (ZIP) 14, and offloaded to iron chaperones for distribution to cellular iron sinks (Bogdan et al., 2016; Eckenroth et al., 2011; Ohgami et al., 2005; Philpott and Jadhav, 2019; Wang and Pantopoulos, 2011). In contrast, proteins that mediate intracellular allocation of internalized Fe3+ and conduct specific biochemical transformations downstream of diatom pTF are unknown. Elucidation of these pathways will advance our understanding of diatom adaptation to survival in iron-deficient ocean waters and provide insight into how climate change will impact diatom physiology and fitness.
To investigate the proximal proteomic neighborhood of P. tricornutum pTF, we employed APEX2, an engineered soybean ascorbate peroxidase that functions as a dual probe for electron microscopy and proximity proteomics (Hung et al., 2016; Lam et al., 2015; Martell et al., 2017). When fused with a protein of interest (bait protein), the ~27 kDa APEX2 enzyme permits spatially resolved proteomic mapping by oxidizing biotin-phenol to short lived (<1 ms) phenoxyl radicals which can covalently react with electron-rich amino acid residues—primarily tyrosine—on the surface of nearby endogenous proteins. Tagged proteins can then be isolated by purification with streptavidin beads and analyzed using mass spectrometry (MS). While the ‘biotinylation radius’ in APEX2 experiments is estimated to be in the 10–20 nm range, it should rather be seen as a ‘probability gradient’ where the likelihood of tagging decreases with distance away from the APEX2-tagged protein of interest (Gingras et al., 2019; Hung et al., 2016; Lam et al., 2015). Proteomic hits can therefore include vicinal proteins which interact with the APEX2-tagged bait protein indirectly in addition to direct interactors of varying strength (Lundberg and Borner, 2019). Tha ability to capture weakly and transiently interacting proteins in a biological process of interest is a major advantage APEX2 offers over traditional pull-down approaches (Gingras et al., 2019). Swaping biotin-phenol for diaminobendizine (DAB)—another APEX2 substrate—enables one to determine fusion protein localization via high-resolution electron microscopy making APEX2 a powerful bifunctional probe (Martell et al., 2017; Martell et al., 2012).
APEX2 and related enzyme-based chemical biology tools have over the past decade been applied to investigate cellular compartments and processes in diverse model systems including mammalian cell culture, yeast, protists, and plant protoplasts (Gingras et al., 2019; Trinkle-Mulcahy, 2019). APEX2 was recently used to map proteomes associated with post-Golgi vesicle trafficking and endocytosis (Del Olmo et al., 2019; Otsuka et al., 2019) as well as receptor-mediated signaling (Paek et al., 2017), successes that are all directly pertinent to the goals of our study.
Here, we provide further evidence for P. tricornutum phytotransferrin pTF endocytosis, demonstrate correct subcellular localization and activity of its APEX2 fusion, and identify nearly 40 likely proximal proteins through a proximity-dependent proteomic mapping experiment conducted in quintuplicate. To the best of our knowledge, this study represents the first application of APEX2 in a marine microbial model system (Figure 1). We then provide initial characterization, using dual fluorophore protein tagging fluorescence microscopy, phylogenetic, domain, and biochemical analyses, for three proteins expressed from a chromosomal gene cluster (pTF.CREG1, pTF.CatCh1, pTF.ap1) believed to act in association with endocytosed pTF downstream of ferric iron binding at the cell surface. Bacterially expressed recombinant pTF.CREG1 was found to be catalytically dead indicating it may have a regulatory, rather than enzymatic, role in pTF endocytosis, similar to functions proposed for human CREG1. Meanwhile, bioinformatic interrogation of pTF.CatCh1 and its localization pattern suggest it might be a chloroplast-associated iron-binding protein. Finally, we introduce a model for an intracellular ferric iron allocation pathway in P. tricornutum based on our proximity proteomics hits.
Results
P. tricornutum phytotransferrin (pTF) is localized to intracellular vesicles
To confirm pTF association with intracellular vesicles as suggested by McQuaid et al., 2018, a pTF-mCherry encoding episome was conjugated into ΔpTF P. tricornutum cells (Diner et al., 2016; Karas et al., 2015; McQuaid et al., 2018). After labeling the fluorescent transconjugant cell line with 100 μM of the membrane dye MDY-64 for 10 min at room temperature, colocalized mCherry and MDY-64 signals were observed on the cell surface and within intracellular vesicles (Figure 2A).
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/52770/elife-52770-fig2-data1-v3.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/52770/elife-52770-fig2-data2-v3.zip
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/52770/elife-52770-fig2-data3-v3.zip
APEX2 fused with P. tricornutum phytotransferrin (pTF) is enzymatically active in vivo
To generate pTF-APEX2 expressing diatom cell lines, an episome encoding pTF with C-terminal APEX2 was assembled and conjugated into WT P. tricornutum cells (Figure 2B). Considering the predicted pTF domains (Figure 2—figure supplement 1A), APEX2 was likely facing the cytosol at the cell surface and once internalized into vesicles. Western blot with pTF-specific antibodies (McQuaid et al., 2018) confirmed the ~84.7 kDa fusion protein expression in four out of five tested transconjugant cell lines. The protein was present in the insoluble cell lysate fraction further indicative of its membrane localization (Figure 2C). Amplex UltraRed, a highly sensitive APEX2 substrate (Hung et al., 2016), was used to assay live pTF-APEX2 expressing cells. A resorufin (reaction product of Amplex UltraRed and APEX2) signal up to 4- and 40-fold above WT background was observed in experiments performed at room temperature (data not shown) and on ice (Figure 2D; Figure 2—source data 1), respectively, indicating active APEX2 with incorporated heme. Resorufin was also directly visualized by confocal microscopy and a strong cytosolic signal not tightly localized to the expected site of origin, similar to previous reports (Martell et al., 2012), was observed (Figure 2—figure supplement 1B), perhaps indirectly supporting our predicted APEX2 orientation. APEX2 was active only in the presence of both Amplex UltraRed and hydrogen peroxide (Figure 2—figure supplement 1C) implying that endogenous H2O2 levels are not sufficient to drive APEX2-catalyzed reactions, and that the overall cell surface and intracellular milieu in P. tricornutum is permissive to the APEX2 catalytic cycle.
pTF-APEX2 is localized to the cell membrane and intracellular vesicles
To confirm pTF-APEX2 is localized to the cell surface and intracellular vesicles similar to the pTF-mCherry fusion protein, pTF-APEX2 expressing cells (cell line s2, Figure 2C–D) were treated with 25 mM 3,3’-diaminobenzidine (DAB) in the presence of 3 mM H2O2 for 15 min on ice (Figure 3A). This reaction leads to DAB polymerization and local precipitation around APEX2 that can be stained with osmium tetroxide and visualized with an electron microscope. Cells were embedded in a 3% agar matrix to prevent losses during numerous washing steps (Figure 3—figure supplement 1A). Tightly localized signal was observed on the cell membrane and in intracellular vesicles indicating that pTF-APEX2 fusion trafficked to the correct subcellular sites (Figure 3B). Additionally, we observed mitochondrial signal in WT and pTF-APEX2 cell lines subjected to the DAB reaction (Figure 3B; Figure 3—figure supplement 1B). Analysis of the P. tricornutum proteome revealed eight peroxidases with Arg38, His42, His163, and Asp208; catalytic amino acid residues that are conserved across all identified ascorbate peroxidases (APX) including soybean APX (Raven, 2003) and its derivative APEX2 (Figure 3—figure supplement 1C). Three of them contain proline in place of alanine at position 134 relative to APEX2, a substrate-binding loop mutation rendering APEX2 highly active (Lam et al., 2015), and another one is predicted to localize to mitochondria (Figure 3—figure supplement 1D–E). It is plausible that one or more of these endogenous APEX2-like peroxidases are responsible for the observed mitochondrial signal. APEX2 was inactivated in iron-deplete conditions (40 nM total Fe) and could be reactivated by supplementing fixed cells with 10 μM hemin chloride for 3 hr (data not shown), but this substantially increased background in the Amplex UltraRed assay. Therefore, subsequent proximity labeling experiments were carried out on cells growing in iron-replete conditions.
Identification of the proximal phytotransferrin proteome with biotin-phenol labeling
To identify proteins proximal to pTF, quintuplicate cultures from one WT and one pTF-APEX2 expressing cell line (cell line s2, Figure 2C–D) were grown to mid- to late-exponential phase and supplemented with 2.5 mM biotin-phenol and 1 mM hydrogen peroxide, following the APEX2 proximity labeling protocol developed for yeast (Hwang and Espenshade, 2016; Hwang et al., 2016). These labeling reaction steps were performed at 4°C and the incubation with hydrogen peroxide was extended to 20 min to mirror Amplex UltraRed assay and DAB reaction conditions (Figure 4A). Increasing biotin-phenol concentration from 0.5 mM, usually used in mammalian cells, to 2.5 mM, and exposing cells to osmotic stress with 1.2 M sorbitol—both of which are critical for efficient labeling in yeast (Hwang and Espenshade, 2016)—was necessary to detect enrichment of biotinylated proteins in experimental samples (Figure 4B). This result confirmed our hypothesis that the lack of heavy silicification in P. tricornutum (Francius et al., 2008; Tesson et al., 2009), making its cell wall composition and cell membrane permeability likely similar to that of yeast, would permit labeling. Streptavidin pull-downs were then performed with clarified cell lysates followed by tandem mass tag (TMT)-based quantitative proteomics. WT and pTF-APEX2 proteomic replicates, with the exception of one WT and one pTF-APEX2 sample, formed two distinct clusters (Figure 4—figure supplement 1; Figure 4—figure supplement 1—source data 1) indicating minimal technical variability. Thirty-eight statistically significant proteins (p value ≤ 0.05) with APEX2/WT ratios of at least 1.5 were identified (Figure 4B; Figure 4—source data 1). These ratios were obtained from average total peptide counts across quintuplicates (Figure 4—source data 1). Endogenous biotinylated proteins were also detected and had APEX2/WT ratios close to one, thus acting as an intrinsic pull-down control (Supplementary file 1—Table S1; Figure 4—source data 1). Some background enrichment in WT cells was likely due to endogenous APEX2-like peroxidases and one would therefore expect mitochondrial proteins with APEX2/WT ratios close to one to be present in our MS dataset. Indeed, at least three were detected: mitochondrial chaperonin CPN60 (Phatr3_J24820, UniProt ID: B7FQ72, APEX2/WT = 0.99), mitochondrial import receptor subunit TOM70 (Phatr3_J47492, UniProt ID: B7G3J4, APEX2/WT = 0.91), and acetyl-CoA dehydrogenase (Phatr3_J11014, UniProt ID: B7FTR6, APEX2/WT = 1.23) (Figure 4—source data 1). Fourteen proteins with an APEX2/WT ratio of at least two were detected. Of these 14 proteins, nine are known to be transcriptionally sensitive to iron availability; this is also the case for an additional five proteins with APEX2/WT ratio of at least 1.5 (Figure 5—figure supplement 1—source data 1; Smith et al., 2016). We note that pTF was present, but not enriched, in our pTF-APEX2 proteomics samples (Figure 4—source data 1). One possible explanation for this result is that pTF may not have surface exposed amino acid residues that are permissive to biotinylation by APEX2-generated phenoxyl radicals and so its presence in all of our pull-down samples could be due to unspecific binding to streptavidin beads.
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/52770/elife-52770-fig4-data1-v3.xlsx
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/52770/elife-52770-fig4-data2-v3.zip
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/52770/elife-52770-fig4-data3-v3.zip
Some of the most apparent biologically interesting hits are summarized in Figure 4C. Considering iron is critically involved in all of the major photosynthetic complexes—photosystem II, photosystem I, and cytochrome b6f (Rochaix, 2011)—we asked if some of the proteins are perhaps predicted to be targeted to the chloroplast. Indeed, five are, albeit with low confidence, by ASAFind (Gruber et al., 2015), including two (Phatr3_J51183 and Phatr3_J54986) that are part of a gene cluster on chromosome 20 and two (Phatr3_J41423 and Phatr3_J55031) that are known to associate with the P. tricornutum chloroplast (Figure 4C; Allen et al., 2012; Kazamia et al., 2019; Kazamia et al., 2018). We elaborate on the potential role for these proteins in intracellular iron trafficking in the Discussion.
Proteins encoded by an iron-sensitive gene cluster on chromosome 20 colocalize with pTF
We focused on three proteins with APEX2/WT ratio of at least two—Phatr3 IDs: J51183 (hereafter pTF.CREG1), J52498 (hereafter pTF.CatCh1), and J54986 (hereafter pTF.ap1)—that are expressed from a previously identified iron- and silicon-sensitive gene cluster on chromosome 20 (Figure 5A; Allen et al., 2008; Sapriel et al., 2009). One additional protein, ISIP2b (Phatr3_J54987, UniProt ID: B7G9B1), which is also transcriptionally sensitive to iron and silicon and that we did not detect in our proximity proteomics experiment, is also encoded by this uncharacterized locus. We note that pTF is not co-located with these genes, but instead lies on chromosome 7 (genomic location 1,000,053–1,001,833, forward strand). All three genes exhibit a transcriptional profile similar to pTF (Figure 5—figure supplement 1), two proteins (pTF.CREG1 and pTF.ap1) are predicted to go to the chloroplast (with low confidence as determined by ASAFind), and two (pTF.CatCh1 and pTF.ap1) contain a C-terminal transmembrane domain (Figure 4C; Supplementary file 1—Table S2). To test whether they colocalize with pTF, co-expression episomes were assembled and conjugated into WT P. tricornutum cells which resulted in diatom cell lines expressing pTF-mCherry and MS hit-EYFP fusion proteins (Figure 5—figure supplement 2A–B). Imaging conditions were optimized with mCherry and Venus (yellow fluorescent protein with spectral properties very similar to those of EYFP) expressing P. tricornutum cell lines for minimal cross-channel bleed-through (Figure 5—figure supplement 2C). pTF.CREG1 and pTF.CatCh1 colocalized with pTF in the chloroplast vicinity and on the chloroplast margin, respectively (Figure 5B). Both proteins consistently exhibited these distinct colocalization patterns as demonstrated by additional three mCherry- and EYFP-positive cells from corresponding co-expression transconjugant cell lines (Figure 5—figure supplement 3A–B). Colocalization of pTF.ap1 with pTF close to the chloroplast was also evident, although somewhat less precise (Figure 5—figure supplement 4).
pTF.CREG1 has homologs across the tree of life and may possess a non-enzymatic role similar to human CREG1
To shed light on possible functions of the uncharacterized proteins colocalizing with pTF and to examine their occurrence in other diatoms and marine phytoplankton beyond Phaeodactylum tricornutum, phylogenetic analysis was performed. pTF.CREG1 homologs were identified across the tree of life (Figure 6; Figure 6—source data 1). Further investigation of the protein sequence alignment underlying the phylogenetic tree revealed 12 amino acid residues that are at least 90% conserved across all homologs (Figure 6—figure supplement 1). Diatom homologs can be seen in a crown group with proteins from other complex plastid-containing algae such as cryptophytes, haptophytes, pelagophytes, chlorarachniophytes, and dinoflagellates (Figure 6; Archibald, 2009; Füssy and Oborník, 2018). pTF.CREG1 and its iron-insensitive paralog Phatr3_J10972—which phylogenetically cluster together—have ~50% amino acid sequence similarity and belong to the flavin mononucleotide (FMN)-binding split barrel fold protein superfamily (SSF50475) according to their Ensembl Protists profiles.
We analyzed pTF.CREG1 amino acid sequence (UniProt ID: B7G9B3, 234 AAs, ~26.4 kDa) and compared it to the human CREG1 ortholog (UniProt ID: O75629, 220 AAs, ~24 kDa). While human CREG1 is evolutionarily and structurally closely related to flavin mononucleotide (FMN)-coordinating dimeric oxidases and oxidoreductases such as pyridoxine 5’-phosphate (PNP) oxidases, it is unable to bind FMN and catalyze NADPH oxidation (Ghobrial et al., 2018; Sacher et al., 2005). pTF.CREG1 has a predicted N-terminal signal peptide ending in SGA that may be recognized by an endoplasmic reticulum-localized signal peptidase (Benham, 2012; Tuteja, 2005), and is likely modified by both N- and O-glycosylation, in line with human CREG1 (Figure 6—figure supplement 2A; Sacher et al., 2005). These features suggest that pTF.CREG1 enters and travels along the secretory pathway (Bard and Chia, 2016; Benham, 2012). Furthermore, the amino acid sequence alignment between pTF.CREG1 and human CREG1 revealed a similar tetrapeptide motif DP(Q/E)S which is located on the loop occluding the FMN-binding pocket at the human CREG1 dimer interface (Figure 6—figure supplement 1; Figure 6—figure supplement 2A; Sacher et al., 2005). This short amino acid sequence—while thought to contribute to the lack of catalytic activity in human CREG1—was found to be critical for mediating cell growth inhibition in human teratocarcinoma cell line NTERA-2 (Sacher et al., 2005). Structural comparison between a pTF.CREG1 homology model created using Phyre2 web portal (Kelley et al., 2015) and the human CREG1 monomer (PDB ID: 1XHN; Sacher et al., 2005) supports our hypothesis that pTF.CREG1 also contains a loop preventing FMN binding (Figure 6—figure supplement 2B). Two arginine amino acid residues involved in binding the terminal FMN phosphate in S. cerevisiae pyridoxine 5′-phosphate oxidase (UniProt ID: P38075, PDB ID: 1CI0) are mutated to aspartic acid and methionine amino acid residues in human CREG1 (Sacher et al., 2005), and one of these two key arginine residues is absent from pTF.CREG1 (Figure 6—figure supplement 3A). Finally, we note that a putative homodimerization region in pTF.CREG1 homology model contains many hydrophobic amino acid residues and overall follows the compositional character of the human CREG1 dimerization interface suggesting pTF.CREG1 is acting as a dimer in vivo (Figure 6—figure supplement 3B).
Given the apparent structural similarity of the P. tricornutum pTF.CREG1 to the catalytically inactive human ortholog, we opted to test its activity, or lack thereof, to delineate it from the evolutionarily related FMN-binding enzymes. Various truncations of His6-tagged pTF.CREG1 were expressed in Escherichia coli and purified using cobalt, Co2+, immobilized metal affinity chromatography (IMAC) to test solubility. A construct encoding pTF.CREG1 without the 31 N-terminal amino acid residues and leading to a soluble protein was selected for larger-scale purification (Figure 6—figure supplement 4A). The Co2+-bound protein fraction (Figure 6—figure supplement 4B) exhibited no flavin reductase activity as measured by the oxidation of NADPH in the presence of potential flavin substrates flavin mononucleotide (FMN) and riboflavin (Figure 6—figure supplement 2C; Figure 6—figure supplement 2—source data 1).
In conclusion, our homology modeling analyses and enzymatic assays indicate that pTF.CREG1 does not bind FMN and does not reduce flavins in an NADPH-dependent manner, respectively, both in agreement with our current knowledge about human CREG1. As such, pTF.CREG1 is an enticing candidate for further study as a key participant in the endocytic pTF receptor-mediated iron uptake pathway in diatoms. We put forth a hypothesis for its function in the Discussion.
pTF.CatCh1 is a heterokont-specific protein with conserved CX(X)C motifs
In contrast to pTF.CREG1, pTF.CatCh1 homologs were identified only in diatoms and other single-celled heterokonts (Figure 7A; Figure 7—source data 1). pTF.CatCh1 is a ~37.4 kDa (349 AAs) large protein with a predicted signal peptide (AAs 1–26) ending in ASA, two N-glycosylation NX(S/T) sequons (i.e. asparagines in the identified NST and NLS motifs may be glycosylated), and 13 predicted O-glycosylation sites (data not shown; Rao and Bernd, 2010). This protein is paralogous to the iron starvation induced protein 2b (ISIP2b/Phatr3_J54987) whose gene lies adjacently to pTF.CatCh1 in the same chromosome 20 gene cluster (Figure 5A). They share three CXC and one CXXC motifs typical of metal-binding, redox-active, and iron–sulfur cluster-containing proteins (Figure 7B; Blaby-Haas et al., 2014; Fomenko and Gladyshev, 2003; Poole, 2015; Przybyla-Toscano et al., 2018). One of the CXC motifs and the CXXC motif are conserved across all but two homologs in our phylogenetic analysis (Figure 7B; Figure 7—figure supplement 1A), and the latter is predicted to be located in one of the two disordered C-terminal regions in pTF.CatCh1 (AAs 196–292, AAs 321–349; Figure 7—figure supplement 1B). With additional 12 cysteines conserved between pTF.CatCh1 and ISIP2b, pTF.CatCh1 contains a total of 20 cysteine amino acid residues (all between AAs 34 and 226). Six of these 20 are 100% conserved across all homologs (Figure 7B). The predicted transmembrane domain (AAs 296–315) is flanked by a polyserine ([Ser]6) and a short arginine-rich stretch (RKL[R]3). Flexible polyserine linkers are found in modular, multidomain proteins (Uversky, 2015a), and positively charged amino acid residue tracts in the vicinity of transmembrane regions were shown to be orientation determinants in outer chloroplast membrane proteins (May and Soll, 1998). These observations are synthesized in a domain organization schematic shown in Figure 7—figure supplement 1C. We elaborate on possible roles of the conserved CX(X)C motifs, as well as other identified motifs and domains in pTF.CatCh1 in the Discussion.
pTF.ap1 homologs are present in marine microeukaryotes including diatoms
Phylogenetic characterization was performed with pTF.ap1 as well and numerous homologs with at least five highly conserved motifs were identified in diatoms and other marine microeukaryotes, including chlorophytes, cryptophytes, and haptophytes (Figure 7—figure supplement 2A–B; Figure 7—figure supplement 2—source data 1). These pTF.ap1 homologs contain five 100% conserved cysteine amino acid residues, two of which are shown in Figure 7—figure supplement 2B, as well as 100% conserved tyrosine, histidine, and aspartic acid residues, indicating this protein may also be redox-active and/or have a metal ion processing role (Dokmanić et al., 2008).
Discussion
In this study, we used the iron receptor protein phytotransferrin (pTF) to implement APEX2-mediated proximity-dependent proteomic mapping in the model marine diatom Phaeodactylum tricornutum. The resulting data provide several advances regarding the function of pTF and its putative protein interaction partners downstream of iron binding at the cell surface.
pTF (ISIP2a; iron starvation induced protein 2a) was identified in P. tricornutum as a marker for iron (Fe) limitation (Allen et al., 2008). It was subsequently shown to be involved in Fe acquisition (Morrissey et al., 2015) and to be widespread in marine and freshwater microeukaryotic phytoplankton communities (Bertrand et al., 2015; McQuaid et al., 2018; Tara Oceans Coordinators et al., 2018). In a breakthrough study, McQuaid et al., 2018, demonstrated that ISIP2a is a type of transferrin, phytotransferrin (pTF), that is crucial for acquisition of dissolved labile ferric, Fe3+, iron, a critical micronutrient for cellular biochemistry. McQuaid et al., 2018, showed that pTF is essential for high-affinity iron uptake in P. tricornutum, and that carbonate, CO32-, and Fe3+ interact synergistically to control the iron uptake rate. The finding that carbonate anions are required for the activity of a key diatom iron transport system is noteworthy and suggests that ocean acidification might inhibit microbial iron uptake. The occurrence of transferrin in marine diatoms and other single-celled marine and freshwater phytoplankton raises significant questions concerning its evolutionary origin and cellular role.
McQuaid et al., 2018 observed pTF in intracellular puncta, showed that the clathrin-mediated endocytosis inhibitor Pitstop 2 reduces iron (Fe) uptake rates, visualized vesicles after adding iron to an iron-limited WT P. tricornutum culture stained with the membrane dye FM 1–43, but our MDY-64 labeling data provide the first direct evidence that diatom pTF is indeed associated with intracellular membranous compartments. Further work is needed to elucidate pTF dynamic upon cellular internalization—especially after iron addition to iron-deplete cells which mimicks oceanic iron repletion events—and its trafficking in and potential exit from the endomembrane system in relation to cell cycle phases and light-dark cycles (Mayle et al., 2012; Naslavsky and Caplan, 2018; Smith et al., 2016).
While APEX2 was expected to be active in P. tricornutum given previous heterologous expression of functional horseradish peroxidase (HRP) in the model marine diatom Thalassiosira pseudonana (Sheppard et al., 2012), it was unclear whether endogenous peroxidases would contribute to the background signal in various APEX2 assays. In our Amplex UltraRed assay experiments, an order of magnitude higher signal-to-noise ratio (i.e. resurofin signal in pTF-APEX2 expressing versus WT P. tricornutum cells) was detected when live cells were reacted on ice as opposed to on room temperature, suggesting endogenous P. tricornutum peroxidases, but not exogenous APEX2, are largely inhibited on ice. This is notable as it indicates that performing APEX2 assays on ice may represent a generalizable experimental strategy for mesophilic microeukaryotic phytoplankton model systems.
Identifying Amplex UltraRed assay conditions which minimized WT background was promising, but not a guarantee that some endogenous enzymes would not be able to process other APEX2 substrates, such as those crucial for electron microscopy (3,3’-diaminobenzidine/DAB) and proximity proteomics (biotin-phenol). Indeed, mitochondrial background signal was detected in our electron microscopy experiments both in WT and pTF-APEX2 expressing cells, and our P. tricornutum proteome analysis revealed candidate endogenous APEX2-like peroxidases that could explain this result. Nevertheless, C-terminal tagging of pTF with APEX2 resulted in the expected cell surface and what appeared to be vesicular fusion protein localization away from mitochondria indicating no spatial overlap between endogenous APEX2-like enzymes and pTF-APEX2. Inconsistent size and shape of vesicle-like structures observed in our electron micrographs were conceivably due to a combination of arbitrary cell orientation inside agar blocks (where diatom cells were embedded) and the exact z position of ultramicrotome cut sites. On the whole, APEX2-based imaging is complementary to super-resolution microscopy in diatoms and represents a basis for electron tomography applications (Gröger et al., 2016; Sengupta et al., 2019). Some off-target biotinylation in our proximity labeling experiments was likely due to background mitochondrial peroxidase activity as indicated by the presence of mitochondrial proteins in our mass spectrometry data. Their APEX2/WT ratios close to one give us more confidence that those with elevated ratios above 1.5 were enriched due to APEX2 activity. Many of the latter have functional annotations for proteins we might expect to play a role in endosomal iron assimilation, and many have previously been shown to have iron-sensitive transcripts, both of which further support this notion. Future efforts to express pTF tagged with N-terminal APEX2, which should lead to intravesicular proximity labeling according to our pTF orientation prediction, optimizing conditions for APEX2 reactivation in iron-deplete conditions with hemin chloride, and more cell wall-permeant APEX2 probes (Li et al., 2020) will allow for identification of additional candidate proteins. Considering iron metabolism genes in P. tricornutum, including pTF, are highly expressed at night, proteomic maps with improved spatial and temporal resolution may emerge from proximity-dependent proteomic mapping experiments conducted with synchronized cell lines sampled at different timepoints throughout diel cycle (Huysman et al., 2014; Lundberg and Borner, 2019; Smith et al., 2016). The successful implementation of TurboID—an engineered biotin ligase for proximity proteomics (Branon et al., 2018)—in model plants Arabidopsis thaliana (Kim et al., 2019; Mair et al., 2019) and Nicotiana benthamiana (Zhang et al., 2019) suggests its adoption in existing and emerging single-celled marine phytoplankton model systems, including diatoms, is likely (Faktorová et al., 2020; Falciatore et al., 2020). TurboID is an iron-independent protein and could thus be particuarly valuable to further advance iron metabolism studies in P. tricornutum. These proposed improvements and alternative experimental approaches should cross-validate the proteins identified in this study as well as illuminate new ones.
Functionally uncharacterized gene cluster on chromosome 20 we zoomed in on is transcriptionally upregulated in iron-deplete (Allen et al., 2008; Smith et al., 2016) and silicic-acid-replete (Sapriel et al., 2009) conditions corroborating a known link between iron and silicon metabolism in diatoms (Brzezinski et al., 2015; Durkin et al., 2012; Hutchins and Bruland, 1998; Leynaert et al., 2004). Physically clustered functionally related genes are a hallmark of prokaryotic biology and operon-like genomic elements have been described in fungi and plants (Nützmann et al., 2018; Osbourn and Field, 2009); however, identifying and characterizing them in marine diatom genomes is in its relative infancy. For example, genes encoding light harvesting fucoxanthin-chlorophyll-binding proteins (FCPs) are co-located on chromosome 2 in P. tricornutum (Bhaya and Grossman, 1993), and a compact ~7 kbp gene cluster encoding enzymes for domoic acid neurotoxin production was characterized in a cosmopolitan diatom Pseudo-nitzschia multiseries (Brunson et al., 2018). Given that diatom genomes have likely been extensively shaped by horizontal gene transfer (Bowler et al., 2008; Diner et al., 2017), it is interesting to note that some yeast species acquired iron acquisition-enabling operon-like genomic elements via horizontal operon transfer (HOT) from an ancient bacterium (Kominek et al., 2019). Our phylogenetic analyses suggest all three proteins—pTF.CREG1, pTF.CatCh1, pTF.ap1—encoded by this locus are present in diatoms and other (marine) microeukaryotes, while pTF.CREG1 homologs are found across the tree of life. They, alongside the rest of our proteomic data, paint an emerging view of intracellular ferric, Fe3+, iron processing events downstream of its phytotransferrin-mediated binding at the cell surface.
In humans, endosome acidification causes structural rearrangements in the transferrin protein and protonation of synergistically coordinated carbonate anion leading to ferric iron release (Cheng et al., 2004). This event is concurrent with or followed by endosomal reductase STEAP3-catalyzed ferric iron reduction and ferrous, Fe2+, iron export through divalent metal transporter DMT1/NRAMP2 or Zrt- and Irt-like protein (ZIP) 14 (Bogdan et al., 2016). Two putative proton, H+, pumps were present in our mass spectrometry data: Phatr3_J23497—predicted as a P-type ATPase—and Phatr3_J27923—a vacuolar-type H+-ATPase subunit A. Multisubunit V-ATPases regulate pH homeostasis in virtually all eukaryotes and are known to function as endosome acidifying protein complexes (Finbow and Harrison, 1997; Maxson and Grinstein, 2014), which makes Phatr3_J27923 a candidate protein involved in acidification of phytotransferrin-rich endosomes in P. tricornutum. Notably, V-ATPase transcripts were overrepresented in diatoms following iron enrichment, possibly to account for an increase in the number of endosomes associated with pTF-driven iron acquisition (Marchetti et al., 2012), and V-ATPase was shown to be involved in diatom silica deposition vesicle (SDV) acidification (Hildebrand et al., 2018; Yee et al., 2020). Endosome acidification, perhaps via V-ATPase, may facilitate phytotransferrin-bound ferric, Fe3+, iron release, possibly for processing by an iron reducing enzyme similar to the human metalloreductase STEAP3.
Human cellular repressor of E1A-stimulated genes 1 (CREG1) is a glycoprotein involved in embryonic development, growth, differentiation, and senescence as part of the endosomal-lysosomal system (Ghobrial et al., 2018; Schähs et al., 2008). Contrary to evolutionarily related flavin mononucleotide (FMN)-binding split barrel fold oxidases and oxidoreductases such as pyridoxine 5′-phosphate (PNP) oxidase, an enzyme involved in vitamin B6 metabolism (Musayev et al., 2003), CREG1 forms a homodimer unable to bind FMN in vitro suggesting it does not act as an enzyme in vivo (Sacher et al., 2005). The lack of FMN in human CREG1 is due to a short loop protruding into the FMN-binding pocket at the dimer interface, and aspartic acid and methionine amino acid residues replacing the corresponding two yeast PNP oxidase (PDB ID: 1CI0) arginines responsible for binding the terminal FMN phosphate (Sacher et al., 2005).
CREG1-like protein pTF.CREG1 was proposed to be central to sustained growth of iron-limited P. tricornutum cells (Allen et al., 2008) or to serve as a cell surface iron receptor (Lommer et al., 2012). CREG transcripts were more common in diatom transcriptomes from the Southern Ocean as opposed to non-Southern Ocean regions indicating these proteins are important in coping with iron limitation in this major high-nutrient, low-chlorophyll (HNLC) zone (Moreno et al., 2018). This may also explain why more diatom homologs are not present in our phylogenetic tree as transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP) database used in our phylogenetic analysis were not obtained from iron-limited cultures (Caron et al., 2017; Keeling et al., 2014). Our protein feature analysis revealed a signal peptide and multiple putative glycosylation sites which suggests that pTF.CREG1 enters the secretory pathway (Bard and Chia, 2016). While the nature of pTF.CREG1 recruitment to and colocalization with pTF-rich endosomes in the vicinity of chloroplast periphery—as suggested by our fluorescence microscopy results—is unclear, it is possible that pTF.CREG1-containing vesicles intersect with pTF-rich endosomes en route through the secretory pathway. Vesicular trafficking is highly dynamic and bidirectional (Naslavsky and Caplan, 2018; Progida and Bakke, 2016), and early endosomal traffic consists of sequential vesicle fusion events (Brenner, 2012), so this is at least in principle plausible. Punctate pTF.CREG1 colocalization with pTF suggests that this protein acts in an endosomal-lysosomal system, and this localization pattern is consistent with the observation that proteins with low confidence ASAFind chloroplast prediction (such as pTF.CREG1) can be associated with ‘blob-like structures’ (BLBs) adjacent to, but not inside, the chloroplast (Gruber et al., 2015; Kilian and Kroth, 2005). Given that our transgenic colocalization P. tricornutum cell lines were neither synchronized nor grown in iron-deplete conditions, it is plausible the protein also acts closer to the cell surface. Likely, but not confirmed, APEX2 orientation away from the cell surface and endosome interior in our proximity labeling experiments, and the reported inability of APEX2-generated phenoxyl radicals to penetrate endomembranes (Rhee et al., 2013), would suggest that pTF.CREG1—at some point during its recruitment to pTF-rich endosomes—has to be exposed to the cytosol to promote efficient biotinylation. Alternatively, it is impossible to rule out the possibility that lipid composition of P. tricornutum endomembranes is permissive to phenoxyl radical diffusion (Zulu et al., 2018).
Domain annotation in pTF.CREG1 (UniProt ID: B7G9B3) encouraged us to perform amino acid sequence and structural alignments with human CREG1 (UniProt ID: O75629, PDB ID: 1XHN). These analyses revealed that pTF.CREG1 lacks at least one arginine amino acid residue required for terminal FMN phosphate binding and likely contains an FMN binding pocket-occluding loop, two key features delineating human CREG1 from the yeast PNP oxidase (Sacher et al., 2005). Our hypothesis that pTF.CREG1 is unable to coordinate FMN was initially supported by our protein purifications which were not yellow as expected from a typical flavoprotein (Theorell, 1935), and was subsequently strengthened by their catalytic inactivity upon supplementation with biologically relevant flavin cofactors. The demonstrated lack of enzymatic activity in pTF.CREG1 and localization data from in vivo fluorescent reporter experiments represent an important extension of the human CREG1 literature to the single-celled eukaryotic phytoplankton field.
Secreted glycosylated human CREG1 was shown to delay G1/S cell cycle phase transition via direct interaction with mannose 6-phosphate/insulin-like growth factor II receptor (M6P/IGF2R), a receptor internalized via clathrin-mediated endocytosis (Di Bacco and Gill, 2003; Ghosh et al., 2003; Han et al., 2009). This cellular growth inhibition was demonstrated to, at least partially, depend on a tetrapeptide motif within an FMN binding pocket-occluding loop (Sacher et al., 2005). In conjunction with our amino acid sequence alignments and structural prediction which point to the conservation of these two features in pTF.CREG1, these studies lend additional support to our data suggesting colocalization and potential direct physical association between phytotransferrin and pTF.CREG1 is non-enzymatic in nature. Internalization of the human transferrin-transferrin receptor complex is not a constitutive, but rather a regulated kinase-activity-dependent process (Cao et al., 2016). We currently hypothesize that pTF.CREG1 may serve as a ‘checkpoint’ to keep phytotransferrin endocytosis at a low level in iron-replete conditions, but when iron is scarce instead positively regulate this pathway. Indeed, transcriptomics experiments show that a second transcriptional peak for pTF.CREG1 occurs in iron-deplete medium approximately 8 hr prior to the pTF peak, a transcriptional dynamic which would be consistent with our hypothesis (Figure 5—figure supplement 1; Smith et al., 2016). Given that glycosylation in human CREG1 appears to be important for interaction with its cognate receptor (Han et al., 2011; Sacher et al., 2005; Schähs et al., 2008), posttranslational modifications, perhaps glycosylation(s), of pTF.CREG1 may further drive this regulatory activity. The presence of Phatr3_J10972, an iron-insensitive pTF.CREG1 paralog, in the P. tricornutum genome suggests that CREG-like proteins are implicated in diverse endocytosis pathways in single-celled eukaryotic phytoplankton. However, the association between iron-sensitive pTF.CREG1 and phytotransferrin in P. tricornutum implies a specific pTF.CREG1 function in non-reductive iron uptake. It will be very interesting to define the exact pTF.CREG1 role in further cellular and mechanistic studies and compare it to that of its homologs in multicellular organisms. Conceivably, an ancient protein like pTF.CREG1 could have remained at the core of linking cell cycle progression, cellular growth, and nutrient homeostasis across the tree of life.
In contrast to pTF.CREG1, no functional annotation could be assigned to and no structural homology models were generated for pTF.CatCh1 in the Pfam database (El-Gebali et al., 2019) and by Phyre2, respectively. This protein has a predicted N-terminal signal peptide and numerous N- and O-glycosylation sites suggesting it enters and trafficks along the secretory pathway. The outermost complex plastid membrane in diatoms corresponds to the chloroplast endoplasmic reticulum (cER) and is continuous with the nuclear envelope (Gibbs, 1981; Prihoda et al., 2012). Nuclear-encoded P. tricornutum proteins destined for the periplastidial compartment (vastly reduced endosymbiont cytoplasm between the inner two and the outer two complex plastid membranes) are first cotranslationally imported into the cER lumen (Grosche et al., 2014; Peschke et al., 2013). Our microscopy results show that pTF.CatCh1 is localized on the chloroplast margin, perhaps anchored to the outer chloroplast (cER) membrane with its C-terminal transmembrane domain, thus raising a possibility that pTF.CatCh1 gets to its final destination via cotranslational membrane insertion.
The cysteine-rich domain contains one CXC and one CXXC motifs which are nearly invariably present in pTF.CatCh1 homologs. Two cysteines spaced by two non-cysteine amino acids (the CXXC motif) are commonly found in metal-binding proteins, including those that bind iron (Fomenko et al., 2008; Gladyshev et al., 2004). This motif in pTF.CatCh1 is preceded by an aspartic acid residue resulting in DCXXC which is known to be particularly permissive to binding divalent metal ions (Figure 7—figure supplement 1A; Banci et al., 2002; Banci et al., 2006). Histidine, tyrosine, methionine, aspartic acid, and glutamic acid residues can also be involved in iron cation coordination, but none of these are conserved across pTF.CatCh1 homologs (Figure 7—source data 1; Dokmanić et al., 2008). Given the nearly completely conserved C(X)XC motifs and additional six 100% conserved cysteine residues, it is conceivable they are critical to pTF.CatCh1 activity.
Based on the fluorescence microscopy data, the conserved cysteine amino acid residues and cysteine-containing motifs, and given that disordered protein domains—two are predicted to flank the pTF.CatCh1 transmembrane domain—are mediators of protein-protein interactions and play an important role in the formation of protein complexes (Uversky, 2015b; Uversky, 2016), we are considering two hypotheses for in vivo pTF.CatCh1 function.
First, pTF.CatCh1 may serve as a metallochaperone that binds ferrous, Fe2+, iron downstream of the endosomal ferric, Fe3+, iron reduction step, and directs it to the chloroplast interior. Metallochaperones are intracellular proteins responsible for protecting and guiding metal ions on their way to correct metalloenzyme sinks (O'Halloran and Culotta, 2000). They engage in protein-protein interactions and help prevent metal-induced toxicity caused by free radical-generating Fenton reactions which are especially common with iron and copper (Philpott and Jadhav, 2019; Rosenzweig, 2002; Valko et al., 2005).
Second, pTF.CatCh1 may recruit pTF to the chloroplast periphery via one of its two C-terminal disordered regions flanking the predicted transmembrane domain, and thus serve as a docking site for iron-laden pTF. It could achieve this via the ‘fly-casting mechanism’ which involves binding-induced protein folding (Shoemaker et al., 2000). The observed colocalization of pTF.CatCh1 with pTF may indicate that Fe3+ can be offloaded directly from pTF to pTF.CatCh1 which would imply an additional, non-endosomal, intracellular Fe3+ reduction step. Notably, a chloroplast-associated Fe3+ reduction pathway mediated by Fe3+ chelate reductase FRO7 is present in Arabidopsis thaliana (Jeong et al., 2008). Abberant cellular localization of heterologously expressed fluorescent protein fusions is a widely observed and acknowledged phenomenon (Jensen, 2012), but we note that pTF localization on the chloroplast margin insdistinguishable from that in this study was reported by McQuaid et al., 2018 suggesting it is biologically relevant. Similarly to pTF.CREG1, pTF.CatCh1 exhibits a second transcriptional peak in iron-deplete medium approx. 16 hr prior to the pTF peak (Figure 5—figure supplement 1). It would make sense to synthesize a protein critical for iron-laden pTF recruitment to the chloroplast margin prior to significant increase in basal pTF endocytosis.
All in all, our data suggest pTF.CatCh1 is a chloroplast margin-localized, pTF-associated, P. tricornutum protein with a possible role in chloroplast iron homeostasis. The two hypotheses are put forth here to help frame its further genetic, cellular, biochemical, and structural dissection.
Chloroplasts are the most iron-rich system in plant cells (López-Millán et al., 2016). Iron—in particular in the form of iron–sulfur (Fe–S) clusters—serves as an essential cofactor for the photosynthetic electron transfer chain, chlorophyll biosynthesis, and chloroplast protein import (Balk and Schaedler, 2014; Marchand et al., 2018; Soll and Schleiff, 2004). A complex vesicle-based system possibly in perpetual dynamic exchange with cytosolic vesicles is present in chloroplasts (Hertle et al., 2020; Lindquist and Aronsson, 2018), an endoplasmic reticulum to Golgi to chloroplast protein trafficking pathway exists (Radhamony and Theg, 2006; Villarejo et al., 2005), and additional hypotheses for such specialized vesicle-mediated plastid targeting of proteins are established (Baslam et al., 2016). Therefore, a direct link between endocytic internalization of phytotransferrin-bound iron and its offloading to or adjacent to chloroplasts as suggested by our results seems plausible. ‘Kiss and run’ mechanism that allows erythroid cell mitochondria to access iron via direct interaction with transferrin-rich endosomes further underscores this idea (Hamdi et al., 2016). Analysis of our mass spectrometry (MS) hits that were not co-expressed with pTF-mCherry reveals three additional proteins—ISIP1 (Phatr3_J55031), FBP1 (Phatr3_J46929), and FBAC5 (Phatr3_J41423)—that support the idea of a direct cell surface-to-chloroplast iron trafficking pathway in Phaeodactylum tricornutum.
ISIP1 (iron starvation induced protein 1; Phatr3_J55031) gene expression in P. tricornutum and other marine microeukaryotes is induced by iron limitation (Allen et al., 2008; Kazamia et al., 2018; Marchetti et al., 2012). P. tricornutum ISIP1—the second most enriched protein in our MS dataset—has a role in siderophore-bound iron uptake and colocalizes with phytotransferrin (pTF) adjacently to the chloroplast in an ‘iron-processing compartment’ (Kazamia et al., 2019; Kazamia et al., 2018). FBP1 (Phatr3_J46929)—iron-sensitive siderophore ferrichrome-binding protein with punctate intracellular localization and central to the carbonate-independent iron acquisition pathway in P. tricornutum (Allen et al., 2008; Coale et al., 2019)—was another highly enriched protein in our dataset indicating the internalization pathways for inorganic (i.e. phytotransferrin-bound) and organic (i.e. siderophore-bound) iron in P. tricornutum may be physically coupled intracellularly and directed to iron-accumulating vesicles close to the chloroplast. Metal homeostasis by such discrete membranous compartments is widespread in nature (Blaby-Haas and Merchant, 2014). Copper- and zinc-storing compartments—cuprosomes and zincosomes, respectively—are present in the model chlorophyte microalga Chlamydomonas reinhardtii (Aron et al., 2015; Hong-Hermesdorf et al., 2014; Merchant, 2019), and lipid-bound iron-accumulating ferrosome compartments exist in diverse bacteria and archaea (Grant and Komeili, 2020).
Finally, the most enriched protein in our MS dataset was FBAC5 (Phatr3_J41423)—iron-sensitive iron-independent class I fructose-biophosphate aldolase—which was shown to be localized to the P. tricornutum pyrenoid—proteinaceous RuBisCO-containing subchloroplastic organelle—and proposed to link Calvin-Benson-Bassham cycle activity with the CO2 concentrating mechanism (Allen et al., 2012; Matsuda et al., 2017). FBAC5 is one of the most reliable markers for diatom iron stress as demonstrated by laboratory culture (Cohen et al., 2018; Lommer et al., 2012), microcosm incubation (Cohen et al., 2017a; Cohen et al., 2017b), and field studies (Bertrand et al., 2015); however, its exact role has remained elusive. Fructose-biphosphate aldolases (FBAs) are glycolytic enzymes known for their ability to moonlight (Boukouris et al., 2016), a term for proteins with multiple physiologically relevant functions encoded by one polypeptide chain (Jeffery, 2018). Moonlighting FBAs have roles in actin cytoskeleton dynamic, cellular growth arrest, and apoptosis (in human cancer cell lines; Gizak et al., 2019), transcriptional regulation (in yeast; Cieśla et al., 2014), and endosome acidification via direct interaction with vacuolar-type ATPase (in a mouse cell line; Merkulova et al., 2011). These reported non-catalytic FBA activities suggest that FBAC5 acting outside of the P. tricornutum chloroplast, perhaps in the endosomal-lysosomal system as our data would imply, is possible.
Our initial characterization of pTF.CREG1 and pTF.CatCh1, together with the discussion on the four other mass spectrometry hits, enable proposition of a model connecting cell surface ferric iron binding, internalization, and intracellular trafficking in Phaeodactylum tricornutum with many outstanding questions ripe for future investigation (Figure 8). Importantly, it remains to be determined which protein is responsible for intracellular ferric, Fe3+, iron reduction in P. tricornutum, activity that would be consistent with the non-reductive cell surface ferric, Fe3+, iron uptake proposed for phytotransferrin (pTF) (McQuaid et al., 2018; Morrissey et al., 2015). Perhaps notably, neither of the five proteins with clear ferric reductase annotation (FRE1–5, Phatr3 IDs: J54486, J46928, J54940, J54409, J54982) nor additional five hypothetical ferric reductases (Coale et al., 2019; Smith et al., 2016) were present in our mass spectrometry data. Further proximity-dependent proteomic mapping experiments will be essential for determining the full set of proteins in the non-reductive iron acquisition pathway in P. tricornutum as enzymatic Fe3+ reduction may happen later in the pathway and/or in a separate subcellular compartment.
In conclusion, the identified proteins and the associated model we introduce in this work represent an important advance for the marine phytoplankton field, given that diatoms are central to the study of cellular iron homeostasis and its link to fluctuating iron levels in the ocean. The presented insights are vital as our understanding of intracellular iron trafficking is relatively incomplete even in significantly more established model systems such as mammalian cells, and our knowledge of chloroplast iron uptake is in its infancy (Blaby-Haas and Merchant, 2012; Philpott and Jadhav, 2019). We anticipate further molecular dissection of our proteomics dataset will not only provide details on the transformation dynamic of acquired ferric iron and its cellular sinks but will thereby also illuminate the extent to which phytotransferrin pathway mirrors the metazoan transferrin cycle.
Implementing the use of APEX2 in a model marine diatom is itself a significant and timely advance as the field is expanding and moving from broad -omics-based surveys to detailed molecular studies (Faktorová et al., 2020; Falciatore et al., 2020). APEX2 and related molecular tools hold great promise to dissect additional key subcellular compartments in diatoms such as pyrenoids and silica deposition vesicles (Figure 8—figure supplement 1; Barrett et al., 2021; Hildebrand et al., 2018; Matsuda et al., 2017; Wang and Jonikas, 2020), two exciting future avenues for these single-celled heterokonts uniquely positioned at the nexus of evolutionary cell biology and blue biotechnology.
Materials and methods
Key resources table
Request a detailed protocolThe quantitative mass spectrometry proteomics data used to generate Figure 4—source data 1 and Figure 4—figure supplement 1—source data 1 have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2019) with the dataset identifier PXD018022.
The amino acid sequence alignment files Figure 6—source data 1, Figure 7—source data 1, and Figure 7—figure supplement 2—source data 1 can be read and interrogated with AlignmentViewer available at https://alignmentviewer.org/.
Vector cloning
pJT_NR_pTF-AP2
Request a detailed protocolGibson Assembly (Gibson et al., 2009) was performed with three DNA fragments: (1) linearized pPtPBR1 episome backbone (Addgene, Cambridge, MA, plasmid pPtPBR1, Catalog #80388) opened ~280 bp downstream of the tetracycline resistance gene, (2) amplicon with nitrate reductase gene (NR) promoter, 5’ NR untranslated region (UTR), and pTF, and (3) gBlocks Gene Fragment (Integrated DNA Technologies (IDT), Coralville, IA) with linker sequence, codon optimized APEX2 (using IDT Codon Optimization Tool and selecting ‘Thalassiosira pseudonana’ from the ‘Organism’ drop-down menu), 3’ NR UTR, and NR terminator.
pJT_native_pTF-mCherry
Request a detailed protocolGibson Assembly was performed with three DNA amplicons: (1) pPtPBR1 episome backbone split into two fragments and (2) expression cassette including pTF, mCherry, and pTF promoter and terminator.
pJT_pTF-mCherry_MS hit-EYFP
Request a detailed protocolpJT_native_pTF-mCherry was split into two fragments via PCR stitching (keeping ampicillin resistance gene split). Genes corresponding to MS hits were amplified using cDNA from iron-starved WT P. tricornutum cells and assembled with EYFP and Phatr3_J23658 (flavodoxin-encoding gene) promoter and terminator through two PCR stitching rounds into a single fragment. Gibson Assembly was then performed to combine all three final amplicons. Detailed assembly scheme is presented in Figure 5—figure supplement 2A.
pJT_Δ31_pTF.CREG1-His6
Request a detailed protocolCloning of gene fragments into the E. coli protein expression vector PtpBAD-CTHF was performed as described previously (Brunson et al., 2018). Briefly, PtpBAD-CTHF was linearized by digestion with XhoI (New England Biolabs (NEB), Ipswich, MA, Catalog #R0146S) and the resulting DNA was column purified. Δ18_pTF.CREG1 gene was obtained by PCR from P. tricornutum gDNA with PrimeSTAR GXL DNA Polymerase (Takara Bio, Kusatsu, Japan, Catalog #R050B) and primer set JT31/JT32 to incorporate the appropriate Gibson Assembly overhangs and remove the predicted N-terminal signal peptide (amino acid residues 1–18). Insertion of truncated pTF.CREG1 into linearized PtpBAD-CTHF was performed using Gibson Assembly Master Mix (NEB, Catalog #E2611S), 1 µL of the Gibson Assembly reaction mixture was transformed via heat shock into chemically competent NEB 5-alpha cells (NEB, Catalog #C2988J), and cells were incubated on lysogeny broth with 10 μg/mL tetracycline (LB-Tet10) 1% agar plates overnight at 37°C. Transformants were screened by colony PCR using the primer set JT37/JT38 and Sapphire polymerase (Takara Bio), and positive clones were selected for outgrowth. Isolated plasmids were sequence-validated by Sanger sequencing (Eurofins, Luxembourg, Luxembourg). A sequence-validated clone was designated as PtpBAD-Δ18_pTF.CREG1-CTHF and transformed into chemically competent BL21 E. coli cells (NEB, Catalog #C2530H) which were spread on LB-Tet10 1% agar plates. The resulting transformants were used for subsequent Δ18_pTF.CREG1 expression experiments. Following unsuccessful expression testing of the Δ18 N-terminal truncation construct, an additional set of genes encoding pTF.CREG1 with N-terminal truncations (Δ26, Δ31, Δ39, Δ43) was generated by PCR using PrimeStar GXL DNA Polymerase (Takara Bio), primer sets JT33/JT32, JT34/JT32, JT35/JT32, JT36/JT32, and PtpBAD-Δ18_pTF.CREG1-CTHF as a template. Assembly, colony PCR screening, plasmid isolation, sequencing, and transformation into chemically competent BL21 E. coli cells of all additional expression vectors, including PtpBAD-Δ31_pTF.CREG1-CTHF (=pJT_Δ31_pTF.CREG1-His6), was performed as above. An additional construct encoding pTF.CREG1 with an N-terminal His6 and Δ39 N-terminal truncation was built (PtpBAD-NTH-Δ39_pTF.CREG1) using the vector NTH-PtpBAD (constructed similarly to PtpBAD-CTHF; see Brunson et al., 2018 and Savitsky et al., 2010), but was not pursued beyond initial expression testing.
pTF (Phatr3_J54465) in all episomes was in its native form (3 exons, 2 introns). Molecular cloning primers are listed in Supplementary file 1—Table S3. Vector details, further amplicon information, and fusion protein sequences are catalogued in Supplementary file 1—Table S4.
Related resource
Request a detailed protocolTurnšek J, Gholami P. 2017. Guidelines for highly efficient construction of diatom episomes using Gibson Assembly. protocols.io. DOI: dx.doi.org/10.17504/protocols.io.jy7cpzn.
Diatom culturing, conjugation, and genotyping
Culturing
Request a detailed protocolSequenced Phaeodactylum tricornutum strain CCMP632 (synonymous to CCMP2561 and CCAP 1055/1; National Center for Marine Algae and Microbiota (NCMA), East Boothbay, ME) was used throughout the study and grown at 18°C, 300 μmol quanta m−2 s−1, and a 10 hr:14 hr dark:light cycle in biotin-free L1 medium prepared by mixing 1 L Aquil salts, 2 mL nitrate and phosphate (NP) nutrient stock, 1 mL trace metal stock, and 1 mL thiamine hydrochloride and cyanocobalamin (TC) stock unless otherwise noted. Preparation of Aquil salts: 0.5 L anhydrous salts (0.5 L Milli-Q, 24.5 g NaCl, 4.09 g Na2SO4, 0.7 g KCl, 0.2 g NaHCO3, 0.1 g KBr, 900 μL 33.3 mg/mL H3BO3 stock, 300 μL 10 mg/mL NaF stock) and 0.5 L hydrous salts (0.5 L Milli-Q, 11.1 g MgCl2 × 6H2O, 1.54 g CaCl2 × 2H2O, 100 μL 170 mg/mL SrCl2 × 6H2O stock) were prepared separately, combined, filter sterilized (0.2 μm), and stored at room temperature. Preparation of NP nutrient stock: 37.5 g NaNO3 and 2.5 g NaH2PO4 were dissolved in 100 mL Milli-Q, filter sterilized (0.2 μm), and stored at 4°C. Preparation of trace metal stock (for 1 L 1000x stock): 3.15 g FeCl3 × 6H2O, 4.36 g Na2EDTA x 2H2O, 0.25 mL 9.8 g/L CuSO4 × 5H2O, 3.0 mL 6.3 g/L Na2MoO4 × 2H2O, 1.0 mL 22 g/L ZnSO4 × 7H2O, 1.0 mL 10 g/L CoCl2 × 6H2O, 1.0 mL 180 g/L MnCl2 × 4H2O, 1.0 mL 1.3 g/L H2SeO3, 1.0 mL 2.7 g/L NiSO4 × 6H2O, 1.0 mL 1.84 g/L Na3VO4, 1.0 mL 1.94 g/L K2CrO4, and Milli-Q up to 1 L were combined, filter sterilized (0.2 μm), and kept at 4°C. Preparation of TC stock: 20 mg thiamine hydrochloride and 0.1 mL 1 g/L cyanocobalamin stock were mixed in 100 mL Milli-Q. The resulting solution was stored at 4°C. ΔpTF P. tricornutum cell line (maintained in the Allen Lab; McQuaid et al., 2018) was additionally supplemented with 200 μg/mL nourseothricin (GoldBio, Saint Louis, MO, Catalog #N-500–1). All transconjugant P. tricornutum cell lines were supplemented with 50 or 100 μg/mL phleomycin (InvivoGen, San Diego, CA, Catalog #ant-ph-10p).
Conjugation
Request a detailed protocol(1) Bacterial donor preparation: Chemically competent pTA-Mob-containing TransforMax EPI300 E. coli cells (Strand et al., 2014; Lucigen, Middleton, WI, Catalog #EC300110) were transformed via heat shock with sequence-verified pPtPBR1 episomes. Transformants were selected on gentamycin-, carbenicillin-, and tetracycline-containing LB 1% agar plates. 3 mL overnight LB cultures supplemented with antibiotics were inoculated from glycerol stocks. (2) P. tricornutum preparation: ~2×108 P. tricornutum cells (in 200 μL) in mid- to late-exponential phase were spread on pre-dried (i.e. plates with lid half open and kept in the laminar flow hood for at least 90 min) ½ L1 1% agar plates with 5% LB and left growing for 1 or 2 days. Plates were additionally supplemented with 200 μg/mL nourseothricin for conjugation of ΔpTF P. tricornutum cells. (3) Conjugation: Overnight donor bacterial cultures were diluted 1:50 in 25 mL LB supplemented with antibiotics, grown at 37°C until OD600 0.8–1, spun down, resuspended in 150 µL Super Optimal broth with Catabolite repression (SOC) medium, and spread as evenly as possible on top of a P. tricornutum lawn. Plates with donor-P. tricornutum co-culture were first left in dark and 30°C for 90 min, then for 1 or 2 days at standard growth conditions. (4) Selection: Co-culture lawn was scraped off of plates with 1 mL fresh L1 medium, transferred to a microcentrifuge tube, and 200 μL spread on pre-dried (see above) ½ L1 1% agar plates with 50 or 100 μg/mL phleomycin. Plates were additionally supplemented with 200 μg/mL nourseothricin for conjugation of ΔpTF P. tricornutum cells. Porous adhesive tape was used to seal the plates and transconjugants emerged after ~10 days of incubation under standard growth conditions. Please see Karas et al., 2015 and Diner et al., 2016 for further description of diatom conjugation.
Genotyping
Request a detailed protocolCandidate transconjugant colonies were inoculated in 300 μL L1 medium supplemented with 50 or 100 μg/mL phleomycin (and 200 μg/mL nourseothricin in case of ΔpTF P. tricornutum cells) and typically grown for ~1 week. 0.5 μL liquid culture was then genotyped using either Phire Plant Direct PCR Master Mix (Thermo Fischer Scientific, Waltham, MA, Catalog #F160L) or Phire Plant Direct PCR Kit (Thermo Fischer Scientific, Catalog #F130WH). 200 μL of each genotype-positive cell line was passaged in 30 mL L1 medium supplemented with 50 or 100 μg/mL phleomycin (and 200 μg/mL nourseothricin in case of ΔpTF P. tricornutum cells).
Related resource
Request a detailed protocolTurnšek J. 2017. Simple and rapid genotyping of marine microeukaryotes. protocols.io. DOI: dx.doi.org/10.17504/protocols.io.jcdcis6.
RNA extraction and cDNA synthesis
Diatom culture
Request a detailed protocol100 mL WT P. tricornutum culture was grown in iron-deplete conditions (L1 medium with 7.5 nM total iron) for two weeks. Cells were then centrifuged, supernatants discarded, pellets flash frozen in liquid nitrogen, and stored at −80°C.
RNA extraction
Request a detailed protocolDirect-zol RNA Miniprep Plus RNA extraction kit was used (Zymo Research, Irvine, CA, Catalog #R2070). Briefly, cell pellets were resuspended in 800 μL Trizol, equal amount of 100% ethanol, and centrifuged in spin columns. Columns were then washed with 400 μL RNA Wash Buffer followed by on-column DNA digestion with 5 μL DNase I in 75 μL DNA Digestion Buffer for 15 min at RT, washing twice with 400 μL Direct-zol RNA PreWash, and once with 700 μL RNA Wash Buffer. RNA was eluted with 50 μL DNase/RNase-Free Water, RNA integrity number (RIN) evaluated with 2200 TapeStation (Agilent Technologies, Santa Clara, CA; measured RIN was 7.0), concentration estimated with Qubit 2.0 Fluorometer (Thermo Fischer Scientific; measured concentration was 27.6 ng/μL), and samples stored at −80°C.
cDNA synthesis
Request a detailed protocolcDNA synthesis kit from Thermo Fischer Scientific was used (SuperScript III First-Strand Synthesis System, Catalog #18080–051). Briefly, 1 μL total RNA was combined with 0.5 μL Oligo(dT)20 Primer, 0.5 μL 10 mM dNTP Mix, and 3 μL nuclease-free water followed by 5 min and 1 min incubation at 65°C and on ice, respectively. After 5 μL cDNA Synthesis Mix (1 μL 10x Reverse Transcription buffer, 2 μL 25 mM MgCl2, 1 μL 0.1 M DTT, 0.5 μL RNase OUT, 0.5 μL SuperScript III Reverse Transcriptase) was added, samples were incubated for 50 min and 5 min at 50°C and 85°C, respectively, then chilled on ice. Finally, 1 μL RNase H was added and sample kept for 20 min at 37°C. cDNA was stored at −20°C until use.
MDY-64 labeling and imaging
Request a detailed protocolTen mL of a pTF-mCherry expressing ΔpTF P. tricornutum cell line grown in iron-deplete conditions (L1 medium with 7.5 nM total iron) were spun down, supernatant discarded, and pellet resuspended in 50 μL phosphate-buffered saline (PBS) (pH 7.4). 0.5 μL 10 mM MDY-64 stock (in DMSO; Thermo Fischer Scientific, Catalog #Y7536) was added, cells incubated for 10 min at RT, pelleted, and resuspended in 50 μL fresh PBS (pH 7.4). Cell suspension was prepared for imaging as follows: 5 μL was placed between a 1.5 mm microscope slide and a cover slip (this setup applies to all imaging experiments in the study). Imaging conditions: Leica TCS SP5 confocal microscope (Leica Microsystems, Wetzlar, Germany), argon laser strength set to 30%, 458 nm laser line at 50% maximum strength, emission window set to 477–517 nm (for visualizing MDY-64; MDY-64 excitation and emission maxima are 451 nm and 497 nm, respectively), 514 nm laser line at 50% maximum strength, emission window set to 620–640 nm (for visualizing mCherry; mCherry excitation and emission maxima are 587 nm and 610 nm, respectively).
Protein expression analyses
pTF-APEX2 detection
Request a detailed protocolCell pellets from 8 mL mid- to late-exponential phase WT or transconjugant P. tricornutum cultures were resuspended in 150 μL cell lysis buffer (50 mM Tris-HCl, 200 mM NaCl, 1 mM DTT, 1 mM PMSF, pH 8.5) and sonicated for 5 min (30 s on, 1 min off) with Bioruptor UCD-200TM (Diagenode, Liège, Belgium). The resulting cell lysates were centrifuged, total protein content in supernatants measured with Bradford Assay Kit (Thermo Fischer Scientific, Catalog #23236), and insoluble fractions resuspended in 150 μL cell lysis buffer. One μg of each soluble protein sample and 0.5 μL of each resuspended insoluble protein fraction were resolved on a NuPage 4–12% Bis-Tris 1.5 mm gel (Thermo Fischer Scientific, Catalog #NP0335BOX), wet transferred to polyvinylidene difluoride (PVDF) membranes (Thermo Fischer Scientific, Catalog #LC2005), and visualized with WesternBreeze Chemiluminescent Kit, anti-rabbit (Thermo Fischer Scientific, Catalog #WB7106).
pTF-mCherry detection
Request a detailed protocolCell pellets from 400 μL mid- to late-exponential phase WT or transconjugant P. tricornutum cultures were resuspended in 50 μL cell lysis buffer (50 mM Tris-HCl, 200 mM NaCl, 1 mM DTT, 1 mM PMSF, pH 8.5) and sonicated 15 min (30 s on, 1 min off) with Bioruptor UCD-200TM. One μL of each whole cell lysate was resolved on a NuPage 4–12% Bis-Tris 1.5 mm gel, wet transferred to PVDF membranes, and visualized with WesternBreeze Chemiluminescent Kit, anti-rabbit.
pTF antibody details
Request a detailed protocolAmino acid residues 32–223 served as the immunogen (N-terminal pTF region just downstream of the signal peptide). The antibody was produced in a rabbit. 1:10,000 dilution of 1.14 mg/mL antibody stock was used in this study for background-free results.
Protein ladder
Request a detailed protocolMagicMark XP Western Protein Standard (Thermo Fischer Scientific, Catalog #LC5602).
Related resource
Request a detailed protocolTurnšek J. 2017. HA tag enables highly efficient detection of heterologous proteins in Phaeodactylum tricornutum (Pt) exconjugants. protocols.io. DOI: dx.doi.org/10.17504/protocols.io.j7ncrme.
Amplex UltraRed assay and resorufin imaging
Amplex UltraRed assay
Request a detailed protocolFive mL of triplicate WT or pTF-APEX2 expressing P. tricornutum cultures in mid- to late-exponential phase were incubated on ice for 5 min, spun down, supernatant discarded, pellet resuspended in 500 μL ice-cold PBS (pH 7.4), and transferred to microcentrifuge tubes. Cells were spun down again, resuspended in 200 μL ice-cold reaction buffer (50 μM Amplex UltraRed [AUR; Thermo Fischer Scientific, Catalog #A36006], 2 mM H2O2, in PBS [pH 7.4]), and incubated on ice for 15 min unless otherwise noted. Fifty μL supernatant was mixed with 50 μL PBS (pH 7.4) and resorufin fluorescence measured in a black microtiter plate with black bottom using Flexstation 3 microtiter plate reader (Molecular Devices, San Jose, CA; excitation: 544 nm, emission: 590 nm; resorufin excitation and emission maxima are 568 nm and 581 nm, respectively). Horseradish peroxidase (HRP) was always included as a positive assay control. Fluorescence was normalized to OD750 of experimental P. tricornutum cultures. Amplex UltraRed was prepared as a 10 mM stock in DMSO and stored in 20 μL aliquots at −20°C. Three percent (w/w) H2O2 stock (Sigma-Aldrich, Saint Louis, MO, Catalog #323381–25 ML) was stored in 100 μL aliquots at −20°C.
Resorufin imaging
Request a detailed protocolWT and pTF-APEX2 expressing P. tricornutum cells after performing the Amplex UltraRed assay were imaged with Leica TCS SP5 confocal microscope using the following parameters: argon laser strength at 30%, 514 nm laser line at 50% maximum strength, resorufin emission window: 575–605 nm, autofluorescence emission window: 700–750 nm.
Transmission electron microscopy (TEM)
Part 1: Labeling
Request a detailed protocolFive mL of triplicate WT or pTF-APEX2 expressing P. tricornutum cultures in mid- to late-exponential phase were spun down (4000 rpm, 4°C, 10 min) and fixed in 5 mL ice-cold 2% (w/v) paraformaldehyde (PFA) and 2% (v/v) glutaraldehyde in 0.15 M sodium cacodylate buffer (pH 7.4) for 30 min on ice. Cells were rinsed in 5 mL 0.15 M sodium cacodylate buffer (pH 7.4) five times for 3 min on ice, then once again in 5 mL 0.15 M sodium cacodylate buffer (pH 7.4) with 10 mM glycine for 3 min on ice. Cells were then treated with 25 mM 3,3’-diaminobenzidine (DAB; Sigma-Aldrich, Catalog #D8001) as follows: 5.36 mg DAB was dissolved in 1 mL 0.1 N HCl and sonicated for 45 min. Five mL of 0.3 M sodium cacodylate buffer (pH 7.4) was added to dissolved DAB, final volume adjusted to 10 mL with ddH2O, solution filtered through a 0.22 μm syringe filter, and 3 μL 30% (w/w) H2O2 added for 3 mM final concentration. Cells were incubated in 5 mL of this reaction buffer for 15 min on ice, rinsed five times for 3 min in 5 mL 0.15 M sodium cacodylate buffer (pH 7.4), post-fixed in 2 mL 1% osmium tetroxide (OsO4) (Electron Microscopy Sciences, Hatfield, PA, Catalog #19150) in 0.15 M sodium cacodylate buffer (pH 7.4) for 30 min on ice, and rinsed in 5 mL ice-cold ddH2O five times for 3 min. Cells were then resuspended in 300 μL melted 3% agar, poured onto a glass slide sitting on ice, left to solidify, and cut into small ~3×3×1 mm pieces which were transferred into scintillation vials with 10 mL ice-cold ddH2O. All shorter spin down steps were done at 4000 rpm and 4°C for 1.5 min. All buffers were used ice-cold.
Part 2: TEM sample preparation
Request a detailed protocolAgar blocks with embedded P. tricornutum cells were fixed in 10 mL 2% (v/v) glutaraldehyde in ddH2O to crosslink agar, rinsed five times for 2 min in 5 mL ice-cold ddH2O, and left incubating in 5 mL 2% uranyl acetate (Electron Microscopy Sciences, Catalog #22400) overnight at 4°C. Next morning, agar blocks were first dehydrated in the following ethanol series: 20%, 50%, 70%, 90%, 100% (on ice, 10 mL), 100%, 100% (at RT, 10 mL), then infiltrated with 10 mL 50% epoxy resin for ~1 hr. To prepare 20 mL 100% resin, Durcupan ACM mixture (Electron Microscopy Sciences, Catalog #14040) components were combined as follows: 11.4 g A (epoxy resin), 10 g B (964 hardener), 0.3 g C (964 accelerator), and 0.1 g D (dibutyl phthalate) in this exact order (for 50% resin 1 part 100% resin was combined with 1 part 100% ethanol). Agar blocks were transferred into 10 mL 100% resin for ~4 hr, then into fresh 10 mL 100% resin overnight, again into fresh 10 mL 100% resin for 4 hr the following morning before finally being poured into aluminum boats and left to polymerize in 60°C oven for at least 48 hr (over the weekend). Polymerized resins were detached from aluminum boats, agar blocks dense with cells cut out and glued to ‘dummy’ blocks by incubation in a 60°C oven for at least 15 min. Approximately 500-nm-thick sections were cut with an ultramicrotome (Leica Ultracut), stained with 1% toluidine blue, and observed under light microscope to make sure embedded cells were exposed. Approximately 100-nm-thick sections were then cut, placed on TEM grids, labeled, and saved until imaging. Imaging was performed with JEOL JEM-1200 (Japan Electron Optics Laboratory, Akishima, Tokyo, Japan) transmission electron microscope at 80 keV. Backscatter scanning electron microscopy was performed with Zeiss Merlin (Oberkochen, Germany) scanning electron microscope (SEM) at two keV by placing ~80-nm-thick sections on a silicon wafer and imaged with inverted contrast which gives a TEM-like image.
Proximity-dependent proteomic mapping
Part 1: Labeling
Request a detailed protocolTwenty-five mL of quintuplicate WT and pTF-APEX2 expressing P. tricornutum cultures in mid- to late-exponential phase—cell density of all cultures just prior to harvest was ~2×107 cells/mL which corresponds to OD750 ~0.4—were cooled on ice for 10 min and pelleted (4000 rpm, 4°C, 10 min). Supernatants were discarded, pellets resuspended in 0.5 mL ice-cold PBS (pH 7.4), transferred to microcentrifuge tubes, and spun down (4000 rpm, 4°C, 10 min). Cells were then resuspended in 0.5 mL ice-cold 1.2 M D-sorbitol in PBS (pH 7.4), supplemented with 2.5 mM biotin-phenol (Berry and Associates, Dexter, MI, Catalog #BT 1015), incubated on a tube rotator at 4°C for 90 min, supplemented with 1 mM H2O2, and incubated on a tube rotator at 4°C for another 20 min. Labeling reaction was quenched by washing cells twice (4000 rpm, 4°C, 5 min) with 0.5 mL ice-cold quenching solution (10 mM sodium ascorbate [VWR International, Radnor, PA, Catalog #95035–692], 5 mM Trolox [Sigma-Aldrich, Catalog #238813–5G], 10 mM sodium azide [VWR International, Catalog #AA14314-22] in PBS [pH 7.4]). Of quenched cell suspension, 50 μL was saved for a streptavidin blot. The remaining 450 μL was spun down (4000 rpm, 4°C, 10 min) and lysed in 250 μL cell lysis buffer (50 mM Tris-HCl, 200 mM NaCl, 1 mM DTT, 1 mM PMSF, pH 8.5) by sonication for 15 min (30 s on, 1 min off). Cell lysates were spun down at 4000 rpm and 4°C for 45 min and protein concentration in supernatants measured using Bradford Assay Kit. For streptavidin blot, 25 μL of saved quenched cells were first sonicated for 15 min (30 s on, 1 min off). A total of 2.5 μL of whole cell lysates and 1 ng biotinylated BSA positive control (Thermo Fischer Scientific, Catalog #29130) were then resolved on a NuPage 4–12% Bis-Tris 1.5 mm gel and wet transferred to a PVDF membrane. Membrane was washed twice with 15 mL PBST (PBS [pH 7.4] with 0.1% [v/v] Tween-20) for 5 min, left blocking overnight at RT and gentle shaking in PBST supplemented with 5% BSA (Sigma-Aldrich, Catalog #A9647-100G), and washed once the next morning with 15 mL PBST for 5 min. It was then incubated for 1 hr at RT and gentle shaking in 15 mL PBST supplemented with 5% BSA and 1:15,000 HRP-conjugated streptavidin (Thermo Fischer Scientific, Catalog #S911). Finally, membrane was washed three times with 15 mL PBST supplemented with 5% BSA for 5 min and once with 15 mL PBST for 10 min after which it was visualized with SuperSignal West Dura Extended Duration reagent (Thermo Fischer Scientific, Catalog #34075) using C-DiGit Blot Scanner (Li-COR, Lincoln, NE).
Protein ladder
Request a detailed protocolWesternSure Pre-stained Chemiluminescent Protein Ladder (HRP-conjugated protein ladder) (Li-COR, Catalog #926–98000).
Part 2: Pull-down and mass spectrometry
Request a detailed protocolPull-down
Request a detailed protocolFifty μL streptavidin beads (Thermo Fischer Scientific, Catalog #88816) transferred to microcentrifuge tubes were pelleted in a magnetic rack (tubes were left standing for ~5 min to pellet fully). Supernatants were discarded and beads were washed twice with 1 mL ice-cold cell lysis buffer (50 mM Tris-HCl, 200 mM NaCl, 1 mM DTT, 1 mM PMSF, pH 8.5). A total of 360 μg proteins in 500 μL total volume (x μL supernatant from cell lysis with 360 μg proteins and 500-x μL cell lysis buffer) were incubated on a tube rotator overnight at 4°C. Streptavidin beads were washed to remove nonspecific binders the next morning as follows: 2 × 1 mL cell lysis buffer, 1 × 1 mL 1 M KCl, 1 × 1 mL 0.1 M Na2CO3 (pH 11.5), 1 × 1 mL 2 M urea (pH 8.0) in 10 mM Tris-HCl, 2 × 1 mL cell lysis buffer, 2 × 1 mL PBS (pH 7.4). PBS after the final washing step was removed before storing beads at −80°C. Notes: All streptavidin beads collection steps were 5 min long. All the washing solutions were ice-cold. Microcentrifuge tubes were either very briefly vortexed (~2 s) or tapped by hand between each washing step to promote bead resuspension.
On-bead digestion and TMT labeling
Request a detailed protocolSamples were prepared as previously described (Kalocsay, 2019). Liquid reagents used were HPLC quality grade. Washed beads were resuspended in 50 μL of 200 mM EPPS (4-(2-hydroxyethyl)−1-piperazinepropanesulfonic acid) buffer (pH 8.5) and 2% (v/v) acetonitrile with 1 μL of 2 mg/mL lysil endoproteinase Lys-C stock solution (FUJIFILM Wako Pure Chemical Corporation, Richmond, VA, Catalog #125–05061), briefly vortexed, and incubated at 37°C for 3 hr. Fifty μL of trypsin stock (Promega, Madison, WI, Catalog #V5111) diluted 1:100 (v/v) in 200 mM EPPS (pH 8.5) was then added. After mixing, digests were incubated at 37°C overnight and beads were magnetically removed. Peptides were then directly labeled as follows: acetonitrile was added to 30% (v/v) concentration and peptides were labeled with TMT 10-plex reagent (Thermo Fisher Scientific, Catalog #90406) for 1 hr. Labeling reactions were quenched with hydroxylamine at a final concentration of 0.3% (v/v) for 15 min and 1% of labeled peptides was analyzed for label incorporation efficiency by mass spectrometry. After quenching, peptide solutions were first acidified with formic acid, trifluoroacetic acid (TFA) was then added to a concentration of 0.1% (v/v), and peptides were desalted by acidic C18 solid phase extraction (StageTip). Labeled peptides were finally resuspended in 1% (v/v) formic acid and 3% (v/v) acetonitrile.
Mass spectrometry
Request a detailed protocolData were collected with a MultiNotch MS3 TMT method (McAlister et al., 2014) using an Orbitrap Lumos mass spectrometer coupled to a Proxeon EASY-nLC 1200 Liquid Chromatography (LC) system (both Thermo Fisher Scientific). The used capillary column was packed with C18 resin (35 cm length, 100 μm inner diameter, 2.6 μm Accucore matrix [Thermo Fisher Scientific]). Peptides were separated for 3 or 4 hr over acidic acetonitrile gradients by LC prior to mass spectrometry (MS) analysis. Data from two 4 hr runs and one 3 hr run were recorded and combined. After an initial MS1 scan (Orbitrap analysis; resolution 120,000; mass range 400–1400 Th), MS2 analysis used collision-induced dissociation (CID, CE = 35) with a maximum ion injection time of 150–300 ms and an isolation window of 0.5 m/z. In order to obtain quantitative information, MS3 precursors were fragmented by high-energy collision-induced dissociation (HCD) and analyzed in the Orbitrap at a resolution of 50,000 at 200 Th. Further details on LC and MS parameters and settings were described recently (Paulo et al., 2016). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2019) with the dataset identifier PXD018022.
pTF-mCherry and MS hit-EYFP co-expression and imaging
Request a detailed protocolSequence-validated episomes with pTF-mCherry and MS hit-EYFP expression cassettes were conjugated into WT P. tricornutum cells, resulting transconjugants genotyped, and screened for fluorescence. Five μL of fluorescent transconjugant cell lines in mid- to late-exponential phase were imaged with settings that minimized cross-channel bleed-through: argon laser strength at 30%, 514 nm laser line at 50% maximum strength, mCherry emission window: 620–640 nm, EYFP emission window: 520–540 nm, autofluorescence emission window: 700–750 nm.
pTF.CREG1 protein expression conditions
Request a detailed protocolSmall-scale expression testing for the Δ31_pTF.CREG1-His6 (and all other pTF.CREG1 truncations not described here) was performed as follows: an overnight culture of BL21 E. coli cells carrying the PtpBAD-Δ31_pTF.CREG1-CTHF (=pJT_Δ31_pTF.CREG1-His6) expression vector was used to inoculate 50 mL of Terrific Broth (TB) supplemented with tetracycline (10 µg/mL). Cultures were grown in a shaking incubator (37°C, 200 rpm) until OD600 of 0.4–0.6 was reached. Incubator temperature was then lowered to 18°C and flasks were allowed to adjust to this temperature for about 30 min before arabinose was added to a final concentration of 0.5% (w/v). Growth at 18°C was continued overnight (12–18 hr) after which 10 mL of cultures were harvested by centrifugation at 6000 g for 10 min. Large-scale Δ31_pTF.CREG1-His6 expression was performed as above, but at 2 L total volume. Following arabinose induction and overnight growth of the 2 L culture, 500 mL was harvested for further processing.
pTF.CREG1 purification conditions
Request a detailed protocolAll purifications were performed using cobalt, Co2+, TALON Metal Affinity Resin (Takara Bio, Catalog #635502) which binds His6-tagged proteins with high affinity. Pellets from small scale expression testing (10 mL culture) were resuspended in 800 µL of lysis buffer (50 mM NaH2PO4 [pH 7.5], 500 mM NaCl, 0.1% [v/v] Triton X-100, 10 mM imidazole, 1 mg/mL lysozyme, and 10 μM β-mercaptoethanol) and subjected to microtip sonication on ice until lysis was complete. Lysates were then clarified by centrifugation (10 min, 15,000 g, 4°C) and supernatants (700 µL) were set aside. For each purification, 25 µL of TALON resin was equilibrated by washing three times with 10 volumes of lysis buffer and pelleting by centrifugation after each wash (30 s, 3000 g, RT). After third wash, 25 µL of lysis buffer was added to the resin to make a 50 µL slurry. Supernatants from previously clarified cell lysates were then added to the TALON resin slurry and incubated for 1 hr at RT with end-over-end mixing. Following the 1 hr incubation, the resin was pelleted by centrifugation (30 s, 3000 g, RT) and washed three times with 10 volumes of wash buffer (50 mM NaH2PO4 [pH 7.5], 500 mM NaCl, 30 mM imidazole). Proteins were eluted with 50 uL of elution buffer (50 mM NaH2PO4 [pH 7.5], 500 mM NaCl, 250 mM imidazole, 10% glycerol) and were either subjected to immediate SDS-PAGE electrophoresis or stored at −80°C.
Pellets from large-scale expression (500 mL harvested bacterial culture) were processed with a few modifications. For lysis, cell pellets were resuspended in 20 mL of lysis buffer. For purification, 1 mL of TALON resin (i.e. 2 mL of equilibrated slurry) was used in combination with approximately 19 mL of clarified cell lysate which was allowed to incubate overnight at 8°C with end-over-end mixing. Final protein elution was performed using 5 mL of elution buffer. Purified protein was subjected to SDS-PAGE electrophoresis and concentrated using an Amicon Ultra-15 10 kDa cutoff concentrator (Millipore Sigma, Burlington, MA, Catalog #UFC901024). Following the concentration step, buffer exchange into 20 mM HEPES (pH 8.0), 10% glycerol, 300 mM KCl was performed using a PD-10 desalting column (GE Healthcare, Chicago, IL, Catalog #17-0851-01). Protein was then concentrated again and total protein content was determined by the Bradford method using the Protein Assay Dye Reagent (Bio-Rad, Hercules, CA, Catalog #5000006).
Flavin reduction assays
Request a detailed protocolTo test for flavin reductase activity, we set up the following enzyme assay conditions in 100 uL final volume: 50 mM Tris-HCl (pH 7.5), 100 mM KCl, 10% glycerol, 30 μM flavin mononucleotide (FMN), 30 μM riboflavin, and 250 μM NADPH (Coves and Fontecave, 1993). Assays were started by the addition of 50 μg of pTF.CREG1, BSA, or Milli-Q water for no enzyme controls. Oxidation of NADPH was measured by the decrease in absorbance at 320 nm on a Flexstation 3 microtiter plate reader and the signal for each reaction was tracked over a period of 60 min. All assays were conducted in triplicate.
Bioinformatic and data analyses
Mass spectrometry data analysis
Request a detailed protocolPeptide-spectrum matches used a SEQUEST (v. 28, rev. 12) algorithm (Eng et al., 1994). Data were searched against a size-sorted forward and reverse database of the Phaeodactylum tricornutum proteome (strain CCAP 1055/1, UniProt reference proteome UP000000759) with added common contaminant proteins and the pTF-APEX2 fusion protein sequence. Spectra were first converted to mzXML and searches were then performed using a mass tolerance of 50 ppm for precursors and a fragment ion tolerance of 0.9 Da. For the searches, maximally two missed cleavages per peptide were allowed. We searched dynamically for oxidized methionine residues (+15.9949 Da) and applied a target decoy database strategy. A false discovery rate (FDR) of 1% was set for peptide-spectrum matches following filtering by linear discriminant analysis (LDA) (Beausoleil et al., 2006; Huttlin et al., 2010). The FDR for final collapsed proteins was 1%. Quantitative peptide information was derived from MS3 scans. Quant tables were generated with the following filter criteria: MS2 isolation specificity of >70% for each peptide and a sum of TMT signal-to-noise (s/n) of >200 over all channels per peptide. Quant tables were exported to Excel and further processed therein. Details of the TMT intensity quantification method and additional applied search parameters were described previously (Paulo et al., 2016). Scaled proteomic data were subjected to two-way hierarchical clustering (Ward’s method) using JMP software package (SAS Institute). Volcano plot with log2-transformed average APEX2/WT ratios and associated p values was made in R using ggplot2 data visualization package (R Development Core Team, 2013). Gene IDs corresponding to protein hits were inferred using Ensembl Protists Phatr3 P. tricornutum genomic database.
Protein feature identification
Request a detailed protocolProtein lengths, molecular weights, and isoelectric points (pI) were determined with ProtParam (Gasteiger et al., 2005). Signal peptides and transmembrane regions were identified with SignalP 4.1 (Petersen et al., 2011) and TMHMM Server v. 2.0 (Krogh et al., 2001), respectively. Protein localizations were predicted with a combination of tools: TargetP 1.1 (Emanuelsson et al., 2000), SignalP 4.1, and ASAFind (Gruber et al., 2015) version 1.1.7. All putative chloroplastic localizations in the study mean that a protein was predicted to be chloroplastic by ASAFind (state-of-the-art plastidial protein localization prediction tool for diatoms); chloroplast localization prediction confidences are noted. Low confidence prediction by ASAFind means that a protein satisfies the following filtering criteria: (1) it contains a signal peptide as detected by SignalP 4.1, (2) +1 position of ASAFind predicted cleavage site is phenylalanine (F), tryptophan (W), tyrosine (Y), or leucine (L) (making the protein ‘potentially plastid targeted’; it was tyrosine in pTF.CREG1), and (3) one or both of the following is false: the ASAFind predicted cleavage site coincides with the SignalP 4.1 prediction and the transit peptide score is higher than 2 (the latter is false for pTF.CREG1: transit peptide score was ~0.81). Peroxidase class prediction was done in RedoxiBase (Savelli et al., 2019). N-glycosylation and GalNAc-type O-glycosylation predictions were performed with NetNGlyc 1.0 (Gupta et al., 2004) and NetOGlyc 4.0 (Steentoft et al., 2013), respectively.
Amino acid sequence alignments
Request a detailed protocolAPEX2-like P. tricornutum peroxidases were aligned with Clustal Omega (Sievers et al., 2011) using the default settings and CLC Sequence Viewer 7.7 (QUIAGEN) using the following parameters: Gap open cost: 10.0, Gap extension cost: 1.0, End gap cost: As any other, Alignment: Very accurate (slow). Putative substrate-binding loops in these peroxidases were evaluated using the Clustal Omega alignment. Conserved motifs in amino acid sequence alignments underlying the phylogenetic trees were displayed in CLC Sequence Viewer 7.7.
Phylogenetic analyses
Request a detailed protocolHomologs of respective P. tricornutum proteins were retrieved from the non-redundant National Center for Biotechnology Information (NCBI) and the Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP) databases using the BLASTP algorithm (e-value threshold set to e−15) (Altschul et al., 1990). The BLAST search retrieved only a handful of homologs for pTF.ap1 and no homologs for pTF.CatCh1 even among closely related diatoms. Therefore, more sensitive HMMER (Eddy, 1998) searches were employed to extend the datasets, which were afterward aligned using the Localpair algorithm as implemented in MAFFT (Katoh et al., 2002). Ambiguously aligned regions, regions composed mostly of gaps, and short fragments were manually removed in SeaView 4 (Gouy et al., 2010). For each alignment, the maximum likelihood analysis was carried out in IQ-TREE (Nguyen et al., 2015) under the best-fitting substitution matrix as inferred by the built-in model finder. Branching support was estimated using ‘thorough’ non-parametric bootstrap analysis from 500 replicates in IQ-TREE.
Identification of disordered protein regions
Request a detailed protocolpTF.CatCh1 amino acid sequence was analyzed for the presence of disordered protein regions with PONDR (Predictor of Naturally Disordered Regions) VSL2 predictor (Peng et al., 2006) and IUPred2a long disorder prediction type (Mészáros et al., 2018).
Homology modeling of pTF.CREG1
Request a detailed protocolFull-length pTF.CREG1 amino acid sequence was used as an input for Phyre2, an online protein structure prediction server (Kelley et al., 2015). Normal modeling mode was selected. Human CREG1 (PDB ID: 1XHN) was identified as the best modeling template (human CREG1 is ~37% identical to pTF.CREG1 as inferred from their Clustal Omega alignment). pTF.CREG1 was modeled with >90% confidence across 182 modeled amino acid residues (Pro43–Lys224) indicating the predicted structure is representative of the actual structure. pTF.CREG1 homology model was aligned with human CREG1 monomer in UCSF Chimera version 1.11.1 using MatchMaker tool with default settings (Pettersen et al., 2004).
Other data analyses
Request a detailed protocolAmplex UltraRed assay and transcriptomic data were plotted in R using ggplot2 data visualization package (R Development Core Team, 2013). Flavin reduction assay data were plotted in Excel.
Appendix 1
Glossary
AA/AAs: amino acid residue/amino acid residues
APEX: engineered ascorbate peroxidase (APX)
CREG: cellular repressor of E1A-stimulated genes
DAB: 3,3’-diaminobenzidine
EYFP: enhanced yellow fluorescent protein
Fe’: dissolved labile iron (all unchelated iron species)
FMN: flavin mononucleotide
HNLC: high-nutrient, low-chlorophyll
ISIP: iron starvation induced protein
MS: mass spectrometry
PBS: phosphate-buffered saline
pTF: phytotransferrin
RT: room temperature
TEM: transmission electron microscopy
TM: transmembrane
TMT: tandem mass tag
UTR: untranslated region
V-ATPase: vacuolar-type H+-ATPase
WT: wild type
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files.
-
PRIDEID PXD018022. Quantitative proximity proteomics suggests a phytotransferrin-mediated cell surface-to-chloroplast iron trafficking axis in marine diatoms.
-
UniProtID UP000000759. Phaeodactylum tricornutum proteome.
-
EnsemblID ASM15095v2. Phaeodactylum tricornutum genome annotation 3 (Phatr3).
-
NCBI BioProjectID PRJNA231566. Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP) database.
References
-
A manganese-rich environment supports superoxide dismutase activity in a lyme disease pathogen, Borrelia burgdorferiJournal of Biological Chemistry 288:8468–8478.https://doi.org/10.1074/jbc.M112.433540
-
Evolution and functional diversification of fructose bisphosphate aldolase genes in photosynthetic marine diatomsMolecular Biology and Evolution 29:367–379.https://doi.org/10.1093/molbev/msr223
-
Basic local alignment search toolJournal of Molecular Biology 215:403–410.https://doi.org/10.1016/S0022-2836(05)80360-2
-
The puzzle of plastid evolutionCurrent Biology 19:R81–R88.https://doi.org/10.1016/j.cub.2008.11.067
-
Recognition- and reactivity-based fluorescent probes for studying transition metal signaling in living systemsAccounts of Chemical Research 48:2434–2442.https://doi.org/10.1021/acs.accounts.5b00221
-
Iron cofactor assembly in plantsAnnual Review of Plant Biology 65:125–153.https://doi.org/10.1146/annurev-arplant-050213-035759
-
Structural basis for metal binding specificity: the N-terminal cadmium binding domain of the P1-type ATPase CadAJournal of Molecular Biology 356:638–650.https://doi.org/10.1016/j.jmb.2005.11.055
-
Cracking the glycome encoder: signaling, trafficking, and glycosylationTrends in Cell Biology 26:379–388.https://doi.org/10.1016/j.tcb.2015.12.004
-
Pyrenoids: CO2-fixing phase separated liquid organellesBiochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1868:118949.https://doi.org/10.1016/j.bbamcr.2021.118949
-
Iron uptake proteins in algae and the role of iron Starvation-Induced proteins (ISIPs)European Journal of Phycology 55:339–360.https://doi.org/10.1080/09670262.2020.1744039
-
Protein secretion and the endoplasmic reticulumCold Spring Harbor Perspectives in Biology 4:a012872.https://doi.org/10.1101/cshperspect.a012872
-
The evolution of diatoms and their biogeochemical functionsPhilosophical Transactions of the Royal Society B: Biological Sciences 372:20160397.https://doi.org/10.1098/rstb.2016.0397
-
The ins and outs of algal metal transportBiochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1823:1531–1552.https://doi.org/10.1016/j.bbamcr.2012.04.010
-
Lysosome-related organelles as mediators of metal homeostasisJournal of Biological Chemistry 289:28129–28136.https://doi.org/10.1074/jbc.R114.592618
-
Regulators of iron homeostasis: new players in metabolism, cell death, and diseaseTrends in Biochemical Sciences 41:274–286.https://doi.org/10.1016/j.tibs.2015.11.012
-
UV-light-driven prebiotic synthesis of iron-sulfur clustersNature Chemistry 9:1229–1234.https://doi.org/10.1038/nchem.2817
-
Metabolic enzymes moonlighting in the nucleus: metabolic regulation of gene transcriptionTrends in Biochemical Sciences 41:712–730.https://doi.org/10.1016/j.tibs.2016.05.013
-
Efficient proximity labeling in living cells and organisms with TurboIDNature Biotechnology 36:880–887.https://doi.org/10.1038/nbt.4201
-
Endocytic traffic: vesicle fusion cascade in the early endosomesCurrent Biology 22:R597–R598.https://doi.org/10.1016/j.cub.2012.06.009
-
Enhanced silica ballasting from iron stress sustains carbon export in a frontal zone within the California currentJournal of Geophysical Research: Oceans 120:4654–4669.https://doi.org/10.1002/2015JC010829
-
Photoferrotrophy: remains of an ancient photosynthesis in modern environmentsFrontiers in Microbiology 8:323.https://doi.org/10.3389/fmicb.2017.00323
-
The endocytic fate of the transferrin receptor is regulated by c-Abl kinaseJournal of Biological Chemistry 291:16424–16437.https://doi.org/10.1074/jbc.M116.724997
-
Probing the evolution, ecology and physiology of marine protists using transcriptomicsNature Reviews Microbiology 15:6–20.https://doi.org/10.1038/nrmicro.2016.160
-
Fructose bisphosphate aldolase is involved in the control of RNA polymerase III-directed transcriptionBiochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1843:1103–1110.https://doi.org/10.1016/j.bbamcr.2014.02.007
-
Transcriptomic and proteomic responses of the oceanic diatom Pseudo-nitzschia granii to iron limitationEnvironmental Microbiology 20:3109–3126.https://doi.org/10.1111/1462-2920.14386
-
Reduction and mobilization of iron by a NAD(P)H:flavin oxidoreductase from Escherichia coliEuropean Journal of Biochemistry 211:635–641.https://doi.org/10.1111/j.1432-1033.1993.tb17591.x
-
BookIron Metabolism: From Molecular Mechanisms to Clinical ConsequencesHoboken, NJ, USA: John Wiley & Sons, Inc.https://doi.org/10.1002/9781118925645
-
Synthesis of iron fertilization experiments: from the iron age in the age of enlightenmentJournal of Geophysical Research 110:C09S16.https://doi.org/10.1029/2004JC002601
-
Refinement of the diatom episome maintenance sequence and improvement of Conjugation-Based DNA delivery methodsFrontiers in Bioengineering and Biotechnology 4:65.https://doi.org/10.3389/fbioe.2016.00065
-
Metals in proteins: correlation between the metal-ion type, coordination number and the amino-acid residues involved in the coordinationActa Crystallographica Section D Biological Crystallography 64:257–263.https://doi.org/10.1107/S090744490706595X
-
Profile hidden markov modelsBioinformatics 14:755–763.https://doi.org/10.1093/bioinformatics/14.9.755
-
The pfam protein families database in 2019Nucleic Acids Research 47:D427–D432.https://doi.org/10.1093/nar/gky995
-
Predicting subcellular localization of proteins based on their N-terminal amino acid sequenceJournal of Molecular Biology 300:1005–1016.https://doi.org/10.1006/jmbi.2000.3903
-
An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein databaseJournal of the American Society for Mass Spectrometry 5:976–989.https://doi.org/10.1016/1044-0305(94)80016-2
-
The vacuolar H+-ATPase: a universal proton pump of eukaryotesBiochemical Journal 324 (Pt 3):697–712.https://doi.org/10.1042/bj3240697
-
Functional diversity of cysteine residues in proteins and unique features of catalytic Redox-Active cysteines in thiol oxidoreductasesMolecules and Cells 26:228–235.
-
Identity and functions of CxxC-derived motifsBiochemistry 42:11214–11225.https://doi.org/10.1021/bi034459s
-
Nanostructure and nanomechanics of live Phaeodactylum tricornutum morphotypesEnvironmental Microbiology 10:1344–1356.https://doi.org/10.1111/j.1462-2920.2007.01551.x
-
BookComplex Endosymbioses I: From Primary to Complex Plastids, Multiple IndependentEvents.In: Maréchal E, editors. Plastids. Methods in Molecular Biology, 1829. Totowa, United States: Humana Press, Inc. pp. 17–35.https://doi.org/10.1007/978-1-4939-8654-5_2
-
BookProtein Identification and Analysis Tools on the ExPASy ServerIn: Walker JM, editors. The Proteomics Protocols Handbook. Totowa, United States: Humana Press, Inc. pp. 571–607.https://doi.org/10.1385/1-59259-890-0:571
-
The structure and biological function of CREGFrontiers in Cell and Developmental Biology 6:136.https://doi.org/10.3389/fcell.2018.00136
-
Mannose 6-phosphate receptors: new twists in the taleNature Reviews Molecular Cell Biology 4:202–213.https://doi.org/10.1038/nrm1050
-
BookThe Chloroplast Endoplasmic Reticulum: Structure, Function, and Evolutionary SignificanceIn: Bourne G. H, Danielli J. F, Jeon K. W, editors. International Review of Cytology. New York, United States: Academic Press, Inc. pp. 49–99.https://doi.org/10.1016/S0074-7696(08)61194-8
-
Getting to know the neighborhood: using proximity-dependent biotinylation to characterize protein complexes and map organellesCurrent Opinion in Chemical Biology 48:44–54.https://doi.org/10.1016/j.cbpa.2018.10.017
-
Identification of trace element-containing proteins in genomic databasesAnnual Review of Nutrition 24:579–596.https://doi.org/10.1146/annurev.nutr.24.012003.132241
-
SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree buildingMolecular Biology and Evolution 27:221–224.https://doi.org/10.1093/molbev/msp259
-
Establishing super-resolution imaging for proteins in diatom biosilicaScientific Reports 6:36824.https://doi.org/10.1038/srep36824
-
The periplastidal compartment: a naturally minimized eukaryotic cytoplasmCurrent Opinion in Microbiology 22:88–93.https://doi.org/10.1016/j.mib.2014.09.017
-
Erythroid cell mitochondria receive endosomal iron by a "kiss-and-run" mechanismBiochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1863:2859–2867.https://doi.org/10.1016/j.bbamcr.2016.09.008
-
CREG inhibits migration of human vascular smooth muscle cells by mediating IGF-II endocytosisExperimental Cell Research 315:3301–3311.https://doi.org/10.1016/j.yexcr.2009.09.013
-
Glycosylation-independent binding to extracellular domains 11-13 of mannose-6-phosphate/insulin-like growth factor-2 receptor mediates the effects of soluble CREG on the phenotypic modulation of vascular smooth muscle cellsJournal of Molecular and Cellular Cardiology 50:723–730.https://doi.org/10.1016/j.yjmcc.2010.12.013
-
Understanding diatom cell wall Silicification—Moving ForwardFrontiers in Marine Science 5:125.https://doi.org/10.3389/fmars.2018.00125
-
Subcellular metal imaging identifies dynamic sites of cu accumulation in ChlamydomonasNature Chemical Biology 10:1034–1042.https://doi.org/10.1038/nchembio.1662
-
Marine phytoplankton and the changing ocean iron cycleNature Climate Change 6:1072–1079.https://doi.org/10.1038/nclimate3147
-
Molecular regulation of the diatom cell cycleJournal of Experimental Botany 65:2573–2584.https://doi.org/10.1093/jxb/ert387
-
A Golgi rhomboid protease Rbd2 recruits Cdc48 to cleave yeast SREBPThe EMBO Journal 35:2332–2349.https://doi.org/10.15252/embj.201693923
-
Proximity-dependent biotin labelling in yeast using the engineered ascorbate peroxidase APEX2Biochemical Journal 473:2463–2469.https://doi.org/10.1042/BCJ20160106
-
Protein moonlighting: what is it, and why is it important?Philosophical Transactions of the Royal Society B: Biological Sciences 373:20160523.https://doi.org/10.1098/rstb.2016.0523
-
Use of fluorescent probes: their effect on cell biology and limitationsThe Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology 295:2031–2036.https://doi.org/10.1002/ar.22602
-
Catalysis of Template-Directed nonenzymatic RNA copying by iron(II)Journal of the American Chemical Society 140:15016–15021.https://doi.org/10.1021/jacs.8b09617
-
BookAPEX Peroxidase-Catalyzed Proximity Labeling and Multiplexed Quantitative ProteomicsIn: Sunbul M, Jäschke A, editors. Proximity Labeling. TotawaUnited States: Human Press, Inc. pp. 41–55.https://doi.org/10.1007/978-1-4939-9537-0_4
-
Designer diatom episomes delivered by bacterial conjugationNature Communications 6:6925.https://doi.org/10.1038/ncomms7925
-
MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transformNucleic Acids Research 30:3059–3066.https://doi.org/10.1093/nar/gkf436
-
ConferenceMonitoring the localization of three iron starvation induced proteins (ISIPs) in a single transformant of Phaeodactylum tricornutum under a range of iron supplementation regimesThe Molecular Life of Diatoms Programme and Abstract Book. pp. 95–96.
-
The Phyre2 web portal for protein modeling, prediction and analysisNature Protocols 10:845–858.https://doi.org/10.1038/nprot.2015.053
-
Life: the first two billion yearsPhilosophical Transactions of the Royal Society B: Biological Sciences 371:20150493.https://doi.org/10.1098/rstb.2015.0493
-
Predicting transmembrane protein topology with a hidden markov model: application to complete genomesJournal of Molecular Biology 305:567–580.https://doi.org/10.1006/jmbi.2000.4315
-
Effect of iron deficiency on diatom cell size and silicic acid uptake kineticsLimnology and Oceanography 49:1134–1143.https://doi.org/10.4319/lo.2004.49.4.1134
-
A clickable APEX probe for Proximity-Dependent proteomic profiling in yeastCell Chemical Biology 27:858–865.https://doi.org/10.1016/j.chembiol.2020.05.006
-
Chloroplast vesicle transportPhotosynthesis Research 138:361–371.https://doi.org/10.1007/s11120-018-0566-0
-
Chloroplast iron transport proteins - Function and impact on plant physiologyFrontiers in Plant Science 7:178.https://doi.org/10.3389/fpls.2016.00178
-
Spatial proteomics: a powerful discovery tool for cell biologyNature Reviews Molecular Cell Biology 20:285–302.https://doi.org/10.1038/s41580-018-0094-y
-
Ion and metabolite transport in the chloroplast of algae: lessons from land plantsCellular and Molecular Life Sciences 75:2153–2176.https://doi.org/10.1007/s00018-018-2793-0
-
Engineered ascorbate peroxidase as a genetically encoded reporter for electron microscopyNature Biotechnology 30:1143–1148.https://doi.org/10.1038/nbt.2375
-
Mechanisms of carbon dioxide acquisition and CO2sensing in marine diatoms: a gateway to carbon metabolismPhilosophical Transactions of the Royal Society B: Biological Sciences 372:20160403.https://doi.org/10.1098/rstb.2016.0403
-
The vacuolar-type H⁺-ATPase at a glance - more than a proton pumpJournal of Cell Science 127:4987–4993.https://doi.org/10.1242/jcs.158550
-
The intracellular trafficking pathway of transferrinBiochimica Et Biophysica Acta (BBA) - General Subjects 1820:264–281.https://doi.org/10.1016/j.bbagen.2011.09.009
-
WebsiteNutritional copper signaling and homeostasisGrantome Research Proposal Abstract. Accessed April 24, 2020.
-
Aldolase directly interacts with ARNO and modulates cell morphology and acidic vesicle distributionAmerican Journal of Physiology-Cell Physiology 300:C1442–C1455.https://doi.org/10.1152/ajpcell.00076.2010
-
The role of unchelated fe in the iron nutrition of phytoplanktonLimnology and Oceanography 53:400–404.https://doi.org/10.4319/lo.2008.53.1.0400
-
Structure and properties of recombinant human pyridoxine 5'-phosphate oxidaseProtein Science 12:1455–1463.https://doi.org/10.1110/ps.0356203
-
The enigmatic endosome - sorting the ins and outs of endocytic traffickingJournal of Cell Science 131:jcs216499.https://doi.org/10.1242/jcs.216499
-
IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogeniesMolecular Biology and Evolution 32:268–274.https://doi.org/10.1093/molbev/msu300
-
Metabolic gene clusters in eukaryotesAnnual Review of Genetics 52:159–183.https://doi.org/10.1146/annurev-genet-120417-031237
-
Metallochaperones, an intracellular shuttle service for metal ionsJournal of Biological Chemistry 275:25057–25060.https://doi.org/10.1074/jbc.R000006200
-
OperonsCellular and Molecular Life Sciences 66:3755–3775.https://doi.org/10.1007/s00018-009-0114-3
-
Parcas is the predominant Rab11-GEF for rhodopsin transport in Drosophila photoreceptorsJournal of Cell Science 132:jcs231431.https://doi.org/10.1242/jcs.231431
-
Length-dependent prediction of protein intrinsic disorderBMC Bioinformatics 7:208.https://doi.org/10.1186/1471-2105-7-208
-
The PRIDE database and related tools and resources in 2019: improving support for quantification dataNucleic Acids Research 47:D442–D450.https://doi.org/10.1093/nar/gky1106
-
UCSF chimera--a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
The ins and outs of iron: escorting iron through the mammalian cytosolFree Radical Biology and Medicine 133:112–117.https://doi.org/10.1016/j.freeradbiomed.2018.10.411
-
The basics of thiols and cysteines in redox biology and chemistryFree Radical Biology and Medicine 80:148–157.https://doi.org/10.1016/j.freeradbiomed.2014.11.013
-
Chloroplast-mitochondria cross-talk in diatomsJournal of Experimental Botany 63:1543–1557.https://doi.org/10.1093/jxb/err441
-
Bidirectional traffic between the golgi and the endosomes - machineries and regulationJournal of Cell Science 129:jcs.185702–3982.https://doi.org/10.1242/jcs.185702
-
Roles and maturation of iron-sulfur proteins in plastidsJBIC Journal of Biological Inorganic Chemistry 23:545–566.https://doi.org/10.1007/s00775-018-1532-1
-
SoftwareR: A Language and Environment for Statistical ComputingR Foundation for Statistical Computing, Vienna, Austria.
-
Evidence for an ER to Golgi to chloroplast protein transport pathwayTrends in Cell Biology 16:385–387.https://doi.org/10.1016/j.tcb.2006.06.003
-
Do N-glycoproteins have preference for specific sequons?Bioinformation 5:208–212.https://doi.org/10.6026/97320630005208
-
Spotlight on protein N-terminal acetylationExperimental & Molecular Medicine 50:1–13.https://doi.org/10.1038/s12276-018-0116-z
-
Regulation of photosynthetic electron transportBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1807:375–383.https://doi.org/10.1016/j.bbabio.2010.11.010
-
Metallochaperones: bind and deliverChemistry & Biology 9:673–677.https://doi.org/10.1016/s1074-5521(02)00156-4
-
High-throughput production of human proteins for crystallization: the SGC experienceJournal of Structural Biology 172:3–13.https://doi.org/10.1016/j.jsb.2010.06.008
-
CryoAPEX - an electron tomography tool for subcellular localization of membrane proteinsJournal of Cell Science 132:jcs222315.https://doi.org/10.1242/jcs.222315
-
Live diatom silica immobilization of multimeric and redox-active enzymesApplied and Environmental Microbiology 78:211–218.https://doi.org/10.1128/AEM.06698-11
-
Protein import into chloroplastsNature Reviews Molecular Cell Biology 5:198–208.https://doi.org/10.1038/nrm1333
-
A global ocean atlas of eukaryotic genesNature Communications 9:373.https://doi.org/10.1038/s41467-017-02342-1
-
Surface chemical composition of diatomsChemBioChem 10:2011–2024.https://doi.org/10.1002/cbic.200800811
-
Preparation in pure state of the effect group of yellow enzymesBiochemische Zeitschrift 275:344–346.
-
Type I signal peptidase: an overviewArchives of Biochemistry and Biophysics 441:107–111.https://doi.org/10.1016/j.abb.2005.07.013
-
The intrinsic disorder alphabet. III. dual personality of serineIntrinsically Disordered Proteins 3:e1027032.https://doi.org/10.1080/21690707.2015.1027032
-
Dancing protein clouds: the strange biology and chaotic physics of intrinsically disordered proteinsJournal of Biological Chemistry 291:6681–6688.https://doi.org/10.1074/jbc.R115.685859
-
Metals, toxicity and oxidative stressCurrent Medicinal Chemistry 12:1161–1208.https://doi.org/10.2174/0929867053764635
-
Regulation of cellular iron metabolismBiochemical Journal 434:365–381.https://doi.org/10.1042/BJ20101825
-
Current trends to comprehend lipid metabolism in diatomsProgress in Lipid Research 70:1–16.https://doi.org/10.1016/j.plipres.2018.03.001
Article and author information
Author details
Funding
Gordon and Betty Moore Foundation (GBMF3828)
- Andrew Ellis Allen
Gordon and Betty Moore Foundation (GBMF5006)
- Andrew Ellis Allen
National Science Foundation (NSF-OCE-1756884)
- Andrew Ellis Allen
National Science Foundation (NSF-MCB-1818390)
- Andrew Ellis Allen
Biological and Environmental Research (DE-SC0018344)
- Andrew Ellis Allen
Gordon and Betty Moore Foundation (GBMF4958)
- Jernej Turnšek
National Institutes of Health (1F31ES030613-01)
- John K Brunson
Czech Science Foundation (21-03224S)
- Miroslav Oborník
European Regional Development Fund (CZ.02.1.01/0.0/0.0/16_019/0000759)
- Miroslav Oborník
National Institutes of Health (F31 1F31ES030613-01)
- John K Brunson
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Andrea Thor and Mason Mackey for help with electron microscopy sample preparation and imaging; Marian Kalocsay for assistance with mass spectrometry experiments and data analyses; Pardis Gholami and Hong Zheng for laboratory assistance; Jeffrey D Martell, Jeffrey B McQuaid, Tyler H Coale, and Sarah R Smith for fruitful discussions. This project was supported by the Gordon and Betty Moore Foundation grants GBMF4958 (JT), GBMF3828 (AEA), and GBMF5006 (AEA), the National Science Foundation grants NSF-OCE-1756884 (AEA) and NSF-MCB-1818390 (AEA), the United States Department of Energy Genomics Science program grant DE-SC0018344 (AEA), the Czech Science Foundation grant 21-03224S (MO), and ERDF/ESF Centre for Research of Pathogenicity and Virulence of Parasites grant No.CZ.02.1.01/0.0/0.0/16_019/0000759 (MO). JKB was supported by the National Institutes of Health (NIH) Ruth L Kirschstein Predoctoral Individual National Research Service Award F31 1F31ES030613-01. National Center for Microscopy and Imaging Research (NCMIR) is supported by the National Institute for General Medical Sciences (NIGMS) of the National Institutes of Health (NIH).
Copyright
© 2021, Turnšek et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,790
- views
-
- 377
- downloads
-
- 18
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Microbiology and Infectious Disease
Phage-derived peptidoglycan hydrolases (i.e. lysins) are considered promising alternatives to conventional antibiotics due to their direct peptidoglycan degradation activity and low risk of resistance development. The discovery of these enzymes is often hampered by the limited availability of phage genomes. Herein, we report a new strategy to mine active peptidoglycan hydrolases from bacterial proteomes by lysin-derived antimicrobial peptide-primed screening. As a proof-of-concept, five peptidoglycan hydrolases from the Acinetobacter baumannii proteome (PHAb7-PHAb11) were identified using PlyF307 lysin-derived peptide as a template. Among them, PHAb10 and PHAb11 showed potent bactericidal activity against multiple pathogens even after treatment at 100°C for 1 hr, while the other three were thermosensitive. We solved the crystal structures of PHAb8, PHAb10, and PHAb11 and unveiled that hyper-thermostable PHAb10 underwent a unique folding-refolding thermodynamic scheme mediated by a dimer-monomer transition, while thermosensitive PHAb8 formed a monomer. Two mouse models of bacterial infection further demonstrated the safety and efficacy of PHAb10. In conclusion, our antimicrobial peptide-primed strategy provides new clues for the discovery of promising antimicrobial drugs.
-
- Ecology
- Microbiology and Infectious Disease
Predicting how species diversity changes along environmental gradients is an enduring problem in ecology. In microbes, current theories tend to invoke energy availability and enzyme kinetics as the main drivers of temperature-richness relationships. Here, we derive a general empirically-grounded theory that can explain this phenomenon by linking microbial species richness in competitive communities to variation in the temperature-dependence of their interaction and growth rates. Specifically, the shape of the microbial community temperature-richness relationship depends on how rapidly the strength of effective competition between species pairs changes with temperature relative to the variance of their growth rates. Furthermore, it predicts that a thermal specialist-generalist tradeoff in growth rates alters coexistence by shifting this balance, causing richness to peak at relatively higher temperatures. Finally, we show that the observed patterns of variation in thermal performance curves of metabolic traits across extant bacterial taxa is indeed sufficient to generate the variety of community-level temperature-richness responses observed in the real world. Our results provide a new and general mechanism that can help explain temperature-diversity gradients in microbial communities, and provide a quantitative framework for interlinking variation in the thermal physiology of microbial species to their community-level diversity.