Next Article in Journal
Novel Treatment for Mitigating Condensate Bank Using a Newly Synthesized Gemini Surfactant
Next Article in Special Issue
Optimisation of Sporosori Purification and Protein Extraction Techniques for the Biotrophic Protozoan Plant Pathogen Spongospora subterranea
Previous Article in Journal
In Vitro Interaction of Organophosphono- and Organophosphorothioates with Human Acetylcholinesterase
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 25
Issue 13
10.3390/molecules25133028
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomic Profiling of Emiliania huxleyi Using a Three-Dimensional Separation Method Combined with Tandem Mass Spectrometry

by
Goyeun Yun
1,†,
Jong-Moon Park
1,2,†,
Van-An Duong
1,†,
Jeong-Hun Mok
1,
Jongho Jeon
1,
Onyou Nam
3,
Joonwon Lee
4,
EonSeon Jin
3,* and
Hookeun Lee
1,*
1
College of Pharmacy, Gachon University, Incheon 21936, Korea
2
Basilbiotech, Seoul 06621, Korea
3
Department of Life Science, Hanyang University, Seoul 04763, Korea
4
College of Letters and Science, University of California Los Angeles, Los Angeles, CA 90095, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Molecules 2020, 25(13), 3028; https://doi.org/10.3390/molecules25133028
Submission received: 29 May 2020 / Revised: 30 June 2020 / Accepted: 30 June 2020 / Published: 2 July 2020
(This article belongs to the Special Issue Advancements in Analytical Techniques for Proteomics)
Browse Figures
Versions Notes

Abstract

:
Emiliania huxleyi is one of the most abundant marine planktons, and it has a crucial feature in the carbon cycle. However, proteomic analyses of Emiliania huxleyi have not been done extensively. In this study, a three-dimensional liquid chromatography (3D-LC) system consisting of strong cation exchange, high- and low-pH reversed-phase liquid chromatography was established for in-depth proteomic profiling of Emiliania huxleyi. From tryptic proteome digest, 70 fractions were generated and analyzed using liquid chromatography-tandem mass spectrometry. In total, more than 84,000 unique peptides and 10,000 proteins groups were identified with a false discovery rate of ≤0.01. The physicochemical properties of the identified peptides were evaluated. Using ClueGO, approximately 700 gene ontology terms and 15 pathways were defined from the identified protein groups with p-value ≤0.05, covering a wide range of biological processes, cellular components, and molecular functions. Many biological processes associated with CO2 fixation, photosynthesis, biosynthesis, and metabolic process were identified. Various molecular functions relating to protein binding and enzyme activities were also found. The 3D-LC strategy is a powerful approach for comparative proteomic studies on Emiliania huxleyi to reveal changes in its protein level and related mechanism.

Graphical Abstract

1. Introduction

Emiliania huxleyi (E. huxleyi), a member of the coccolithophores group (calcareous alga), is a unicellular marine phytoplankton that can be found throughout the ocean with a unique morphology [1]. It has been known to play a critical role in global biogeochemistry through the uptake of CO2 by photosynthesis and calcification [2]. Calcification is a biogenic process of calcium carbonate production in E. huxleyi and other coccolithophores. During calcification, E. huxleyi produces mineral plates consisting of calcium carbonate on exoskeletons known as coccoliths [3]. The sedimentation of these compartments in the ocean contributes negatively to ocean acidification [4]. E. huxleyi is able to form massive blooms covering the ocean surface up to 250,000 km2, which are sometimes detected by satellite. It reflects the amount of CO2 absorbed by the ocean and is closely related to climate changes on Earth [5]. As the biogeochemical and ecological importance of E. huxleyi has been recognized, various physiological, biochemical, and genetic studies have been extensively carried out with this species [6,7,8,9,10]. The genome sequence database of E. huxleyi (strain CCMP 1516) was previously constructed in 2013, which consisted of 30,569 protein-coding genes [11]. Nevertheless, proteomic analysis of this species has just been reported in a few studies with a small number of identified proteins. Jones et al. identified 99 proteins from E. huxleyi (strain NZEH) using one-dimensional SDS-PAGE and liquid chromatography-tandem mass spectrometry (LC–MS/MS) [12]. This group later utilized two-dimensional liquid chromatography (2D-LC) and identified 115 homologous protein groups from the same strain [13]. Another group used LC-MS/MS to identify 346 to 500 proteins from E. huxleyi (strain CCMP 1516) [14,15]. Thus, it is necessary to perform a proteomic profiling study on E. huxleyi to identify a large dataset of its protein identification.
Proteomics is the study of the entire proteins (proteome) in a sample. The analysis of highly complex samples has been a main technical task in proteomic research aiming at genome-wide analysis with the identification of low-abundance proteins [16]. Various approaches have attempted to achieve complete proteome coverage of complex samples, and yet reducing sample complexity remains a bottleneck against reaching a fundamental goal in proteomics [17]. In proteomic studies, separations in protein or peptide levels are frequently used to reduce sample complexity prior to mass spectrometric analysis [18,19,20,21]. Many separation techniques have been widely used in proteomic studies, including two-dimensional electrophoresis (2-DE) [22], reversed-phase liquid chromatography (RPLC) [23], isoelectric focusing [24], and capillary zone electrophoresis (CZE) [25]. To improve the identification of peptides and proteins, several multi-dimensional separation methods have been developed and evaluated [26,27,28]. Multi-dimensional separation is currently considered as one of the most powerful approaches to increase peak capacity [29,30,31]. The concept of multi-dimensional separation was described by Giddings in 1984, which stated that two or more independent separation methods could be coupled based on the orthogonality in elution mechanisms to resolve complex mixtures [32]. Based on that, the first Multi-dimensional Protein Identification Technology, an online proteomic technique packing strong cation exchange (SCX) and RP resins into a single capillary, was successfully developed [23,33]. Recently, an integrated spintip consisting of SCX beads and RP disk in one pipet tip was developed for deep proteomic profiling [34]. Multi-dimensional separation is classified into three main strategies [35,36]. The classic approach starts with the separation of proteins by 2-DE or LC prior to enzymatic digestion and LC-MS/MS analysis of the peptide digest [37]. The second strategy is applied in top-down proteomics, which allows the multi-dimensional separation in protein level, followed by MS/MS analysis [38]. The third method is widely used in bottom-up proteomics studies, which first digests the proteins into peptides prior to multi-dimensional separation and MS/MS analysis [39]. In bottom-up proteomics, SCX-RPLC has been widely used due to their high orthogonality, resulting in better separation and resolving power [40,41,42]. Besides, the combination of two RPLC under extremely different pH conditions showed the highest peak capacity among various chromatographic combinations [43,44]. The effectiveness of two-dimensional (2D) separation with high-pH (HpH) and low-pH (LpH) RPLC was confirmed through the increases in peptide and protein identification [45,46]. Some three-dimensional LC (3D-LC) platforms have been demonstrated. Wei et al. developed a 3D (RPLC-SCX-RPLC) system to analyze proteins in yeast and identified 5954 unique peptides and 1457 proteins [47]. Similar strategies were used for proteomic profiling of monkey brain tissue [48], human hepatocellular carcinoma tissues [49], and plasma of patients with sepsis and systemic inflammatory response syndrome [50]. Betancourt et al. combined SCX-HpH RPLC-LpH RPLC to identify more than 5000 proteins in mouse embryonic fibroblast cells [51]. Besides, several combinations of 3D separation methods have been established, including isoelectric focusing-SCX-RPLC [52], RPLC-strong anion exchange-RPLC [53], electrostatic repulsion hydrophilic interaction chromatography-HpH RPLC-LpH RPLC [54], and SCX-RPLC-CZE [55]. Recently, Spicer et al. developed a 3D system consisting of three consecutive RPLC, which identified more than 14,000 proteins across 126 fractions [56].
In this study, in-depth proteomic profiling of E. huxleyi (CCMP371) was performed using a 3D-LC system. The 3D-LC strategy consisted of SCX and HpH RPLC fractionation, followed by LpH RPLC separation and tandem mass spectrometry (MS/MS) analysis. Seventy SCX-HpH RPLC fractions were generated from E. huxleyi proteome digest. Peptide and protein identification was performed using Trans-Proteomics Pipeline (TPP). The physicochemical properties of the identified peptides and proteins were evaluated. In addition, the identified proteins were used to define functional classifications based on gene ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway through ClueGO.

2. Results

2.1. Design of an Off-line 3D-LC (SCX-HpH RPLC-LpH RPLC) System

In this study, an off-line 3D-LC separation system was coupled with the Q Exactive™ Hybrid Quadrupole-Orbitrap mass spectrometer to identify a large number of proteins from E. huxleyi cell lysates (Figure 1). The proteome digest of E. huxleyi cell lysates was first fractionated into 14 fractions (C to P) using SCX. After that, each of them was further separated into 16 fractions using HpH RPLC. A set of 16 fractions generated from one SCX fraction was concatenated into five final fractions (1-6-11, 2-7-12, 3-8-13, 4-9-14, and 5-10-15-16). Combination of fractions with different hydrophobicity in the HpH RP dimension could reduce the LpH RPLC–MS/MS analysis time without a significant reduction in protein identification [43]. Overall, 70 fractions were generated from the original proteome digest of E. huxleyi cell lysates for LC–MS/MS analyses.

2.2. Proteomic Analysis of E. huxleyi Using 3D-LC System

2.2.1. Identification of Peptide and Protein Group

The proteome of E. huxleyi was searched against a Uniprot database using TPP. After MS/MS search with Comet, the data were further analyzed with PeptideProphet and ProteinProphet. The list of peptides identified by PeptideProphet at a false discovery rate (FDR) ≤0.01 is shown in Table S1. The term “protein” means a distinct protein, whereas the “protein group” indicates indistinguishable proteins (with shared peptides) grouped by ProteinProphet [57]. For simplification, “protein group” was used for both proteins and protein groups. Figure 2 presents the results of E. huxleyi proteomic profiling using the developed 3D-LC system. In total, 84,753 unique peptides and 15,331 protein groups (single hit included) were identified (Figure 2a). From 35,707 protein entries in the database, 43% of protein coverage was obtained (single hit included). The identification of peptides and protein groups gradually increased with the number of fractions. However, the increase rate reduced in the last fractions (~50–70), most likely due to the fraction-to-fraction overlap of the peptides, although this overlap was relatively small (Figure 2c). Figure 2b exhibits the distribution of protein groups by the number of unique peptides. When single hits were excluded, the number of protein groups was 10,145. It can be observed that over 90% of the identified peptides constrained within one to two fractions, and approximately 96% of them spanned in one to three fractions (Figure 2c). This finding demonstrates minimal fraction-to-fraction overlap, as well as the sufficient power of the SCX-HpH RPLC fractionation. The effectiveness of tryptic digestion was demonstrated in Figure 2d,e. Nearly 90% of the peptides were identified with two tryptic termini, whereas approximately 7% of the peptides contained one to two missing cleavages.

2.2.2. Physicochemical Properties of Identified Peptides

The physicochemical properties of the identified peptides from the 3D-LC system were also evaluated. Retention time (Rt), molecular weight (MW), and isoelectric point (pI) of peptides were obtained using TPP. The pI value is the pH where the peptide carries no net charge. At a pH below pI, the peptide has a positive charge and vice versa. The grand average of hydropathy (GRAVY) values of peptides were calculated using the GRAVY calculator (http://gravy-calculator.de/) [55]. Negative GRAVY values signify hydrophilicity, and positive values indicate hydrophobicity. Figure 3 presents the cumulative distribution of MW, pI, and GRAVY of identified peptides. The majority of the identified peptides have MW in the range 800–2000 Da (86%) and GRAVY <0 (72%). Approximately 78% of the peptides have pI in the range 3.5–7.0. The distribution of pI is similar to a previous report using SCX-LpH RPLC [42]. We constructed various two-dimensional graphs to elucidate the relationship among different parameters of the identified peptides (Figure S1). We can observe that the majority of high MW peptides with MW >2000 Da was eluted at a later time in LpH RPLC (Rt of 60 to 130 min) (Figure S1a, Supplementary Materials). They mostly had GRAVY of −1.5~+1 (Figure S1d) and pI of 3 to 7 (Figure S1e). Acidic peptides seemed to appear throughout LpH RPLC running time, whereas basic peptides, particularly those with pI >11, mostly presented at 30 to 90 min (Figure S1b). As shown in Figure S1c, highly hydrophilic peptides were likely associated with a short Rt, whereas hydrophobic peptides were eluted throughout the LpH RPLC running time. GRAVY and pI do not show any distinct relationship (Figure S1f).

2.2.3. Effects of Strong Cation Exchange (SCX) on Peptide Separation

The first separation step (SCX) generated 14 fractions (which were named as C–P) corresponding to the increase in the KCl concentration of the elution buffer. The identified peptides were grouped by SCX fractions (overlap peptides were kept in each SCX fraction), and various box-plots were constructed to elucidate the effect of SCX (i.e., KCl concentration) on the separation of peptides (Figure S2). It is evident that SCX did not considerably separate the peptides regarding their MW and GRAVY (Figure S2a,d). However, as shown in Figure S2b, peptides with low pI values were eluted at a low salt concentration and vice versa. This finding is in agreement with the theoretical characteristics of SCX [41]. In SCX, electrostatic interaction mainly affects the retention of peptides, whereas hydrophobic interactions only play a modest role [58]. In this study, SCX was operated at a low pH (2.7) to reduce the dissociation of carboxylic groups of peptides and thereby promote interactions between the sulfonate groups of SCX resins and the protonated basic amino acid residues [40]. Under pH 2.7, peptides with pI >2.7 were positively charged and bound to the aliphatic sulfonic acid groups of SCX resins. The peptides with higher pI values bound stronger to the stationary phase than those with lower pI values [41]. A salt gradient (KCl) was used to separate the peptides depending on their charge. Peptides with pI values closer to 2.7 would be eluted first when applying a lower ionic strength. As the salt concentration was higher, peptides with higher pI would be eluted. Besides, Figure S2c shows that peptides corresponding to a low salt strength in SCX (eluted at a low Rt) seemed to have higher Rt in LpH RPLC. It suggests that SCX and LpH RPLC have a high orthogonality.

2.2.4. Protein Identification Using ProteinProphet

The list of protein groups identified by ProteinProphet at an FDR ≤0.01 is shown in Table S2. A total of 10,448 protein groups were identified with at least two unique peptides. Figure 4 exhibits the distribution of the protein groups by their length and protein sequence coverage (%). Most of the protein groups had <600 amino acids (80%). Additionally, about 79% of the protein group had a sequence coverage percentage of 10 to 50%. Among more than 10,000 identified proteins groups, only 1770 protein groups (~17%) have been characterized, whereas the majority were defined as uncharacterized proteins in the Uniprot database.

2.3. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG)

ClueGO add-in (Cytoscape) was used to categorize the identified proteins/protein groups based on GO and KEGG pathways. The database of E. huxleyi contained 2802 biological processes (9885 genes), 509 cellular components (8596 genes), 1452 molecular functions (12,470 genes), and 104 KEGG pathways (3259 genes), which was updated on 31 October, 2019. In total, 1454 biological processes (51.9%), 235 cellular components (46.2%), 678 molecular functions (46.7%), and 100 KEGG pathways (96.2%) were found. After filtering with a p-value of ≤0.05, 15 KEGG pathways were defined and listed in Table 1. The percentages of associated genes vary from 21.1 to 83.9%. Some KEGG pathways were identified with high percentages of associated genes, such as basal transcription factors, citrate cycle (TCA cycle), proteasome, spliceosome, mRNA surveillance pathway, protein processing in endoplasmic reticulum, and amino sugar and nucleotide sugar metabolism (≥67%). In addition, after filtering with a p-value of ≤0.05, the identified protein groups are involved in 395 biological processes, 110 cellular components, and 181 molecular functions (Table S3). Top 10 GO terms with highest −log10(p-value) were presented in Figure 5a–c. The top 10 biological processes are biosynthetic processes, metabolic processes, translation, protein modification, and protein ubiquitination. The top 10 cellular components are components of membrane, nucleus, ribosome, and different organelles. The top 10 molecular functions are transferase, ligase, and kinase activity, and different binding functions. Figure 5d classifies GO terms regarding the percentage of associated genes. It is evident that most of the GO terms have 50 to 80% associated genes.

3. Discussion

Various 2D-LC and few 3D-LC separation systems have been applied to proteomic researches for decades to solve the under-sampling issue (i.e., the complexity of peptide mixtures is too high to be fully detected by mass spectrometry). The developed 3D-LC system can minimize this limitation by increasing peptide and protein identification. This system is based on the high orthogonality between SCX and RPLCs as well as the efficacy of two RPLC with an extreme pH difference [33,44,45]. SCX was performed as the first-dimensional separation to separate peptides based on their charges [55,59]. Then, peptides were fractionated by HpH RPLC with respect to their hydrophobicity [60]. Lastly, LC-MS/MS with a low-pH solvent was performed as the third separation. The orthogonality between HpH RPLC and LpH RPLC was due to the pH of mobile phases, as previously reported. The pH affected the net charge of peptides according to the pKa values of their ionizable amino acid groups. Thus, basic peptides were more retained at high pH, while acidic peptides were retained at low pH [43,44]. The developed sequential 3D-LC platform was performed in an off-line manner to overcome the limitation of small loading sample amounts, an obvious drawback of online multi-dimensional separation approaches [29]. Besides, an off-line fractionation method could dramatically improve separation power, since it is highly flexible and easy to manipulate between separations, and thus each dimensional separation can be optimized by varying buffer and elution conditions [29,48,61]. Another improvement of our 3D-LC system is the arrangement of the three different LC systems. SCX has a high sample loading capacity serving as the first separation. The desalting step is often required after SCX to remove salts in samples, which can reduce sensitivity or cause blockage in the LC-MS/MS column [62]. However, the C18 column of HpH RPLC following SCX can serve as the desalting column for SCX-fractions. Therefore, this 3D-LC system reduces additional steps after SCX, as previously reported by Betancourt et al. [51]. This procedure, however, was complicated due to additional steps (acylation of amino groups before SCX fractionation and regeneration of amino groups after that). In addition, the concatenation of fractions in the second separation dimension reduces the total number of SCX-HpH RPLC from 224 to 70, which is quite moderate for LpH RPLC–MS/MS analysis [43]. Seventy fractions were analyzed by the LpH RPLC–MS/MS in 175 h (150 min per analysis). It is relatively long (~1 week), as an inevitable result of a large fraction number in a 3D-LC approach, but it is necessary for a deep proteomic profiling study [56]. Depending on the biological question of each study, the number of SCX-HpH RPLC fractions may be adjusted, and thus, the MS/MS analysis time can also be altered. In comparative proteomics, due to the increase in the number of samples, the analysis time may increase substantially with this 3D-LC method. However, a labeling method can allow the simultaneous process of many samples, and thereby, the MS/MS analysis time remains unchanged [26].
Some previous studies investigated the proteome of E. huxleyi. A shotgun proteomic analysis identified 99 proteins from E. huxleyi (NZEH) [12]. Another study used tandem mass-spectrometry with isobaric tagging (iTRAQ) to reveal changes in protein abundance between normal and high CO2 conditions. However, due to the limited number of quantified proteins, acclimation mechanisms tolerating high CO2 amount remained unrevealed [13]. A study of McKew et al. investigated acclimation of E. huxleyi to nutrient limitation and found the marked increases in the abundance of some proteins belonging to inorganic nutrient transport and the internal remobilization of N and P compounds [15]. In this study, utilizing the 3D-LC system for the in-depth profiling of E. huxleyi (CCMP371), we identified more than 84,000 unique peptides and above 10,000 protein groups (with at least two unique peptides) at an FDR of ≤0.01. The 3D-LC system was not evaluated regarding the reproducibility. It could be carried out by dividing the proteome digest into two to three parts after digestion and simultaneously performing the peptide fractionation and LC-MS/MS analysis. Generally, 2D- and 3D-LC systems suffer low reproducibility in quantification. However, these issues have been solved in recent studies using different labeling methods, which allow the simultaneous handling of many samples [26]. This 3D-LC approach may also present a low reproducibility. Nevertheless, it could be used in comparative proteomic studies to reveal changes in protein levels and define related mechanisms under different conditions. Some labeling methods may necessarily be incorporated to allow the implementation of many samples under identical conditions, such as metabolic labeling [63], stable isotope labeling using amino acids in cell culture (SILAC) [64], tandem mass tags (TMT) [65], and iTRAQ [66].
Some of the identified protein groups have been defined as playing critical roles in various processes of E. huxleyi. For example, proteins involved in carbon transport, pH homeostasis, and biomineralization were found, such as carbonic anhydrases, vacuolar H+ ATPase, plasma membrane ATPase, H+/Ca2+ exchanger, Na+/Ca2+-K+ exchanger, Ca2+ ATPase, and cation/Ca2+ exchanger. The identification of these proteins in a proteomic study could allow a direct measurement of their abundance or activity, giving a better insight into the cellular changes under different conditions compared to the measurement of transcripts [67]. Carbonic anhydrases are able to catalyze the interconversion water and CO2 with H+ and bicarbonate [12]. The regulation of carbonic anhydrases is linked to biomineralization in coccolithophores [68]. Vacuolar H+ ATPase and clathrin are related to calcifying vesicle membranes. Some of the ATPase subunits identified in this study may be a part of a complex in E. huxleyi membranes, which maintains an alkali pH together with the formation of CaCO3 [12]. Clathrin-coated vesicles contain vacuolar H+ ATPase as previously reported [69]. These vesicles participate in protein and lipid transportation, as well as some membrane and trans-Golgi network trafficking pathways [70]. In E. huxleyi, coccolith formation occurs in Golgi-derived coccolith vesicles, and clathrin-coated vesicles containing vacuolar H+ ATPase play a role in the alkalization for coccolith formation [12,71]. Under ocean acidification, calcifying marine organisms tend to shift the interconversion to increase H+ and bicarbonate concentration. A proteomic approach may help to reveal proteome changes and related mechanisms. In E. huxleyi (strain NZEH), under high CO2 level, abundances of ribosomal proteins (30S ribosomal protein S7) and histones (H2A, H3, and H4) decreased, indicating a reduction in DNA and chromatin synthesis [13]. Besides calcification, biomineralization also includes silicification. In some coccolithophores, diatom-like silicon transporters are present for silicification, as Si also plays a crucial role in the formation of calcite coccoliths. However, E. huxleyi is one of the species that lack the requirement for Si [72]. The proteomic results in our study also confirmed the absence of any Si-related proteins in E. huxleyi. In cellular component, adaptor protein (AP)-type membrane coat adaptor complex was identified, specifically subunits of AP1 and AP2 complexes. Five different AP complexes (AP1-5) relate to biomineralization by transporting vesicles intracellularly. It was reported that a partial or complete loss of AP3 and AP5 in the gene-level has occurred in three coccolithophores (Isochrysis galbana, Gephyrocapsa oceanica, and E. huxleyi strain Van 556, 92A, EH2, and CCMP1516) as a consequence of evolutionary events [73]. Similarly, AP3 and AP5 were not identified in E. huxleyi (strain CCMP371) in this study, suggesting their loss or low abundances.
Various identified proteins groups are involved in photosynthesis (77), photosynthesis, light-harvesting (60), and photosynthesis, light reaction (68). In addition, chloroplast (78 protein groups), chloroplast thylakoid (62 protein groups), and chloroplast thylakoid membrane (62 protein groups) were defined, which are essential cellular components that mediate the effect of light on calcification in E. huxleyi [74]. A previous proteomic study found 49 light-harvesting complex proteins and 12 photosynthetic electron transfer chain proteins. The abundances of these proteins varied depending on suboptimal or supra-optimal light. Proteins in photosystem I and II increased their abundance in low light, whereas light-harvesting fucoxanthin–chlorophyll proteins and photoprotective LI818 proteins up-regulated in high light [14]. Some transport processes are involved in the transport of protein, peptide, amine, nitrogen compound, organic substance, Golgi vesicle, xenobiotic, cation, anion, ammonium, carbohydrate, metal ion, RNA, and nucleic acid. Upon changes in environmental conditions, the relevant transport proteins may show abundance changes. For example, during periods of phosphate starvation, proteins involved in phosphorus transport are more abundant and vice versa [15,75]. Similarly, upon the limitation of nitrogen-containing compounds, various nitrogen transporters are up-regulated [15]. In addition, the nitrogen starvation also leads to the up-regulation of glutamine synthetase, a protein frequently used by marine bacterioplankton to dissolve organic or inorganic nitrogen [76,77]. In this study, glutamine synthetase was also detected, which related to the assimilation of ammonium into amino acids.
Many metabolic processes relating to peptide, amide, carboxylic acid, organic acid, amino acid, cellular nitrogen compound, cellular carbohydrate, RNA, ribose phosphate, ribonucleotide, nucleoside phosphate, heme, and glutathione were identified. Biosynthesis of nucleoside, ribonucleoside, peptidyl-diphthamide, DNA, glycosyl compound, quinone, and inositol phosphate were also listed in biological processes. These results can favor further studies to identify some molecules in E. huxleyi, such as primary/secondary metabolites and bioactive molecules that can be potentially applied in pharmaceutical, cosmeceutical, and nutraceutical, like other marine organisms [78]. The glycosyl compound-quinone biosynthetic process and glutathione metabolic process were defined. They hold potential for the discovery of biopharmaceuticals and antioxidant cosmeceuticals. Previously, several secondary metabolite molecules from marine cyanobacteria were defined as potential anticancer biomolecules [79]. Phlorotannins from brown algae have hypopigmentation effects, which can be applied in the cosmetic industry. Docosahexaenoic acid-rich oils derived from Thraustochytrids are currently available on the market as a dietary supplement [80]. In addition, various molecular functions relating to binding (proteins, l-ascorbic acid, vitamin B6, and iron ion) and enzyme activities (transferase, kinase, peptidase, lipase, and phosphatase) were identified. Vitamins and minerals in E. huxleyi can possibly be exploited as a useful source of food ingredients or nutraceutical supplementary. This study only revealed the proteome of E. huxleyi, whereas the metabolome was not accessed, and thus, many critical features relating to metabolites in this species remained unknown. Several studies have performed an integrated metabolomics-proteomics, which could show a better landscape of the marine plankton [81,82]. Integrating “-omics” studies (e.g., genomics, proteomics, transcriptomics, metabolomics, lipidomics, and glycomics) will enable the investigation of the impacts of ecological interactions on organism physiology as well as evaluation of physiological and biochemical properties relating to human health for future applications [83,84].

4. Materials and Methods

4.1. Materials

E. huxleyi strain CCMP371 was purchased from the Provasoli-Guillard National Center for Marine Algae and Microbiota (NCMA, East Boothbay, ME, USA). Tris (2-carboxyethyl)phosphine (TCEP) was supplied by Thermo Fisher Scientific (Rockford, IL, USA). Formic acid, iodoacetamide (IAA), ammonium bicarbonate (ABC), and potassium chloride were obtained from Sigma-Aldrich (St. Louis, MO, USA). Trypsin was received from Promega (Madison, WI, USA). HPLC-grade water and acetonitrile were purchased from J.T. Baker (Phillipsburg, NJ, USA). All the chemicals were used as received without further purification.

4.2. Cell Culture

E. huxleyi (CCMP371) cells were cultured in 500 mL-Erlenmeyer flasks under 12-h light/dark cycles at 20 °C. The media contained sterile artificial seawater enriched with nitrates, phosphates, trace metals, vitamins at f/2 concentration, and selenium [85,86]. The irradiance of fluorescent lamps was set at 50 μmol m−2 s−1.

4.3. Sample Preparation

4.3.1. Protein Extraction

The cell pellets were placed in Maintainor® Tissue cards and stabilized at 95 °C in StabilizorTM T1 (Denator, Gothenburg, Sweden) for denaturing. The samples were then loaded into pre-chilled TT1 tissue TUBETM (Covaris, Woburn, MA, USA), frozen in liquid nitrogen followed by pulverization using CryoPrep® (Covaris). Proteins were extracted in lysis buffer (8 M urea and 0.1 M Tris-HCl, pH 8.5) under sonication (12 min, 18 °C, Covaris). Four volumes of acetone (−20 °C) were added to the sample, followed by incubation at −20 °C overnight for protein precipitation. After centrifugation (4000 rpm, 4 °C, 10 min, Centrifuge 5810 R, Eppendorf, Hamburg, Germany), the supernatant was discarded, and the protein pellets were dried in a ScanSpeed 40 centrifugal evaporator (1800 rpm, 3 h, Labogene, Lillerød, Denmark). The purified proteins were re-suspended in the lysis buffer followed by protein quantification using the Pierce BCA Protein Assay kit (Thermo Fisher Scientific, Rockford, IL, USA).

4.3.2. Automated Filter-Aided Sample Preparation (FASP)

A liquid handling robotic system (Agilent Technologies, Santa Clara, CA, USA) controlled by VWorks software was used to perform the automated FASP [87,88]. Multiscreen Vacuum Manifold™ (Millipore, Billerica, MA, USA) was used to supply vacuum during sample preparation. In brief, protein samples were first loaded onto 50 wells of a 96-well plate (AcroPrep Advance Filter Plates, 350 µL, Omega 30K MWCO) with 100 µg protein/well. Proteins were reduced with 100 µL of 5 mM TCEP (37 °C, 30 min), alkylated with 100 µL of 50 mM IAA (25 °C, 30 min, in the dark), and digested with trypsin at an enzyme: protein ratio of 1:50 (w/w) (37 °C, 18 h). The digestion was stopped by reducing the pH of the samples to 2–3 with formic acid. The number of samples was reduced from 50 to 5 by pooling before desalting.

4.3.3. Desalting with Reverse-Phase Solid-Phase Extraction

Sample desalting was conducted using Sep-Pak® Vac 1cc C18 cartridge (Waters, Milford, MA, USA). The cartridge was washed with 1 mL of methanol followed by 1 mL of solvent A1 (0.1% formic acid in water), equilibrated with 1 mL of solvent B1 (0.1% formic acid in acetonitrile: water (80:20, v/v)), and then washed with 1 mL of solvent A1. The digested samples were loaded onto the column and washed with 500 μL of solvent A1. The peptides were eluted with 1 mL of solvent B1 and subsequently dried in the ScanSpeed 40 centrifugal evaporator (1800 rpm, 3 h). Samples were dissolved in solvent A1 and pooled to one prior to fractionation.

4.3.4. SCX Fractionation

The peptide mixtures were fractionated by SCX chromatography using an HPLC system consisting of a vacuum degasser, a binary pump, an autosampler (Agilent series 1200, Agilent Technologies), and a UV detector (Agilent series 1100, Agilent Technologies). Sample separation was carried out on a Polysulfoethyl ATM column (4.6 × 250 mm, 5 μm particles, 100-Å pores, PolyLC, Columbia, MD, USA) using a flow rate of 0.5 mL/min. The mobile phase solvent consisted of A2 (10 mM KH2PO4 in acetonitrile: water (20:80, v/v), pH 2.7) and B2 (10 mM KH2PO4 and 0.6 M KCl in acetonitrile: water (20:80, v/v), pH 2.7). The column was equilibrated with solvent A2 for 30 min. The gradient was applied as follows for solvent B2: 0–10% for 2 min, 10–80% for 68 min, 80–100% for 5 min, and holding at 100% for 15 min. The column was re-equilibrated using solvent A2 for 30 min and stored at −4 °C. The eluting peptides were monitored at 214 nm and collected automatically into tubes at 3 min intervals. The eluted samples were dried in the ScanSpeed 40 centrifugal evaporator (1800 rpm, 3 h). According to the peptide concentration of each fraction, some fractions were unchanged, whereas several two to three adjacent fractions were pooled to make 14 SCX-fractions (C to P).

4.3.5. HpH RPLC Fractionation

HpH RPLC was carried out using Sep-Pak® Vac 1cc tC18 cartridges (Waters, Milford, MA, USA). The column was washed in the following steps: 1 mL of 100% methanol, 1 mL of solvent A1 (0.1% formic acid in water), 1 mL of solvent B1 (0.1% formic acid in acetonitrile: water (80:20, v/v)), and 1 mL of solvent A1. Peptide mixtures in each SCX fraction were loaded onto the column and washed with 0.5 mL of solvent A1. Peptides were eluted with 16 eluting solvents containing 0.1% trimethylamine in different acetonitrile: water mixtures. The amount of acetonitrile was varied as follows: 2.5, 5, 6.25, 7.5, 8.75, 10, 11.25, 12.5, 13.75, 15, 16.25, 17.5, 18.75, 20, 35, and 50% (v/v). Sixteen fractions (numbered from 1 to 16) were pooled to make 5 HpH RPLC fractions (1-6-11, 2-7-12, 3-8-13, 4-9-14, and 5-10-15-16) [43]. The obtained fractions were labeled, e.g., C1, C2, C3, C4, and C5 were the five HpH RPLC fractions derived from SCX-fraction C. The samples were then dried in the ScanSpeed 40 centrifugal evaporator (1800 rpm, 3 h).

4.4. LC–MS/MS Analysis

The dried peptides were reconstituted in 50 μL of solvent A1 (0.1% formic acid in water) and injected (1 μg) for analysis. The LC-MS/MS system consisted of a Dionex Ultimate 3000 HPLC coupled with a Q Exactive™ Hybrid Quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific, Rockford, IL, USA). The peptides were loaded onto an Acclaim™ PepMap™ 100 C18 nano-trap column (75 μm × 2 cm, 3 μm particles, 100-Å pores, Thermo Fisher Scientific) using solvent A1 (0.1% formic acid in water) at a flow rate of 2.5 μL/min for 5 min. Then, the peptides were separated in an Acclaim™ PepMap™ C18 100A RSLC nano-column (75 μm × 50 cm, 2 μm particles, 100-Å pores, Thermo Fisher Scientific) at a flow rate of 300 nL/min. The mobile phase solvent consisted of A3 (0.1% formic acid in water) and B3 (0.1% formic acid in acetonitrile: water (90:10, v/v)). The gradient was set up as follows for solvent B3: equilibration at 4% for 10 min, 4–40% for 110 min, 40–96% for 0.1 min, holding at 96% for 9.9 min, 96–4% for 0.1 min, and holding at 4% for 19.9 min for re-equilibration of the column. The following parameters were used: spray voltage, 2.2 kV; capillary temperature, 320 °C; isolation width, ± 2 m/z; scan range, 400–2000 m/z; resolution in full-MS scans, 70,000; and resolution in MS/MS scans at 200 m/z, 17,500. The MS was operated using a data-dependent acquisition method. Top ten precursor ions with the highest intensity were isolated in the quadrupole and fragmented by the higher-energy collisional dissociation with 27% normalized collisional energy. Dynamic exclusion was set at 20 s to minimize the repeated analyses of the same abundant precursor ions.

4.5. Data Analysis

The proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository [89], with the dataset identifier PXD018511. Raw MS/MS data files were converted to. mzXML format by MSConvert. Comet (version 2017.01 rev.0) was used to search MS/MS spectra against a database of E. huxleyi (CCMP371) obtained from Uniprot. The draft release of E. huxleyi genome assembly version 1.0 includes a total of 39,126 predicted gene models and functional annotations using the JGI annotation pipeline [11]. The best gene model at each locus, including diploid alleles, was the filtered set selected by the JGI annotation pipeline. The following parameters were set for the search: maximum of two missed cleavages with trypsin; semitryptic cleavage; 10 ppm and 0.02 Da tolerances of precursor ion masses and fragment ion mass, respectively; static carbamidomethylation of cysteine; and variable modifications including methionine oxidation (+15.995 Da) and carbamylation of protein in N-term (+43.0006 Da). The search result files in pepXML format were transported to the TPP version 5.1.0 [90], and PeptideProphet [91] and ProteinProphet [57] were operated. Peptides and proteins were filtered at an FDR of ≤0.01. The data were processed and visualized using MS. Excel 2016. The boxplots were generated using R version 3.6.1. GO [92] and KEGG pathway [93] were categorized using Cytoscape version 3.7.1 (National Institute of General Medical Sciences, MD, USA) via ClueGO version 2.5.4 (Cordeliers Research Center, France) [94]. The GO terms and KEGG pathways were filtered at a p-value of ≤0.05.

5. Conclusions

In this study, a 3D-LC separation method consisting of SCX-HpH RPLC-LpH RPLC was developed for proteomic profiling of E. huxleyi. More than 84,000 unique peptides and 10,000 protein groups were identified from 70 SCX-HpH RPLC fractions. Approximately 700 GO terms and 15 KEGG pathways were defined from the identified protein groups, which relate to various important biological processes, cellular components, and molecular functions of E. huxleyi. The identification of E. huxleyi proteins in this study will facilitate further studies on this species, particularly those that aim to reveal proteome changes of E. huxleyi under different conditions as well as related mechanisms.

Supplementary Materials

The following are available online, Figure S1: Relationship between different parameters of identified peptides, Table S1: List of identified peptides (PeptideProphet) of Emiliania huxleyi (CCMP371) in three-dimensional separation method, Table S2: List of identified proteins/protein groups (ProteinProphet) of Emiliania huxleyi (CCMP371) in three-dimensional separation method, Table S3: Gene ontology analysis of identified proteins/protein groups from Emiliania huxleyi (CCMP371) using ClueGO. The proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD018511. Reviewers may access it using username [email protected] and password wyy3eSj5).

Author Contributions

Conceptualization, E.J. and H.L.; methodology, J.-M.P., E.J. and H.L.; validation, J.-M.P. and V.-A.D.; formal analysis, G.Y. and V.-A.D.; investigation, G.Y., J.-H.M., J.J., O.N. and J.L.; resources, J.-M.P. and J.-H.M.; data curation, G.Y. and V.-A.D.; writing—original draft preparation, G.Y.; writing—review and editing, J.-M.P. and V.-A.D.; visualization, V.-A.D.; supervision, J.-M.P.; project administration, E.J. and H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant, funded by the Korean government (MSIT) (No. NRF-2017M3D9A1073784, NRF-2016M1A5A1027461).

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Cubillos, J.; Wright, S.; Nash, G.; De Salas, M.; Griffiths, B.; Tilbrook, B.; Poisson, A.; Hallegraeff, G. Calcification morphotypes of the coccolithophorid Emiliania huxleyi in the Southern Ocean: changes in 2001 to 2006 compared to historical data. Marine Ecol. Progress Series 2007, 348, 47–54. [Google Scholar] [CrossRef]
  2. Iglesias-Rodríguez, M.D.; Brown, C.W.; Doney, S.C.; Kleypas, J.; Kolber, D.; Kolber, Z.; Hayes, P.K.; Falkowski, P.G. Representing key phytoplankton functional groups in ocean carbon cycle models: Coccolithophorids. Global Biogeochem. Cy. 2002, 16. [Google Scholar] [CrossRef]
  3. Obata, T.; Shiraiwa, Y. A novel eukaryotic selenoprotein in the haptophyte alga Emiliania huxleyi. J. Biol. Chem. 2005, 280, 18462–18468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tyrrell, T.; Merico, A. Emiliania huxleyi: Bloom observations and the conditions that induce them. In Coccolithophores: From Molecular Processes to Global Impact; Thierstein, H.R., Young, J.R., Eds.; Springer: Berlin/Heidelberg, Germany, 2004; pp. 75–97. [Google Scholar]
  5. Hallegraeff, G.M. Ocean climate change, phytoplankton community responses, and harmful algal blooms: A formidable predictive challenge. J. Phycol. 2010, 46, 220–235. [Google Scholar] [CrossRef]
  6. Medlin, L.K.; Barker, G.L.A.; Campbell, L.; Green, J.C.; Hayes, P.K.; Marie, D.; Wrieden, S.; Vaulot, D. Genetic characterisation of Emiliania huxleyi (Haptophyta). J. Marine Sys. 1996, 9, 13–31. [Google Scholar] [CrossRef] [Green Version]
  7. Wahlund, T.M.; Hadaegh, A.R.; Clark, R.; Nguyen, B.; Fanelli, M.; Read, B.A. Analysis of Expressed Sequence Tags from Calcifying Cells of Marine Coccolithophorid (Emiliania huxleyi). Marine Biotechnol. 2004, 6, 278–290. [Google Scholar] [CrossRef]
  8. Dyhrman, S.T.; Haley, S.T.; Birkeland, S.R.; Wurch, L.L.; Cipriano, M.J.; McArthur, A.G. Long Serial Analysis of Gene Expression for Gene Discovery and Transcriptome Profiling in the Widespread Marine Coccolithophore Emiliania huxleyi. Appl. Environ. Microbiol. 2006, 72, 252–260. [Google Scholar] [CrossRef] [Green Version]
  9. Quinn, P.; Bowers, R.M.; Zhang, X.; Wahlund, T.M.; Fanelli, M.A.; Olszova, D.; Read, B.A. cDNA Microarrays as a Tool for Identification of Biomineralization Proteins in the Coccolithophorid Emiliania huxleyi (Haptophyta). Appl. Environmen. Microbiol. 2006, 72, 5512–5526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Richier, S.; Kerros, M.-E.; de Vargas, C.; Haramaty, L.; Falkowski, P.G.; Gattuso, J.-P. Light-Dependent Transcriptional Regulation of Genes of Biogeochemical Interest in the Diploid and Haploid Life Cycle Stages of Emiliania huxleyi. Appl. Environ. Microbiol. 2009, 75, 3366–3369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Read, B.A.; Kegel, J.; Klute, M.J.; Kuo, A.; Lefebvre, S.C.; Maumus, F.; Mayer, C.; Miller, J.; Monier, A.; Salamov, A.; et al. Pan genome of the phytoplankton Emiliania underpins its global distribution. Nature 2013, 499, 209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Jones, B.M.; Edwards, R.J.; Skipp, P.J.; O’Connor, C.D.; Iglesias-Rodriguez, M.D. Shotgun Proteomic Analysis of Emiliania huxleyi, a Marine Phytoplankton Species of Major Biogeochemical Importance. Marine Biotechnol. 2011, 13, 496–504. [Google Scholar] [CrossRef] [PubMed]
  13. Jones, B.M.; Iglesias-Rodriguez, M.D.; Skipp, P.J.; Edwards, R.J.; Greaves, M.J.; Young, J.R.; Elderfield, H.; O’Connor, C.D. Responses of the Emiliania huxleyi Proteome to Ocean Acidification. PLoS ONE 2013, 8, e61868. [Google Scholar] [CrossRef] [PubMed]
  14. McKew, B.A.; Lefebvre, S.C.; Achterberg, E.P.; Metodieva, G.; Raines, C.A.; Metodiev, M.V.; Geider, R.J. Plasticity in the proteome of Emiliania huxleyi CCMP 1516 to extremes of light is highly targeted. New Phytol. 2013, 200, 61–73. [Google Scholar] [CrossRef] [PubMed]
  15. McKew, B.A.; Metodieva, G.; Raines, C.A.; Metodiev, M.V.; Geider, R.J. Acclimation of Emiliania huxleyi (1516) to nutrient limitation involves precise modification of the proteome to scavenge alternative sources of N and P. Environmen. Microbiol. 2015, 17, 4050–4062. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Mamontova, T.; Afonin, A.M.; Ihling, C.; Soboleva, A.; Lukasheva, E.; Sulima, A.S.; Shtark, O.Y.; Akhtemova, G.A.; Povydysh, M.N.; Sinz, A.; et al. Profiling of Seed Proteome in Pea (Pisum sativum L.) Lines Characterized with High and Low Responsivity to Combined Inoculation with Nodule Bacteria and Arbuscular Mycorrhizal Fungi. Molecules 2019, 24, 1603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Xie, F.; Smith, R.D.; Shen, Y. Advanced proteomic liquid chromatography. J. Chromatogr. A 2012, 1261, 78–90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Van-An, D.; Jeeyun, A.; Na-Young, H.; Jong-Moon, P.; Jeong-Hun, M.; Tae Wan, K.; Hookeun, L. Proteomic Analysis of the Vitreous Body in Proliferative and Non-Proliferative Diabetic Retinopathy. Curr. Proteomics 2020, 17. [Google Scholar] [CrossRef]
  19. Ferrari, E.; Wittig, A.; Basilico, F.; Rossi, R.; De Palma, A.; Di Silvestre, D.; Sauerwein, W.A.G.; Mauri, P.L. Urinary Proteomics Profiles Are Useful for Detection of Cancer Biomarkers and Changes Induced by Therapeutic Procedures. Molecules 2019, 24, 794. [Google Scholar] [CrossRef] [Green Version]
  20. Guan, W.; Liu, Y.; Liu, Y.; Wang, Q.; Ye, H.-L.; Cheng, Y.-G.; Kuang, H.-X.; Jiang, X.-C.; Yang, B.-Y. Proteomics Research on the Protective Effect of Mangiferin on H9C2 Cell Injury Induced by H2O2. Molecules 2019, 24, 1911. [Google Scholar] [CrossRef] [Green Version]
  21. Zimny-Arndt, U.; Schmid, M.; Ackermann, R.; Jungblut, P.R. Classical Proteomics: Two-Dimensional Electrophoresis/MALDI Mass Spectrometry. In Mass Spectrometry of Proteins and Peptides: Methods and Protocols; Lipton, M.S., Paša-Tolic, L., Eds.; Humana Press: Totowa, NJ, USA, 2009; pp. 65–91. [Google Scholar]
  22. O’Farrell, P.H. High resolution two-dimensional electrophoresis of proteins. J. Biol. Chem. 1975, 250, 4007–4021. [Google Scholar]
  23. Link, A.J.; Eng, J.; Schieltz, D.M.; Carmack, E.; Mize, G.J.; Morris, D.R.; Garvik, B.M.; Yates, J.R. Direct analysis of protein complexes using mass spectrometry. Nat. Biotechnol. 1999, 17, 676–682. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, J.; Balgley, B.M.; DeVoe, D.L.; Lee, C.S. Capillary Isoelectric Focusing-Based Multidimensional Concentration/Separation Platform for Proteome Analysis. Anal. Chem. 2003, 75, 3145–3152. [Google Scholar] [CrossRef] [PubMed]
  25. Kašička, V. Capillary electrophoresis of peptides. Electrophoresis 1999, 20, 3084–3105. [Google Scholar] [CrossRef]
  26. Wu, Q.; Yuan, H.; Zhang, L.; Zhang, Y. Recent advances on multidimensional liquid chromatography–mass spectrometry for proteomics: From qualitative to quantitative analysis—A review. Anal. Chim. Acta 2012, 731, 1–10. [Google Scholar] [CrossRef] [PubMed]
  27. Stoll, D.R.; Lhotka, H.R.; Harmes, D.C.; Madigan, B.; Hsiao, J.J.; Staples, G.O. High resolution two-dimensional liquid chromatography coupled with mass spectrometry for robust and sensitive characterization of therapeutic antibodies at the peptide level. J. Chromatogr. B 2019, 1134–1135, 121832. [Google Scholar] [CrossRef]
  28. Capriotti, A.L.; Cavaliere, C.; Cavazzini, A.; Gasparrini, F.; Pierri, G.; Piovesana, S.; Laganà, A. A multidimensional liquid chromatography–tandem mass spectrometry platform to improve protein identification in high-throughput shotgun proteomics. J. Chromatogr. A 2017, 1498, 176–182. [Google Scholar] [CrossRef]
  29. Fournier, M.L.; Gilmore, J.M.; Martin-Brown, S.A.; Washburn, M.P. Multidimensional Separations-Based Shotgun Proteomics. Chem. Rev. 2007, 107, 3654–3686. [Google Scholar] [CrossRef]
  30. Duong, V.-A.; Park, J.-M.; Lee, H. Review of Three-Dimensional Liquid Chromatography Platforms for Bottom-Up Proteomics. Int. J. Mol. Sci. 2020, 21, 1524. [Google Scholar] [CrossRef] [Green Version]
  31. Di Palma, S.; Hennrich, M.L.; Heck, A.J.R.; Mohammed, S. Recent advances in peptide separation by multidimensional liquid chromatography for proteome analysis. J. Proteomics 2012, 75, 3791–3813. [Google Scholar] [CrossRef]
  32. Giddings, J.C. Two-dimensional separations: concept and promise. Anal. Chem. 1984, 56, 1258A–1270A. [Google Scholar] [CrossRef]
  33. Washburn, M.P.; Wolters, D.; Yates, J.R. Large-scale analysis of the yeast proteome by multidimensional protein identification technology. Nat. Biotechnol. 2001, 19, 242–247. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, W.; Wang, S.; Adhikari, S.; Deng, Z.; Wang, L.; Chen, L.; Ke, M.; Yang, P.; Tian, R. Simple and Integrated Spintip-Based Technology Applied for Deep Proteome Profiling. Anal. Chem. 2016, 88, 4864–4871. [Google Scholar] [CrossRef]
  35. Jungblut, P.R.; Holzhütter, H.G.; Apweiler, R.; Schlüter, H. The speciation of the proteome. Chem. Central J. 2008, 2, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Schlüter, H.; Apweiler, R.; Holzhütter, H.-G.; Jungblut, P.R. Finding one’s way in proteomics: a protein species nomenclature. Chem. Central J. 2009, 3, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Thiede, B.; Koehler, C.J.; Strozynski, M.; Treumann, A.; Stein, R.; Zimny-Arndt, U.; Schmid, M.; Jungblut, P.R. High resolution quantitative proteomics of HeLa cells protein species using stable isotope labeling with amino acids in cell culture(SILAC), two-dimensional gel electrophoresis(2DE) and nano-liquid chromatograpohy coupled to an LTQ-OrbitrapMass spectrometer. Mol. Cell. Proteomics 2013, 12, 529–538. [Google Scholar] [PubMed] [Green Version]
  38. Takemori, A.; Butcher, D.S.; Harman, V.M.; Brownridge, P.J.; Shima, K.; Higo, D.; Ishizaki, J.; Hasegawa, H.; Suzuki, J.; Yamashita, M.; et al. PEPPI-MS: polyacrylamide gel-based prefractionation for analysis of intact proteoforms and protein complexes by mass spectrometry. J. Proteome Res. 2020. [Google Scholar] [CrossRef]
  39. Zhao, P.; Schulz, T.C.; Sherrer, E.S.; Weatherly, D.B.; Robins, A.J.; Wells, L. The human embryonic stem cell proteome revealed by multidimensional fractionation followed by tandem mass spectrometry. Proteomics 2015, 15, 554–566. [Google Scholar] [CrossRef] [Green Version]
  40. Motoyama, A.; Xu, T.; Ruse, C.I.; Wohlschlegel, J.A.; Yates, J.R. Anion and Cation Mixed-Bed Ion Exchange for Enhanced Multidimensional Separations of Peptides and Phosphopeptides. Anal. Chem. 2007, 79, 3623–3634. [Google Scholar] [CrossRef] [Green Version]
  41. Wolters, D.A.; Washburn, M.P.; Yates, J.R. An Automated Multidimensional Protein Identification Technology for Shotgun Proteomics. Anal. Chem. 2001, 73, 5683–5690. [Google Scholar] [CrossRef]
  42. Dai, J.; Jin, W.-H.; Sheng, Q.-H.; Shieh, C.-H.; Wu, J.-R.; Zeng, R. Protein Phosphorylation and Expression Profiling by Yin-Yang Multidimensional Liquid Chromatography (Yin-Yang MDLC) Mass Spectrometry. J. Proteome Res. 2007, 6, 250–262. [Google Scholar] [CrossRef] [Green Version]
  43. Dwivedi, R.C.; Spicer, V.; Harder, M.; Antonovici, M.; Ens, W.; Standing, K.G.; Wilkins, J.A.; Krokhin, O.V. Practical Implementation of 2D HPLC Scheme with Accurate Peptide Retention Prediction in Both Dimensions for High-Throughput Bottom-Up Proteomics. Anal. Chem. 2008, 80, 7036–7042. [Google Scholar] [CrossRef]
  44. Gilar, M.; Olivova, P.; Daly, A.E.; Gebler, J.C. Orthogonality of Separation in Two-Dimensional Liquid Chromatography. Anal. Chem. 2005, 77, 6426–6434. [Google Scholar] [CrossRef]
  45. Wang, Y.; Yang, F.; Gritsenko, M.A.; Wang, Y.; Clauss, T.; Liu, T.; Shen, Y.; Monroe, M.E.; Lopez-Ferrer, D.; Reno, T.; et al. Reversed-phase chromatography with multiple fraction concatenation strategy for proteome profiling of human MCF10A cells. Proteomics 2011, 11, 2019–2026. [Google Scholar] [CrossRef] [PubMed]
  46. Yang, F.; Shen, Y.; Camp, D.G.; Smith, R.D. High-pH reversed-phase chromatography with fraction concatenation for 2D proteomic analysis. Expert Rev. Proteomics 2012, 9, 129–134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Wei, J.; Sun, J.; Yu, W.; Jones, A.; Oeller, P.; Keller, M.; Woodnutt, G.; Short, J.M. Global Proteome Discovery Using an Online Three-Dimensional LC−MS/MS. J. Proteome Res. 2005, 4, 801–808. [Google Scholar] [CrossRef] [PubMed]
  48. Law, H.C.H.; Kong, R.P.W.; Szeto, S.S.W.; Zhao, Y.; Zhang, Z.; Wang, Y.; Li, G.; Quan, Q.; Lee, S.M.Y.; Lam, H.C.; et al. A versatile reversed phase-strong cation exchange-reversed phase (RP–SCX–RP) multidimensional liquid chromatography platform for qualitative and quantitative shotgun proteomics. Analyst 2015, 140, 1237–1252. [Google Scholar] [CrossRef] [PubMed]
  49. Xu, B.; Wang, F.; Song, C.; Sun, Z.; Cheng, K.; Tan, Y.; Wang, H.; Zou, H. Large-Scale Proteome Quantification of Hepatocellular Carcinoma Tissues by a Three-Dimensional Liquid Chromatography Strategy Integrated with Sample Preparation. J. Proteome Res. 2014, 13, 3645–3654. [Google Scholar] [CrossRef]
  50. Shen, Z.; Want, E.J.; Chen, W.; Keating, W.; Nussbaumer, W.; Moore, R.; Gentle, T.M.; Siuzdak, G. Sepsis Plasma Protein Profiling with Immunodepletion, Three-Dimensional Liquid Chromatography Tandem Mass Spectrometry, and Spectrum Counting. J. Proteome Res. 2006, 5, 3154–3160. [Google Scholar] [CrossRef]
  51. Betancourt, L.H.; De Bock, P.-J.; Staes, A.; Timmerman, E.; Perez-Riverol, Y.; Sanchez, A.; Besada, V.; Gonzalez, L.J.; Vandekerckhove, J.; Gevaert, K. SCX charge state selective separation of tryptic peptides combined with 2D-RP-HPLC allows for detailed proteome mapping. J. Proteomics 2013, 91, 164–171. [Google Scholar] [CrossRef]
  52. Bandhakavi, S.; Stone, M.D.; Onsongo, G.; Van Riper, S.K.; Griffin, T.J. A Dynamic Range Compression and Three-Dimensional Peptide Fractionation Analysis Platform Expands Proteome Coverage and the Diagnostic Potential of Whole Saliva. J. Proteome Res. 2009, 8, 5590–5600. [Google Scholar] [CrossRef] [Green Version]
  53. Zhou, F.; Lu, Y.; Ficarro, S.B.; Adelmant, G.; Jiang, W.; Luckey, C.J.; Marto, J.A. Genome-scale proteome quantification by DEEP SEQ mass spectrometry. Nature Commun. 2013, 4, 2171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Boichenko, A.P.; Govorukhina, N.; van der Zee, A.G.J.; Bischoff, R. Multidimensional separation of tryptic peptides from human serum proteins using reversed-phase, strong cation exchange, weak anion exchange, and fused-core fluorinated stationary phases. J. Separ. Science 2013, 36, 3463–3470. [Google Scholar] [CrossRef] [PubMed]
  55. Chen, D.; Shen, X.; Sun, L. Strong cation exchange-reversed phase liquid chromatography-capillary zone electrophoresis-tandem mass spectrometry platform with high peak capacity for deep bottom-up proteomics. Anal. Chim. Acta 2018, 1012, 1–9. [Google Scholar] [CrossRef] [PubMed]
  56. Spicer, V.; Ezzati, P.; Neustaeter, H.; Beavis, R.C.; Wilkins, J.A.; Krokhin, O.V. 3D HPLC-MS with Reversed-Phase Separation Functionality in All Three Dimensions for Large-Scale Bottom-Up Proteomics and Peptide Retention Data Collection. Anal. Chem. 2016, 88, 2847–2855. [Google Scholar] [CrossRef]
  57. Nesvizhskii, A.I.; Keller, A.; Kolker, E.; Aebersold, R. A Statistical Model for Identifying Proteins by Tandem Mass Spectrometry. Anal. Chem. 2003, 75, 4646–4658. [Google Scholar] [CrossRef]
  58. Gussakovsky, D.; Neustaeter, H.; Spicer, V.; Krokhin, O.V. Sequence-Specific Model for Peptide Retention Time Prediction in Strong Cation Exchange Chromatography. Anal. Chem. 2017, 89, 11795–11802. [Google Scholar] [CrossRef]
  59. Lorne Burke, T.W.; Mant, C.T.; Black, J.A.; Hodges, R.S. Strong cation-exchange high-performance liquid chromatography of peptides: Effect of non-specific hydrophobic interactions and linearization of peptide retention behaviour. J. Chromatogr. A 1989, 476, 377–389. [Google Scholar] [CrossRef]
  60. Wang, Z.; Ma, H.; Smith, K.; Wu, S. Two-dimensional separation using high-pH and low-pH reversed phase liquid chromatography for top-down proteomics. Int. J. Mass Spectr. 2018, 427, 43–51. [Google Scholar] [CrossRef]
  61. Motoyama, A.; Yates, J.R. Multidimensional LC Separations in Shotgun Proteomics. Anal. Chem. 2008, 80, 7187–7193. [Google Scholar] [CrossRef] [Green Version]
  62. Metwally, H.; McAllister, R.G.; Konermann, L. Exploring the Mechanism of Salt-Induced Signal Suppression in Protein Electrospray Mass Spectrometry Using Experiments and Molecular Dynamics Simulations. Anal. Chem. 2015, 87, 2434–2442. [Google Scholar] [CrossRef] [PubMed]
  63. Oda, Y.; Huang, K.; Cross, F.R.; Cowburn, D.; Chait, B.T. Accurate quantitation of protein expression and site-specific phosphorylation. Proc. Natl. Acad. Sci. USA 1999, 96, 6591–6596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Ong, S.-E.; Blagoev, B.; Kratchmarova, I.; Kristensen, D.B.; Steen, H.; Pandey, A.; Mann, M. Stable Isotope Labeling by Amino Acids in Cell Culture, SILAC, as a Simple and Accurate Approach to Expression Proteomics. Mol. Cell. Proteom. 2002, 1, 376–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Thompson, A.; Schäfer, J.; Kuhn, K.; Kienle, S.; Schwarz, J.; Schmidt, G.; Neumann, T.; Hamon, C. Tandem Mass Tags:  A Novel Quantification Strategy for Comparative Analysis of Complex Protein Mixtures by MS/MS. Anal. Chem. 2003, 75, 1895–1904. [Google Scholar] [CrossRef] [PubMed]
  66. Ross, P.L.; Huang, Y.N.; Marchese, J.N.; Williamson, B.; Parker, K.; Hattan, S.; Khainovski, N.; Pillai, S.; Dey, S.; Daniels, S.; et al. Multiplexed Protein Quantitation in Saccharomyces cerevisiae Using Amine-reactive Isobaric Tagging Reagents. Mol. Cell. Proteom. 2004, 3, 1154–1169. [Google Scholar] [CrossRef] [Green Version]
  67. Bach, L.T.; Mackinder, L.C.M.; Schulz, K.G.; Wheeler, G.; Schroeder, D.C.; Brownlee, C.; Riebesell, U. Dissecting the impact of CO2 and pH on the mechanisms of photosynthesis and calcification in the coccolithophore Emiliania huxleyi. New Phytol. 2013, 199, 121–134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Isenberg, H.D.; Lavine, L.S.; Weissfellner, H. The Suppression of Mineralization in a Coccolithophorid By an Inhibitor of Carbonic Anhydrase. J. Protozool. 1963, 10, 477–479. [Google Scholar] [CrossRef]
  69. Forgac, M. Structure and Properties of the Clathrin-Coated Vesicle and Yeast Vacuolar V-ATPases. J. Bioenerget. Biomembr. 1999, 31, 57–65. [Google Scholar] [CrossRef]
  70. Kirchhausen, T. Clathrin. Ann. Rev. Biochem. 2000, 69, 699–727. [Google Scholar] [CrossRef]
  71. Young, J.R.; Henriksen, K. Biomineralization Within Vesicles: The Calcite of Coccoliths. Rev. Mineral. Geochem. 2003, 54, 189–215. [Google Scholar] [CrossRef]
  72. Durak, G.M.; Taylor, A.R.; Walker, C.E.; Probert, I.; de Vargas, C.; Audic, S.; Schroeder, D.; Brownlee, C.; Wheeler, G.L. A role for diatom-like silicon transporters in calcifying coccolithophores. Nat. Commun. 2016, 7, 10543. [Google Scholar] [CrossRef] [Green Version]
  73. Lee, L.J.Y.; Klute, M.J.; Herman, E.K.; Read, B.; Dacks, J.B. Losses, Expansions, and Novel Subunit Discovery of Adaptor Protein Complexes in Haptophyte Algae. Protist 2015, 166, 585–597. [Google Scholar] [CrossRef] [PubMed]
  74. Paasche, E. A review of the coccolithophorid Emiliania huxleyi (Prymnesiophyceae), with particular reference to growth, coccolith formation, and calcification-photosynthesis interactions. Phycologia 2001, 40, 503–529. [Google Scholar] [CrossRef]
  75. Martiny, A.C.; Coleman, M.L.; Chisholm, S.W. Phosphate acquisition genes in Prochlorococcus ecotypes: Evidence for genome-wide adaptation. Proc. Natl. Acad. Sci. USA 2006, 103, 12552–12557. [Google Scholar] [CrossRef] [Green Version]
  76. Hoch, M.P.; Snyder, R.A.; Jeffrey, W.H.; Dillon, K.S.; Coffin, R.B. Expression of glutamine synthetase and glutamate dehydrogenase by marine bacterioplankton: Assay optimizations and efficacy for assessing nitrogen to carbon metabolic balance in situ. Limnol. Oceanogr. Methods 2006, 4, 308–328. [Google Scholar] [CrossRef] [Green Version]
  77. Sowell, S.M.; Abraham, P.E.; Shah, M.; Verberkmoes, N.C.; Smith, D.P.; Barofsky, D.F.; Giovannoni, S.J. Environmental proteomics of microbial plankton in a highly productive coastal upwelling system. ISME J. 2011, 5, 856–865. [Google Scholar] [CrossRef] [PubMed]
  78. Abida, H.; Ruchaud, S.; Rios, L.; Humeau, A.; Probert, I.; De Vargas, C.; Bach, S.; Bowler, C. Bioprospecting marine plankton. Mar. Drugs 2013, 11, 4594–4611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Tan, L.T.; Phyo, M.Y. Marine Cyanobacteria: A Source of Lead Compounds and their Clinically-Relevant Molecular Targets. Molecules 2020, 25, 2197. [Google Scholar] [CrossRef]
  80. Azam, M.S.; Choi, J.; Lee, M.-S.; Kim, H.-R. Hypopigmenting Effects of Brown Algae-Derived Phytochemicals: A Review on Molecular Mechanisms. Mar. Drugs 2017, 15, 297. [Google Scholar] [CrossRef] [Green Version]
  81. Favre, L.; Ortalo-Magné, A.; Pichereaux, C.; Gargaros, A.; Burlet-Schiltz, O.; Cotelle, V.; Culioli, G. Metabolome and proteome changes between biofilm and planktonic phenotypes of the marine bacterium Pseudoalteromonas lipolytica TC8. Biofouling 2018, 34, 132–148. [Google Scholar] [CrossRef]
  82. Favre, L.; Ortalo-Magné, A.; Kerloch, L.; Pichereaux, C.; Misson, B.; Briand, J.-F.; Garnier, C.; Culioli, G. Metabolomic and proteomic changes induced by growth inhibitory concentrations of copper in the biofilm-forming marine bacterium Pseudoalteromonas lipolytica. Metallomics 2019, 11, 1887–1899. [Google Scholar] [CrossRef] [Green Version]
  83. Maghembe, R.; Damian, D.; Makaranga, A.; Nyandoro, S.S.; Lyantagaye, S.L.; Kusari, S.; Hatti-Kaul, R. Omics for Bioprospecting and Drug Discovery from Bacteria and Microalgae. Antibiotics 2020, 9, 229. [Google Scholar] [CrossRef] [PubMed]
  84. Mishra, A.; Medhi, K.; Malaviya, P.; Thakur, I.S. Omics approaches for microalgal applications: Prospects and challenges. Biores. Technol. 2019, 291, 121890. [Google Scholar] [CrossRef]
  85. Guillard, R.R.L. Culture of Phytoplankton for Feeding Marine Invertebrates. In Culture of Marine Invertebrate Animals: Proceedings—1st Conference on Culture of Marine Invertebrate Animals Greenport; Smith, W.L., Chanley, M.H., Eds.; Springer: Boston, MA, USA, 1975; pp. 29–60. [Google Scholar]
  86. Araie, H.; Sakamoto, K.; Suzuki, I.; Shiraiwa, Y. Characterization of the Selenite Uptake Mechanism in the Coccolithophore Emiliania huxleyi (Haptophyta). Plant Cell Physiol. 2011, 52, 1204–1210. [Google Scholar] [CrossRef] [Green Version]
  87. Joo, M.; Park, J.-M.; Duong, V.-A.; Kwon, D.; Jeon, J.; Han, M.; Cho, B.-K.; Choi, H.-K.; Lee, C.-G.; Kang, H.-G.; et al. An automated high-throughput sample preparation method using double-filtration for serum metabolite LC-MS analysis. Anal. Methods 2019, 11, 4060–4065. [Google Scholar] [CrossRef]
  88. Arul, A.-B.; Byambadorj, M.; Han, N.-Y.; Park, J.M.; Lee, H. Development of an Automated, High-throughput Sample Preparation Protocol for Proteomics Analysis. Bull. Korean Chem. Soc. 2015, 36, 1791–1798. [Google Scholar] [CrossRef]
  89. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2018, 47, D442–D450. [Google Scholar] [CrossRef] [PubMed]
  90. Deutsch, E.W.; Mendoza, L.; Shteynberg, D.; Slagel, J.; Sun, Z.; Moritz, R.L. Trans-Proteomic Pipeline, a standardized data processing pipeline for large-scale reproducible proteomics informatics. Proteomics 2015, 9, 745–754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Keller, A.; Nesvizhskii, A.I.; Kolker, E.; Aebersold, R. Empirical Statistical Model To Estimate the Accuracy of Peptide Identifications Made by MS/MS and Database Search. Anal. Chem. 2002, 74, 5383–5392. [Google Scholar] [CrossRef] [PubMed]
  92. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene Ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [Green Version]
  93. Kanehisa, M.; Goto, S.; Kawashima, S.; Nakaya, A. The KEGG databases at GenomeNet. Nucleic Acids Res. 2002, 30, 42–46. [Google Scholar] [CrossRef]
  94. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.-H.; Pagès, F.; Trajanoski, Z.; Galon, J. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [Green Version]
Sample Availability: Samples of the Emiliania huxleyi cells (strain CCMP371) are available from the authors.
Molecules 25 03028 g001 550
Figure 1. Workflow of proteomic analysis of Emiliania huxleyi (E. huxleyi) using three-dimensional liquid chromatography (3D-LC) system. The proteome digest of E. huxleyi (CCMP371) cell lysates was fractionated into 14 strong cation exchange (SCX) fractions. Each of them was further separated into 16 fractions using high-pH reversed-phase liquid chromatography (HpH RPLC). A set of 16 fractions generated from one SCX fraction was concatenated into five final fractions (1-6-11, 2-7-12, 3-8-13, 4-9-14, and 5-10-15-16). Overall, 70 fractions were generated from the original proteome digest of E. huxleyi cell lysates for liquid chromatography -tandem mass spectrometry (LC-MS/MS) analyses.
Figure 1. Workflow of proteomic analysis of Emiliania huxleyi (E. huxleyi) using three-dimensional liquid chromatography (3D-LC) system. The proteome digest of E. huxleyi (CCMP371) cell lysates was fractionated into 14 strong cation exchange (SCX) fractions. Each of them was further separated into 16 fractions using high-pH reversed-phase liquid chromatography (HpH RPLC). A set of 16 fractions generated from one SCX fraction was concatenated into five final fractions (1-6-11, 2-7-12, 3-8-13, 4-9-14, and 5-10-15-16). Overall, 70 fractions were generated from the original proteome digest of E. huxleyi cell lysates for liquid chromatography -tandem mass spectrometry (LC-MS/MS) analyses.
Molecules 25 03028 g001
Molecules 25 03028 g002a 550Molecules 25 03028 g002b 550
Figure 2. Summary of the results from proteomic profiling of E. huxleyi using 3D-LC system. (a) The accumulated numbers of identified peptides and protein groups (single hit included) versus the number of fractions. (b) Number of unique peptides per protein group. (c) Distribution of peptides by frequencies in 70 fractions. (d) Number of tryptic termini in identified peptides. (e) Number of missed cleavages in identified peptides.
Figure 2. Summary of the results from proteomic profiling of E. huxleyi using 3D-LC system. (a) The accumulated numbers of identified peptides and protein groups (single hit included) versus the number of fractions. (b) Number of unique peptides per protein group. (c) Distribution of peptides by frequencies in 70 fractions. (d) Number of tryptic termini in identified peptides. (e) Number of missed cleavages in identified peptides.
Molecules 25 03028 g002aMolecules 25 03028 g002b
Molecules 25 03028 g003 550
Figure 3. Cumulative distribution of (a) molecular weight, (b) isoelectric point (pI), and (c) the grand average of hydropathy (GRAVY) value of identified peptides. GRAVY values of peptides were calculated using the GRAVY calculator (http://gravy-calculator.de/); negative GRAVY values signify hydrophilicity, and positive values indicate hydrophobicity.
Figure 3. Cumulative distribution of (a) molecular weight, (b) isoelectric point (pI), and (c) the grand average of hydropathy (GRAVY) value of identified peptides. GRAVY values of peptides were calculated using the GRAVY calculator (http://gravy-calculator.de/); negative GRAVY values signify hydrophilicity, and positive values indicate hydrophobicity.
Molecules 25 03028 g003
Molecules 25 03028 g004 550
Figure 4. Distribution of identified proteins groups according to (a) protein length and (b) protein sequence coverage.
Figure 4. Distribution of identified proteins groups according to (a) protein length and (b) protein sequence coverage.
Molecules 25 03028 g004
Molecules 25 03028 g005a 550Molecules 25 03028 g005b 550
Figure 5. Gene ontology (GO) analysis of identified protein groups from E. huxleyi (CCMP371) using ClueGO. Top 10 of (a) biological processes, (b) cellular components, and (c) molecular functions regarding −log10(p-value). (d) Summary of GO term number according to % associated genes.
Figure 5. Gene ontology (GO) analysis of identified protein groups from E. huxleyi (CCMP371) using ClueGO. Top 10 of (a) biological processes, (b) cellular components, and (c) molecular functions regarding −log10(p-value). (d) Summary of GO term number according to % associated genes.
Molecules 25 03028 g005aMolecules 25 03028 g005b
Table
Table 1. Kyoto Encyclopedia of Gene and Genomes (KEGG) pathways involved in identified protein groups from E. huxleyi (CCMP371) using ClueGO.
Table 1. Kyoto Encyclopedia of Gene and Genomes (KEGG) pathways involved in identified protein groups from E. huxleyi (CCMP371) using ClueGO.
IDDescription% Associated GenesNumber of GenesTerm p-Value−log10(p-Value)
KEGG:03040Spliceosome69.41200.00053.3052
KEGG:00531Glycosaminoglycan degradation31.0130.00093.0608
KEGG:04145Phagosome38.0400.00103.0222
KEGG:03022Basal transcription factors83.9260.00172.7735
KEGG:01040Biosynthesis of unsaturated fatty acids21.140.00202.6914
KEGG:04141Protein processing in endoplasmic reticulum67.71050.00462.3355
KEGG:00020Citrate cycle (TCA cycle)78.8260.01231.9104
KEGG:00330Arginine and proline metabolism41.9260.02011.6972
KEGG:03050Proteasome71.9410.02131.6717
KEGG:00430Taurine and hypotaurine metabolism23.130.02171.6635
KEGG:04933AGE-RAGE signaling pathway in diabetic complications33.380.02361.6280
KEGG:04070Phosphatidylinositol signaling system45.3390.03621.4411
KEGG:03410Base excision repair43.5270.03891.4097
KEGG:03015mRNA surveillance pathway67.9550.04091.3885
KEGG:00520Amino sugar and nucleotide sugar metabolism67.0590.04951.3052

Share and Cite

MDPI and ACS Style

Yun, G.; Park, J.-M.; Duong, V.-A.; Mok, J.-H.; Jeon, J.; Nam, O.; Lee, J.; Jin, E.; Lee, H. Proteomic Profiling of Emiliania huxleyi Using a Three-Dimensional Separation Method Combined with Tandem Mass Spectrometry. Molecules 2020, 25, 3028. https://doi.org/10.3390/molecules25133028

AMA Style

Yun G, Park J-M, Duong V-A, Mok J-H, Jeon J, Nam O, Lee J, Jin E, Lee H. Proteomic Profiling of Emiliania huxleyi Using a Three-Dimensional Separation Method Combined with Tandem Mass Spectrometry. Molecules. 2020; 25(13):3028. https://doi.org/10.3390/molecules25133028

Chicago/Turabian Style

Yun, Goyeun, Jong-Moon Park, Van-An Duong, Jeong-Hun Mok, Jongho Jeon, Onyou Nam, Joonwon Lee, EonSeon Jin, and Hookeun Lee. 2020. "Proteomic Profiling of Emiliania huxleyi Using a Three-Dimensional Separation Method Combined with Tandem Mass Spectrometry" Molecules 25, no. 13: 3028. https://doi.org/10.3390/molecules25133028

Article Metrics

Back to TopTop