Next Article in Journal
A New Protocol of Computer-Assisted Image Analysis Highlights the Presence of Hemocytes in the Regenerating Cephalic Tentacles of Adult Pomacea canaliculata
Previous Article in Journal
Broad-Spectrum Bactericidal Activity of a Synthetic Random Copolymer Based on 2-Methoxy-6-(4-Vinylbenzyloxy)-Benzylammonium Hydrochloride
Previous Article in Special Issue
LPIAT, a lyso-Phosphatidylinositol Acyltransferase, Modulates Seed Germination in Arabidopsis thaliana through PIP Signalling Pathways and is Involved in Hyperosmotic Response
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Membrane Profiling by Free Flow Electrophoresis and SWATH-MS to Characterize Subcellular Compartment Proteomes in Mesembryanthemum crystallinum

1
Southern Cross Plant Science, Faculty of Science and Engineering, Southern Cross University, Lismore, NSW 2480, Australia
2
Department of Biochemistry and Molecular Biology, University of Nevada, Reno, NV 89557, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(9), 5020; https://doi.org/10.3390/ijms22095020
Submission received: 22 March 2021 / Revised: 4 May 2021 / Accepted: 6 May 2021 / Published: 9 May 2021
(This article belongs to the Special Issue Plant Lipids: From Physiology to Biotechnological Applications)

Abstract

:
The study of subcellular membrane structure and function facilitates investigations into how biological processes are divided within the cell. However, work in this area has been hampered by the limited techniques available to fractionate the different membranes. Free Flow Electrophoresis (FFE) allows for the fractionation of membranes based on their different surface charges, a property made up primarily of their varied lipid and protein compositions. In this study, high-resolution plant membrane fractionation by FFE, combined with mass spectrometry-based proteomics, allowed the simultaneous profiling of multiple cellular membranes from the leaf tissue of the plant Mesembryanthemum crystallinum. Comparisons of the fractionated membranes’ protein profile to that of known markers for specific cellular compartments sheds light on the functions of proteins, as well as provides new evidence for multiple subcellular localization of several proteins, including those involved in lipid metabolism.

1. Introduction

The plant membrane bilayer provides a barrier between cells and organelles and their surroundings, providing protection from a constantly challenging environment. Cellular membranes also allow compartmentalization of biochemical reactions and pathways into defined units, conferring specific organelles with distinct functions within the cell [1,2]. Despite their importance for this, the contribution and complexity of the broad range of lipids and proteins that make up the different membranes are far from fully understood, and especially, the underlying mechanisms that lead to the diversity of lipid composition remain to be explored [3,4]. A comprehensive analysis of the molecular composition of subcellular membranes would help to provide a more detailed understanding of their organization and function. While mass spectrometry-based proteomics has transformed the ability to perform large-scale protein identification, the ability to profile cellular organelles and endomembranes is still hampered by the limited techniques available to fractionate the membranes. Moreover, when doing this, it is important to be able to associate protein composition with lipid profiles when trying to explain the complexity in membrane processes.
In the past, plant researchers have tried to link the biological activity of membrane fractions with their lipid composition using fractionated subcellular membranes. However, data are generally restricted to studying a single organelle rather than providing a global cellular overview. Approaches such as Localization of Organelle Proteins by Isotope Tagging (LOPIT) have shown promise to decipher protein localizations based on subcellular fractionation by centrifugation and isotope tagging [5,6]. However, this method is unable to associate protein features with lipid profiles of subcellular membranes. Therefore, alternative approaches, which rely on distinct physicochemical properties of membranes while maintaining the structural integrity of lipids and proteins, are required.
Free Flow Electrophoresis (FFE) allows for the continuous electrophoretic separation and fractionation of cells, organelles, membranes, and proteins from a wide variety of organisms and cell types in a thin, non-denaturing liquid film based on net surface charge [7]. The surface charge of subcellular membranes is the sum of mainly the protein and lipids’ combined charges [8] (Figure 1). High resolution and extreme reproducibility are the highlights of this technique [9,10,11,12,13,14,15,16,17,18,19,20,21]. FFE technique involves injecting a complex mixture of membranes that have been initially collected by ultra-centrifugation of cell lysates, into a chamber with a laminar flow buffer stream (Figure 1). These membranes are then deflected by a perpendicular electric field, with the degree of deflection directly related to the charge of the membrane. Fractions (up to 96) are then collected at the end of the separation chamber (Figure 1).
In plants, FFE has been used to increase the purity of specific organelle and endomembrane fractions by first isolating/purifying the organelle via conventional methods and then sub-fractionating those samples. In this way, FFE fractionation has been successfully achieved for thylakoid membranes, ER/Golgi membranes, peroxisomes, mitochondria, vacuole, and plasma membranes [22,23,24,25,26,27], but until now, this method has not been used as a method for global profiling across all cellular membranes in plants.
Combining the high-throughput, high-resolution FFE technology for fractionation of cellular membranes with LC-MS/MS proteomic approaches allows for the comparison of the abundance profiles of proteins throughout all separated fractions. These can then be compared to the profile of known marker proteins of specific compartments to yield classes of co-fractionated proteins from the same organelles, and can allow for the identification and characterization of subcellular compartments. Further, fractionated samples with known subcellular localization can be used for lipidomic studies. Additionally, with the knowledge of the locations of lipid metabolism-related proteins, subcellular metabolic maps of plant lipids may be achieved.
As a model plant for the study of salinity in plants, Mesembryanthemum crystallinum employs both structural and metabolic adaptations including the presence of epidermal bladder cells for salt accumulation and shifting from C3 photosynthesis to CAM to increase water use efficiency [28]. Additionally, at the molecular level, salinity has been show to induce changes to membrane lipids and proteins involved in transport, signaling and compatible solute synthesis [29,30,31]. In this study the ability to fractionate subcellular membranes of M. crystallinum and identify proteins and lipids will facilitate our understanding of the importance of membrane regulation and remodeling in salinity tolerance.

2. Results

2.1. M. crystallinum Membrane Proteome

Free Flow Electrophoresis of microsomal membranes from leaf tissue of M. crystallinum plants resulted in the high-resolution fractionation of membranes into 96 samples. Every second fraction was then combined to give a total of 48 protein samples (i.e., sample 1 = FFE fractions 1 and 2, sample 2 = FFE fraction 3 and 4, etc.). Fractions 15 to 70 (combined samples 8 to 35) with positive protein values at O.D.280 (Figure 2B) were subsequently analyzed by IDA for ion library generation and SWATH-MS for protein quantification. A schematic representation of the approach is illustrated in Figure 2.
Spectra acquired by information-dependent acquisition (IDA) from all samples (28 samples × 3 biological replicates) were submitted to ProteinPilot 5.0.2 to generate the reference spectral library, which included the fragmentation patterns and retention time of each peptide that was required for targeted identification in SWATH-MS [32]. Peptides with confidence scores higher than 95% were selected and summarized. The final ion library contained 1917 distinct proteins at a critical false discovery rate of 1%. The identification and quantification of proteins were performed using a SWATH-MS approach, which resulted in the quantitative export of 1462 unique proteins across all the FFE samples.
In order to describe the composition of the M. crystallinum membrane proteome, transmembrane domains were predicted using the transmembrane topology (HMMTOP v2.0 and TMHMM v2.0) and beta-barrel membrane (MCMBB and TMBETADISC-RBF) prediction programs [33,34,35,36]. The results are shown in Table S1. Proteins predicted to possess a transmembrane domain by any of the four prediction programs (Figure 3) were considered as membrane proteins. Over 62% (908 out of 1462) of the M. crystallinum membrane proteome could be assigned to integral membrane proteins with at least one transmembrane domain, while the others were considered as peripheral or membrane-associated proteins.

2.2. Functional Annotation and Subcellular Localization

Functional annotation of the identified proteins was based on the best BLASTP hits found in the National Center for Biotechnology Information (NCBI), UniprotKB/SwissProt A. thaliana database (May 2020), with 10−5 as the e-value cut-off. The majority (1290, ~88%) of the proteins were matched to homologs from A. thaliana (Table S2). The matched A. thaliana genome accessions (AGIs) were then submitted to the Subcellular Localisation Database For Arabidopsis Proteins 4 (SUBA4) [37] for subcellular localization information. According to SUBA4, approximately 66% of the 1290 proteins were assigned to a high-confidence marker (HCM) for subcellular compartments (Table S2), and the percentage of proteins assigned to each compartment is displayed in the pie graph in Figure 4. Disregarding the unassigned category, proteins classified to a chloroplast HCM comprised the largest proportion of the M. crystallinum membrane proteome (22.3%), followed by cytosol (15.8%), plasma membrane (10.4%), and mitochondria (6.7%) (Table S2). The secretory pathway (from the ER to the TGN) together accounted for 4.3%, while the remaining categories were less than 3%.
Although abundant experimental data from fluorescent protein tagging or mass spectrometry are available for the A. thaliana proteome, only approximately 30% of the proteome is covered by high-confidence location data, whereas the remaining 70% has only been computationally predicted [38]. Using a naive Bayes classifier, the SUBcellular Arabidopsis consensus (SUBAcon) algorithm integrates 22 computational prediction algorithms, for experimental GFP and MS localization, protein–protein interactions, and co-expression data to derive a consensus location of proteins [38]. The SUBAcon predictions of the submitted proteome are listed in Table S2. Based on these predictions, 31.5% of the M. crystallinum leaf proteome were categorized as chloroplast proteins, while 20.7 and 13.1% were from cytosol and plasma membrane, respectively (Figure 4). Proteins assigned to other locations accounted for less than 10%, with 8.7, 8.4, and 5.3% classified to the secretory pathway, mitochondria, and nuclear proteins, respectively. The remaining categories comprised less than 5%, with vacuolar and peroxisome proteins representing 3.3 and 1.5%, respectively. Notably, 5.7% of the identified proteome was predicted to be present in multiple subcellular locations (Figure 4, Table S2).

2.3. FFE Profile of Representative Subcellular Membrane Markers

The averaged abundance of proteins from three biological replicates across all analyzed FFE samples was then carried out to obtain “digital westerns” (Figure 2D), whereby peak area for representative well-characterized marker proteins present in the samples was aligned and visualized in Table 1, demonstrating the positioning of specific subcellular membranes. These markers have been previously verified experimentally with high confidence using cell biology approaches and are accepted to be representative of their particular compartment [37].
Vacuole membrane (tonoplast)-associated proteins appeared in two distinct populations (Table 1 and Table 2). The first population eluting at the farthest left of the FFE chamber showed the highest mobility moving towards the anode electrode from samples 13 to 19; the second population, which showed reduced mobility, was observed as a smaller population with a peak at sample 24. The FFE profiles of the tonoplast H+-ATPase subunits exhibited isoform- and subunit-specific patterns. For example, VHA-B2, VHA-E3, and VHA-F were present evenly in both populations, while the abundance of VHA-a3, VHA-c2, and VHA-d2 were much higher in the anodic population (Table 2). Surprisingly, VHA-A and VHA-D were mainly identified in the second population with a peak at sample 24 and were barely detected in samples 13 to 19. Similar to VHA-a3, VHA-c2, and VHA-d2, tonoplast aquaporins and ion channels, including delta-TIP (TIP2-1) and gamma-TIP (TIP1-1), two-pore calcium channel protein 1 (TPC1), chloride channel protein B (CLC-B), and CLC-C were also present in both populations (Table 2), despite the peak of the second population being much smaller compared to their anodic populations (Table 2).
ER membrane markers, represented by cytochrome b5 isoform B (CYTB5-B) and very-long-chain 3-oxoacyl-CoA reductase 1 (KCR1) (Table 1), were present in two populations from samples 18 to 20, and 23 to 25 (Table 3). Markers of the Golgi apparatus overlapped with ER markers, although the protein abundance was relatively greater in the second population representing samples 23 and 24, as shown by the profiles for polygalacturonate 4-alpha-galacturonosyltransferase (GAUT1) and probable methyltransferase PMT8 (At1g04430) in Table 3. Notably, proteins that are known to frequently exchange between the ER and Golgi were in very low abundance in the anodic population. As an example, Coat Protein Complex II (COPII)-associated Sec23/Sec24, and COPI-associated coatomer subunits, were found almost exclusively in samples 23 and 24 (Table 3). In addition, post-Golgi trafficking markers, shown as AP-1 complex subunit gamma-2 and clathrin heavy chain 1, were mostly concentrated in a single FFE sample, 24, suggesting a slightly different profile of these proteins from those exchanged between ER and Golgi.
Proteins localized to chloroplast envelope membranes, such as inner/outer envelope membrane protein translocon at the inner envelope membrane of chloroplasts 110 (TIC110) and translocon at the outer envelope membrane of chloroplasts 34 (TOC34), were present in a single population fractionating to the middle of the FFE chamber from samples 23 to 25 (Table 1). Nevertheless, proteins localized to thylakoid additionally presented in high abundance in sample 20, as was observed for ATP synthase gamma chain 1 (ATPC1), ATP synthase subunit delta (ATPD), and light reaction-related proteins, with only two exceptions found for photosystem I (PSI) reaction center subunit II-2 (PSAD2) and subunit psaK (PSAK) which presented almost exclusively in FFE samples 21 and 22 (Table 4).
Proteins with known localization to mitochondria were also found in sample 20 (Table 5), although their relative abundance in sample 20 was dramatically lower than that of thylakoid markers (Table 4), as is shown for the mitochondrial outer membrane protein porin 1 (VDAC1) and ATP synthase subunit beta-3 (mitochondrial inner membrane localized) (Table 1 and Table 5). It is also interesting to find that alternative NAD(P)H-ubiquinone oxidoreductase C1 (NDC1), which was previously shown to have dual subcellular locations in both chloroplast and mitochondria [39], has a slightly different profile than either chloroplast or mitochondrial markers in the FFE fractions, with a higher abundance than mitochondrial markers in a population with a peak at sample 20, and also a wider distribution (samples 23–26) than that of the cathodic population of the chloroplast markers which were present mostly in samples 23 and 24 (Table 5).
The PM was the least negatively charged subcellular compartment, fractionating closest to the cathode during FFE. The FFE profile of the PM was determined by the presence of plasma membrane marker proteins in FFE samples 24 to 30 with a peak of abundance at sample 25, for instance, P-type H+-ATPase including AHA2, AHA4, and AHA5, as well as specific PM water channel proteins such as MIP-A, MIP-D, MIP-H, and PIP2;8, with only the profile of PIP2;5 slightly moved towards the anode electrode with a peak of abundance at sample 24 (Table 1 and Table 6). PM-localized transporters, such as sucrose transport protein SUC2 and polyol transporter 5 (PLT5), also showed a typical protein distribution in line with other PM markers. Worth mentioning, as previously reported, is PM nanodomain-localized proteins, fasciclin-like arabinogalactan proteins 2 and 6 (FLA2 and FLA6), and early nodulin-like protein 2 [40], distributed within a relatively narrower population, were only identified in samples 24 to 26 (Table 6).
Table 7 summarizes the distribution and corresponding abundance of several proteins which show unexpected location profiles in the samples, which do not match the published reports. These include sodium/hydrogen exchanger 7 (NHX7/SOS1) [41], sodium/calcium exchanger (NCL) [42], ammonium transporter 1 member 1 (AMT1-1) [43], ABC transporter G family member 40 (ABCG40) [44], ABC transporter C family member 4 (ABCC4) [45], and equilibrative nucleotide transporter 1 (ENT1) [46]. According to previous experimental evidence and SUBAcon, these proteins are suggested to be PM localized. However, in this study, they all fractionated into anodic FFE samples 13-19, showing the greatest overlap with tonoplast markers. Indeed, there are previous reports that have experimentally provided evidence for the dual location of two of these, ABCC4 and ENT1, which were both shown to localize to TP and PM [45,46,47,48].

2.4. FFE Profiles of Lipid Metabolism-Related Proteins

With the help of “digital westerns” constructed using well-characterized membrane marker proteins that revealed the separation and distribution of subcellular compartments across the FFE fractions, a map of the distribution of lipid metabolism-related proteins can be generated. Protein profiles of these proteins are shown in Table 8. The predicted localization of the listed proteins indicates, as expected, that lipid metabolism occurs mainly in Golgi, ER, PM, and chloroplast (Table 8), and these proteins can be assigned to pathways including fatty acid synthesis and elongation, galactolipid and sulfolipid synthesis, phospholipid synthesis and signaling, triacylglycerol synthesis, sterol synthesis, and suberin synthesis.
A good overlap of the protein profiles of lipid metabolism-related proteins with marker proteins from the expected subcellular compartments was observed from the results. For example, lipid-related proteins suggested as chloroplast membrane localized, such as esterase lipase thioesterase family protein, trigalactosyldiacyl-glycerol 2 and 4 (TGD2 and TGD4), fatty acid desaturase 6 and 8 (FAD6 and FAD8), alpha beta-hydrolases superfamily protein (DALL3), single hybrid motif superfamily protein (BCCP2), and long-chain acyl-CoA synthetase 9 (LACS9), were present mainly in the single fraction sample 20 and in the population fractionating to the samples 23 and 24, which match the protein profiles of chloroplast markers present in Table 4. However, the protein abundance of these proteins in the single chloroplast fraction at sample 20 was much lower than was observed for fractions 23 and 24. Moreover, phospholipid glycerol acyltransferase family protein (ATS1/GPAT), AMP-dependent synthetase, and ligase family protein (AAE16) were completely absent from that fraction. The protein distribution for acetyl CoA carboxylase carboxyltransferase alpha subunit (CAC3), however, was different from the other chloroplast-localized proteins, as it was also identified in FFE fraction samples 15 to 19.
Proteins identified as ER proteins, such as beta-ketoacyl reductase 1 (KCR1), cinnamate-4-hydroxylase (CYP73A5), cytochrome b5 (CYTB5-B), phospholipid: diacylglycerol acyltransferase 1 (PDAT1), membrane-associated progesterone-binding protein 3 (MSBP2), 7-dehydrocholesterol reductase (DWF5), and very-long-chain enoyl-CoA reductase (ECR) all presented primarily in two distinct populations (from samples 18 to 20, and 23 to 25), matching what was observed for ER marker proteins (Table 8). Nevertheless, an exception was observed for synaptotagmin-1 (SYT1), which was previously classified as both an ER- and PM-localized protein [49,50]. In this study, the profile of SYT1 matched that of PM markers, being identified in samples 24 to 29 (Table 1 and Table 6), similar to other PM-localized lipid metabolism-related proteins (Table 8), including non-specific phospholipase C4 (NPC4), PLC-like phosphodiesterases superfamily protein, phospholipase D delta (PLDDELTA), non-specific lipid transfer protein GPI-anchored 2 (LTPG2), and phosphatidylinositol-specific phospholipase C4 (PLC4). Additionally, phosphatidylinositol 3- and 4-kinase family protein (PI4KA1), AMP-dependent synthetase, ligase family protein (LACS4), 3-ketoacyl-CoA synthase 6 (CUT1), ABC-2 type transporter ABCG15, and ABC-2 type transporter ABCG22 were also suggested as PM-localized lipid-related proteins by SUBAcon. However, these latter proteins were also identified in samples 13 to 19 overlapping with the tonoplast and the ER profiles, suggesting there may be uncharacterized subcellular locations and biological functions of these proteins.

3. Discussion

The ability to fractionate subcellular membranes with high resolution can provide important information for understanding the subcellular location and biological function of a protein [52]. Combining high-resolution membrane fractionation by Free Flow Electrophoresis (FFE) with the unbiased quantification method of SWATH-MS in this study allowed the characterization of protein profiles for multiple membrane compartments. Subcellular localization of numerous membrane proteins from M. crystallinum leaf tissue was confirmed by their FFE profiles and shown to overlap with known marker proteins. These results also enabled the identification of novel localizations for proteins outside their characterized resident membrane, providing insight into poorly characterized compartments.

3.1. V-ATPase VHA Subunit Localization

The vacuolar proton-pumping V-ATPase consists of 14 subunits divided between two distinct domains; the membrane integral V0 domain and the membrane peripheral V1 domain. In this study, we identified 11 subunits (VHA-A, VHA-B, VHA-C, VHA-D, VHA-E3, VHA-F, VHA-G1, VHA-H, VHA-a3, VHA-c2, and VHA-d2-). The majority of these had similar FFE profiles, showing two distinct populations (Table 2). The more negative anodic population showed close overlap with other tonoplast proteins, including pyrophosphate-energized vacuolar membrane proton pump 1 (AVP1), two tonoplast intrinsic proteins, aquaporins TIP1-1, TIP2-1, and the TPC1 calcium channel. However, there were several VHA subunits that only presented in one of these two populations.
Evidence has shown that V-ATPase subunits are localized not only to the tonoplast but also to the TGN in plants [53]. Previous proteomics studies have identified VHA-A subunit isoforms on Golgi in A. thaliana cell culture proteomic studies [54]. Moreover, pollen grains from an A. thaliana vha-A mutant exhibited severe alterations in the morphology of Golgi stacks and Golgi-derived vesicles, while vacuole morphology remained unaffected [55]. In a vha-a2/vha-a3 double mutant, which showed increased sodium sensitivity, reduction in V-ATPase activity in the trans-Golgi network/early endosome (TGN/EE) was observed, suggesting a role for the V-ATPase in salt tolerance in these membranes [53]. In this study, VHA-A and VHA-D appeared to be absent from the anodic population (samples 13–19), presenting almost exclusively in sample 24 (Table 2). The single population of these two V-ATPase subunit isoforms had almost identical FFE profiles to that of Golgi-specific markers, such as GAUT1, GAUT9, and PMT8 (Table 1, Table 2 and Table 3), supporting their Golgi localization and suggesting that these subunits may have a specific function in this membrane.
In contrast, the VHA-a3 subunit presented mostly in the anodic population from samples 13 to 19, suggesting it was solely tonoplast localized (Table 2), and indicated a different localization to that observed for the isoform in A. thaliana [56]. Additionally, the relatively higher protein abundance (by peak area) of the cathodic population of subunits VHA-B2, VHA-E3, VHA-F, and VHA-H suggested these may also play a role in the Golgi/TGN V-ATPase holoenzyme organization in M. crystallinum [57]. Previous MS/MS-based experiments have also shown their alternative subcellular locations on Golgi-related compartments [54,58].

3.2. Proteins with Unexpected FFE Profiles

The FFE profiles of several PM-classified transporter proteins were unexpected, with proteins exhibiting a fractionation profile overlapping with vacuolar markers, including sodium/hydrogen exchanger 7 (SOS1/NHX7), sodium/calcium exchanger (NCL), ammonium transporter 1-1 (AMT1-1), sulfite exporter TauE/SafE family protein 4 (TauE), ABC transporters (ABCG40 and ABCC4), and equilibrative nucleotide transporter 1 (ENT1) (Table 7). Among these proteins, SOS1, a well-studied member of the Na+/H+ exchanger (NHX7) family, shown to have an important role in plant salt tolerance, showed a profile similar to vesicle trafficking-related proteins, such as PRA1F2, RABB1C, RABG3A, and RABG3F (Table S3). While this protein has been characterized as functioning at the PM in plants [41,59], a study characterizing a sos1 mutant in A. thaliana demonstrated multiple consequences in the cell due to the lack of this transporter, distinct from a role solely in plasma membrane sodium transport; these included the inhibition of vesicle trafficking and endocytosis in the mutant as well as a disintegrated tonoplast and a decrease in vacuolar pH [60]. Further evidence for a non-PM role was shown in another study which demonstrated that the function of the vacuolar-type H+-PPase in A. thaliana under salt stress required the regulation of SOS1, as the pyrophosphate-dependent proton pump activity of the vacuolar-type H+-PPase was significantly lower in the AVP1/sos1 line compared with AVP1 overexpression line, suggesting SOS1 is epistatic to AVP1 [61]. However, evidence for direct interaction between these two proteins was not presented.
In addition, while NCL and AMT1-1 are accepted as PM-localized proteins involved in the transport of sodium/calcium and ammonium, respectively [42,43,62], evidence for non-PM localization has been presented, including vacuolar localization of these proteins as shown by fluorescent-tagging and MS/MS-based approaches [48,63,64], although the biological function for their role on the tonoplast has not been determined. Similar to what was shown for SOS1 and mentioned above, FFE fractionation may have also captured these proteins in an endomembrane trafficking compartment as they share a similar profile to a number of small GTPases (Table S3).
Our study also captured a number of mitochondrial proteins that appeared to be co-localized with plasma membrane fractions, such as cytochrome b-c1 complex subunit Rieske-1 (UCR1-1), mitochondrial-processing peptidase subunit beta (MPPbeta), and mitochondrial carnitine/acylcarnitine carrier-like protein (BOU) (Table 5 and Table 6). Mitochondrial membranes have been shown to be tethered to the PM in yeast [65] and have been identified in plasma membrane proteomics studies in plants, including a recent study using highly purified plasma membrane fractionated in tandem by two-phase partitioning, followed by FFE [26]. This suggests that rather than contaminants, they may represent mitochondrial proteins associated with the plasma membrane through specific contact points. As FFE is a non-denaturing technique, it is highly likely that membrane interactions remain intact during the fractionation.

3.3. A Snapshot of the Components of Electrochemical Reactions

The FFE fractionation of photosynthetic membranes was also observed in this study. Photosynthesis occurs on the thylakoid membrane inside the lumen of the chloroplast, enclosed by a double-membrane envelope composed of the outer envelope membrane (OEM) and the inner envelope membrane (IEM) [66]. The separation of OEM and IEM was not observed in this study, as the OEM marker TOC34 and IEM markers, TIC55 and TIC110, were all mostly present in samples 23 and 24.
However, while the chloroplast inner and outer membranes fractionated together in samples 23 and 24, a particularly large number of thylakoid membrane proteins were identified in sample 20 (Table 4). The thylakoid membrane is a lipoprotein system that has an inner and outer surface with the inner face being more negatively charged [67], and might explain the two distinct populations of thylakoid membrane markers in this study (Table 4). Barber [67] has suggested that the thylakoid membrane is not homogeneous and is likely to be derived mainly from protein rather than lipid components. Interestingly, the protein profile of the primary electron acceptor of photosystem I, ferredoxin-1 [68] was observed to be highly negatively charged in this study (Table 9).
Ferredoxin-1, which based on sequence analysis does not possess any transmembrane domains (tr_41050, Table S1), is likely associated with the membrane via protein/protein interactions with the PS1 trimer. In support of this, a recent study using affinity chromatography and nuclear magnetic resonance demonstrated that the trimeric PSI complex binds three ferredoxins, but in a non-equivalent manner, with one ferredoxin showing stronger binding to PSI than the other two [69]. This might explain our results with the negatively charged population of ferredoxins in the chamber (samples 8 to 14) representing unbound protein that has been released from its protein interaction with PS1 due to the FFE conditions and is more strongly attracted to the anodic electrode during FFE, while the small population found in samples 23 and 24 represents protein still firmly bound to the PSI on the thylakoid membrane, as shown by the small population of this protein which overlapped with thylakoid profiles (Table 9). Interestingly, this is the only protein that shows this profile, and no other proteins were detected in samples 8 to 14 in this study.

3.4. Subcellular Mapping of Lipid Metabolism

Characterization of the membrane distribution of lipid metabolism-related proteins may give a better understanding of lipid signaling, as well as the compositional variance in lipids within subcellular membranes, and how these change under certain developmental and environmental conditions. While previous efforts have demonstrated lipid metabolism pathways based on the types of lipids produced and their subcellular locations [51], the exact localization of some of the enzymes in the biosynthetic pathways remained under question, or in some cases were identified in different subcellular compartments by different studies [70]. The comparative analysis of the protein profiles of lipid metabolism-related proteins with well-known markers in this study would help provide a means for more precise localization of these proteins.
Lipid biosynthesis enzymes are often promiscuous, utilizing a broad range of more or less similar substrates, and the production of these substrates is usually strictly spatially compartmentalized within cells [3], making it easier for us to observe the biological importance of each pathway. To give an example, the plastid is the main site for fatty acid synthesis and glycerolipid production in plants and algae [70,71,72]. Our results agree with this, as lipid biosynthesis proteins suggested as plastid localized by SUBAcon, present mainly in sample 20 and in the population fractionating in samples 23 to 25 (Table 8), matching the profiles of chloroplast markers shown in Table 1 and Table 4. However, unlike what was observed for thylakoid lumenal protein, chlorophyll a-b binding proteins, and photosystem I/II reaction center subunits (Table 4), lipid metabolism-related proteins (chloroplast localized) had much lower abundance in sample 20 (Table 8). This is likely due to the fact that biosynthesis of plastid glycerolipids takes place in the envelope membranes, which is the site of fatty acid assembly [71,73,74], whereas the build-up of thylakoid membranes requires the transport of lipids from envelope membranes during plastid biogenesis [75,76]. Notably, different from the other chloroplast-localized proteins, acetyl CoA carboxylase carboxyltransferase alpha subunit was localized to samples 15-19 (Table 8). Although this protein is a high-confidence marker for chloroplasts, it was previously observed in Golgi [58], vacuole [77], and plasma membrane [78] by MS/MS approaches, indicating this protein may have uncharacterized subcellular locations and biological functions.
In addition to the plastid (prokaryotic pathway), the ER is the other main site for fatty acid elongation, acyl editing, and lipid assembly through the eukaryotic pathway. It was suggested that in A. thaliana, around 2/3 of the fatty acids are thought to be exported to the eukaryotic pathway, although half of which will be returned for plastid lipid assembly [79]. Our data show that proteins suggested as being involved in the fatty acid synthesis and elongation distributed to the same FFE fractions as the ER markers (Table 3 and Table 8), fractionating in samples 18 to 20, and from 23 to 25.
The ER-type fractionation profiles were also seen for proteins that participate in other ER-localized lipid metabolism pathways, such as triacylglycerol biosynthesis (phospholipid: diacylglycerol acyltransferase) and steroid metabolism (membrane-associated progesterone-binding protein 3 and 7-dehydrocholesterol reductase). Notably, the protein abundance in samples 25 to 28 for phospholipid: diacylglycerol acyltransferase, membrane-associated progesterone-binding protein 3,7-dehydrocholesterol reductase, and very-long-chain enoyl-CoA reductase family protein was relatively high compared to the other ER-localized proteins. Interestingly, these proteins were also identified as abundant plasma membrane-localized proteins in previous studies using MS/MS proteomic approaches [26,78,80,81], indicating their subcellular location may need to be redefined and re-validated. By contrast, concrete evidence by fluorescent protein methods has shown the dual localization at both ER and PM of the lipid-binding protein synaptotagmin-1 [49,50], whereas its protein profile is more like that of PM markers in this study (Table 6 and Table 8).
Lipid metabolism-related proteins, including non-specific phospholipase C4, PLC-like phosphodiesterases superfamily protein, phospholipase D delta, non-specific lipid transfer protein GPI-anchored 2, and phosphatidylinositol-specific phospholipase C4, were suggested as PM-localized proteins by SUBAcon and in this study presented in samples 24 to 29 with a peak at 25, which matched the protein profiles of PM markers (Table 6 and Table 8). Nevertheless, several proteins, such as phosphatidylinositol 3- and 4-kinase family protein, AMP-dependent synthetase and ligase family protein, lipid phosphate phosphatase 2, 3-ketoacyl-CoA synthase 6, and ABC-2 type transporter family proteins ABCG15 and ABCG22 were present in samples 13 to 19. Notably, even though 3-ketoacyl-CoA synthase 6 was suggested as PM localized by SUBAcon, it has been hypothesized that this protein is likely involved in the transport of lipids (wax) from the ER to the PM [82], while further verification of its ER localisation has been provided by fluorescent microscopy [83]. In agreement with this, our data also showed the ER marker-like profile of this protein, demonstrating the sensitivity of the method used in this study to characterize protein subcellular location. Furthermore, evidence of non-PM locations was also available for other proteins (suggested as PM localized by SUBAcon). For example, phosphatidylinositol 4-kinase was observed in the cytosol by fluorescent protein assay [84] and AMP-dependent synthetase and ligase family protein was detected in vacuole and trans-Golgi network by MS/MS methods [58], suggesting this approach can also provide valuable information on proteins located in multiple cellular compartments.

4. Materials and Methods

4.1. M. crystallinum Plant Material

Mesembryanthemum crystallinum L. seeds were sown on commercial soil mix supplemented with dolomite, slow-release fertilizer (Scotts Osmocote®, Bella Vista, NSW, Australia), and micronutrients (Micromax®, Unanderra, NSW, Australia) in a seedling propagation tray. Three-week-old seedlings of uniform size were transplanted to pots containing the same soil mixture, with two plants per 15 cm diameter pot. Plants were watered daily and supplemented with one-half strength Hoagland’s solution weekly. Once the plants had reached eight weeks old, the second pair of mature leaves were collected for downstream experimental analysis. Plants were grown in a glasshouse with evaporative cooling to maintain the temperature between 15 and 32 °C and under natural irradiation and photoperiod (photon flux density maximum of 1300 mmol m−2 s−1).

4.2. Microsomal Membrane Extraction

Leaf material (60 g) from M. crystallinum was harvested, cut into small pieces, and immediately placed into 300 mL of ice-cold homogenization medium consisting of 400 mM mannitol, 10% (w/v) glycerol, 5% (w/v) polyvinylpyrrolidine-10, 0.5% (w/v) bovine serum albumin, 1 mM phenylmethylsulphonyl fluoride, 2 mM dithiothreitol, 30 mM trisaminomethane, 5 mM MgSO4, 5 mM EGTA, 0.5 mM butylated hydroxytoluene, 0.25 mM dibucaine, 1 mM benzamidine, and 26 mM K+-metabisulfite, adjusted to pH 8.0 with H2SO4. Once homogenized in a commercial blender, leaf tissue was filtered through cheesecloth and concentrated by centrifugation in a Beckman L8-M ultracentrifuge at 100,000× g for 50 min at 4 °C. The supernatant was removed, and microsomal pellets were resuspended using a Teflon pestle in suspension buffer containing 400 mM mannitol, 10% (w/v) glycerol, 6 mM Tris/Mes (pH 8.0), and 2 mM dithiothreitol, then immediately frozen in liquid nitrogen for storage at −80 °C.

4.3. Free Flow Electrophoresis (FFE)

Microsomal membranes (3 biological replicates) from leaf tissue of M. crystallinum were fractionated using a Free Flow Electrophoresis System (FFE Service GmbH). Samples were diluted 2:1 (v:v) in separation medium containing 250 mM sucrose, 10 mM acetic acid, 10 mM triethanolamine, and 2 mM KCl, and centrifuged at 14,000× g for 20 min at 4 °C prior to injection into the FFE chamber. Diluted samples (3 mg/mL protein) were injected continuously through a peristaltic pump at a rate of 1.2 mL/h via the anodic sample inlet while media inlet composition was as follows: inlets 1 and 9, stabilization medium comprising 180 mM sucrose, 40 mM acetic acid, 40 mM triethylamine, and 8 mM KCl; inlets 2 to 6, separation medium comprising 250 mM sucrose, 10 mM acetic acid, 10 mM triethylamine, and 2 mM KCl: cathodic and anodic circuit electrolyte solutions consisted of 100 mM triethylamine, 100 mM acetic acid, and 10 mM KCl adjusted to pH 7.4 with NaOH. In order to avoid the loss of chloride by anodic oxidation, formaldehyde (0.4%) was added to the anodic solution. The counter flow medium was the same as the separation medium for inlets C1, C2, and C3. FFE was conducted in a horizontal mode at a constant voltage of 750 V (~135 mA) with a media and counterflow rate of 250 mL/h. The temperature of the chamber was maintained at 5 °C throughout the whole FFE process by the continual flow of coolant from a circulating water bath. Each injected sample was fractionated into 96 fractions, which were then continually collected into 15 mL polypropylene tubes placed on ice. Fractionation by FFE was monitored by collecting microtiter plates (200 mL/well) at the beginning and end of each run and protein in each well measured by taking the O.D.280 reading using a microplate scanning spectrophotometer (Victor X4 2030 Multilabel Reader, PerkinElmer, Turku, Finland). Every second fraction was then combined to give a total of 48 protein samples (i.e., fractions 1 and 2, 3 and 4, 5 and 6, etc.). Samples from numbers 8 to 35, which gave positive protein values at O.D.280, were used in subsequent applications.

4.4. TCA Precipitation and Trypsin Digestion of Membrane Proteins

Two hundred microliters of 10X TE, 0.3% (w/v) sodium deoxycholate, and 72% (w/v) trichloroacetic acid were added to 1 mL of each sample in a centrifuge tube in sequence and vortexed between every step. Samples were incubated on ice for 1 h, and the supernatant was discarded by aspiration after centrifugation in a Sigma 4K15 laboratory centrifuge at 14,000× g for 20 min at 4 °C. Pellets were resuspended with 90% methanol then incubated overnight at −20 °C. Following the repeated steps of centrifugation and aspiration, protein pellets were dried in a fume cabinet on ice.
Samples were made up to 1 mg/mL concentrations using 2 M urea in 50 mM ammonium bicarbonate (pH 8.0). One hundred microliters of each sample were digested with 20 µg trypsin and incubated in an ultrasonic bath for 10 min. Tryptic digestion was carried out overnight at 37 °C followed by 4 min in a microwave on the lowest power setting. Samples were dried down using a Heto vacuum centrifuge at 45 °C for 2 h or until dry and dissolved in 100 µL 1% trifluoroacetic acid in milliQ water, then 50 µL samples were aliquoted into autosampler vials for injecting into the nanoHPLC/MS MS/MS system.

4.5. Information-Dependent Acquisition (IDA)

The tryptic peptide samples were analyzed using Ekspert nano LC400 uHPLC (SCIEX, Concord, ON, Canada) coupled with a TripleTOF® 6600 quadrupole time-of-flight (QTOF) mass analyzer (SCIEX, Canada) equipped with a PicoView nanoflow (New Objective, Woburn, MA, USA) ion source. A trap column (5 mm × 300 µm, C18 3 µm, SGE, Ringwood, VIC, Australia) and an analytical column (75 µm × 150 mm ChromXP C18 CL 3 µm, SCIEX, Canada) were used. A linear gradient consisting of 0.1% formic acid in water (Solvent A) and 0.1% formic acid in acetonitrile (Solvent B) was employed with Solvent B from 2 to 40% over 60 min, 40 to 90% for 5 min, 90% for 5 min, and equilibration at 2% between samples. The mobile phase flow rate was 400 nL/min, and the column temperature was set at 45 °C. The peptide samples (5 µL each) were loaded onto the trap column (10 µL/min, 5 min), and eluted onto the analytical column at 400 nL/min. The ionization parameters were as follows: interface heater 150 °C, ionspray voltage 2600 V, declustering potential 80 V, and both nebulizer gas 1 (GS1) and curtain gas flow values set at 30. The protein identification for peptide library production was achieved using the information-dependent acquisition (IDA) mode of the QTOF. During the analysis, a full scan TOF (50 ms, 350–1800 m/z) was first acquired, and then the 30 most intense ions (with minimum 200 counts and a charge state of +2 to +5) of the TOF scan were selected to generate MS/MS spectra (100 ms each, with a rolling collision energy, product ion scan 100–1500 m/z) using the QTOF. Analyst TF 1.7 software (SCIEX, Canada) was used for data acquisition and analysis.

4.6. Sequential Window Acquisition of All Theoretical Mass Spectra (SWATH-MS)

Peptide and protein quantification data were acquired under the same HPLC and mass spectrometry experimental conditions as described for IDA, with the following exceptions. Data were acquired using a SWATH, product ion ms/ms all approach, which continuously acquires all fragment ion spectra (MS2) in an unbiased fashion for every scan cycle for the complete gradient The SWATH experiment was set to acquire 100 product ion spectra from m/z 350 to 1500 per scan cycle with a product ion window set to 6 Da and collision energy from 16 to 60 V with an energy spread of 5 V. The TOF-MS scan acquisition time was set to 50 ms and each product ion scan to 25 ms. The data were acquired and processed using Analyst TF 1.7 software (SCIEX, Canada).

4.7. Ion Library Generation

For ion library generation, all mass spectrometry files were searched in unison in ProteinPilot™ software (Version 5.0.2, SCIEX) [85,86] using the paragon algorithm. Samples were uploaded as unlabeled samples using the following parameters: protein identification, digestion with trypsin, and no special factors (for sample preparation). The search was performed by searching a M. crystallinum L. genomic database of all M. crystallinum ORFs [87]. Peptides with confidence higher than 95% were selected, and the distinct proteins in the final ion library were filtered at a critical false discovery rate of 1%.

4.8. Peptide and Protein Quantification

Spectral alignment and targeted data extraction of SWATH-MS data were conducted in the SWATH Processing Micro App in PeakView (Version 1.2, SCIEX) using the spectral library generated above. All raw files were loaded using an extraction window of 15 min with the following parameters: 6 peptides, 5 transitions and peptide confidence of >95%, excluded shared peptides, and extracted ion chromatogram (XIC) width set at 75 ppm. The peak area of proteins from all samples was normalized to protein concentration measured at O.D.280.

4.9. Trans Membrane Domain, Functional Annotation and Subcellular Localisation

Transmembrane domains of identified proteins were predicted using the transmembrane topology (HMMTOP v2.0 and TMHMM v2.0) and beta-barrel membrane (MCMBB and TMBETADISC-RBF) prediction programs to describe the composition of the M. crystallinum membrane proteome [33,34,35,36]. Functional annotation of the identified M. crystallinum proteome was based on the best BLASTP hits found in the National Center for Biotechnology Information (NCBI), UniprotKB/SwissProt Arabidopsis thaliana database (May 2020), with 10−5 as the e-value cut-off. Arabidopsis thaliana genome accessions (AGIs) of the matched proteins were then submitted to Subcellular Localisation Database For Arabidopsis Proteins 4 (SUBA4) [37] for subcellular location annotation.

5. Conclusions and Perspectives

The future of plant biology studies is the progression from the mere identification of biomolecules in a given organelle to the functional characterization of biomolecular properties in complex cellular environments, that is, associating the results of multiple “omics” studies to help explain the physiological mechanism in plants under different conditions. Relying on the distinct surface charges instead of buoyant densities, this study highlights the advantages of FFE separation technology to provide a means to simultaneously isolate multiple subcellular membranes at high resolution without the need for initial organelle purification, to give insight into protein function and provide valuable information on proteins located in multiple cellular compartments. Due to the non-denaturing feature of the FFE technique, fractionated microsomal membrane samples assigned to known subcellular compartments can also be used for lipidomic studies in the future, and with the knowledge of the abundance and locations of lipid metabolism-related proteins characterized in this study, a sophisticated model for plant membrane metabolism under various developmental and environmental conditions may be achieved.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22095020/s1, Table S1: Transmembrane predictions of identified proteins by four different programs, Table S2: Functional annotation and subcellular localization of identified proteins, Table S3: Digital western FFE profiles of Ras-related proteins.

Author Contributions

Conceptualization, B.J.B.; methodology, B.J.B. and Q.G.; formal analysis, Q.G.; resources, W.C.Y. and J.C.C.; writing—original draft preparation, Q.G.; writing—review and editing, Q.G., B.J.B., L.L., and J.C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Southern Cross Plant Science HDR scholarship for Q.G., and J.C.C acknowledges support from the U.S. Department of Energy, Office of Science, Genomic Science Program (Biosystems Design to Enable Next-generation Biofuels) under award number DE-SC0008834.

Data Availability Statement

Raw spectral data were submitted to the Southern Cross University research portal (https://researchportal.scu.edu.au/, accessed on 4 May 2021) and are available under DOI:10.25918/data.134.

Acknowledgments

The authors thank Rosario Vera-Estrella (Instituto de Biotecnología, Universidad Nacional Autónoma de México) for initial technical work on a pilot project and Alun Jones (Institute for Molecular Bioscience, The University of Queensland, Australia) for mass spectrometry support and technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maron, L. Divide and conquer: Giving organelle context to genomic data. Plant J. 2019, 100, 436–437. [Google Scholar] [CrossRef] [PubMed]
  2. Sweetlove, L.J.; Fernie, A.R. The spatial organization of metabolism within the plant cell. Annu. Rev. Plant Biol. 2013, 64, 723–746. [Google Scholar] [CrossRef] [PubMed]
  3. Harayama, T.; Riezman, H. Understanding the diversity of membrane lipid composition. Nat. Rev. Mol. Cell Biol. 2018, 19, 281–296. [Google Scholar] [CrossRef]
  4. Casares, D.; Escriba, P.V.; Rossello, C.A. Membrane lipid composition: Effect on membrane and organelle structure, function and compartmentalization and therapeutic avenues. Int. J. Mol. Sci. 2019, 20, 2167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Dunkley, T.P.J.; Hester, S.; Shadforth, I.P.; Runions, J.; Weimar, T.; Hanton, S.L.; Griffin, J.L.; Bessant, C.; Brandizzi, F.; Hawes, C.; et al. Mapping the Arabidopsis organelle proteome. Proc. Natl. Acad. Sci. USA 2006, 103, 6518. [Google Scholar] [CrossRef] [Green Version]
  6. Geladaki, A.; Kočevar Britovšek, N.; Breckels, L.M.; Smith, T.S.; Vennard, O.L.; Mulvey, C.M.; Crook, O.M.; Gatto, L.; Lilley, K.S. Combining LOPIT with differential ultracentrifugation for high-resolution spatial proteomics. Nat. Commun. 2019, 10, 331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Roman, M.C.; Brown, P.R. Free Flow Electrophoresis. Anal. Chem. 1994, 66, 86–94. [Google Scholar] [CrossRef]
  8. Burry, R.W.; Wood, J.G. Contributions of lipids and proteins to the surface charge of membranes. An electron microscopy study with cationized and anionized ferritin. J. Cell Biol. 1979, 82, 726–741. [Google Scholar] [CrossRef] [Green Version]
  9. Harms, E.; Kern, H.; Schneider, J.A. Human lysosomes can be purified from diploid skin fibroblasts by free-flow electrophoresis. Proc. Natl. Acad. Sci. USA 1980, 77, 6139. [Google Scholar] [CrossRef] [Green Version]
  10. Barth, F.; Gruetter, M.G.; Kessler, R.; Manz, H.-J. The use of free flow electrophoresis in the purification of recombinant human tissue plasminogen activator expressed in yeast. Electrophoresis 1986, 7, 372–375. [Google Scholar] [CrossRef]
  11. Marsh, M.; Schmid, S.; Kern, H.; Harms, E.; Male, P.; Mellman, I.; Helenius, A. Rapid analytical and preparative isolation of functional endosomes by free flow electrophoresis. J. Cell Biol. 1987, 104, 875–886. [Google Scholar] [CrossRef] [Green Version]
  12. Froshauer, S.; Kartenbeck, J.; Helenius, A. Alphavirus RNA replicase is located on the cytoplasmic surface of endosomes and lysosomes. J. Cell Biol. 1988, 107, 2075–2086. [Google Scholar] [CrossRef] [Green Version]
  13. Schreiber, G.A.; Spiegelberg, A.; Butkowskyj-Walkiw, T.; Bögl, K.W. Assessment of electrophoretic mobility changes of human erythrocytes by free flow electrophoresis after in vitro and in vivo irradiation. Int. J. Radiat. Appl. Instrum. Part B. Nucl. Med. Biol. 1990, 17, 377–380. [Google Scholar] [CrossRef]
  14. Van Der Sluijs, P.; Hull, M.; Zahraoui, A.; Tavitian, A.; Goud, B.; Mellman, I. The small GTP-binding protein rab4 is associated with early endosomes. Proc. Natl. Acad. Sci. USA 1991, 88, 6313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Morré, D.J.; Morré, D.M.; Van Alstine, J.M. Separation of endosomes by aqueous two-phase partition and free-flow electrophoresis. J. Chromatogr. B Biomed. Sci. Appl. 1998, 711, 203–215. [Google Scholar] [CrossRef]
  16. Kronenberger, P.; Schober, D.; Prchla, E.; Ofori-Anyinam, O.; Vrijsen, R.; Rombaut, B.; Blaas, D.; Fuchs, R.; Boeyé, A. Uptake of poliovirus into the endosomal system of HeLa cells. Arch. Virol. 1998, 143, 1417–1424. [Google Scholar] [CrossRef]
  17. Völkl, A.; Mohr, H.; Fahimi, H.D. Peroxisome subpopulations of the rat liver: Isolation by immune Free Flow Electrophoresis. J. Histochem. Cytochem. 1999, 47, 1111–1117. [Google Scholar] [CrossRef] [Green Version]
  18. Xie, H.; Rhodus, N.L.; Griffin, R.J.; Carlis, J.V.; Griffin, T.J. A catalogue of human saliva proteins identified by Free Flow Electrophoresis-based peptide separation and tandem mass spectrometry. Mol. Cell Proteom. 2005, 4, 1826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Wase, N.V.; Ow, S.Y.; Salim, M.; Nissum, M.; Whalley, M. A Free Flow Electrophoresis separation strategy for segregation of high abundant phycobilisomes from cyanobacterium Nostoc punctiforme PCC 73102. J. Proteom. Bioinform. 2012, 5. [Google Scholar] [CrossRef] [Green Version]
  20. Murillo, A.; Vera-Estrella, R.; Barkla, B.J.; Méndez, E.; Arias, C.F. Identification of host cell factors associated with astrovirus replication in Caco-2 Cells. J. Virol. 2015, 89, 10359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Barkla, B.J.; Vera-Estrella, R.; Hernandez-Coronado, M.; Pantoja, O. Quantitative proteomics of the tonoplast reveals a role for glycolytic enzymes in salt tolerance. Plant Cell 2009, 21, 4044–4058. [Google Scholar] [CrossRef] [Green Version]
  22. Eubel, H.; Meyer, E.H.; Taylor, N.L.; Bussell, J.D.; O’Toole, N.; Heazlewood, J.L.; Castleden, I.; Small, I.D.; Smith, S.M.; Millar, A.H. Novel proteins, putative membrane transporters, and an integrated metabolic network are revealed by quantitative proteomic analysis of Arabidopsis cell culture peroxisomes. Plant Physiol. 2008, 148, 1809–1829. [Google Scholar] [CrossRef] [Green Version]
  23. Eubel, H.; Lee, C.P.; Kuo, J.; Meyer, E.H.; Taylor, N.L.; Millar, A.H. TECHNICAL ADVANCE: Free-flow electrophoresis for purification of plant mitochondria by surface charge. Plant J. 2007, 52, 583–594. [Google Scholar] [CrossRef] [PubMed]
  24. Parsons, H.T.; Heazlewood, J.L. Beyond the Western front: Targeted proteomics and organelle abundance profiling. Front. Plant Sci. 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Eichacker, L.A.; Weber, G.; Sukop-Köppel, U.; Wildgruber, R. Free Flow Electrophoresis for separation of native membrane protein complexes. In Proteomic Profiling: Methods and Protocols; Posch, A., Ed.; Springer: New York, NY, USA, 2015; pp. 415–425. [Google Scholar] [CrossRef]
  26. De Michele, R.; McFarlane, H.E.; Parsons, H.T.; Meents, M.J.; Lao, J.; Gonzalez Fernandez-Nino, S.M.; Petzold, C.J.; Frommer, W.B.; Samuels, A.L.; Heazlewood, J.L. Free-Flow Electrophoresis of plasma membrane vesicles enriched by two-phase partitioning enhances the quality of the proteome from Arabidopsis seedlings. J. Proteome Res. 2016, 15, 900–913. [Google Scholar] [CrossRef] [Green Version]
  27. Parsons, H.T.; Stevens, T.J.; McFarlane, H.E.; Vidal-Melgosa, S.; Griss, J.; Lawrence, N.; Butler, R.; Sousa, M.M.L.; Salemi, M.; Willats, W.G.T.; et al. Separating Golgi proteins from cis to trans reveals underlying properties of cisternal localization. Plant Cell 2019, 31, 2010–2034. [Google Scholar] [CrossRef] [Green Version]
  28. Guan, Q.; Kong, W.; Zhu, D.; Zhu, W.; Dufresne, C.; Tian, J.; Chen, S. Comparative proteomics of Mesembryanthemum crystallinum guard cells and mesophyll cells in transition from C3 to CAM. J. Proteom. 2021, 231, 104019. [Google Scholar] [CrossRef]
  29. Barkla, B.J.; Garibay-Hernandez, A.; Melzer, M.; Rupasinghe, T.W.T.; Roessner, U. Single cell-type analysis of cellular lipid remodelling in response to salinity in the epidermal bladder cells of the model halophyte Mesembryanthemum crystallinum. Plant Cell Environ. 2018, 41, 2390–2403. [Google Scholar] [CrossRef] [PubMed]
  30. Klink, R.; Haschke, H.-P.; Kramer, D.; Lüttge, U. Membrane Particles, Proteins and ATPase Activity of Tonoplast Vesicles of Mesembryanthemum crystallinum in the C-3 and CAM State*. Bot. Acta 1990, 103, 24–31. [Google Scholar] [CrossRef]
  31. Barkla, B.J.; Zingarelli, L.; Blumwald, E.; Smith, J.A.C. Tonoplast Na+/H+ antiport activity and tts energization by the bacuolar H+-ATPase in the halophytic plant Mesembryanthemum crystallinum. Plant Physiol. 1995, 109, 549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Schubert, O.T.; Gillet, L.C.; Collins, B.C.; Navarro, P.; Rosenberger, G.; Wolski, W.E.; Lam, H.; Amodei, D.; Mallick, P.; MacLean, B.; et al. Building high-quality assay libraries for targeted analysis of SWATH MS data. Nat. Protoc. 2015, 10, 426–441. [Google Scholar] [CrossRef]
  33. Tusnády, G.E.; Simon, I. The HMMTOP transmembrane topology prediction server. Bioinformatics 2001, 17, 849–850. [Google Scholar] [CrossRef]
  34. Sonnhammer, E.L.L.; Heijne, G.v.; Krogh, A. A hidden markov model for predicting transmembrane helices in protein sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1998, 6, 175–182. [Google Scholar] [PubMed]
  35. Bagos, P.G.; Liakopoulos, T.D.; Spyropoulos, I.C.; Hamodrakas, S.J. A Hidden Markov Model method, capable of predicting and discriminating β-barrel outer membrane proteins. BMC Bioinform. 2004, 5, 29. [Google Scholar] [CrossRef] [Green Version]
  36. Ou, Y.Y.; Gromiha, M.M.; Chen, S.A.; Suwa, M. TMBETADISC-RBF: Discrimination of beta-barrel membrane proteins using RBF networks and PSSM profiles. Comput. Biol. Chem. 2008, 32, 227–231. [Google Scholar] [CrossRef]
  37. Hooper, C.M.; Castleden, I.R.; Tanz, S.K.; Aryamanesh, N.; Millar, A.H. SUBA4: The interactive data analysis centre for Arabidopsis subcellular protein locations. Nucleic Acids Res. 2017, 45, D1064–D1074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Hooper, C.M.; Tanz, S.K.; Castleden, I.R.; Vacher, M.A.; Small, I.D.; Millar, A.H. SUBAcon: A consensus algorithm for unifying the subcellular localization data of the Arabidopsis proteome. Bioinformatics 2014, 30, 3356–3364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Carrie, C.; Murcha, M.W.; Kuehn, K.; Duncan, O.; Barthet, M.; Smith, P.M.; Eubel, H.; Meyer, E.; Day, D.A.; Millar, A.H.; et al. Type II NAD(P)H dehydrogenases are targeted to mitochondria and chloroplasts or peroxisomes in Arabidopsis thaliana. FEBS Lett. 2008, 582, 3073–3079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Tapken, W.; Murphy, A.S. Membrane nanodomains in plants: Capturing form, function, and movement. J. Exp. Bot. 2015, 66, 1573–1586. [Google Scholar] [CrossRef] [Green Version]
  41. Qiu, Q.S.; Barkla, B.J.; Vera-Estrella, R.; Zhu, J.K.; Schumaker, K.S. Na+/H+ exchange activity in the plasma membrane of Arabidopsis. Plant Physiol. 2003, 132, 1041–1052. [Google Scholar] [CrossRef] [Green Version]
  42. Wang, P.; Li, Z.; Wei, J.; Zhao, Z.; Sun, D.; Cui, S. A Na+/Ca2+ exchanger-like protein (AtNCL) involved in salt stress in Arabidopsis. J. Biol. Chem. 2012, 287, 44062–44070. [Google Scholar] [CrossRef] [Green Version]
  43. Loqué, D.; Yuan, L.; Kojima, S.; Gojon, A.; Wirth, J.; Gazzarrini, S.; Ishiyama, K.; Takahashi, H.; Von Wirén, N. Additive contribution of AMT1;1 and AMT1;3 to high-affinity ammonium uptake across the plasma membrane of nitrogen-deficient Arabidopsis roots. Plant J. 2006, 48, 522–534. [Google Scholar] [CrossRef]
  44. Kang, J.; Hwang, J.-U.; Lee, M.; Kim, Y.-Y.; Assmann, S.M.; Martinoia, E.; Lee, Y. PDR-type ABC transporter mediates cellular uptake of the phytohormone abscisic acid. Proc. Natl. Acad. Sci. USA 2010, 107, 2355–2360. [Google Scholar] [CrossRef] [Green Version]
  45. Klein, M.; Geisler, M.; Suh, S.J.; Kolukisaoglu, H.Ü.; Azevedo, L.; Plaza, S.; Curtis, M.D.; Richter, A.; Weder, B.; Schulz, B.; et al. Disruption of AtMRP4, a guard cell plasma membrane ABCC-type ABC transporter, leads to deregulation of stomatal opening and increased drought susceptibility. Plant J. 2004, 39, 219–236. [Google Scholar] [CrossRef]
  46. Li, J.; Wang, D. Cloning and in vitro expression of the cDNA encoding a putative nucleoside transporter from Arabidopsis thaliana. Plant Sci. 2000, 157, 23–32. [Google Scholar] [CrossRef]
  47. Bernard, C.; Traub, M.; Kunz, H.H.; Hach, S.; Trentmann, O.; Mohlmann, T. Equilibrative nucleoside transporter 1 (ENT1) is critical for pollen germination and vegetative growth in Arabidopsis. J. Exp. Bot. 2011, 62, 4627–4637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Li, P.; Zhang, G.; Gonzales, N.; Guo, Y.; Hu, H.; Park, S.; Zhao, J. Ca2+-regulated and diurnal rhythm-regulated Na+/Ca2+ exchanger AtNCL affects flowering time and auxin signalling in Arabidopsis. Plant Cell Environ. 2016, 39, 377–392. [Google Scholar] [CrossRef] [Green Version]
  49. Pérez-Sancho, J.; Vanneste, S.; Lee, E.; McFarlane, H.E.; Esteban Del Valle, A.; Valpuesta, V.; Friml, J.; Botella, M.A.; Rosado, A. The Arabidopsis synaptotagmin1 is enriched in endoplasmic reticulum-plasma membrane contact sites and confers cellular resistance to mechanical stresses. Plant Physiol. 2015, 168, 132–143. [Google Scholar] [CrossRef] [Green Version]
  50. Yamazaki, T.; Takata, N.; Uemura, M.; Kawamura, Y. Arabidopsis synaptotagmin SYT1, a type I signal-anchor protein, requires tandem C2 domains for delivery to the plasma membrane. J. Biol. Chem. 2010, 285, 23165–23176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Li-Beisson, Y.; Shorrosh, B.; Beisson, F.; Andersson, M.X.; Arondel, V.; Bates, P.D.; Baud, S.; Bird, D.; DeBono, A.; Durrett, T.P.; et al. Acyl-Lipid Metabolism. Arab. Book 2013, e0161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Heazlewood, J.L.; Tonti-Filippini, J.; Verboom, R.E.; Millar, A.H. Combining experimental and predicted datasets for determination of the subcellular location of proteins in Arabidopsis. Plant Physiol. 2005, 139, 598–609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Krebs, M.; Beyhl, D.; Gorlich, E.; Al-Rasheid, K.A.; Marten, I.; Stierhof, Y.D.; Hedrich, R.; Schumacher, K. Arabidopsis V-ATPase activity at the tonoplast is required for efficient nutrient storage but not for sodium accumulation. Proc. Natl. Acad. Sci. USA 2010, 107, 3251–3256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Parsons, H.T.; Christiansen, K.; Knierim, B.; Carroll, A.; Ito, J.; Batth, T.S.; Smith-Moritz, A.M.; Morrison, S.; McInerney, P.; Hadi, M.Z.; et al. Isolation and proteomic characterization of the Arabidopsis Golgi defines functional and novel components involved in plant cell wall biosynthesis. Plant Physiol. 2012, 159, 12–26. [Google Scholar] [CrossRef] [Green Version]
  55. Dettmer, J.; Schubert, D.; Calvo-Weimar, O.; Stierhof, Y.D.; Schmidt, R.; Schumacher, K. Essential role of the V-ATPase in male gametophyte development. Plant J. 2005, 41, 117–124. [Google Scholar] [CrossRef]
  56. Dettmer, J.; Hong-Hermesdorf, A.; Stierhof, Y.D.; Schumacher, K. Vacuolar H+-ATPase activity is required for endocytic and secretory trafficking in Arabidopsis. Plant Cell 2006, 18, 715–730. [Google Scholar] [CrossRef] [Green Version]
  57. Strompen, G.; Dettmer, J.; Stierhof, Y.D.; Schumacher, K.; Jurgens, G.; Mayer, U. Arabidopsis vacuolar H+-ATPase subunit E isoform 1 is required for Golgi organization and vacuole function in embryogenesis. Plant J. 2005, 41, 125–132. [Google Scholar] [CrossRef]
  58. Heard, W.; Sklenář, J.; Tomé, D.F.A.; Robatzek, S.; Jones, A.M.E. Identification of regulatory and cargo proteins of endosomal and secretory pathways in Arabidopsis thaliana by proteomic dissection. Mol. Cell Proteom. 2015, 14, 1796–1813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Shi, H.; Quintero, F.J.; Pardo, J.M.; Zhu, J.-K. The putative plasma membrane Na+/H+ antiporter SOS1 controls long-distance Na+ transport in plants. Plant Cell 2002, 14, 465. [Google Scholar] [CrossRef] [Green Version]
  60. Oh, D.-H.; Lee, S.Y.; Bressan, R.A.; Yun, D.-J.; Bohnert, H.J. Intracellular consequences of SOS1 deficiency during salt stress. J. Exp. Bot. 2010, 61, 1205–1213. [Google Scholar] [CrossRef] [Green Version]
  61. Undurraga, S.F.; Santos, M.P.; Paez-Valencia, J.; Yang, H.; Hepler, P.K.; Facanha, A.R.; Hirschi, K.D.; Gaxiola, R.A. Arabidopsis sodium dependent and independent phenotypes triggered by H+-PPase up-regulation are SOS1 dependent. Plant Sci. 2012, 183, 96–105. [Google Scholar] [CrossRef]
  62. Nühse, T.S.; Stensballe, A.; Jensen, O.N.; Peck, S.C. Phosphoproteomics of the Arabidopsis plasma membrane and a new phosphorylation site database. Plant Cell 2004, 16, 2394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Jaquinod, M.; Villiers, F.; Kieffer-Jaquinod, S.; Hugouvieux, V.; Bruley, C.; Garin, J.; Bourguignon, J. A Proteomics dissection of Arabidopsis thaliana vacuoles isolated from cell culture. Mol. Cell Proteom. 2007, 6, 394. [Google Scholar] [CrossRef] [Green Version]
  64. Whiteman, S.A.; Serazetdinova, L.; Jones, A.M.; Sanders, D.; Rathjen, J.; Peck, S.C.; Maathuis, F.J. Identification of novel proteins and phosphorylation sites in a tonoplast enriched membrane fraction of Arabidopsis thaliana. Proteomics 2008, 8, 3536–3547. [Google Scholar] [CrossRef]
  65. Ping, H.A.; Kraft, L.M.; Chen, W.; Nilles, A.E.; Lackner, L.L. Num1 anchors mitochondria to the plasma membrane via two domains with different lipid binding specificities. J. Cell Biol. 2016, 213, 513–524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Staehelin, L.A. Chloroplast structure and supramolecular organization of photosynthetic membranes. In Photosynthesis III: Photosynthetic Membranes and Light Harvesting Systems; Staehelin, L.A., Arntzen, C.J., Eds.; Springer: Berlin/Heidelberg, Germany, 1986. [Google Scholar]
  67. Barber, J. Membrane surface charges and potentials in relation to photosynthesis. Biochim. Biophys. Acta BBA Rev. Bioenergy 1980, 594, 253–308. [Google Scholar] [CrossRef]
  68. Hanke, G.; Mulo, P. Plant type ferredoxins and ferredoxin-dependent metabolism. Plant Cell Environ. 2013, 36, 1071–1084. [Google Scholar] [CrossRef]
  69. Kubota-Kawai, H.; Mutoh, R.; Shinmura, K.; Setif, P.; Nowaczyk, M.M.; Rogner, M.; Ikegami, T.; Tanaka, H.; Kurisu, G. X-ray structure of an asymmetrical trimeric ferredoxin-photosystem I complex. Nat. Plants 2018, 4, 218–224. [Google Scholar] [CrossRef] [PubMed]
  70. Joyard, J.; Ferro, M.; Masselon, C.; Seigneurin-Berny, D.; Salvi, D.; Garin, J.; Rolland, N. Chloroplast proteomics highlights the subcellular compartmentation of lipid metabolism. Prog. Lipid Res. 2010, 49, 128–158. [Google Scholar] [CrossRef] [PubMed]
  71. Kobayashi, K.; Kondo, M.; Fukuda, H.; Nishimura, M.; Ohta, H. Galactolipid synthesis in chloroplast inner envelope is essential for proper thylakoid biogenesis, photosynthesis, and embryogenesis. Proc. Natl. Acad. Sci. USA 2007, 104, 17216–17221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Li-Beisson, Y.; Beisson, F.; Riekhof, W. Metabolism of acyl-lipids in Chlamydomonas reinhardtii. Plant J. 2015, 82, 504–522. [Google Scholar] [CrossRef]
  73. Kobayashi, K. Role of membrane glycerolipids in photosynthesis, thylakoid biogenesis and chloroplast development. J. Plant Res. 2016, 129, 565–580. [Google Scholar] [CrossRef] [Green Version]
  74. Hölzl, G.; Dörmann, P. Chloroplast lipids and their biosynthesis. Annu. Rev. Plant Biol. 2019, 70, 51–81. [Google Scholar] [CrossRef]
  75. Andersson, M.X.; Dörmann, P. Chloroplast Membrane Lipid Biosynthesis and Transport. In The Chloroplast Plant Cell Monographs; Sandelius, A.S., Aronsson, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; Volume 13. [Google Scholar]
  76. Andersson, M.X.; Kjellberg, J.M.; Sandelius, A.S. Chloroplast biogenesis. Regulation of lipid transport to the thylakoid in chloroplasts isolated from expanding and fully expanded leaves of pea. Plant Physiol. 2001, 127, 184–193. [Google Scholar] [CrossRef] [Green Version]
  77. Yoshida, K.; Ohnishi, M.; Fukao, Y.; Okazaki, Y.; Fujiwara, M.; Song, C.; Nakanishi, Y.; Saito, K.; Shimmen, T.; Suzaki, T.; et al. Studies on vacuolar membrane microdomains isolated from Arabidopsis suspension-cultured cells: Local distribution of vacuolar membrane proteins. Plant Cell Physiol. 2013, 54, 1571–1584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Mitra, S.K.; Walters, B.T.; Clouse, S.D.; Goshe, M.B. An efficient organic solvent based extraction method for the proteomic analysis of Arabidopsis plasma membranes. J. Proteome Res. 2009, 8, 2752–2767. [Google Scholar] [CrossRef] [PubMed]
  79. Browse, J.; Warwick, N.; Somerville, C.R.; Slack, C.R. Fluxes through the prokaryotic and eukaryotic pathways of lipid synthesis in the ‘16:3′ plant Arabidopsis thaliana. Biochem. J. 1986, 235, 25–31. [Google Scholar] [CrossRef] [Green Version]
  80. Zargar, S.M.; Kurata, R.; Inaba, S.; Oikawa, A.; Fukui, R.; Ogata, Y.; Agrawal, G.K.; Rakwal, R.; Fukao, Y. Quantitative proteomics of Arabidopsis shoot microsomal proteins reveals a cross-talk between excess zinc and iron deficiency. Proteomics 2015, 15, 1196–1201. [Google Scholar] [CrossRef] [PubMed]
  81. Benschop, J.J.; Mohammed, S.; O’Flaherty, M.; Heck, A.J.; Slijper, M.; Menke, F.L. Quantitative phosphoproteomics of early elicitor signaling in Arabidopsis. Mol. Cell Proteom. 2007, 6, 1198–1214. [Google Scholar] [CrossRef] [Green Version]
  82. Kunst, L.; Samuels, A.L. Biosynthesis and secretion of plant cuticular wax. Prog. Lipid Res. 2003, 42, 51–80. [Google Scholar] [CrossRef]
  83. Joubès, J.; Raffaele, S.; Bourdenx, B.; Garcia, C.; Laroche-Traineau, J.; Moreau, P.; Domergue, F.; Lessire, R. The VLCFA elongase gene family in Arabidopsis thaliana: Phylogenetic analysis, 3D modelling and expression profiling. Plant Mol. Biol. 2008, 67, 547–566. [Google Scholar] [CrossRef]
  84. Okazaki, K.; Miyagishima, S.Y.; Wada, H. Phosphatidylinositol 4-phosphate negatively regulates chloroplast division in Arabidopsis. Plant Cell 2015, 27, 663–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Tang, W.H.; Shilov, I.V.; Seymour, S.L. Nonlinear fitting method for determining local false discovery rates from decoy database searches. J. Proteome Res. 2008, 7, 3661–3667. [Google Scholar] [CrossRef] [PubMed]
  86. Shilov, I.V.; Seymour, S.L.; Patel, A.A.; Loboda, A.; Tang, W.H.; Keating, S.P.; Hunter, C.L.; Nuwaysir, L.M.; Schaeffer, D.A. The Paragon Algorithm, a next generation search engine that uses sequence temperature balues and feature probabilities to identify peptides from tandem mass spectra. Mol. Cell Proteomics 2007, 6, 1638–1655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Yim, W.C.; Cushman, J.C. Caryophyllales genomes provide insights into the evolution of Crassulacean acid metabolism (CAM) and halophytism. In Proceedings of the 26th Plant and Animal Genome Conference, San Diego, CA, USA, 12–17 January 2018. [Google Scholar]
Figure 1. Schematic diagram of the working principle of Free Flow Electrophoresis (FFE) for plant subcellular membrane fractionation. Subcellular membranes with different origins (organelle/subcellular compartment) possess different overall surface charges mainly due to the chemical and compositional diversity of membrane lipids and proteins. In general, membranes exhibit negative surface charges around a neutral pH. Subcellular membranes form vesicles during the microsomal extraction procedure as described in “Materials and Methods”. Microsomal membranes are injected into the chamber and subjected to a high-voltage electric field. Membranes separate according to their distinct surface charge.
Figure 1. Schematic diagram of the working principle of Free Flow Electrophoresis (FFE) for plant subcellular membrane fractionation. Subcellular membranes with different origins (organelle/subcellular compartment) possess different overall surface charges mainly due to the chemical and compositional diversity of membrane lipids and proteins. In general, membranes exhibit negative surface charges around a neutral pH. Subcellular membranes form vesicles during the microsomal extraction procedure as described in “Materials and Methods”. Microsomal membranes are injected into the chamber and subjected to a high-voltage electric field. Membranes separate according to their distinct surface charge.
Ijms 22 05020 g001
Figure 2. Schematic overview of sample fractionation by Free Flow Electrophoresis (FFE) and downstream analysis by mass spectrometry. (A) Microsomal membranes from M. crystallinum leaf tissue were fractionated by FFE into 96 fractions based on their different net surface charge; (B) The absorbance of FFE fractions was determined at 280 nm to identify the range of fractions with positive protein values. Fractions 15 to 70 with positive protein values at O.D.280 were used for subsequent protein identification; (C) Proteins were analyzed by SCIEX TripleTOF 6600 in IDA and SWATH-MS modes for reference library generation, protein identification, and quantification; (D) The average protein abundance of three biological replicates of each sample was used for the “digital western” FFE profile generation. Here we use the term digital western to refer to the detection of proteins by MS/MS in specific FFE fractions, similar to a traditional Western where proteins in the wells of a gel are detected by antibodies. Different colored columns indicate different subcellular membrane origins based on marker protein profiles.
Figure 2. Schematic overview of sample fractionation by Free Flow Electrophoresis (FFE) and downstream analysis by mass spectrometry. (A) Microsomal membranes from M. crystallinum leaf tissue were fractionated by FFE into 96 fractions based on their different net surface charge; (B) The absorbance of FFE fractions was determined at 280 nm to identify the range of fractions with positive protein values. Fractions 15 to 70 with positive protein values at O.D.280 were used for subsequent protein identification; (C) Proteins were analyzed by SCIEX TripleTOF 6600 in IDA and SWATH-MS modes for reference library generation, protein identification, and quantification; (D) The average protein abundance of three biological replicates of each sample was used for the “digital western” FFE profile generation. Here we use the term digital western to refer to the detection of proteins by MS/MS in specific FFE fractions, similar to a traditional Western where proteins in the wells of a gel are detected by antibodies. Different colored columns indicate different subcellular membrane origins based on marker protein profiles.
Ijms 22 05020 g002
Figure 3. Number of transmembrane domains predicted in the M. crystallinum membrane proteome. (A) Transmembrane helix prediction by HMMTOP and TMHMM, and (B) beta-barrel predictions by TMBETADISC-RBF and MSMBB transmembrane prediction data are summarized in Table S1.
Figure 3. Number of transmembrane domains predicted in the M. crystallinum membrane proteome. (A) Transmembrane helix prediction by HMMTOP and TMHMM, and (B) beta-barrel predictions by TMBETADISC-RBF and MSMBB transmembrane prediction data are summarized in Table S1.
Ijms 22 05020 g003
Figure 4. Percentage of the number of distinct proteins for each compartment by (inner pie chart) high-confidence markers (HCM) assignment and (outer pie chart) SUBcellular Arabidopsis consensus (SUBAcon) predictions.
Figure 4. Percentage of the number of distinct proteins for each compartment by (inner pie chart) high-confidence markers (HCM) assignment and (outer pie chart) SUBcellular Arabidopsis consensus (SUBAcon) predictions.
Ijms 22 05020 g004
Table 1. Digital western FFE profiles of selected membrane markers from different subcellular compartments a.
Table 1. Digital western FFE profiles of selected membrane markers from different subcellular compartments a.
Uniprot Recommended Name (Gene) bLocation cDigital Western FFE Profiles d
Ijms 22 05020 i001
V-type proton ATPase subunit a3
(VHA-a3)
TP Ijms 22 05020 i002
aquaporin TIP1-1/gamma-TIP
(TIP1-1)
TP Ijms 22 05020 i003
very-long-chain 3-oxoacyl-CoA reductase 1
(KCR1)
ER Ijms 22 05020 i004
cytochrome b5 isoform B
(CYTB5-B)
ER Ijms 22 05020 i005
polygalacturonate 4-alpha-galacturonosyltransferase
(GAUT1)
GA Ijms 22 05020 i006
galacturonosyltransferase 9
(GAUT9)
GA Ijms 22 05020 i007
translocase of chloroplast 110
(TIC110)
Chl Ijms 22 05020 i008
chlorophyll a-b binding protein 1
(LHCB1.3)
Chl Ijms 22 05020 i009
mitochondrial outer membrane protein porin 2 (VDAC2)MT Ijms 22 05020 i010
prohibitin-3 (PHB3)MT Ijms 22 05020 i011
ATPase 2, plasma membrane-type
(AHA2)
PM Ijms 22 05020 i012
MipA
(mipA)
PM Ijms 22 05020 i013
a Digital western FFE profiles were generated based on the peak area from SWATH-MS quantification; b Sequences of M. crystallinum proteins were BLAST in NCBI against the UniprotKB/Swissprot A. thaliana database; descriptions are based on UniprotKB/Swissprot A. thaliana annotations; c Subcellular location predictions by SUBAcon. TP: tonoplast, ER: endoplasmic reticulum, GA: Golgi apparatus, Chl: chloroplast, MT: mitochondria, PM: plasma membrane; d The x-axis of FFE profiles indicates FFE samples (8–35), the y-axis indicates the average relative abundance (peak area) of three biological replicates for the proteins in the corresponding FFE sample.
Table 2. Digital western FFE profiles of tonoplast marker proteins a.
Table 2. Digital western FFE profiles of tonoplast marker proteins a.
Uniprot Recommended Name (Gene) bDigital Western FFE Profiles c
Ijms 22 05020 i014
V-type proton ATPase subunit A
(VHA-A)
Ijms 22 05020 i015
V-type proton ATPase subunit B2
(VHA-B2)
Ijms 22 05020 i016
V-type proton ATPase subunit C
(VHA-C)
Ijms 22 05020 i017
V-type proton ATPase subunit D
(VHA-D)
Ijms 22 05020 i018
V-type proton ATPase subunit E3
(VHA-E3)
Ijms 22 05020 i019
V-type proton ATPase subunit F
(VHA-F)
Ijms 22 05020 i020
V-type proton ATPase subunit G1
(VHA-G1)
Ijms 22 05020 i021
V-type proton ATPase subunit H
(VHA-H)
Ijms 22 05020 i022
V-type proton ATPase subunit a3
(VHA-a3)
Ijms 22 05020 i023
V-type proton ATPase subunit c2
(VHA-c2)
Ijms 22 05020 i024
V-type proton ATPase subunit d2
(VHA-d2)
Ijms 22 05020 i025
pyrophosphate-energized vacuolar membrane proton pump 1
(AVP1)
Ijms 22 05020 i026
aquaporin TIP2-1/Delta-TIP
(TIP2-1)
Ijms 22 05020 i027
aquaporin TIP1-1/gamma-TIP
(TIP1-1)
Ijms 22 05020 i028
two-pore calcium channel protein 1
(TPC1)
Ijms 22 05020 i029
chloride channel protein CLC-b
(CLC-B)
Ijms 22 05020 i030
chloride channel protein CLC-c
(CLC-C)
Ijms 22 05020 i031
a,b,c Footnotes are the same as for Table 1.
Table 3. Digital western FFE profiles of markers of cargo transport-involved subcellular compartments a.
Table 3. Digital western FFE profiles of markers of cargo transport-involved subcellular compartments a.
Uniprot Recommended Name (Gene) bLocation cDigital Western FFE Profiles d
Ijms 22 05020 i032
calreticulin-1
(CRT1)
ER Ijms 22 05020 i033
SEC12-like protein 2
(STL2P)
ER Ijms 22 05020 i034
protein glycosyltransferase subunit 1B
(OST1B)
ER Ijms 22 05020 i035
translocon-associated protein subunit alpha
(At2g21160)
ER Ijms 22 05020 i036
heat shock 70 kDa protein BIP2
(BIP2)
ER Ijms 22 05020 i037
ADP-ribosylation factor 1
(ARF1)
ER-GA Ijms 22 05020 i038
coatomer subunit alpha-1
(At1g62020)
ER-GA Ijms 22 05020 i039
coatomer subunit beta-1
(At4g31480)
ER-GA Ijms 22 05020 i040
coatomer subunit beta’-2
(At1g52360)
ER-GA Ijms 22 05020 i041
coatomer subunit gamma
(At4g34450)
ER-GA Ijms 22 05020 i042
transport protein Sec23/24-like
(At3g07100)
ER-GA Ijms 22 05020 i043
apyrase 2
(APY2)
GA Ijms 22 05020 i044
polygalacturonate 4-alpha-galacturonosyltransferase
(GAUT1)
GA Ijms 22 05020 i045
galacturonosyltransferase 9
(GAUT9)
GA Ijms 22 05020 i046
probable methyltransferase PMT8
(At1g04430)
GA Ijms 22 05020 i047
AP-1 complex subunit gamma-2
(At1g60070)
TGN Ijms 22 05020 i048
clathrin heavy chain 1
(CHC1)
TGN Ijms 22 05020 i049
a,b,d Footnotes are the same as for Table 1; c Subcellular location abbreviations: ER, endoplasmic reticulum; ER-GA, protein exchange between ER and Golgi; GA, Golgi apparatus; TGN, trans-Golgi network.
Table 4. Digital western FFE profiles of markers localized to chloroplast envelope and thylakoid membranes a.
Table 4. Digital western FFE profiles of markers localized to chloroplast envelope and thylakoid membranes a.
Uniprot Recommended Name (Gene) bLocation cDigital Western FFE Profiles d
Ijms 22 05020 i050
translocase of chloroplast 55
(TIC55)
IME Ijms 22 05020 i051
translocase of chloroplast 110
(TIC110)
IME Ijms 22 05020 i052
translocase of chloroplast 34
(TOC34)
OME Ijms 22 05020 i053
thylakoid lumenal 16.5 kDa protein
(At4g02530)
TLK Ijms 22 05020 i054
ATP synthase gamma chain 1
(ATPC1)
TLK Ijms 22 05020 i055
ATP synthase subunit delta
(ATPD)
TLK Ijms 22 05020 i056
chlorophyll a-b binding protein 1
(LHCB1.3)
TLK Ijms 22 05020 i057
chlorophyll a-b binding protein 6
(LHCA1)
TLK Ijms 22 05020 i058
photosystem I reaction center subunit II-2
(PSAD2)
TLK Ijms 22 05020 i059
photosystem I reaction center subunit III
(PSAF)
TLK Ijms 22 05020 i060
photosystem I reaction center subunit VI-2
(PSAH2)
TLK Ijms 22 05020 i061
photosystem I reaction center subunit XI
(PSAL)
TLK Ijms 22 05020 i062
photosystem I reaction center subunit psaK
(PSAK)
TLK Ijms 22 05020 i063
photosystem II stability/assembly factor HCF136
(HCF136)
TLK Ijms 22 05020 i064
photosystem II 10 kDa polypeptide
(PSBR)
TLK Ijms 22 05020 i065
a,b,d Footnotes are the same as for Table 1; c Subcellular location abbreviations: IME/OME, chloroplast inner/outer envelope; TLK, thylakoid.
Table 5. Digital western FFE profiles of mitochondrial markers a.
Table 5. Digital western FFE profiles of mitochondrial markers a.
Uniprot Recommended Name (Gene) bLocation cDigital Western FFE Profiles d
Ijms 22 05020 i066
mitochondrial outer membrane protein porin 1
(VDAC1)
MOM Ijms 22 05020 i067
mitochondrial outer membrane protein porin 2
(VDAC2)
MOM Ijms 22 05020 i068
mitochondrial outer membrane protein porin 4
(VDAC4)
MOM Ijms 22 05020 i069
mitochondrial import receptor subunit TOM9-2
(TOM9-2)
MOM Ijms 22 05020 i070
prohibitin-2
(PHB2)
MIM Ijms 22 05020 i071
prohibitin-3
(PHB3)
MIM Ijms 22 05020 i072
ATP synthase subunit beta-3
(At5g08680)
MIM Ijms 22 05020 i073
cytochrome b-c1 complex subunit rieske-1
(UCR1-1)
MIM Ijms 22 05020 i074
cytochrome c oxidase subunit 5b-2
(COX5B-2)
MIM Ijms 22 05020 i075
mitochondrial-processing peptidase subunit beta
(MPPbeta)
MIM Ijms 22 05020 i076
mitochondrial carnitine/acylcarnitine carrier-like protein
(BOU)
MIM Ijms 22 05020 i077
mitochondrial phosphate carrier protein
(MPT3)
MIM Ijms 22 05020 i078
mitochondrial dicarboxylate/tricarboxylate transporter DTC
(DTC)
MIM Ijms 22 05020 i079
alternative NAD(P)H-ubiquinone oxidoreductase C1
(NDC1)
DL Ijms 22 05020 i080
a,b,d Footnotes are the same as for Table 1; c Subcellular location abbreviations: MIM/MOM, mitochondrial inner/outer membrane; DL, protein dual-localized to chloroplast and mitochondrial membranes.
Table 6. Digital western FFE profiles of plasma membrane markers a.
Table 6. Digital western FFE profiles of plasma membrane markers a.
Uniprot Recommended Name (Gene) bDigital Western FFE Profiles c
Ijms 22 05020 i081
ATPase 2, plasma membrane-type
(AHA2)
Ijms 22 05020 i082
ATPase 4, plasma membrane-type
(AHA4)
Ijms 22 05020 i083
ATPase 5, plasma membrane-type
(AHA5)
Ijms 22 05020 i084
MipA
(mipA)
Ijms 22 05020 i085
MipD
(mipD)
Ijms 22 05020 i086
MipH
(mipH)
Ijms 22 05020 i087
aquaporin pip2;5
(PIP2-5)
Ijms 22 05020 i088
aquaporin pip2;8
(PIP2-8)
Ijms 22 05020 i089
CSC1-like protein ERD4
(ERD4)
Ijms 22 05020 i090
sucrose transport protein SUC2
(SUC2)
Ijms 22 05020 i091
polyol transporter 5
(PLT5)
Ijms 22 05020 i092
fasciclin-like arabinogalactan protein 2
(FLA2)
Ijms 22 05020 i093
fasciclin-like arabinogalactan protein 6
(FLA6)
Ijms 22 05020 i094
early nodulin-like protein 2
(At4g27520)
Ijms 22 05020 i095
a,b,c Footnotes are the same as for Table 1.
Table 7. Unexpected FFE profiles of proteins known to localize to plasma membrane a.
Table 7. Unexpected FFE profiles of proteins known to localize to plasma membrane a.
Uniprot Recommended Name (Gene) bDigital Western FFE Profiles c
Ijms 22 05020 i096
sodium/hydrogen exchanger 7
(NHX7)
Ijms 22 05020 i097
sodium/calcium exchanger
(NCL)
Ijms 22 05020 i098
ammonium transporter 1 member 1
(AMT1-1)
Ijms 22 05020 i099
sulfite exporter TauE/SafE family protein 4
(TauE)
Ijms 22 05020 i100
ABC transporter G family member 40
(ABCG40)
Ijms 22 05020 i101
ABC transporter C family member 4
(ABCC4)
Ijms 22 05020 i102
equilibrative nucleotide transporter 1
(ENT1)
Ijms 22 05020 i103
a,b,c Footnotes are the same as for Table 1.
Table 8. Digital western FFE profiles of lipid metabolism-related proteins a.
Table 8. Digital western FFE profiles of lipid metabolism-related proteins a.
Uniprot Recommended Name (Gene) bLipid Metabolism Pathway cLocation dDigital Western FFE Profiles e
Ijms 22 05020 i104
sterol methyltransferase 1
(SMT1)
sterol synthesisGA Ijms 22 05020 i105
glycerol-3-phosphate acyltransferase 4
(GPAT4)
suberin synthesis and
transport
GA Ijms 22 05020 i106
3-ketoacyl-CoA synthase 10
(FDH)
fatty acid synthesisGA Ijms 22 05020 i107
very-long-chain 3-oxoacyl-CoA reductase 1
(KCR1)
fatty acid elongation and
wax biosynthesis
ER Ijms 22 05020 i108
cinnamate-4-hydroxylase
(CYP73A5)
fatty acid elongation and
wax biosynthesis
ER Ijms 22 05020 i109
cytochrome b5 isoform B
(CYTB5-B)
fatty acid synthesisER Ijms 22 05020 i110
phospholipid: diacylglycerol acyltransferase
(PDAT1)
triacylglycerol biosynthesisER Ijms 22 05020 i111
membrane-associated progesterone-binding protein 3
(MSBP2)
steroid bindingER Ijms 22 05020 i112
7-dehydrocholesterol reductase
(DWF5)
sterol synthesisER Ijms 22 05020 i113
very-long-chain enoyl-CoA
reductase
(ECR)
fatty acid elongation and
wax biosynthesis
ER Ijms 22 05020 i114
synaptotagmin-1
(SYT1)
lipid bindingER,PM Ijms 22 05020 i115
mitochondrial acyl carrier
protein 2
(ACP)
mitochondrial fatty acid and
lipoic acid synthesis
MT Ijms 22 05020 i116
non-specific phospholipase C4
(NPC4)
eukaryotic galactolipid and
sulfolipid synthesis
PM Ijms 22 05020 i117
PLC-like phosphodiesterases superfamily protein
(At5g67130)
lipid signalingPM Ijms 22 05020 i118
phospholipase D delta
(PLDDELTA)
phospholipid signalingPM Ijms 22 05020 i119
non-specific lipid transfer protein GPI-anchored 2
(LTPG2)
fatty acid elongation and
wax biosynthesis
PM Ijms 22 05020 i120
phosphatidylinositol-specific phospholipase C4
(PLC4)
phospholipid signalingPM Ijms 22 05020 i121
phosphatidylinositol 3- and 4-kinase family protein
(PI4KA1)
phospholipid signalingPM Ijms 22 05020 i122
AMP-dependent synthetase and ligase family protein
(LACS4)
eukaryotic phospholipid synthesis and editingPM Ijms 22 05020 i123
delta(14)-sterol reductase
(FK)
sterol synthesisPM Ijms 22 05020 i124
3-ketoacyl-CoA synthase 6
(CUT1)
fatty acid elongation and
wax biosynthesis
PM Ijms 22 05020 i125
ABC-2 type transporter family protein
(ABCG15)
lipid transportPM Ijms 22 05020 i126
ABC-2 type transporter family protein
(ABCG22)
lipid transportPM Ijms 22 05020 i127
esterase lipase thioesterase family protein
(At1g54570)
fatty acid elongation and
wax biosynthesis
Chl Ijms 22 05020 i128
trigalactosyldiacyl-glycerol 2
(TGD2)
eukaryotic galactolipid and sulfolipid synthesisChl Ijms 22 05020 i129
trigalactosyldiacyl-glycerol 4 (TGD4)eukaryotic galactolipid and sulfolipid synthesisChl Ijms 22 05020 i130
translocon at the outer envelope membrane of chloroplasts 159
(TOC159)
galactolipid, sulfolipid, and phospholipid synthesisChl Ijms 22 05020 i131
alpha beta-hydrolases superfamily protein
(DALL3)
oxylipin metabolismChl Ijms 22 05020 i132
single hybrid motif superfamily protein
(BCCP2)
fatty acid synthesisChl Ijms 22 05020 i133
long-chain acyl-CoA synthetase 9
(LACS9)
fatty acid synthesis
galactolipid, sulfolipid, and phospholipid synthesis
Chl Ijms 22 05020 i134
fatty acid desaturase 6
(FAD6)
fatty acid desaturation
galactolipid, sulfolipid, and phospholipid synthesis
Chl Ijms 22 05020 i135
fatty acid desaturase 8
(FAD8)
fatty acid desaturation
galactolipid, sulfolipid, and phospholipid synthesis
Chl Ijms 22 05020 i136
phospholipid glycerol acyltransferase family protein
(ATS1/GPAT)
galactolipid, sulfolipid, and phospholipid synthesisChl Ijms 22 05020 i137
AMP-dependent synthetase and ligase family protein
(AAE16)
fatty acid metabolismChl Ijms 22 05020 i138
acetyl CoA carboxylase carboxyltransferase alpha subunit
(CAC3)
fatty acid synthesisChl Ijms 22 05020 i139
a,b,d,e Footnotes are the same as for Table 1; c Lipid metabolism pathways were assigned based on the A. thaliana Acyl-Lipid Metabolism Database [51] and Uniprot annotations.
Table 9. Protein profile of ferredoxin compared with three other thylakoid markers a.
Table 9. Protein profile of ferredoxin compared with three other thylakoid markers a.
Uniprot Recommended Name bDigital Western FFE Profiles c
Ijms 22 05020 i140
ferredoxin-1
(FD1)
Ijms 22 05020 i141
ATP synthase gamma chain 1
(ATPC1)
Ijms 22 05020 i142
ATP synthase subunit delta
(ATPD)
Ijms 22 05020 i143
photosystem I reaction center subunit VI-2
(PSAH2)
Ijms 22 05020 i144
a,b,c Footnotes are the same as for Table 4.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, Q.; Liu, L.; Yim, W.C.; Cushman, J.C.; Barkla, B.J. Membrane Profiling by Free Flow Electrophoresis and SWATH-MS to Characterize Subcellular Compartment Proteomes in Mesembryanthemum crystallinum. Int. J. Mol. Sci. 2021, 22, 5020. https://doi.org/10.3390/ijms22095020

AMA Style

Guo Q, Liu L, Yim WC, Cushman JC, Barkla BJ. Membrane Profiling by Free Flow Electrophoresis and SWATH-MS to Characterize Subcellular Compartment Proteomes in Mesembryanthemum crystallinum. International Journal of Molecular Sciences. 2021; 22(9):5020. https://doi.org/10.3390/ijms22095020

Chicago/Turabian Style

Guo, Qi, Lei Liu, Won C. Yim, John C. Cushman, and Bronwyn J. Barkla. 2021. "Membrane Profiling by Free Flow Electrophoresis and SWATH-MS to Characterize Subcellular Compartment Proteomes in Mesembryanthemum crystallinum" International Journal of Molecular Sciences 22, no. 9: 5020. https://doi.org/10.3390/ijms22095020

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop