- Research article
- Open Access
Analysis of the lipid body proteome of the oleaginous alga Lobosphaera incisa
BMC Plant Biology volume 17, Article number: 98 (2017)
Lobosphaera incisa (L. incisa) is an oleaginous microalga that stores triacylglycerol (TAG) rich in arachidonic acid in lipid bodies (LBs). This organelle is gaining attention in algal research, since evidence is accumulating that proteins attached to its surface fulfill important functions in TAG storage and metabolism.
Here, the composition of the LB proteome in L incisa was investigated by comparing different cell fractions in a semiquantitative proteomics approach. After applying stringent filters to the proteomics data in order to remove contaminating proteins from the list of possible LB proteins (LBPs), heterologous expression of candidate proteins in tobacco pollen tubes, allowed us to confirm 3 true LBPs: A member of the algal Major Lipid Droplet Protein family, a small protein of unknown function and a putative lipase. In addition, a TAG lipase that belongs to the SUGAR DEPENDENT 1 family of TAG lipases known from oilseed plants was identified. Its activity was verified by functional complementation of an Arabidopsis thaliana mutant lacking the major seed TAG lipases.
Here we describe 3 LBPs as well as a TAG lipase from the oleaginous microalga L. incisa and discuss their possible involvement in LB metabolism. This study highlights the importance of filtering LB proteome datasets and verifying the subcellular localization one by one, so that contaminating proteins can be recognized as such. Our dataset can serve as a valuable resource in the identification of additional LBPs, shedding more light on the intriguing roles of LBs in microalgae.
Microalgae are a highly diverse group of organisms with regards to both their evolutionary and ecological background. Some species of microalgae accumulate the neutral lipid triacylglycerol (TAG) within the cells as a form of carbon and energy storage, a process that can be manipulated by exposure to different forms of abiotic stress. Among these, nitrogen starvation has been shown to result in increased TAG synthesis , a phenomenon that is accompanied by growth arrest [2,3,4,5], alterations in the proteome reducing the nitrogen content [6, 7], reduced photosynthesis on a transcription level [6, 8] and chloroplast degradation [3, 4, 6, 9, 10]. The fact that these changes are reversible makes oleaginous algae particularly interesting for studies of lipid metabolism.
TAG needs to be stored within the cell in a way that permits mobilization when needed. This function is fulfilled by cytosolic organelles the so-called lipid bodies (LBs) . Apart from a core of neutral lipids, they consist of a monolayer of polar lipids and proteins that can be directly or indirectly attached to its surface . A small number of structural elements is known that makes it possible for proteins to bind to LBs directly: A hydrophobic domain that protrudes into the LB , a “proline knot”  or “proline knob”  within this domain, a β-barrel  or amphipathic helices [17, 18]. However, a conserved amino acid sequence that targets proteins to the LB has so far not been identified.
In many organisms, LB proteins (LBPs) have been found to be involved in lipid metabolism or in maintaining LB structure [19, 20]. The former includes TAG lipases, which are required to remove the fatty acids from the backbone of the storage lipid. Mono- and Diacylglycerol (MAG and DAG) lipases may then cleave off the remaining two acyl moieties, which can subsequently be broken down by β-oxidation to access the stored energy. LB-associated TAG lipases have been described in yeast [21, 22], castor bean , rapeseed and mustard seed  as well as Arabidopsis thaliana (A. thaliana) seeds [25, 26]. Until recently, no TAG lipase activity has been established in any alga, with the exception of Phaeodactylum tricornutum TAG LIPASE1. It has been demonstrated that this enzyme is capable of hydrolyzing the TAG analog para-nitrophenyl butyrate in vitro and that suppressed expression of the gene leads to TAG accumulation in vivo , but a direct connection to TAG (rather than DAG or MAG) degradation has not been established to date.
The second category of LBPs comprises proteins that maintain the structural integrity of LBs. The most intensely studied family of structural LBPs was first described in maize seeds , has since been characterized in many species of higher plants and was named oleosins . They have been accredited with functions in preventing LB coalescence, thus maintaining a high surface/volume ratio that is beneficial for rapid degradation when needed  and for freezing tolerance . It has also been shown that oleosins need to be degraded in order for LB breakdown to take place . Similar functions are served by perilipins in mammals [33, 34] and Drosophila  as well as LIPID DROPLET ASSOCIATED PROTEINs (LDAPs) in the vegetative tissue of olive, avocado and oil palm [36, 37]. Microalgae harbor yet another family of structural LBPs, the members of which are mostly called MAJOR LIPID DROPLET PROTEINs (MLDPs) and do not contain any of the known structural elements. The C. reinhardtii MLDP was the first one to be characterized  and it has since become evident that it recruits other proteins, particularly tubulins, to the LB surface . Moreover, C. reinhardtii MLDP transcript abundance has been used as a marker for TAG accumulation . Members of the MLDP family have been characterized in the microalgae Haematococcus pluvialis  and Dunaliella salina , while homologous genes can be found in the genomes of further members of the Volvocales and Chlorellales order . The diatom Phaeodactylum tricornutum  and the heterokont microalga Nannochloropsis  each contain unique structural LB proteins, which, in contrast to MLDPs, contain prominent hydrophobic domains similar to oleosins.
Lobosphaera incisa is an oleaginous alga belonging to the class of Trebouxiophyceae. Here we investigated strain SAG 2468, which was originally isolated on a Japanese glacier . It is unusual in accumulating large amounts of TAG that is rich in arachidonic acid (ARA, 20:4 (n-6)) . Such partitioning of a PUFA into TAG is not common among microalgae  and nitrogen starvation can be used to further push the TAG content from 43%  to 87% of total fatty acids (TFAs), increasing the proportion of ARA at the same time .
The mitochondrial and plastidial genomes of this strain have been published [46, 47] and chemical mutagenesis has been successfully employed, resulting in a clear phenotype . In addition, stable transformation has been achieved , albeit with low efficiency. A much larger array of molecular biology tools and resources is currently available for the model green alga C. reinhardtii, however this organism requires nutrient starvation to exhibit substantial LB formation. In contrast to this, the same process can be readily observed without manipulation in L. incisa.
So far, several algal LB proteomes have been analyzed, but only for very few postulated LBPs the subcellular localization has been verified by fluorescence microscopy [39, 41, 50] or immunolabelling . In order to better understand lipid storage in the oleaginous alga L. incisa, the objective of this study was to identify proteins that associate with LBs and to verify this localization in order to reveal true LBPs. We furthermore used the recently sequenced nuclear genome as well as expression studies to reveal genes of interest for TAG degradation.
As LB-associated proteins have been found in many organisms, the L. incisa genome was initially searched for homologs of these. Neither oleosins, caleosins or steroleosins, which are known from oil seed plants, are encoded in the algal genome, nor could a homolog of plant LDAPs or mammalian perilipins be found. LBs were therefore isolated from L. incisa in order to identify novel LBPs.
Putative LBPs were identified by a proteomics approach
The steps required for LB isolation are outlined in Fig. 1. An L. incisa culture was starved of nitrogen for 3 days prior to cell fractionation, so as to promote LB formation. A large array of mechanical and enzymatic methods of cell disruption (as well as combinations thereof) had proven to be ineffective in breaking the adamant cell wall while leaving LBs intact, whereas grinding in liquid nitrogen yielded sufficient amounts of intact LBs as verified by Nile Red staining (Fig. 1). Repeated cycles of resuspension and centrifugation allowed partial removal of adhering membranes. Samples taken from the total extract, the soluble fraction and total membranes were used as controls in the subsequent identification of true LBPs.
The process of proteome data filtering leading up to this goal is summarized in Table 1. Proteins that had been detected in all three technical replicates of the LB sample were considered for further analysis, if they fulfilled at least one of four criteria: (i) Strong enrichment in the LB fraction compared to one or more control samples, (ii) high abundance in the LB fraction, (iii) a clear induction of gene expression in nitrogen-limiting conditions or (iv) homology with proteins that point to a function in lipid metabolism or LB homeostasis. Those candidates for which the coding sequences could be successfully amplified from cDNA were analyzed for subcellular localization. The genes of interest were fused to the reporter gene mVenus and transiently expressed in Nicotiana tabacum (N. tabacum) pollen tubes. This readily transformable tissue contains a large number of LBs and is therefore an ideal system to confirm the localization of putative LBPs of plant or algal origin .
The 11 proteins that were studied further are presented in more detail in Table 2. For most of the candidates investigated, at least ¼ of the amino acid sequence was covered by peptides detected (see Additional file 1). Enrichment in the LB fraction compared to the three control samples varies greatly, however out of the 5 most abundant proteins, three could be verified as LBPs and an additional protein was found to be partially associated with LBs.
LiMLDP is a small and highly abundant LBP
The protein encoded by gene g555 was identified as an LBP candidate based on its high abundance in the LB protein extract (Table 2) and its similarity to a known algal LBP. A protein BLAST search revealed Haematococcus pluvialis oil globule protein (HpOGP) as the closest characterized homolog in any organism with 30% sequence identity, moreover the two proteins share a striking similarity in the distribution of hydrophobic residues along the amino acid sequence (Fig. 2b). These properties include the protein in the family of algal MLDPs (Fig. 2e) and the resemblance led us to name the L. incisa protein LiMLDP. Heterologous expression in N. tabacum pollen tubes confirmed its localization (Fig. 2a). Analyses of gene expression under conditions of nitrogen starvation and resupply indicate an upregulation during the first 3 days of growth limiting conditions (Fig. 2c and d, Additional file 2). The total fatty acid (TFA) content is shown for comparison (Fig. 2d, Additional file 3).
In contrast to oleosins, MLDPs including LiMLDP do not possess a prominent hydrophobic domain (Fig. 2b) that could directly anchor them at the LB surface. Nevertheless, LiMLDP appears to attach to LBs in some way and a possible function in maintaining LB structural integrity was tested. The gene was placed under the control of the seed specific Brassica napus napin A (napA) promoter and stably expressed in the A. thaliana oleo1 mutant. LBs are dramatically enlarged in imbibed mutant seeds compared to the wildtype, however a functional complementation could not be observed in T2 seeds of 15 independent lines expressing LiMLDP (Fig. 3). See Additional file 4 for confirmation of gene expression.
LiLBP62 localizes to LBs
LiLBP62 is strongly enriched in the LB fraction of L. incisa compared to the control samples (Table 2) and also localizes to the same organelle in tobacco pollen tubes (Fig. 4a). Similarly to LiMLDP, it lacks a distinct hydrophobic domain (Fig. 4b). It contains a domain of unknown function (DUF 4057) [52, 53] that is absolutely specific of plants and green algae. LiLBP62 undergoes no major expression changes over the 7-day course of nitrogen starvation, but the transcript drops down to below detection immediately following nitrogen resupply (Fig. 4c and d, Additional file 2).
Due to its high abundance at the LBs and its small size of 287 amino acids, a potential oleosin-like role for LiLBP62 in the structural integrity of LBs was investigated. Stable expression of the gene in the A. thaliana oleo1 mutant as described above again did not reveal a complementation of LB size in imbibed T2 seeds of 13 independent lines.
LiLBP36 is a putative LB lipase
LiLBP36 is not only highly abundant at L. incisa LBs, but also at least 5× enriched compared to all three control samples (Table 2) and the subcellular localization could be clearly confirmed in tobacco pollen tubes (Fig. 5a). The corresponding transcript levels respond to nitrogen starvation with a sharp increase within 3 days and drop to their initial value upon nitrogen resupply (Fig. 5b and c, Additional file 2).
LiLBP36 was annotated by Interproscan as a class 3 lipase, and the predicted secondary structures of a section of 224 amino acids (see Additional file 5) match very closely the experimentally determined structures of a group of fungal and yeast secreted lipases that act on MAG, DAG and TAG. These enzymes contain a catalytic triad of serine, aspartate and histidine that are found in many lipases , the serine residue being part of a conserved GXSXG motif. Except for the histidine residue, all of these elements are present in the structurally homologous section of LiLBP36 (Fig. 5d).
Based on these observations along with the strong association with LBs, we suspected a TAG lipase activity for this protein. We therefore investigated the ability of LiLBP to complement TAG degradation in an A. thaliana mutant lacking the major TAG lipases that drive postgerminative growth: SUGAR-DEPENDENT1 (SDP1) and SDP1-LIKE . Etiolated sdp1/sdp1-L seedlings require sucrose in order to form equally long hypocotyls as the wildtype, a property that permits an indirect observation of TAG breakdown [25, 26]. We introduced the LiLBP36 gene into the mutant under the control of the napA promoter and measured the hypocotyl length of T2 seedlings after 5 days of germination in the dark in the presence or absence of sucrose in the media (Fig. 6a). For each independent line, the average hypocotyl length with sucrose was then set to 100% the relative length that was reached without sucrose was calculated and then compared to mutant and wildtype seedlings carrying only an empty vector (Fig. 6b). In 4 out of 5 lines, a modest improvement of postgerminative growth could be observed, however a strong functional complementation could not be ascertained (Fig. 6b and c, Additional file 6). See Additional file 4 for confirmation of gene expression.
Several LBP candidates do not localize to LBs when tested in N. tabacum pollen tubes
The search for LBPs yielded 8 more candidates (see Table 2) that were further analyzed. g9582.p1 was selected based on its strong enrichment in the LB fraction and its Interproscan annotation as a class 3 lipase. When transiently expressed in tobacco pollen tubes, the protein localized to some LBs but was also abundant in other parts of the cell. g9864 expression is drastically upregulated in response to nitrogen starvation and encodes a protein of unknown function that accumulates to some extent at the LBs. g13714.p1, a protein with high similarity to A. thaliana fatty acid amide hydrolase, is clearly enriched at LBs. At this point, a possible function for g4703.p1 cannot be predicted, but Interproscan reveals a sterile alpha motif possibly involved in protein binding and the protein is much more abundant in the LB fraction than in total membranes. g13209 encodes a predicted short chain dehydrogenase/reductase and is highly expressed during nitrogen starvation. The putative FAD-binding monooxygenase g12144.p1 is similarly induced, while g13747.p1 and the patatin-related g14373.p1 are characterized by a striking accumulation at the LBs compared to the control samples. These 7 candidates did not localize to the LBs when expressed in tobacco pollen tubes and were thus not investigated further in this study.
LiSDP1 is a TAG lipase
We were surprised not to find a TAG lipase among the proteins that passed our criteria for filtering LB proteomics data and therefore turned to the L. incisa genome in search of promising candidates. This identified a protein that shares 44% identity with A. thaliana SDP1 and that was therefore named LiSDP1. Like its plant homolog, the protein features a patatin domain that is characteristic of a diverse group of lipases . The domain harbors a catalytic dyad consisting of a serine residue within a conserved GXSXG motif as well as an aspartate residue (Fig. 7d). That LiSDP1 was not a LB protein was confirmed by transient expression in N. tabacum pollen tubes, which showed no overlap of the LiSDP1-mVenus and Nile Red signals (Fig. 7a). Gene expression only underwent minor changes during nitrogen starvation (Fig. 7b and c, Additional file 2), whereas it markedly increased within 6 h of nitrogen resupply, a growth phase characterized by substantial fatty acid mobilization (Fig. 7c).
TAG lipase activity was tested by complementing the A. thaliana sdp1/sdp1-L mutant as described above. In 7 out of 8 independent lines carrying the LiSDP1 gene controlled by the napA promoter, a significantly improved postgerminative growth could be observed compared to the mutant (Fig. 8a) and three lines expressing the gene constitutively showed an even stronger effect (Fig. 8b, see Additional file 4 for confirmation of gene expression). The assay was extended for selected lines by quantifying fatty acids in seedlings and relating them to the average seed fatty acid content of the same line (Fig. 8c, Additional file 7). In wildtype seedlings, the majority of fatty acids is metabolized during postgerminative growth, whereas the mutant retains 80% of fatty acids. For most complemented lines, the analysis revealed a significantly reduced amount of residual fatty acids after 5 d of germination without light or sucrose.
LBs are not considered as passive cell components anymore and interest in the functions they harbor is growing . There is a high degree of variation in the amount and identity of LBPs that have been reported in previous studies even within the same algal species. This becomes obvious when comparing recent LB proteome studies in C. reinhardtii: Different strains were all cultivated photoheterotrophically and while one study did not reveal a single LBP , other groups reported 259  and 248 , with an overlap of less than half. Co-immunoprecipitation with C. reinhardtii MLDP revealed 124 additional proteins that had not been found by the same group in the same strain . This high degree of variation may illustrate the strong influence of sample preparation on the detection of LBPs in proteomics studies.
Filtering of proteomics data revealed true LBPs
In this study, highly sensitive liquid chromatography-tandem mass spectrometry (LC-MS/MS) yielded a large number of putative LBPs. LB isolates from microalgae often contain large amounts of proteins originating from other cellular components [57,58,59,60], which can be due to an interaction of other organelle membranes with the LB surface  or contamination as a result of the cell fractionation process. This explains why extensive filtering of proteomics data was necessary in order to identify the three clearly LB-associated proteins described in this study.
A contamination of the LB fraction with plastoglobules could also be feared, but this seems unlikely. To date, the plastoglobule proteome has been characterized in A. thaliana and the green alga Dunaliella bardawil [62, 63] and for both organisms, only very little overlap with the LB proteome was found. In an electron microscopy study of L. incisa, no plastoglobules could be observed  and like other chlorophytes, the L. incisa genome lacks a homolog of plastoglobulin, which is characteristic of these structures .
The proteins that have been identified as LBPs in other algae so far can be grouped into highly abundant structural proteins [3, 10, 38, 40, 41, 66, 67] and enzymes that play a role in lipid metabolism or trafficking [38, 57]. Confirmation of the subcellular localization is essential in investigating the LB proteome and tobacco pollen tubes are a useful tool to do this . In our study, three proteins could be identified as true LBPs.
While oleosins in higher plants are anchored in the LB surface via a hydrophobic domain that contains a proline knot , no such element can be found in any of the LBPs described here. Mammalian perilipins attach by means of amphipathic helices, which are formed by tandem repeats in the amino acid sequence . A search for such repeats in the LiLBPs described here using HHrepID [68, 69] yielded no results. This is also the case for a range of other proteins that are known to associate with LBs including algal MLDPs [25, 38, 70]. These proteins possibly attach to the LB surface through acylation or by interaction with other proteins.
LiMLDP and LBP62 are LB-associated proteins of unknown function
LiMLDP has similarities to other algal MLDPs: some sequence homology, a highly comparable distribution of hydrophobicity along the amino acid sequence and a clear induction of gene expression in early stages of nitrogen starvation. Unlike Nannochloropsis LDSP (lipid droplet surface protein)  however, in our experiment LiMLDP is unable to complement the high surface/volume ratio of LBs in A. thaliana oleo1 embryos. This may be due to the napA promoter driving gene expression at too low a level or at an inappropriate time during seed maturation. A different role for this protein is also conceivable, for instance in recruitment of enzymes to the LB. Similar functions can be proposed for LiLBP62, which bears no features that could point to a particular physiological role. It is likely to be involved in LB formation or homeostasis rather than LB breakdown, since transcription of the LiLBP62 gene rapidly ceases within 6 h of nitrogen resupply.
LBP36 is a LB-localized putative lipase
In spite of its showing only weak sequence similarity to other known protein, LiLBP36 can be classified as a putative lipase based on predicted structural resemblance to a group of yeast and fungal lipases degrading MAG, DAG and TAG. This group of enzymes includes the first TAG lipase for which a crystal structure was solved and a catalytic triad of serine, histidine and aspartate identified in the active site . Modeling LiLBP36 on the experimentally determined structure of this enzyme permits the identification of catalytic residues, except for the catalytic histidine, which is located in a section of the crystal structure that does not exhibit extensive homology with the modeled LiLBP36 structure.
Expression of the LiLBP36 gene appears to contradict a possible involvement in TAG degradation, as transcription is strongly induced during nitrogen starvation, whereas it decreases as soon as nitrogen is resupplied and lipid reserves are degraded. Nevertheless, gene expression could be induced in anticipation of a future need for TAG hydrolysis, possibly resulting in an inactive form of the protein that can be activated as soon as environmental conditions permit growth once again. A similar reasoning has been proposed to explain the expression pattern of the major TAG lipase in A. thaliana seeds, SDP1 .
In spite of these similarities with TAG lipases, LiLBP36 only effected a minor complementation of postgerminative growth in etiolated A. thaliana sdp1/sdp1-L seedlings. This could be due to a slow translation rate for this protein caused by differences in codon usage, considering that the LiLBP36 coding sequence contains almost 62% guanidine and cytosine residues, while the average for A. thaliana coding sequences is 44% (The Arabidopsis Information Resource, TAIR). Another possible explanation lies in the substrate specificity, as it is extremely difficult to predict the exact substrate range solely based on the modeled structure of a lipase. The distribution of molecular TAG species in L. incisa (mostly 18:1 (n-9) and ARA, ) differs from that in A. thaliana seeds (mostly 18:1 (n-9), 18:2 (n-6), 18:3 (n-3) and 20:1 (n-9), ) and a TAG lipase that acts on one TAG pool will not necessarily act with high efficiency on the other. Among neutral lipids, the enzyme may have a higher affinity for MAG or DAG, neither of which were tested in this study. Alternatively, it may have a primary function in hydrolyzing membrane lipids that form the polar monolayer on the LB surface, thereby permitting subsequent TAG degradation by other lipases.
LiSDP1 is a TAG lipase
The proteomics approach employed in this study proved useful in the identification of LBPs and the LB proteome dataset can serve as a resource for the identification of additional LBPs. We were surprised to find that the L. incisa homolog of SDP1 is not present in our LB protein dataset. We were interested in TAG hydrolysis and therefore went on to analyze this gene in more detail. SDP1 and its paralog SDP1-LIKE have been shown to degrade TAG in A. thaliana [25, 26, 74].
Our analysis of LiSDP1 transcript levels revealed that gene expression is rapidly induced as soon as nitrogen is no longer limiting in the culture and lipid mobilization begins. We observed that the enzyme is in fact able to take over a TAG lipase function when introduced into the A. thaliana sdp1/sdp1-L mutant. The fatty acid profile of seedlings and seeds did not differ significantly between WT and transgenic lines, pointing to a broad substrate specificity of L. incisa SDP1 as reported for its A. thaliana homolog . The overall level of TAG lipase activity exerted by the algal enzyme appeared low, possibly owing to a specialization for TAG species that are absent from A. thaliana seeds. The effect on postgerminative growth of etiolated seedlings was most pronounced when expression was controlled by the 35S promoter, which has been shown to be equally effective as the endogenous A. thaliana sdp1 promoter in complementing the sdp1 mutant . The napA promoter effected a weaker complementation, which may be due to a lower level or an unfavorable temporal course of gene expression, even though the promoter is reportedly active not only during seed maturation but also during germination in Brassica oleracea .
The observations that LiSDP1 was not present in our LB isolate from nitrogen starved algae and that it does not localize to the LBs in tobacco pollen tubes suggests that the TAG lipase may be among the proteins that are specifically recruited to the LBs in L. incisa when needed. It has been shown for C. reinhardtii that the LB proteome is subject to vast changes in response to nitrogen resupply . However, recent evidence obtained for A. thaliana seedlings may support our findings. Here the picture emerges that SDP1 may be rather localized at the peroxisomal membrane and is only translocated to the LB membrane upon TAG mobilization [76, 77]. The gene expression pattern may support this notion, as LiSDP1 expression is only induced in response to nitrogen resupply.
L. incisa is an oleaginous microalga that holds great promise for the study of lipid metabolism. In this study, an extensive LB proteome dataset was generated and successfully filtered for true LBPs. N. tabacum pollen tubes were used as an efficient tool for the verification of LB localization in our experimental setup. Two LBPs of unknown function are suspected to play a role in maintaining the structural integrity of the organelle and preventing its premature degradation, while another is a promising candidate for lipase activity that could be required for storage lipid mobilization. Failure to confirm the LB localization of other candidate LBPs that were investigated here emphasizes the necessity of verifying results from cell fractionation. At the same time, the use of a heterologous system such as tobacco pollen possibly entails differences in targeting signal recognition or a lack of certain scaffolding proteins required for LB localization. Nevertheless, the dataset generated here will be a useful resource for the identification of further LBPs, possibly revealing novel aspects of LB metabolism and homeostasis.
Screening the L. incisa genome for further lipase candidates revealed a homolog of the major TAG lipase in A. thaliana seedlings, LiSDP1. The algal gene was able to reconstitute hypocotyl elongation in etiolated seedlings lacking endogenous SDP1 and SDP1-L, verifying TAG lipase activity for this enzyme. The enzyme does not localize to the LBs in nitrogen-starved L. incisa or the tobacco pollen tube system, instead it is likely to be recruited to the sites of TAG storage when these reserves need to be mobilized. The robust cell wall of L. incisa is a major obstacle in the isolation of relatively pure organelles.
L. incisa Genome sequencing
Assembly of the nuclear genome of L. incisa was carried out along with that of the mitochondrial and plastidial genomes as described previously . Briefly, Illumina paired-end and Long Jumping Distance mate pair reads were assembled using the super-read-based assembler MSR-CA (MaSuRCA) . Protein-coding genes were predicted with AUGUSTUS 2.6.1 (PMID:16469098)  and annotated using InterProScan 5.9 (PMID:24451626) , USEARCH/UBLAST 7.0.1001 (Edgar 2010, Bioinformatics 26, 2460–2461) and KEGG (using the KAAS server: Moriya et al. 2007, Nucleic Acids Research, 35, W182–W185). The sequences and annotations are available at https://giavap-genomes.ibpc.fr. Please contact email@example.com for authorizations.
Liquid culture of L. incisa
L. incisa strain SAG 2468 had originally been isolated on a Japanese glacier  and was kindly provided by Dr. Inna Khozin-Goldberg, Ben-Gurion University of the Negev, Israel. It was cultivated in BG11 media  in 300 mL volumes at 25 °C and subject to continuous illumination with 190 μmol photons m−2s−1 as well as aeration with a supplement of 1% (v/v) CO2. Nitrogen starvation was achieved by washing and resuspending the cells in BG11 media which was modified by omission of NaNO3 and replacement of ammonium ferric citrate with ferric citrate . Nitrogen was resupplied by sedimenting cells and once again taking them up in full BG11 media.
Fatty acids in algal material were quantified by GC as previously described  with minor modifications. Following lyophilization of the material, lipids were isolated by Methyl tert-butyl ether (MTBE) extraction according to  and fatty acids were derivatized by acidic methanolysis . The temperature gradient employed in GC measurements comprised the following steps: 150 °C for 1 min, 150 °C to 200 °C at 4 °C/min, 200 °C to 250 °C at 20 °C/min and 250 °C for 3 min.
A. thaliana seeds were dried at 60 °C over night and frozen seedlings (including seed coats) were lyophilized prior to direct acidic methanolysis . Tri-15:0 was included as an internal standard for quantification and samples were mechanically disrupted during the reaction using a metal rod to ensure efficient methanolysis.
Centrifugation steps were generally carried out at 4 °C in an Eppendorf 5810R centrifuge, an Optima LE-80 K ultracentrifuge with an SW40 rotor or an Optima TLX ultracentrifuge with a TLS55 rotor for 50 mL tubes, 12 mL thin wall polypropylene tubes and 2.5 mL thin wall polypropylene tubes, respectively. Ice-cold buffers were used and samples were kept on ice between fractionation steps.
An L. incisa culture was starved of nitrogen for 3 days prior to LB isolation. Cells from 50 mL of culture were sedimented by centrifugation at 2500 x g for 10 min, washed once with distilled water and then ground in liquid nitrogen. LBs were subsequently isolated according to  with some modifications. 60 mL centrifugation buffer A (100 mM Tris-HCl pH 7.5, 3 mM ethylenediaminetetraacetic acid (EDTA), 10 mM Dithiotreitol (DTT), 1% (v/v) Plant Protease Inhibitor Cocktail (Sigma-Aldrich), 0.6 M sucrose) were added to the homogenized material and a 200 μL sample of the total cell extract was taken before removal of cell debris by centrifugation at 10,000 x g for 10 min. The entire supernatant was transferred to 12 mL tubes, each was carefully overlain with 10 mL centrifugation buffer B (buffer A with 0.4 M sucrose) and centrifugation was repeated. LBs floating on top were removed with a spatula and the remaining supernatant was separated into membranes and a soluble fraction by ultracentrifugation at 105,000 x g for 90 min. A 200 μL sample of the soluble fraction was taken and the sedimented membranes were combined. The LBs were transferred to a Potter-Elvehjem tissue grinder and carefully resuspended in 10 mL centrifugation buffer A. Overlaying with 8 mL buffer B and another centrifugation step at 10,000 x g yielded an LB fraction floating on top, which was then washed this way two more times. Finally, LBs were resuspended in 1.5 mL buffer A, transferred to a 2.5 mL ultracentrifugation tube and overlain with buffer B before ultracentrifugation at 100,000 x g for 1 h.
LB protein identification
Proteins were extracted from each cell fraction without further separation as described previously [85, 86] with some modifications. They were precipitated in 90% ethanol at −80 °C for 2 h and sedimented at 20,000 x g and 4 °C for 15 min followed by washing with 80% ethanol three times. Proteins were solubilized in 100 μL denaturing protein solubilization buffer (4% sodium dodecyl sulfate (SDS), 4 mM DTT, 8% (v/v), 80 mM Tris-HCl pH 6.8, 0.02% (w/v) Bromophenolblue, 7 M urea, 2 M thiourea [59, 87] at 37 °C for 2 h before applying each sample in triplicate to SDS-PAGE [87,88,89,90] until it had migrated 1 cm. The gel was stained with Coomassie [91, 92], destained in water and the entire band was excised. The proteins were then subjected to an in-gel tryptic digest as described previously  and peptides were identified by liquid chromatography coupled to tandem mass spectrometry (see Additional file 1 for details). Proteins were then identified by comparison to a list of all L. incisa proteins that had been predicted based on the nuclear and plastidial genomes.
Proteins that were detected in all three technical replicates of the LB sample were considered for further analysis and for each one, abundance was estimated by calculating the Normalized Spectral Abundance Factor (NSAF, [94, 95]). Enrichment compared to the three control samples (total extract, membranes, soluble fraction) was determined by dividing NSAFs.
Protein in silico analyses
Coding sequence amplification and plasmid construction
RNA was isolated from an L. incisa culture that had been deprived of nitrogen for 3 days and cDNA was synthesized (see below). Coding sequences of interest were amplified and restriction sites were added using the primers listed in Additional file 8 before ligation into a subcloning vector using the CloneJET PCR cloning kit (Thermo Fisher Scientific).
The coding sequences without a stop codon were inserted into the pUC-LAT52-mVenus vector [107, 108] for expression in tobacco pollen tubes by restriction cloning using endonucleases obtained from Thermo Fisher Scientific. Similarly, the full coding sequences were inserted into the pEntry-E vector  for further transfer into pCambia plant expression vector 43.0 by Gateway cloning technology. pCambia 43.0 had been as a new expression vector with mCherry as a reporter gene in plants under the control of the seed specific napin promoter from Brassica napus. To generate the vector, a napin::mCherry construct was amplified via PCR, creating restriction sites SacI and XhoI at the 5′ and 3′-end. pCambia33.0G was digested with SacI and XhoI and ligated with the napin::mCherry construct to substitute for the bar gene under the control of the 35S promoter. For vector details see Additional file 9.
Transient gene expression in tobacco pollen tubes
Transient expression in tobacco pollen tubes was used to verify the LB localization of LBP candidates as previously described . N. tabacum (ecotype Samsun NN) pollen grains were transiently transformed by particle bombardment as described previously and allowed to germinate on glass slides for 6 h. Pollen tubes were then fixed with PT media  containing formaldehyde (1% (w/v) final concentration) and stained for LBs with Nile Red (0.17 ng/μL final concentration).
Fluorescence images were obtained using a Zeiss LSM Meta confocal microscope. Nile Red and mVenus were excited at 561 and 488 nm, imaged with HFT 405/488/561 and HFT 405/514/633 nm major beam splitters and fluorescence was detected at 571-593 and 507-539 nm, respectively.
Stable transformation of A. thaliana
The A. thaliana mutants sdp1/sdp1-L ) and oleo1 (SM_3.29875, kindly provided along with the corresponding wildtype by Prof. Dr. Christoph Benning, Michigan State University) were transformed with pCambia vectors containing algal genes along with the MCHERRY reporter gene by the floral dip method using Agrobacterium tumefaciens strain EHA 105 . The mutants as well as the Col-0 wildtype were also transformed with an empty vector to generate lines to be used as negative controls. T1 and T2 seeds expressing the reporter gene were identified by fluorescence microscopy using an M165C stereomicroscope equipped with an EL6000 mercury lamp, a M205FA/M165FC filter set and a DFC3000 G camera (all from Leica Microsystems). Gene expression was confirmed as described below.
Gene expression analyses
For RNAseq, L. incisa cells grown in control conditions or nitrogen-starved for 12 or 72 h, either in low light (75 μmol photons.m−2.s−1) or in high light (150 μmol photons.m−2.s−1), from two biological replicates. RNA was isolated from frozen samples using the SV Total RNA isolation kit (Promega) after breaking the material with iron beads in liquid nitrogen. RNA quality was examined on a 2100 Electrophoresis Bioanalyzer. The transcriptome was sequenced using the Illumina Truseq high-throughput short-read technology (using 50-nt single-end and 100-nt paired-end reads) on a HiSeq 1000 instrument by the transcriptomics platform of the Institut de Biologie de l’École Normale Supérieure (Paris, France). For each growth condition, two biological replicate samples were sequenced, each in four technical replicates. They were mapped onto the transcripts and analyzed for differential expression using the protocol of Haas et al.  using tools from the Trinity (release r20131110), RSEM 1.2.9 , and DESeq , the latter also providing an analysis of deviance (ANODEV). A false discovery rate of 5% was used for cutoff.
For selected genes, RNA-Seq results were confirmed and expanded by means of qRT-PCR. For this purpose, total RNA was isolated from three L. incisa cultures at various time points during nitrogen starvation as well as following nitrogen resupply. Samples were freeze-dried and subsequently ground in liquid nitrogen to ensure efficient cell disruption prior to Trizol extraction . DNase I and RevertAid H Minus (Thermo Fisher Scientific) were then used to remove residual genomic DNA and to synthesize cDNA. Primer3Prefold  and Primer3Plus  were used for primer design (see Additional file 8 for primer sequences) and RIBOSOMAL PROTEIN S21 was used as a reference gene for normalization. qRT-PCR was carried out with an iQ5 qPCR cycler (Biorad Laboratories) and the Takyon No Rox SYBR Core Kit blue dTTP (Eurogentec).
In dry seeds of A. thaliana, gene expression was confirmed by extracting RNA as previously described , generating cDNA as outlined above and detecting transcripts by PCR.
Postgerminative growth assay
Hypocotyl length was used to evaluate postgerminative growth in etiolated seedlings of A. thaliana sdp1/sdp1-L mutant lines expressing putative TAG lipase genes as previously described . Seed-specific or constitutive transgene expression were controlled by the Brassica napus napin A promoter or the Cauliflower Mosaic Virus 35S promoter, respectively. T2 seeds were preselected for a size of 250 – 300 μm, surface sterilized and lined up on ½ MS agar plates with or without 1% sucrose. Reporter gene expression was documented by fluorescence microscopy and germination was synchronized by stratification at 4 °C for 2 days. Seeds were then exposed to light for 30 min to induce germination and afterwards kept at room temperature in the dark for 5 days. At this point seedlings were documented using a scanner and for those that originated from transgenic seeds, hypocotyl length was measured. For each line, 4 batches of at least 10 seedlings were measured. Gas chromatography was used to determine the TFA content for these seedlings (including seed coats) as well as transgenic seeds of the same line, using the configurations as detailed above. For each line, 4 batches of at least 10 seedlings were analyzed and the TFA content in each batch was divided by the number of seedlings. The 4 values for TFA per seedling were then used to calculate the mean value and standard error of the mean for each line. Transgenic seeds were analyzed in the same way.
LB size assay
The impact of algal proteins on the structural integrity of LBs was investigated by expressing the corresponding genes in the A. thaliana mutant oleo1 and observing LB size in embryo hypocotyl cells. T2 seeds were preselected for reporter gene expression, imbibed in 0.1% agar at 4 °C over night and embryos were isolated by gentle squeezing between two glass slides. 100 μL of 0.1% agar containing the embryos was supplemented with 100 ng Nile Red (dissolved in acetone) and incubated in the dark at room temperature for 10 min before confocal microscopy as outlined above. Only embryos that appeared fully intact were observed more closely and the hypocotyl cells of three embryos per line were analyzed. The intensity of fluorescence emitted by the MCHERRY protein was negligible compared to the highly intense Nile Red stain, therefore permitting clear observation of stained LBs in the red part of the spectrum.
- A. thaliana :
Domain of unknown function
Fragments per kilobase of transcript per million mapped reads
- L. incisa :
Lipid body protein
Liquid chromatography-tandem mass spectrometry
Lipid droplet associated protein
Lipid droplet surface protein
Major lipid droplet protein
Methyl tert-butyl ether
- N. tabacum :
Normalized spectral abundance factor
Oil globule protein
Sugar dependent 1
Sodium dodecyl sulfate
Total fatty acids
Zienkiewicz K, Du ZY, Ma W, Vollheyde K, Benning C. Stress-induced neutral lipid biosynthesis in microalgae - molecular, cellular and physiological insights. Biochim Biophys Acta. 2016;1861(9 Pt B):1269–81.
Cheng-Wu Z, Cohen Z, Khozin-Goldberg I, Richmond A. Characterization of growth and arachidonic acid production of Parietochloris incisa comb. Nov (Trebouxiophyceae, Chlorophyta). J Appl Phycol. 2002;14(6):453–60.
Davidi L, Katz A, Pick U. Characterization of major lipid droplet proteins from Dunaliella. Planta. 2012;236(1):19–33.
Msanne J, Xu D, Konda AR, Casas-Mollano JA, Awada T, Cahoon EB, et al. Metabolic and gene expression changes triggered by nitrogen deprivation in the photoautotrophically grown microalgae Chlamydomonas reinhardtii and Coccomyxa sp. C-169. Phytochemistry. 2012;75:50–9.
Tsai CH, Warakanont J, Takeuchi T, Sears BB, Moellering ER, Benning C. The protein compromised hydrolysis of Triacylglycerols 7 (CHT7) acts as a repressor of cellular quiescence in Chlamydomonas. Proc Natl Acad Sci U S A. 2014;111(44):15833–8.
Schmollinger S, Muhlhaus T, Boyle NR, Blaby IK, Casero D, Mettler T, et al. Nitrogen-sparing mechanisms in Chlamydomonas affect the transcriptome, the proteome, and photosynthetic metabolism. Plant Cell. 2014;26(4):1410–35.
Longworth J, Wu D, Huete-Ortega M, Wright PC, Vaidyanathan S. Proteome response of Phaeodactylum tricornutum, during lipid accumulation induced by nitrogen depletion. Algal Res. 2016;18:213–24.
Miller R, Wu G, Deshpande RR, Vieler A, Gartner K, Li X, et al. Changes in transcript abundance in Chlamydomonas reinhardtii following nitrogen deprivation predict diversion of metabolism. Plant Physiol. 2010;154(4):1737–52.
Allen JW, DiRusso CC, Black PN. Triacylglycerol synthesis during nitrogen stress involves the prokaryotic lipid synthesis pathway and acyl chain remodeling in the microalgae Coccomyxa subellipsoidea. Algal Res. 2015;10:110–20.
Peled E, Leu S, Zarka A, Weiss M, Pick U, Khozin-Goldberg I, et al. Isolation of a novel oil globule protein from the green alga Haematococcus pluvialis (Chlorophyceae). Lipids. 2011;46(9):851–61.
Chapman KD, Ohlrogge JB. Compartmentation of triacylglycerol accumulation in plants. J Biol Chem. 2012;287(4):2288–94.
Yatsu LY, Jacks TJ. Spherosome membranes: half unit-membranes. Plant Physiol. 1972;49(6):937–43.
Tzen JTC, Lie GC, Huang AHC. Characterization of the charged components and their topology on the surface of plant seed oil bodies. J Biol Chem. 1992;267(22):15626–34.
Abell BM, Holbrook LA, Abenes M, Murphy DJ, Hills MJ, Moloney MM. Role of the proline knot motif in oleosin endoplasmic reticulum topology and oil body targeting. Plant Cell. 1997;9(8):1481–93.
Lin LJ, Tai SS, Peng CC, Tzen JT. Steroleosin, a sterol-binding dehydrogenase in seed oil bodies. Plant Physiol. 2002;128(4):1200–11.
May C, Hohne M, Gnau P, Schwennesen K, Kindl H. The N-terminal beta-barrel structure of lipid body lipoxygenase mediates its binding to liposomes and lipid bodies. Eur J Biochem. 2000;267(4):1100–9.
Najt CP, Lwande JS, McIntosh AL, Senthivinayagam S, Gupta S, Kuhn LA, et al. Structural and functional assessment of perilipin 2 lipid binding domain(s). Biochemistry. 2014;53(45):7051–66.
Rowe ER, Mimmack ML, Barbosa AD, Haider A, Isaac I, Ouberai MM, et al. Conserved amphipathic helices mediate lipid droplet targeting of perilipins 1-3. J Biol Chem. 2016;291(13):6664–78.
Jolivet P, Acevedo F, Boulard C, d'Andrea S, Faure JD, Kohli A, et al. Crop seed oil bodies: from challenges in protein identification to an emerging picture of the oil body proteome. Proteomics. 2013;13(12-13):1836–49.
Goold H, Beisson F, Peltier G, Li-Beisson Y. Microalgal lipid droplets: composition, diversity, biogenesis and functions. Plant Cell Rep. 2015;34(4):545–55.
Athenstaedt K, Daum G. YMR313c/TGL3 encodes a novel triacylglycerol lipase located in lipid particles of Saccharomyces cerevisiae. J Biol Chem. 2003;278(26):23317–23.
Athenstaedt K, Daum G. Tgl4p and Tgl5p, two triacylglycerol lipases of the yeast Saccharomyces cerevisiae are localized to lipid particles. J Biol Chem. 2005;280(45):37301–9.
Eastmond PJ. Cloning and characterization of the acid lipase from castor beans. J Biol Chem. 2004;279(44):45540–5.
Lin YH, Huang AH. Lipase in lipid bodies of cotyledons of rape and mustard seedlings. Arch Biochem Biophys. 1983;225(1):360–9.
Eastmond PJ. SUGAR-DEPENDENT1 encodes a patatin domain triacylglycerol lipase that initiates storage oil breakdown in germinating Arabidopsis seeds. Plant Cell. 2006;18(3):665–75.
Kelly AA, Quettier AL, Shaw E, Eastmond PJ. Seed storage oil mobilization is important but not essential for germination or seedling establishment in Arabidopsis. Plant Physiol. 2011;157(2):866–75.
Barka F, Angstenberger M, Ahrendt T, Lorenzen W, Bode HB, Büchel C. Identification of a triacylglycerol lipase in the diatom Phaeodactylum tricornutum. Biochim Biophys Acta. 2016;1861(3):239–48.
Vance VB, Huang AH. The major protein from lipid bodies of maize. Characterization and structure based on cDNA cloning. J Biol Chem. 1987;262(23):11275–9.
Frandsen GI, Mundy J, Tzen JTC. Oil bodies and their associated proteins, oleosin and caleosin. Physiol Plantarum. 2001;112(3):301–7.
Siloto RM, Findlay K, Lopez-Villalobos A, Yeung EC, Nykiforuk CL, Moloney MM. The accumulation of oleosins determines the size of seed oilbodies in Arabidopsis. Plant Cell. 2006;18(8):1961–74.
Shimada TL, Shimada T, Takahashi H, Fukao Y, Hara-Nishimura I. A novel role for oleosins in freezing tolerance of oilseeds in Arabidopsis thaliana. Plant J. 2008;55(5):798–809.
Deruyffelaere C, Bouchez I, Morin H, Guillot A, Miquel M, Froissard M, et al. Ubiquitin-mediated proteasomal degradation of oleosins is involved in oil body mobilization during post-germinative seedling growth in Arabidopsis. Plant Cell Physiol. 2015;56(7):1374–87.
Sztalryd C, Kimmel AR. Perilipins: lipid droplet coat proteins adapted for tissue-specific energy storage and utilization, and lipid cytoprotection. Biochimie. 2014;96(1):96–101.
Wang H, Bell M, Sreenivasan U, Hu H, Liu J, Dalen K, et al. Unique regulation of adipose triglyceride lipase (ATGL) by perilipin 5, a lipid droplet-associated protein. J Biol Chem. 2011;286(18):15707–15.
Gronke S, Beller M, Fellert S, Ramakrishnan H, Jäckle H, Kühnlein RP. Control of fat storage by a drosophila PAT domain protein. Curr Biol. 2003;13(7):603–6.
Horn PJ, James CN, Gidda SK, Kilaru A, Dyer JM, Mullen RT, et al. Identification of a new class of lipid droplet-associated proteins in plants. Plant Physiol. 2013;162(4):1926–36.
Gidda SK, Park S, Pyc M, Yurchenko O, Cai Y, Wu P, et al. Lipid droplet-associated proteins (LDAPs) are required for the dynamic regulation of neutral lipid compartmentation in plant cells. Plant Physiol. 2016;170(4):2052–71.
Moellering ER, Benning C. RNA interference silencing of a major lipid droplet protein affects lipid droplet size in Chlamydomonas reinhardtii. Eukaryot Cell. 2010;9(1):97–106.
Tsai CH, Zienkiewicz K, Amstutz CL, Brink BG, Warakanont J, Roston R, et al. Dynamics of protein and polar lipid recruitment during lipid droplet assembly in Chlamydomonas reinhardtii. Plant J. 2015;83(4):650–60.
Yoneda K, Yoshida M, Suzuki I, Watanabe MM. Identification of a major lipid droplet protein in a marine diatom Phaeodactylum tricornutum. Plant Cell Physiol. 2016;57(2):397–406.
Vieler A, Brubaker SB, Vick B, Benning C. A lipid droplet protein of Nannochloropsis with functions partially analogous to plant oleosins. Plant Physiol. 2012;158(4):1562–9.
Watanabe S, Hirabayashii S, Boussiba S, Cohen Z, Vonshak A, Richmond A. Parietochloris incisa comb. Nov. (Trebouxiophyceae, Chlorophyta). Phycol Res. 1996;44(2):107–8.
Bigogno C, Khozin-Goldberg I, Boussiba S, Vonshak A, Cohen Z. Lipid and fatty acid composition of the green oleaginous alga Parietochloris incisa, the richest plant source of arachidonic acid. Phytochemistry. 2002;60(5):497–503.
Hu Q, Sommerfeld M, Jarvis E, Ghirardi M, Posewitz M, Seibert M, et al. Microalgal triacylglycerols as feedstocks for biofuel production: perspectives and advances. Plant J. 2008;54(4):621–39.
Khozin-Goldberg I, Bigogno C, Shrestha P, Cohen Z. Nitrogen starvation induces the accumulation of arachidonic acid in the freshwater green alga Parietochloris incisa (Trebuxiophyceae). J Phycol. 2002;38(5):991–4.
Tourasse NJ, Shtaida N, Khozin-Goldberg I, Boussiba S, Vallon O. The complete mitochondrial genome sequence of the green microalga Lobosphaera (Parietochloris) incisa reveals a new type of palindromic repetitive repeat. BMC Genomics. 2015;16(1):580.
Tourasse NJ, Barbi T, Waterhouse JC, Shtaida N, Leu S, Boussiba S, et al. The complete sequence of the chloroplast genome of the green microalga Lobosphaera (Parietochloris) incisa. Mitochondrial DNA A DNA Mapp Seq Anal. 2016;27(3):2073–5.
Iskandarov U, Khozin-Goldberg I, Cohen Z. Selection of a DGLA-producing mutant of the microalga Parietochloris incisa: I. Identification of mutation site and expression of VLC-PUFA biosynthesis genes. Appl Microbiol Biotechnol. 2011;90(1):249–56.
Zorin B, Grundman O, Khozin-Goldberg I, Leu S, Shapira M, Kaye Y, et al. Development of a nuclear transformation system for oleaginous green alga Lobosphaera (Parietochloris) incisa and genetic complementation of a mutant strain, deficient in arachidonic acid biosynthesis. PLoS One. 2014;9(8):e105223.
Shemesh Z, Leu S, Khozin-Goldberg I, Didi-Cohen S, Zarka A, Boussiba S. Inducible expression of Haematococcus oil globule protein in the diatom Phaeodactylum tricornutum: association with lipid droplets and enhancement of TAG accumulation under nitrogen starvation. Algal Res. 2016;18:321–31.
Müller AO, Blersch KF, Gippert AL, Ischebeck T. Tobacco pollen tubes - a fast and easy tool to study lipid droplet association of plant proteins. Plant J. 2016; doi:10.1111/tpj.13441.
The Pfam web interface [http://pfam.xfam.org/].
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(Database issue):D222–30.
Derewenda ZS, Derewenda U. Relationships among serine hydrolases: evidence for a common structural motif in triacylglyceride lipases and esterases. Biochem Cell Biol. 1991;69(12):842–51.
Kienesberger PC, Oberer M, Lass A, Zechner R. Mammalian patatin domain containing proteins: a family with diverse lipolytic activities involved in multiple biological functions. J Lipid Res. 2009;50(Suppl):S63–8.
Wang ZT, Ullrich N, Joo S, Waffenschmidt S, Goodenough U. Algal lipid bodies: stress induction, purification, and biochemical characterization in wild-type and starchless Chlamydomonas reinhardtii. Eukaryot Cell. 2009;8(12):1856–68.
Nguyen HM, Baudet M, Cuine S, Adriano JM, Barthe D, Billon E, et al. Proteomic profiling of oil bodies isolated from the unicellular green microalga Chlamydomonas reinhardtii: with focus on proteins involved in lipid metabolism. Proteomics. 2011;11(21):4266–73.
Ding Y, Zhang S, Yang L, Na H, Zhang P, Zhang H, et al. Isolating lipid droplets from multiple species. Nat Protoc. 2013;8(1):43–51.
Zienkiewicz A, Rejón JD, Alché JDD, Rodríguez-García MI, Castro AJ. A protocol for protein extraction from lipid-rich plant tissues suitable for electrophoresis. In: Jorrin-Novo JV, Komatsu S, Weckwerth W, Wienkoop S, editors. Plant proteomics: methods and Protocolls, vol. 1072. Totowa: Humana Press; 2014. p. 85–91.
Nojima D, Yoshino T, Maeda Y, Tanaka M, Nemoto M, Tanaka T. Proteomics analysis of oil body-associated proteins in the oleaginous diatom. J Proteome Res. 2013;12(11):5293–301.
Barbosa AD, Savage DB, Siniossoglou S. Lipid droplet-organelle interactions: emerging roles in lipid metabolism. Curr Opin Cell Biol. 2015;35:91–7.
Davidi L, Levin Y, Ben-Dor S, Pick U. Proteome analysis of cytoplasmatic and plastidic beta-carotene lipid droplets in Dunaliella bardawil. Plant Physiol. 2015;167(1):60–79.
Lundquist PK, Poliakov A, Bhuiyan NH, Zybailov B, Sun Q, van Wijk KJ. The functional network of the Arabidopsis plastoglobule proteome based on quantitative proteomics and genome-wide coexpression analysis. Plant Physiol. 2012;158(3):1172–92.
Merzlyak MN, Chivkunova OB, Gorelova OA, Reshetnikova IV, Solovchenko AE, Khozin-Goldberg I, et al. Effect of nitrogen starvation on optical properties, pigments, and arachidonic acid content of the unicellular green alga Parietochloris incisa (Trebouxiophyceae, Chlorophyta). J Phycol. 2007;43(4):833–43.
Brehelin C, Kessler F, van Wijk KJ. Plastoglobules: versatile lipoprotein particles in plastids. Trends Plant Sci. 2007;12(6):260–6.
Sharma KK, Li Y, Schenk PM. Rapid lipid induction in Chlorella sp by UV-C radiation. Bioenerg Res. 2015;8(4):1824–30.
Lin IP, Jiang PL, Chen CS, Tzen JT. A unique caleosin serving as the major integral protein in oil bodies isolated from Chlorella sp. cells cultured with limited nitrogen. Plant Physiol Biochem. 2012;61:80–7.
The HHrepID web interface [http://omictools.com/hhrepid-tool].
Biegert A, Soding J. De novo identification of highly diverged protein repeats by probabilistic consistency. Bioinformatics. 2008;24(6):807–14.
Ploier B, Scharwey M, Koch B, Schmidt C, Schatte J, Rechberger G, et al. Screening for hydrolytic enzymes reveals Ayr1p as a novel triacylglycerol lipase in Saccharomyces cerevisiae. J Biol Chem. 2013;288(50):36061–72.
Brady L, Brzozowski AM, Derewenda ZS, Dodson E, Dodson G, Tolley S, et al. A serine protease triad forms the catalytic centre of a triacylglycerol lipase. Nature. 1990;343(6260):767–70.
Bigogno C, Khozin-Goldberg I, Cohen Z. Accumulation of arachidonic acid-rich triacylglycerols in the microalga Parietochloris incisa (Trebuxiophyceae, Chlorophyta). Phytochemistry. 2002;60(2):135–43.
Taylor DC, Giblin EM, Reed DW, Hogge LR. Stereospecific analysis and mass spectrometry of triacylglycerols from Arabidopsis thaliana (L.) heynh. columbia seed. J Am Oil Chem Soc. 1995;72(3):305–8.
Kelly AA, van Erp H, Quettier AL, Shaw E, Menard G, Kurup S, et al. The SUGAR-DEPENDENT1 lipase limits triacylglycerol accumulation in vegetative tissues of Arabidopsis. Plant Physiol. 2013;162(3):1282–9.
Soeda Y, Konings MC, Vorst O, van Houwelingen AM, Stoopen GM, Maliepaard CA, et al. Gene expression programs during Brassica oleracea seed maturation, osmopriming, and germination are indicators of progression of the germination process and the stress tolerance level. Plant Physiol. 2005;137(1):354–68.
Cui S, Hayashi Y, Otomo M, Mano S, Oikawa K, Hayashi M, et al. Sucrose production mediated by lipid metabolism suppresses the physical interaction of Peroxisomes and oil bodies during germination of Arabidopsis thaliana. J Biol Chem. 2016;291(38):19734–45.
Thazar-Poulot N, Miquel M, Fobis-Loisy I, Gaude T. Peroxisome extensions deliver the Arabidopsis SDP1 lipase to oil bodies. Proc Natl Acad Sci U S A. 2015;112(13):4158–63.
Zimin AV, Marçais G, Puiu D, Roberts M, Salzberg SL, Yorke JA. The MaSuRCA genome assembler. Bioinformatics. 2013;29(21):2669–77.
Stanke M, Steinkamp R, Waack S, Morgenstern B. AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 2004;32(Web Server issue):W309–12.
Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.
Stanier RY, Kunisawa R, Mandel M, Cohen-Bazire G. Purification and properties of unicellular blue-green algae (order Chroococcales). Bacteriol Rev. 1971;35(2):171–205.
Lang I, Hodac L, Friedl T, Feussner I. Fatty acid profiles and their distribution patterns in microalgae: a comprehensive analysis of more than 2000 strains from the SAG culture collection. BMC Plant Biol. 2011;11(1):124.
Matyash V, Liebisch G, Kurzchalia TV, Shevchenko A, Schwudke D. Lipid extraction by methyl-tert-butyl ether for high-throughput lipidomics. J Lipid Res. 2008;49(5):1137–46.
Miquel M, Browse J. Arabidopsis mutants deficient in polyunsaturated fatty acid synthesis. Biochemical and genetic characterization of a plant oleoyl-phosphatidylcholine desaturase. J Biol Chem. 1992;267(3):1502–9.
Feussner I, Kindl H. A lipoxygenase is the main lipid body protein in cucumber and soybean cotyledons during the stage of triglyceride mobilization. FEBS Lett. 1992;298(2-3):223–5.
Valledor L, Weckwerth W: An improved detergent-compatible gel-fractionation LC-LTQ-Orbitrap-MS workflow for plant and microbial proteomics. In. Edited by Jorrin-Novo JV, Komatsu S, Weckwerth W, Wienkoop S, 2 edn; 2014: 85-91.
Laemmli UK. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970;227(5259):680–5.
Kellenberger E. The components of the T4 capsid and of other, capsid-related structures. Virology. 1968;34(3):549–61.
Davis BJ. Disc electrophoresis. II. Method and application to human serum proteins. Ann N Y Acad Sci. 1964;121:404–27.
Ornstein L. Disc electrophoresis. I. Background and theory. Ann N Y Acad Sci. 1964;121:321–49.
Meyer TS, Lamberts BL. Use of coomassie brilliant blue R250 for the electrophoresis of microgram quantities of parotid saliva proteins on acrylamide-gel strips. Biochim Biophys Acta. 1965;107(47):144–5.
Neuhoff V, Arold N, Taube D, Ehrhardt W. Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie brilliant blue G-250 and R-250. Electrophoresis. 1988;9(6):255–62.
Shevchenko A, Tomas H, Havlis J, Olsen JV, Mann M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat Protoc. 2006;1(6):2856–60.
Powell DW, Weaver CM, Jennings JL, McAfee KJ, He Y, Weil PA, et al. Cluster analysis of mass spectrometry data reveals a novel component of SAGA. Mol Cell Biol. 2004;24(16):7249–59.
Zybailov B, Mosley AL, Sardiu ME, Coleman MK, Florens L, Washburn MP. Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiae. J Proteome Res. 2006;5(9):2339–47.
The BLASTP web interface [http://blast.ncbi.nlm.nih.gov/BlastAlign.cgi].
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
The Phyre2 web interface [http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index].
Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJ. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10(6):845–58.
The ExPASy ProtScale web interface [http://web.expasy.org/protscale/].
Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, et al. Protein identification and analysis tools on the ExPASy server. In: Walker JM, editor. The proteomics protocols handbook. New York: Humana Press; 2005. p. 571–607.
Kyte J, Doolittle RF. A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982;157(1):105–32.
The PhyML web interface [http://www.trex.uqam.ca/index.php?action=phyml2].
Boc A, Diallo AB, Makarenkov V. T-REX: a web server for inferring, validating and visualizing phylogenetic trees and networks. Nucleic Acids Res. 2012;40(Web Server issue):W573–9.
The Clustal Omega web interface [https://www.ebi.ac.uk/Tools/msa/clustalo/].
Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.
Mahs A, Steinhorst L, Han JP, Shen LK, Wang Y, Kudla J. The calcineurin B-like Ca2+ sensors CBL1 and CBL9 function in pollen germination and pollen tube growth in Arabidopsis. Mol Plant. 2013;6(4):1149–62.
Twell D, Wing R, Yamaguchi J, McCormick S. Isolation and expression of an anther-specific gene from tomato. Mol Gen Genet. 1989;217(2-3):240–5.
Heilmann M, Iven T, Ahmann K, Hornung E, Stymne S, Feussner I. Production of wax esters in plant seed oils by oleosomal cotargeting of biosynthetic enzymes. J Lipid Res. 2012;53(10):2153–61.
Read SM, Clarke AE, Bacic A. Stimulation of growth of cultured Nicotiana tabacum W 38 pollen tubes by poly(ethylene glycol) and Cu(II) salts. Protoplasma. 1993;177:1–14.
Hood EE, Gelvin SB, Melchers LS, Hoekema A. New Agrobacterium helper plasmids for gene transfer to plants. Transgenic Res. 1993;2(4):208–18.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protocols. 2013;8(8):1494–512.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12(1):323.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987;162(1):156–9.
Markham NR, Zuker M. UNAFold: software for nucleic acid folding and hybridization. Methods Mol Biol. 2008;453(1):3–31.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3-new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):1–12.
Onate-Sanchez L, Vicente-Carbajosa J. DNA-free RNA isolation protocols for Arabidopsis thaliana, including seeds and siliques. BMC Res Notes. 2008;1:93.
We are grateful to Prof. Dr. Volker Lipka and Dr. Hassan Ghareeb (Department for Molecular Biology of Plant-Microbe Interactions) for providing the fluorescence stereomicroscope to select transgenic seeds. We also thank Prof. Dr. Steven Johnsen and Dr. Florian Wegwitz (Clinic for General, Visceral and Pediatric Surgery) for the opportunity to use the confocal microscope as well as Prof. Dr. Jörg Kudla and Dr. Leonie Steinhorst (University of Münster) for generating and providing the pUC-LAT52-mVenus expression vector for tobacco pollen transformation. We are very grateful to Sabine Freitag for help with lipid extraction and GC analyses and to Andrea Nickel as well as Janett Wochnik for support in genotyping and measuring hypocotyls. We would also like to thank Dr. Ellen Hornung and Dr. Amelie Kelly for very helpful discussions.
This research was financially supported by the European Commission’s Seventh Framework Program for Research and Technology Development (FP7), project GIAVAP, grant no. 266401 as well as the Göttingen Graduate School of Neurosciences and Molecular Biology (GGNB). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The Lobosphaera incisa genome sequence and annotation can be retrieved and analyzed at https://giavap-genomes.ibpc.fr. Please contact firstname.lastname@example.org for authorizations.
GenBank accession numbers of coding sequences used in this study:
Illumina RNA-Seq reads can be downloaded from SRA (Bioproject PRJNA262782).
HS cultivated algae, isolated LBs, extracted and digested proteins, analyzed proteomics data, performed in silico analyses, generated cDNA and cloned all coding sequences except for LiSDP1, generated and analyzed transgenic A. thaliana lines and N. tabacum pollen, interpreted the data and drafted the manuscript. OVale performed LC-MS/MS measurements and gave advice on data interpretation. TI gave advice on experimental setup as well as data interpretation for pollen tube assays and performed confocal microscopy of LiSDP1-expressing pollen tubes. JP amplified the LiSDP1 coding sequence from cDNA and performed GC analyses. NJT and OVall designed, sequenced and assembled the L. incisa genome and transcriptome. IK-G designed the genome and transcriptome experiment, co-supervised the lipid body isolation experiments, provided algal samples for both sequencing experiments as well as advice on L. incisa cultivation. GHB participated in study design and revised the manuscript. IF designed the study, analyzed the data and reviewed the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
LC-MS/MS dataset used for LBP identification. Specifications of LC-MS/MS measurements and data analysis are given and the data obtained are presented along with each filtering step in the identification of LBP candidates. (XLSX 7844 kb)
Confirmation of gene expression in transgenic A. thaliana lines. RNA was isolated from dry seeds and gene expression was analyzed by Reverse Transcriptase-PCR (RT-PCR) using primers detailed in Additional file 8. A, A. thaliana oleo1 expressing LiMLDP. B, A. thaliana sdp1/sdp1-L expressing LiLBP36. C, A. thaliana sdp1/sdp1-L expressing LiSDP1. (PDF 1766 kb)
Structural homology of LBP36 with yeast and fungal lipases. Phyre2 results are shown for the prediction of the LiLBP36 secondary structure and its alignment with known crystal structures. Additionally, an alignment of the predicted LiLBP36 structure (modeled on Malassezia globosa LIPASE1 using Phyre2) and Malassezia globosa LIPASE1 is shown. Only the structurally related section of LiLBP36 was aligned and the GXSXG motif as well as the catalytic aspartate are highlighted in green, with the two active residues being depicted as sticks. The 3D depiction was generated using PyMOL (The PyMOL Molecular Graphics System, Version 1.8 Schrödinger, LLC.). (PDF 1255 kb)
Quantification of total fatty acids in complemented A. thaliana lines. GC measurements of seeds and seedlings for Fig. 8c are shown. (XLSX 127 kb)
Primers used in this study. Primers used for RT-PCR and addition of restriction sites, qRT-PCR and confirmation of gene expression in transgenic lines are listed. Restriction sites are indicated capital letters. All primers were obtained from Sigma-Aldrich Chemie GmbH, Steinheim, Germany. CDS = Coding sequence. (DOCX 17 kb)
pCambia 43.0 vector information. A map of vector features is shown in addition to the entire nucleotide sequence. (DOCX 88 kb)
About this article
Cite this article
Siegler, H., Valerius, O., Ischebeck, T. et al. Analysis of the lipid body proteome of the oleaginous alga Lobosphaera incisa . BMC Plant Biol 17, 98 (2017). https://doi.org/10.1186/s12870-017-1042-2
- Lipid bodies
- Lipid droplets
- Lobosphaera incisa
- Nitrogen starvation
- Oil bodies
- TAG lipase