Next Article in Journal
Genomic Features of Open Chromatin Regions (OCRs) in Wild Soybean and Their Effects on Gene Expressions
Next Article in Special Issue
α-Fodrin in Cytoskeletal Organization and the Activity of Certain Key Microtubule Kinesins
Previous Article in Journal
Deletion of the Imprinted Phlda2 Gene Increases Placental Passive Permeability in the Mouse
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 5
10.3390/genes12050638
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Mitochondrial Processes during Early Development of Dictyostelium discoideum: From Bioenergetic to Proteomic Studies

by
Monika Mazur
1,
Daria Wojciechowska
1,2,
Ewa Sitkiewicz
3,
Agata Malinowska
3,
Bianka Świderska
3,
Hanna Kmita
1 and
Małgorzata Wojtkowska
1,*
1
Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, 61-614 Poznan, Poland
2
Department of Macromolecular Physics, Faculty of Physics, Adam Mickiewicz University, 61-614 Poznan, Poland
3
Institute of Biochemistry and Biophysics Polish Academy of Sciences, 02-106 Warszawa, Poland
*
Author to whom correspondence should be addressed.
Genes 2021, 12(5), 638; https://doi.org/10.3390/genes12050638
Submission received: 1 April 2021 / Revised: 18 April 2021 / Accepted: 20 April 2021 / Published: 25 April 2021
(This article belongs to the Special Issue Proteins: Proteomics and Beyond)
Browse Figures
Review Reports Versions Notes

Abstract

:
The slime mold Dictyostelium discoideum’s life cycle includes different unicellular and multicellular stages that provide a convenient model for research concerning intracellular and intercellular mechanisms influencing mitochondria’s structure and function. We aim to determine the differences between the mitochondria isolated from the slime mold regarding its early developmental stages induced by starvation, namely the unicellular (U), aggregation (A) and streams (S) stages, at the bioenergetic and proteome levels. We measured the oxygen consumption of intact cells using the Clarke electrode and observed a distinct decrease in mitochondrial coupling capacity for stage S cells and a decrease in mitochondrial coupling efficiency for stage A and S cells. We also found changes in spare respiratory capacity. We performed a wide comparative proteomic study. During the transition from the unicellular stage to the multicellular stage, important proteomic differences occurred in stages A and S relating to the proteins of the main mitochondrial functional groups, showing characteristic tendencies that could be associated with their ongoing adaptation to starvation following cell reprogramming during the switch to gluconeogenesis. We suggest that the main mitochondrial processes are downregulated during the early developmental stages, although this needs to be verified by extending analogous studies to the next slime mold life cycle stages.

1. Introduction

The unicellular slime mold Dictyostelium discoideum, one of the 150 species of Dictyostelia belonging to the Amoebozoa eukaryotic supergroup, displays the multicellularity form, namely aggregative or sorocarp multicellularity [1,2]. Due to starvation, the slime mold aggregates, which results in an intermediate migrating ‘slug’ stage and cells specializing into spores as well as three somatic cell types forming a stalk and structures to support the stalk and spore mass (Figure 1). This process is mainly based on chemical signals settled on cAMP and cAMP-dependent kinase (PkaC) [3], which induces the expression of the genes responsible for aggregation such as the receptor of cAMP (CarA), adenylate cyclase (AcaA), and cAMP phosphodiesterase (PdsA) [2,4]. cAMP pulses are initially secreted by a few starving cells, then cells move chemotactically toward the local cAMP and collect into molds, following which they form a slug and, finally, the fruiting body [5]. In addition, cAMP pulses upregulate the genes that are required during and after aggregation as tgrB and tgrC. These genes are responsible for cell adhesion and are known as the marker genes of post-aggregative cell differentiation [6].
D. discoideum cells, during the early steps of aggregation when unicellular cells start to move and form waves, are referred to as early developmental stages induced by starvation (Figure 1) [1]. These early developmental stages have been studied previously [7,8]. Accordingly, transcriptomic and cell proteomic data on mutant strains lacking genes that are crucial for cell aggregation are available [7]. The data revealed the downregulation of mitochondria proteins of the Krebs cycle and oxidative phosphorylation (OXPHOS). Another application of comparative MitoProteome studies for mitochondria isolated from different developmental stages including the early ones, is the general downregulation of crucial mitochondrial proteins and the upregulation of alternative oxidation (AOX) [8] during the early stages. The results of these studies were related to the mitochondrial energy status that was detected in vivo with a MitoTracker. However, these studies were not functionally based, which could have provided more precise answers to the question regarding mitochondrial functionality during starvation-induced aggregation.
Undoubtedly, D. discoideum is a popular model organism for studying cellular processes such as cell signaling, cell differentiation, and morphogenesis [9], the evolution of sociality and prey–predator relationships [3]. Due to the wide range of human orthologue genes, slime mold can be used to study human diseases, particularly neurodegeneration [10]. Moreover, D. discoideum is also a perfect model for research concerning intracellular and intercellular mechanisms such as early and late aggregation from unicellularity to multicellularity, in which mitochondria seem to have a crucial function [4].
It has been shown that cell differentiation is under the control of the nutrient glucose [11,12,13]. In our study, we performed an estimation of mitochondria’s functional state in intact cells representing the unicellular (U), aggregation (A), and streams (S) stages and combined the data with comparative proteomic analysis. The obtained results indicate that the observed differences among the studied stages may be associated with ongoing adaptation to starvation following the cell reprogramming that is supported by the metabolic background of starving cells. To our knowledge, this is the first report presenting combined D. discoideum mitochondrial proteomics and bioenergetic studies.

2. Methods

2.1. Cell Culture and Mitochondria Isolation

D. discoideum AX2 cell cultures were performed according to [14]. For the measurement of oxygen uptake using the Clarke electrode (Section 2.3), cells were cultured exponentially until they reached a density of 2.0–3.0 × 106 cells/mL (the doubling time of 10–11 h) in an axenic growth medium (1% w/v protease peptone, 0.5% w/v yeast extract, 0.45% w/v D-glucose, 0.09% w/v Na2HPO4, and 0.23% w/v KH2PO4). Cells were then cultured in 100 mL at 19 °C using a Certomat MO rotary shaker at 90 rpm. The cells of stage U were transferred to sterile 24-well plates at 1 million cells per well and left for 1 h for cell adhesion. Cells at stages A and S were obtained by washing the cell culture twice in the developmental buffer (DB) starvation medium containing 5 mM of Na2HPO4, 5 mM of KH2PO4, 1 mM of CaCl2, and 2 mM of MgCl2, with a pH of 6.5 [15]. The cells were then transferred to the sterile 24-well plates for 23 h for the A stage and 25 h for the S stage cells. The presence of the U, A, and S stages was confirmed by using a Nikon Eclipse TE 2000-U microscope. For mitochondria isolation, cell cultures (850 mL)at a given stage were centrifuged at 600× g for 5 min and then washed twice in phosphate medium A (14.5 mM KH2PO4, 5 mM Na2HPO4, pH 6) by centrifugation at 600× g for 5 min. Next, the cells were homogenized in medium B (0.38 M sucrose, 20 mM Tris-Cl, 0.5 mM EDTA, 1% defatted bovine serum albumin [BSA], pH 7.5) using a glass/teflon homogenizer and centrifuged at 860× g for 5 min. After this step, the supernatant was centrifuged at 10,000× g for 15 min and the obtained pellet was suspended in medium C (0.38 M sucrose, 20 mM TrisCl, 0.5 mM EDTA, 0.6% defatted BSA, pH 7.2), and centrifuged at 860× g for 5 min to remove debris. The supernatant was centrifuged at 10,000× g for 10 min, the mitochondrial pellet was suspended in medium D (0.38 M sucrose; 20 mM, Tris-Cl, pH 7.2), and centrifuged at 10,000× g for 15 min. The mitochondrial pellet was suspended in medium D and used for further analysis. All centrifugation procedures were carried out at 4 °C and media used for mitochondria isolation were ice cold. The calculated mean values of the yield of mitochondria for 1 g of cells were 5 mg, 1.4 mg, and 1.7 mg for stage U, A, and S, respectively.

2.2. The Measurements of Cell Respiration and the Status of Mitochondrial Coupling in Intact D. discoideum Cells

To measure the oxygen consumption rates of D. discoideum, cells from the U, A, and S stages collected from the 24-well plates were centrifuged (11,000× g, 3 min) and resuspended in the DB containing 5 mM of Na2HPO4, 5 mM of KH2PO4, 1 mM of CaCl2, and 2 mM of MgCl2 at a pH of 6.5 [15]. The rate of the cells’ oxygen uptake (2 × 106) was measured in 0.5 mL of the DB in a water-thermostated incubation chamber with a computer-controlled Clark-type O2 electrode at 18 °C (Oxygraph, Hansatech, UK). To estimate the functional states of the mitochondria at the studied stages, mitochondrial coupling capacity and mitochondrial coupling efficiency were calculated based on basal (actual) respiration, state 4 (the resting state), state 3 (the phosphorylating state), and maximal respiration. Basal cell respiration was measured in the DB with supplementation of pyruvate at a final concentration of 10 mM. State 4 was enforced through the addition of 0.5–1.5 µM of tributyltin (TBT), which is an inhibitor of mitochondrial ATP synthase, whereas maximal respiration was induced by using 0.1–0.24 µM of carbonyl cyanide-p-trifluoromethoxyphenylhydrazone, an uncoupler (FCCP). State 3 was calculated by subtracting state 4 from basal respiration. The contribution of state 3 to the basal respiration corresponded to mitochondrial coupling efficiency, while the mitochondrial coupling capacity that corresponded to FCCP uncoupling capacity was estimated as the maximal respiration to state 4 ratio. The TBT and FCCP concentrations were titrated for each stage [16,17,18]. To estimate the ability of the substrate supply and electron transport to respond to an increase in energy demand, mitochondrial spare respiratory capacity was estimated through the determination of the difference between maximal respiration and basal respiration [19].

2.3. Sample Preparation for Mass Spectrometry Analysis

Mass spectrometry experiments were performed at the Mass Spectrometry Laboratory at the Institute of Biochemistry and Biophysics, PAS. A total of 5 µg of isolated mitochondria was suspended in a dissolution buffer (8 M of urea, 2% CHAPS, 10 mM of HEPES, with a pH of 8.5). Samples were processed using the modified single-pot solid-phase-enhanced sample preparation (SP3) method [20]. The exact sample preparation protocol was as follows. Cysteines were reduced through 1 h of incubation with 20 mM of tris (2-carboxyethyl)phosphine (TCEP) at 37 °C, followed by 10 min of incubation at room temperature with 50 mM of methyl methanethiosulfonate (MMTS). The magnetic bead mix was prepared by combining equal parts of Sera-Mag Carboxyl hydrophilic and hydrophobic particles (09-981-121 and 09-981-123, GE Healthcare Lifesciences, Upsala, Sweden). The bead mix was washed three times with MS-grade water and resuspended in a working concentration of 10 µg/µL. Furthermore, 8 µL of the prepared bead mix, along with 5 µL of 10% formic acid and 500 µL of acetonitrile, was added to each sample. The proteins that bound to the beads were washed with 75% ethanol, isopropanol, and acetonitrile, followed by overnight digestion with 1 µg of trypsin/Lys-C mix (Promega). After digestion, the peptides were washed with acetonitrile and eluted from the beads by subsequent incubation with MS-grade water and 2% DMSO, with sonication during each step. The pulled aliquots were dried in a SpeedVac and resuspended in 60 µL of 2% acetonitrile and 0.1% formic acid. Peptide concentrations were measured using the Pierce Quantitative Colorimetric Peptide assay (Thermo Scientific, Rockford, IL, USA).

2.4. Mass Spectrometry

Three µg of digested samples were analyzed using the LC–MS system comprising a high-performance liquid chromatography (UPLC) chromatograph (nanoAcquity, Waters) and a Q Exactive mass spectrometer (Thermo). Peptides were trapped on a C18 pre-column (180 µm × 20 mm, Waters) using a 0.1% water solution of FA as a mobile phase then transferred to a BEH C18 column (75 µm × 250 mm, 1.7 µm, Waters) using an ACN gradient (0–35% ACN in 160 min) in the presence of 0.1% FA at a flow rate of 250 nL/min. The spectrometer was working in a data-dependent mode. To prevent cross-contamination, blank runs were performed between the sample runs.

2.5. Data Analysis

The obtained data were pre-processed with Mascot Distiller software (Matrixscience), and protein identification was performed using Mascot Server 2.5 (Matrixscience) against the D. discoideum protein sequences (13921 sequences) deposited in the NCBInr database (20180903, 167148673 sequences; 60963227986 residues). The parameters were set as follows: enzyme—trypsin, missed cleavages—1, fixed modifications—methylthio (C), variable modifications—GlyGly (K), oxidation (M), and instrument—HCD. To reduce mass errors, peptide and fragment mass tolerance settings were established separately for each file after an off-line mass recalibration [21]. The assessment of confidence was based on the target/decoy database search strategy [22], which provided q-value estimates for each peptide spectrum match. All queries with a q-value of >0.01, subset proteins, and proteins identified with one peptide were discarded from further analysis. The mass recalibration, false discovery rate (FDR) computations, and data filtering were done with Mscan software that was developed in-house (MScan 2.0.3 (accessed on 26 June 2018)).
The lists of identified peptides were merged and overlaid onto 2-D heatmaps generated from the LC-MS spectra, and the volumes were obtained from the assigned peaks (a more detailed description of data extraction procedures can be found in [23]). Quantitative values were then exported into text files for statistical analysis with Diffprot [21] software. The calculated p-values were adjusted for multiple testing using a procedure that controlled for a FDR. Diffprot was run with the following parameters: the number of random peptide sets = 106, clustering of peptide sets—only when 90% were identical, normalization through LOWESS, and quantification based on unique peptides. Only proteins with a q-value below 0.05 or those present in only one of the two compared analytical groups were taken into consideration during further analysis. The tendency in comparative studies as down in A and up in S was calculated from the ratio values of comparative qualitative studies for phase to phase ratio quotient U/A and A/S (Table 1). The expression profile was calculated against the U state as a basal of expression level (1). The ratio value (>1) and (<1) was regarded as upregulated and downregulated, respectively. The phase to phase ratio quotient “[S/U]” was taken for the proteins isocitrate dehydrogenase (NAD+), cytochrome c oxidase subunit IV, and heat shock protein Hsp70 family protein. For other proteins, the [S/U] was calculated from the equation [A/U]/[A/S], where [A/U] and [A/S] are given as the ratio in Section for “Phase to phase ratio quotient [A/U]” and Section for “Phase to phase ratio quotient [A/S]”.

3. Results and Discussion

3.1. The Functional Status of Intact D. discoideum Cells’ Mitochondria at the Early Developmental Stages

The respiratory state of mitochondria is a dynamic state resulting from crosstalk between the two bioenergetic states termed state 3 and state 4 [16,19]. State 4, also referred to as the proton-selective leak, is characterized by low oxygen uptake and higher inner membrane potential, while state 3, also referred to as the phosphorylating state due to resulting ADP phosphorylation to ATP, reflects high oxygen uptake and lower inner membrane potential. Maximal respiration is triggered in the presence of an uncoupler and represents the respiratory chain’s ability to support additional energetic demands [16,19]. All of the states are crucial for the estimation and comparison of the functional status of different intact cells’ mitochondria, as they allow for the calculation of the following parameters: mitochondrial coupling efficiency, mitochondrial coupling capacity, and mitochondrial spare respiratory capacity, as described in the Methods (Section 2.2). The parameter calculation allows for the comparison of cells containing different amounts of mitochondria, which was observed herein (Methods, Section 2.5) and has already been reported for cell in stages U, A, and S [8]. The representative traces registered for cells in stages U, A, and S are presented in the Figure 2A whereas raw data used for the calculation of all bioenergetic parameters are presented in Supplementary file S1. As shown in Figure 2B, intact D. discoideum cells in stages U, A, and S displayed the same rate of state 4 respiration (i.e., the proton leak), but differed in their basal, state 3, and maximal respiration rates. The values of the parameters were higher for cells in the U stage than for cells in the A and S stages, while the cells in stages A and S did not differ in the values of their parameters, except for maximal respiration, which was lower for cells in the S stage. However, according to the yield of mitochondria isolation, it could be assumed that cells in stages A and S contained the comparable mitochondria mass that was lower than in cells in the U stage.
Basal respiration is usually strongly controlled by ATP turnover; hence, it alters in response to ATP demand [19]. The observed differences in basal respiration rates co-occurred with a decrease in the rate of state 3, suggesting that the cells in stages A and S displayed a lower demand for ATP (Figure 2B). Accordingly, mitochondrial coupling efficiency corresponding to state 3 respiration contribution to basal respiration decreased to the same level for the cells in stages A and S when compared to the cells in stage U (Figure 2C. Thus, the cells in stages A and S displayed lower levels of ATP synthesis than those in stage U. However, the cells in stages A and S differed in the uncoupling capacity of FCCP (Figure 2B), which resulted in differences in mitochondrial coupling capacity (Figure 2D), defined as the ratio of maximal respiration to state 4 respiration. The value of the parameter is sensitive to changes in substrate oxidation and proton leakage [19]. Accordingly, FCCP-imposed uncoupling should reflect the maximum activity of electron transport and substrate oxidation, which are achievable by cells under the given conditions. Therefore, we compared the mitochondrial spare respiratory capacity of the studied cells. The parameter corresponds to the difference between maximal respiration and basal respiration and reflects how close to its bioenergetic limit a cell is operating [19]. One of the known factors that influence the extent of mitochondrial spare respiratory capacity is substrate availability, which is based on substrate entry into the Krebs cycle that is synchronized with the electron transport of the respiratory chain. As shown in Figure 2E, the cells in stage S displayed a distinct decrease in mitochondrial spare respiratory capacity when compared to the cells in stages U and A, but the value of the parameter increased for the cells in stage A when compared to the cells in stage U. Interestingly, it has been shown that mitochondrial spare respiratory capacity increases in starvation-resistant cells [24], which may permit rapid adaptation to metabolic changes [25]. On the other hand, cell reprogramming may correlate with a decrease in mitochondrial spare respiratory capacity [26], which could explain the distinct decrease in the parameter value of the cells in stage S [27]. Overall, the analysis of bioenergetic parameters indicated decreased ATP synthesis during stages A and S, which co-occurred in the case of the latter, with the limited access of substrates to the respiratory chain. These changes may be associated with an ongoing adaptation to starvation, followed by cell reprogramming. Therefore, we decided to study the changes in mitochondrial proteome differences among the cells in stages U, A, and S.

3.2. Quantitative Comparative Analysis of the Mitochondrial Proteomes of the Early Developmental Stages

In mitochondria isolated from the cells in stages U, A, and S, and due to the mass spectrometry analysis, we detected 294 proteins overall (Supplementary file S2 and segregated them into nine functional groups of proteins: protein synthesis and degradation, signaling, protein import, substrate transport, mitochondrial gene expression, the Krebs cycle, OXPHOS, other known processes (fatty acid metabolism, redox), and uncharacterized proteins (Figure 3).
The numbers of proteins assigned to the mentioned functional groups were comparable for the studied mitochondria. As presented in Figure 3, the most abundant groups were protein biosynthesis and degradation, OXPHOS, and other known processes. The less-abundant groups were transport, protein import, the Krebs cycle, and uncharacterized proteins. The least-numerous groups of proteins were the signaling and gene expression functional groups. Thus, the mitochondrial proteomes obtained for the studied early developmental stages were comparable. Therefore, we decided to perform a qualitative comparative analysis of the obtained mitochondrial proteome data.

3.3. Qualitative Comparative Analysis of the Mitochondrial Proteomes of the Early Developmental Stages

Comparative qualitative analyses of the detected proteins were performed for the following pairs of stages: A/U, S/U, and A/S regarding the previously described functional groups of proteins. Analyses were done according to the procedure described in the Methods section (Section 2.3, Section 2.4 and Section 2.5).

3.3.1. The Krebs Cycle and OXPHOS

The Krebs cycle and OXPHOS are well-known, crucial mitochondrial processes of energy transformation. The comparative analysis revealed the presence of four proteins of the Krebs cycle (citrate synthase (XP_643860.1), isocitrate dehydrogenase (NAD+) (XP_628920.1), succinate-CoA ligase (XP_636263.1), and malate dehydrogenase (XP_629516.1) and two proteins functionally linked to this cycle (3-oxoacid CoA-transferase and phosphoenolpyruvate carboxykinase) (XP_636911.1) (Table 1). We observed decreased levels in the majority of the mentioned proteins in stage A and increased levels in stage S, termed the ‘down in A and up in S tendency’ (Table 1), which often resulted in down- or upregulation in stage S when compared to stage U (Figure 4A,B).
As shown in Figure 4A, when the detected Krebs cycle proteins were compared between stages S and U, succinate-CoA ligase was upregulated, citrate synthase and isocitrate dehydrogenase (NAD+) were downregulated and malate dehydrogenase returned to stage U level, which denoted no change in its expression profile. The mitochondrial malate dehydrogenase catalyzed the reaction of the malate-to-oxaloacetate conversion with the use of NAD+ reduction to NADH [28,29]. It has been shown that enzyme activity is inhibited by differentiation-inducing-factor 1 (DIF-1), which affects cell energy transformation and leads to the inhibition of cell proliferation. These data suggest that the described inhibition of malate dehydrogenase by DIF-1 could be one of the mechanisms that induce the anti-proliferative effects that remain in agreement with the differentiation processes that start to occur at stage A. However, the steady state of its expression profile seemed not to influence the Krebs cycle.
The tendency of ‘down in A and up in S’ was also observed for the two enzymes functionally linked to the Krebs cycle. The cycle 3-oxoacid CoA-transferase was finally downregulated in stage S when compared to stage U (Figure 4B). Interestingly, the product of the human 3-oxoacid CoA-transferase (OXCT1) was converted to acetyl-CoA and finally fed into the Krebs cycle [30]. Mitochondrial phosphoenolpyruvate carboxykinase (PCK2) (XP_645396.1) [31] is a gluconeogenesis enzyme that catalyzes the phosphorylation reaction of oxaloacetate. PCK2 was the enzyme that was finally upregulated in stage S when compared to stage U (Figure 4A). We assumed that starving D. discoideum cells in stages A and S caused glucose deficiency that could trigger gluconeogenesis [32]. Interestingly, in our study, we observed an increase in phosphoenolpyruvate carboxykinase (PCK2) in stage S. The decrease in Krebs-cycle enzyme expression in stage A may coincide with enhanced gluconeogenesis, which is used to produce glucose during starvation. Unexpectedly, in our comparative analysis, we did not observe a significant change in proteins in the OXPHOS functional group, except for the cytochrome c-oxidase, which was downregulated (Figure 4C).
Our obtained data coincided with the low state 3 contributions to basal respiration in stages A and S, reflecting decreased OXPHOS, decreased spare respiratory capacity in stage S, and possibly enhanced gluconeogenesis. Nevertheless, the obtained data did not reveal an upregulation of any glycolytic enzymes that could have explained the aerobic glycolysis that is regarded as typical for proliferative cells, similar to cancer [33]. This is because the study was focused on mitochondrial proteins. However, glycolytic enzymes could interact with mitochondria, as glycolysis occurs in the cytosol.
The data indicating the downregulation of the Krebs cycle enzymes coincide with the results that were obtained in [7,8]. We concluded that the main function of the mitochondria in the early developmental stages emerged temporarily, rewriting the mitochondrial metabolism from the Krebs cycle into gluconeogenesis. We showed, for the first time, the presence of PCK2 upregulation during D. discoideum starvation. The regulation of PCK2 at the level of the transcription of the encoding gene is under the positive control of cAMP and is crucial for cell aggregation. Indeed, the deficiency of glucose during starvation activates adenylate cyclase, which actives CREB and, finally, PCK2. This could take place when phosphoenolpyruvate (PEP) cannot be converted into pyruvate, which is used in the gluconeogenesis pathway (Figure 4C) [31].

3.3.2. Protein Metabolism

In the case of the protein biosynthesis and degradation group, the obtained data revealed three enzymes that were responsible for the metabolism of mitochondrial amino acids (Table 1). The general tendency concerning their expression level remained ‘down in A and up in S’ and was observed for acetyl glutamate kinase (XP_637813.1), 4-aminobutyrate transaminase (XP_647552.1) (Figure 4A), and branched-chain amino acid aminotransferase (XP_638096.1). The latter did not return to the level of stage U, so the enzyme was classified as downregulated in stage S when compared to stage U (Figure 4B). Acetyl glutamate kinase is a NAD+ binding enzyme engaged in arginine biosynthesis [28]. In D. discoideum, the ArgC enzyme is similar to the bifunctional acetyl glutamate kinase and N-acetyl-γ-glutamyl-phosphate reductase from yeast, which catalyze the second and third steps in the biosynthesis of the arginine precursor ornithine. The detection of this protein revealed the known connection between the urea and Krebs cycles. However, we did not observe any changes in protein expression in the case of the urea cycle. The 4-aminobutyrate transaminase catalyzed the conversion of 4-aminobutanoate and 2-oxoglutarate into succinate semialdehyde and L-glutamate [33]. This process is dedicated to glutamate degradation and is also correlated with the urea and Krebs cycles. The expression profile of these proteins included downregulation in A and an increase in S, resulting in upregulation (Figure 4A). In the case of the branched-chain amino acid aminotransferase, another example of aminotransferase, the transamination of branch-amino acids as leucine, isoleucine, and valine as well as their respective alfa-ketoacids, was reported [28]. For this protein, an opposite decreasing tendency was observed, so we classified the protein as downregulated in stage S compared to stage U (Figure 4B). We speculate that the amino acids’ biosynthesis protein expression profile revealed the metabolic flexibility of D. discoideum cells determined by their ability to reprogram anabolic and catabolic pathways, presumably through altering gene expression programs and intercellular interactions within the glucose-free microenvironment during starvation [34]. This reprogramming could be based on the temporary shift from the Krebs cycle to gluconeogenesis, with amino acids as the fuel for gluconeogenesis (Figure 4C).
The comparative study on the discussed groups of proteins also revealed that the presence of cytosolic chaperone Hsp70 (XP_643155.1) significantly increased in stage A when compared to stage U and slightly decreased in stage S when compared to stage U (Figure 4A). Hsp70 is a chaperone that plays a crucial role in protein folding, disaggregation, and degradation [35]. It was also shown that ribosomes stay in close proximity to mitochondria [36], which implies their co-translational import into mitochondria [37]. Therefore, the presence of the mentioned chaperone could ensure sufficient protein refolding, which is essential for the successful import of protein into mitochondria.
To summarize, the detected changes in the expression levels of the functional protein group suggest that at the early developmental stages, the processes of mitochondrial amino acid metabolism are upregulated, which we showed here for the first time.

3.3.3. The Import of Proteins into Mitochondria and the Transport of Other Molecules

Most mitochondrial proteins are synthesized in the cytosol and have to be delivered to the various mitochondrial sub-locations. This process is regulated by the presence and activity of import apparatus consisting of proteins that, in general, recognize and translocate proteins across mitochondrial membranes (for details, see [37]). In this comparative study, we detected mtHsp70 (XP_629204.1) [38], which showed the ‘down in A and up in S’ expression profile tendency. We estimated that mtHsp70 was upregulated in stage S when compared to stage U (Figure 4A). This mitochondrial chaperone is responsible for the protein assembly of proteins delivered into the mitochondrial matrix and acts as a subunit of the PAM import machinery that is crucial for the proper unidirectional transfer of precursor proteins across the inner mitochondrial membrane. It seems that this process could be upregulated in stage S. The elevated level of cytosolic Hsp70 could be arranged in the same processes as those involved when the chaperone assists in the partial refolding of precursor proteins for import across the protein channel [8]. Moreover, mtHp70 is a chaperone that plays an important role in iron–sulfur cluster (IS) biogenesis [39]. Studies on yeast and human ISC involve gene mutations to show genome instability and the induction of the DNA damage repair pathway [40]. This, in turn, indicates that pathways related to ISC biosynthesis are strongly conserved—more so than for ATP production, for example—among diverse mitochondrial homologues [41]. This also explains ISC’s presence (as only one common module to mitochondria) in mitosomes. This suggests that the essential process of ISC biosynthesis occurs in stage S.
In the case of the functional transport protein group, we detected an uncharacterized mitochondria substrate carrier protein named mt carrier Q (XP_645160.1). Its level significantly increased in stage A when compared to stage U and decreased in stage S when compared to stage U, although its level was higher when in stage U (Figure 4A). We conclude that the organization of the mitochondrial import machineries do not undergo major changes during the early developmental stages. We identified a new, uncharacterized transport protein, which was significantly upregulated in stage S (Figure 4A) and is worth further study to discern its exact function in the process of aggregation.

3.3.4. Other Known Processes

The functional group, termed ‘other known processes’, included enzymes engaged in fatty acid metabolism and redox processes. The oxidation of fatty acids may provide another source for biosynthetic growth under glucose limitation during D. discoideum starvation [42]. Our comparative study provided four enzymes of the fatty acid metabolism acetyl-CoA C-acetyltransferase (XP_645587.1), which showed the ‘down in A and up in S’ tendency, but its level in stage S was lower than in stage U, denoting a final downregulation (Table 1 and Figure 4B). The ‘down in A and up in S’ tendency was also observed for the electron transfer of flavoprotein α and β subunits (XP_635485.1 and XP_642058.1). However, the comparison of the proteins’ levels between stages U and S indicate their upregulation (Figure 4A). Interestingly, acyl-CoA oxidase A (XP_643323.1) engaged in β-oxidation that was elevated in stage A when compared to stage U and strongly decreased in stage S when compared to stage U (Figure 4B). This profile was distinctly different from that observed for the previous groups of proteins. Therefore, we can assume the diminution of β-oxidation.
During β-oxidation, acetyl-CoA is produced, then enters the Krebs cycle and consequently stimulates ATP synthesis. However, during starvation, acetyl-CoA could activate the first enzyme of gluconeogenesis—the pyruvate carboxykinase that converts pyruvate into oxaloacetate. As shown, the detected β-oxidation proteins were mostly downregulated in stage S when compared to stage U. Therefore, D. discoideum cells probably do not use acetyl-CoA during starvation to stimulate gluconeogenesis (the first enzyme of this pathway). Consequently, glycerol is not used as a carbon source for D. discoideum’s gluconeogenesis. Thus, as described, during starvation, gluconeogenesis may start due to PCK2, which has been proven to upregulate in stage S when compared to stage U and uses amino acids as fuel. Thus, amino acids are metabolized into the oxaloacetate, then converted into PEP by mitochondrial PCK2 (Figure 4C) [43].
Within the described group of other known processes, proteins engaged in the redox reaction were also enclosed. For this group, we detected mitochondrial superoxidase dismutase (SOD2) (XP_645815.1). SOD2 followed the ‘down in A and up in S’ pattern. It was finally upregulated in stage S when compared to stage U, which coincides with the available data [8]. This upregulation was not very significant (Figure 4A) but remained in agreement with the data concerning the role of ROSs as signaling molecules during aggregation, where SOD control signaling molecules are essential for proper development [44]. The ROS production could increase as a response of mitochondrial metabolism reprogramming. Thus, additional studies in this field are needed.

3.3.5. Uncharacterized Proteins

Two uncharacterized mitochondrial proteins were detected in the comparative analysis (Table 1). The putative mitochondrial transferase caf17 (XP_639996.1) and an uncharacterized protein (XP_639145.1), named in the D.discoideum database (dictybase) ‘as putative delta-1-pyrroline-5-carboxylate dehydrogenase’, were discovered. Both followed the ‘down in A and up in S’ tendency. The putative mitochondrial transferase caf17 was presumably involved in the incorporation of iron–sulfur clusters (ISC) into mitochondrial aconitase-type proteins. ISC was previously mentioned as being a common but crucial molecule present in mitochondria and related to organelles [45], which is thought to reflect its basic mitochondrial function. The ‘putative delta-1-pyrroline-5-carboxylate dehydrogenase’ is probably involved in synthesis of L-glutamate from L-proline, suggesting its function in the amino acid degradation process (Figure 4C). Both uncharacterized proteins were upregulated in stage S when compared to stage U (Figure 4A) and warrant a deeper investigation.

4. Conclusions

Our study on D. discoideum cells in early developmental stages revealed differences in their mitochondria functional state. The differences correspond to the comparative proteomic analysis results and may be associated with diminution of β oxidation, Krebs cycle, and oxidative phosphorylation, and enhancement of the metabolism of amino acids, being a known fuel for gluconeogenesis during starvation. This in turn may support metabolically ongoing adaptation to starvation and following cell reprogramming.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12050638/s1, Supplementary file S1 show raw data obtained from oxygraphic measurements and Supplementary file S2 show proteins and peptides identified in LC-MS analysis.

Author Contributions

Conceptualization, M.W. and H.K.; Methodology, M.W., E.S. and M.M.; Software, E.S.; Validation, E.S., D.W., M.M. and A.M.; Formal analysis, M.W., D.W., E.S. and A.M.; Investigation, M.M., E.S., D.W. and B.Ś.; Resources, M.M.; Data curation, M.W., D.W. and M.M.; Writing—original draft preparation M.W. and H.K.; Writing—review and editing, M.W. and H.K.; Visualization, M.M., D.W., A.M. and E.S.; Supervision, H.K.; Project administration, M.W.; Funding acquisition, M.W., M.M. and H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the National Science Center (Poland) project (grant no. 2012/05/N/NZ3/00293) and the “KNOW RNA Research Center in Poznan” (grant no. 01/KNOW2/2014). Special Support was provided by the Dean of Faculty of Biology UAM Poznań, P. Wojtaszek.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

LC-MS data are available on request at the Mass Spectrometry Laboratory IBB PAS.

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schilde, C.; Schaap, P. The Amoebozoa. Methods Mol. Biol. 2013, 983, 1–15. [Google Scholar] [PubMed] [Green Version]
  2. Kawabe, Y.; Du, Q.; Schilde, C.; Schaap, P. Evolution of multicellularity in Dictyostelia. Int. J. Dev. Biol. 2019, 63, 359–369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Du, Q.; Kawabe, Y.; Schilde, C.; Zhi-Hui, C.; Schaap, P. The evolution of aggregative multicellularity and cell-cell commu-nication in the Dictyostelia. J. Mol. Biol. 2015, 427, 3722–3733. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Maeda, Y. Cell-cycle checkpoint for transition from cell division to differentiation. Dev. Growth Differ. 2011, 53, 463–481. [Google Scholar] [CrossRef]
  5. Dormann, D.; Kim, J.Y.; Devreotes, P.N.; Weijer, C.J. cAMP receptor affinity controls wave dynamics, geometry and morphogenesis in Dictyostelium. J. Cell Sci. 2001, 114, 2513–2523. [Google Scholar]
  6. Hirose, S.; Chen, G.; Kuspa, A.; Shaulsky, G. The polymorphic proteins TgrB1 and TgrC1 function as a ligand-receptor pair in Dictyostelium allorecognition. J. Cell Sci. 2017, 130, 4002–4012. [Google Scholar] [CrossRef] [Green Version]
  7. González-Velasco, Ó.; De Las Rivas, J.; Lacal, J. Proteomic and transcriptomic profiling identifies early developmentally reg-ulated proteins in Dictyostelium discoideum. Cells 2019, 8, 1187. [Google Scholar] [CrossRef] [Green Version]
  8. Czarna, M.; Mathy, G.; Mac’Cord, A.; Dobson, R.; Jarmuszkiewicz, W.; Sluse-Goffart, C.M.; Leprince, P.; De Pauw, E.; Sluse, F.E. Dynamics of the Dictyostelium discoideum mitochondrial proteome during vegetative growth, starvation and early stages of development. Proteomics 2009, 10, 6–22. [Google Scholar] [CrossRef]
  9. Mathavarajah, S.; Flores, A.; Huber, R.J. Dictyostelium discoideum: A Model System for Cell and Developmental Biology. Curr. Protoc. Essent. Lab. Tech. 2017, 15, 14.1.1–14.1.19. [Google Scholar] [CrossRef]
  10. Pearce, X.G.; Annesley, S.J.; Fisher, P.R. The Dictyostelium model for mitochondrial biology and disease. Int. J. Dev. Biol. 2019, 63, 497–508. [Google Scholar] [CrossRef] [Green Version]
  11. Leach, C.K.; Ashworth, J.M.; Garrod, D.R. Cell sorting out during the differentiation of mixtures of metabolically distinct populations of Dictyostelium discoideum. J. Embryol. Exp. Morphol. 1973, 29, 647–661. [Google Scholar] [PubMed]
  12. Tasaka, M.; Takeuchi, I. Role of Cell Sorting in Pattern Formation in Dictyostelium discoideum. Differentiation 1981, 18, 191–196. [Google Scholar] [CrossRef]
  13. Thompson, C.R.; Kay, R.R. Cell-Fate Choice in Dictyostelium: Intrinsic Biases Modulate Sensitivity to DIF Signaling. Dev. Biol. 2000, 227, 56–64. [Google Scholar] [CrossRef]
  14. Mazur, M.; Wojtkowska, M.; Skalski, M.; Słocińska, M.; Kmita, H. The The TOB/SAM complex composition in mitochondria of Dictyostelium discoideum during progression from unicellularity to multicellularity. Acta Biochim. Pol. 2019, 66, 551–557. [Google Scholar] [CrossRef]
  15. Fey, P.; Kowal, A.S.; Gaudet, P.; Pilcher, K.E.; Chisholm, R.L. Protocols for growth and development of Dictyostelium dis-coideum. Nat. Protoc. 2007, 2, 1307–1316. [Google Scholar] [CrossRef]
  16. Michejda, J.; Guo, X.J.; Lauquin, G.J.-M. The respiration of cells and mitochondria of porin deficient yeast mutants is coupled. Biochem. Biophys. Res. Commun. 1990, 171, 354–361. [Google Scholar] [CrossRef]
  17. Karachitos, A.; Jordán, J.; Kmita, H. Cytoprotective activity of minocycline includes improvement of mitochondrial coupling: The importance of minocycline concentration and the presence of VDAC. J. Bioenerg. Biomembr. 2012, 44, 297–307. [Google Scholar] [CrossRef] [Green Version]
  18. Karachitos, A.; Grobys, D.; Antoniewicz, M.; Jedut, S.; Jordan, J.; Kmita, H. Human VDAC isoforms differ in their capability to interact with minocycline and to contribute to its cytoprotective activity. Mitochondrion 2016, 28, 38–48. [Google Scholar] [CrossRef]
  19. Brand, M.D.; Nicholls, D.G. Assessing mitochondrial dysfunction in cells. Biochem. J. 2011, 435, 297–312. [Google Scholar] [CrossRef] [Green Version]
  20. Hughes, C.S.; Foehr, S.; Garfield, D.A.; Furlong, E.E.; Steinmetz, L.M.; Krijgsveld, J. Ultrasensitive proteome analysis using paramagnetic bead technology. Mol. Syst. Biol. 2014, 10, 757. [Google Scholar] [CrossRef] [PubMed]
  21. Malinowska, A.; Kistowski, M.; Bakun, M.; Rubel, T.; Tkaczyk, M.; Mierzejewska, J.; Dadlez, M. Diffprot—Software for non-parametric statistical analysis of differential proteomics data. J. Proteom. 2012, 75, 4062–4073. [Google Scholar] [CrossRef]
  22. Elias, J.; Haas, W.; Faherty, B.; Gygi, S.P. Comparative evaluation of mass spectrometry platforms used in large scale prote-omics investigations. Nat. Methods 2005, 2, 667–675. [Google Scholar] [CrossRef] [PubMed]
  23. Bakun, M.; Karczmarski, J.; Poznanski, J.; Rubel, T.; Rozga, M.; Malinowska, A.; Sands, D.; Hennig, E.; Oledzki, J.; Ostrowski, J.; et al. An integrated LC-ESI-MS platform for quantitation of serum peptide ladders. Application for colon carcinoma study. Proteom. Clin. Appl. 2009, 3, 932–946. [Google Scholar] [CrossRef] [PubMed]
  24. Isono, T.; Chano, T.; Yonese, J.; Yuasa, T. Therapeutic inhibition of mitochondrial function induces cell death in starva-tion-resistant renal cell carcinomas. Sci. Rep. 2016, 9, 25669. [Google Scholar] [CrossRef] [Green Version]
  25. Keuper, M.; Jastroch, M.; Yi, C.; Fischer-Posovszky, P.; Wabitsch, M.; Tschöp, M.H.; Hofmann, S.M. Spare mitochondrial respiratory capacity permits human adipocytes to maintain ATP homeostasis under hypoglycemic conditions. FASEB J. 2014, 28, 761–770. [Google Scholar] [CrossRef]
  26. Zhou, Y.; Al-Saaidi, R.A.; Fernandez-Guerra, P.; Freude, K.K.; Olsen, R.K.J.; Jensen, U.B.; Gregersen, N.; Hyttel, P.; Bolund, L.; Aagaard, L.; et al. Mitochondrial Spare Respiratory Capacity Is Negatively Correlated with Nuclear Reprogramming Efficiency. Stem Cells Dev. 2017, 26, 166–176. [Google Scholar] [CrossRef] [PubMed]
  27. Pfleger, J.M.; He, M.; Abdellatif, M.M.M. Mitochondrial complex II is a source of the reserve respiratory capacity that is regulated by metabolic sensors and promotes cell survival. Cell Death Dis. 2015, 6, e1835. [Google Scholar] [CrossRef] [Green Version]
  28. Journet, A.; Klein, G.; Brugière, S.; Vandenbrouck, Y.; Chapel, A.; Kieffer, S.; Bruley, C.; Masselon, C.; Aubry, L. Investigating the macropinocytic proteome of Dictyostelium amoebae by high-resolution mass spectrometry. Proteomics 2012, 12, 241–245. [Google Scholar] [CrossRef]
  29. Matsuda, T.; Takahashi-Yanaga, F.; Yoshihara, T.; Maenaka, K.; Watanabe, Y.; Miwa, Y.; Morimoto, S.; Kubohara, Y.; Hirata, M.; Sasaguri, T. Dictyostelium differentiation-inducing factor-1 binds to mitochondrial malate dehydrogenase and inhibits its activity. J. Pharmacol. Sci. 2010, 112, 320–326. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, S.; Xie, C. The role of OXCT1 in the pathogenesis of cancer as a rate-limiting enzyme of ketone body metabolism. Life Sci. 2017, 183, 110–115. [Google Scholar] [CrossRef]
  31. Li, C.-L.F.; Santhanam, B.; Webb, A.N.; Zupan, B.; Shaulsky, G. Gene discovery by chemical mutagenesis and whole-genome sequencing in Dictyostelium. Genome Res. 2016, 26, 1268–1276. [Google Scholar] [CrossRef] [Green Version]
  32. Vincent, E.E.; Sergushichev, A.; Griss, T.; Gingras, M.-C.; Samborska, B.; Ntimbane, T.; Coelho, P.P.; Blagih, J.; Raissi, T.C.; Choinière, L.; et al. Mitochondrial Phosphoenolpyruvate Carboxykinase Regulates Metabolic Adaptation and Enables Glucose-Independent Tumor Growth. Mol. Cell 2015, 60, 195–207. [Google Scholar] [CrossRef] [Green Version]
  33. Wu, Y.; Janetopoulos, C. Systematic Analysis of γ-Aminobutyric Acid (GABA) Metabolism and Function in the Social Amoeba Dictyostelium discoideum. J. Biol. Chem. 2013, 288, 15280–15290. [Google Scholar] [CrossRef] [Green Version]
  34. Pavlova, N.N.; Thompson, C.B. The Emerging Hallmarks of Cancer Metabolism. Cell Metab. 2016, 23, 27–47. [Google Scholar] [CrossRef] [Green Version]
  35. Fernández-Fernández, M.R.; Valpuesta, J.M. Hsp70 chaperone: A master player in protein homeostasis. F1000Research 2018, 7, 1497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Vardi-Oknin, D.; Arava, Y. Characterization of factors involved in localized translation near mitochondria by ribosome-proximity labeling. Front. Cell Dev. Biol. 2019, 7, 305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hansen, K.G.; Herrmann, J.M. Transport of Proteins into Mitochondria. Protein J. 2019, 38, 330–342. [Google Scholar] [CrossRef]
  38. Wojtkowska, M.; Buczek, D.; Suzuki, Y.; Shabardina, V.; Makałowski, W.; Kmita, H. The emerging picture of the mitochon-drial protein import complexes of Amoebozoa supergroup. BMC Genom. 2017, 18, 997. [Google Scholar]
  39. Shan, Y.; Cortopassi, G. Mitochondrial Hspa9/Mortalin regulates erythroid differentiation via iron-sulfur cluster assembly. Mitochondrion 2016, 26, 94–103. [Google Scholar] [CrossRef]
  40. Stehling, O.; Wilbrecht, C.; Lill, R. Mitochondrial iron–sulfur protein biogenesis and human disease. Biochimie 2014, 100, 61–77. [Google Scholar] [CrossRef]
  41. Embley, T.M.; Martin, W. Eukaryotic evolution, changes and challenges. Nat. Cell Biol. 2006, 440, 623–630. [Google Scholar] [CrossRef] [PubMed]
  42. Vacanti, N.M.; Divakaruni, A.S.; Green, C.R.; Parker, S.J.; Henry, R.R.; Ciaraldi, T.P.; Murphy, A.N.; Metallo, C.M. Regula-tion of substrate utilization by the mitochondrial pyruvate carrier. Mol. Cell 2014, 56, 425–435. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Satre, M.; Mattei, S.; Aubry, L.; Gaudet, P.; Pelosi, L.; Brandolin, G.; Klein, G. Mitochondrial carrier family: Repertoire and peculiarities of the cellular slime mould Dictyostelium discoideum. Biochimie 2007, 89, 1058–1069. [Google Scholar] [CrossRef] [PubMed]
  44. Bloomfield, G.; Pears, C. Superoxide signalling required for multicellular development of Dictyostelium. J. Cell Sci. 2003, 116, 3387–3397. [Google Scholar] [CrossRef] [Green Version]
  45. Roger, A.J.; Muñoz-Gómez, S.A.; Kamikawa, R. The Origin and Diversification of Mitochondria. Curr. Biol. 2017, 27, R1177–R1192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 12 00638 g001 550
Figure 1. The life cycle of D. discoideum. The studied stages were marked by bold letters.
Figure 1. The life cycle of D. discoideum. The studied stages were marked by bold letters.
Genes 12 00638 g001
Genes 12 00638 g002 550
Figure 2. The functional state of mitochondria in intact D. discoideum cells in the studied stages. (A) Representative traces for cells in studied stages applied for calculations as well as explanation of the calculated parameters; (B) oxygen consumption rates under the applied oxygen uptake rate measurement condition; (C) coupling efficiency (contribution of state 3 to basal respiration); (D) coupling capacity (maximal respiration to state 4 ratio); (E) spare respiratory capacity (the difference between maximal respiration and basal respiration). The rate of oxygen uptake by cells (2 × 106) was measured in 0.5 mL of DB in a water-thermostated incubation chamber with a computer-controlled Clark-type O2 electrode at 18 °C (Oxygraph, Hansatech, UK). Tributyltin (TBT) was added in a concentration (0.5–1.5 µM) to enforce state 4. The maximal respiration was induced by using 0.1–0.24 µM of carbonyl cyanidep-trifluoromethoxyphenylhydrazone (FCCP); sodium dithionite is a zero oxygen solution powder used for calibration of oxygen sensors at zero oxygen. The data are presented as mean values ± SD of three independent experiments. * p < 0.05; ** p < 0.005, *** p < 0.001. Calculated original data are presented in Supplementary file S1.
Figure 2. The functional state of mitochondria in intact D. discoideum cells in the studied stages. (A) Representative traces for cells in studied stages applied for calculations as well as explanation of the calculated parameters; (B) oxygen consumption rates under the applied oxygen uptake rate measurement condition; (C) coupling efficiency (contribution of state 3 to basal respiration); (D) coupling capacity (maximal respiration to state 4 ratio); (E) spare respiratory capacity (the difference between maximal respiration and basal respiration). The rate of oxygen uptake by cells (2 × 106) was measured in 0.5 mL of DB in a water-thermostated incubation chamber with a computer-controlled Clark-type O2 electrode at 18 °C (Oxygraph, Hansatech, UK). Tributyltin (TBT) was added in a concentration (0.5–1.5 µM) to enforce state 4. The maximal respiration was induced by using 0.1–0.24 µM of carbonyl cyanidep-trifluoromethoxyphenylhydrazone (FCCP); sodium dithionite is a zero oxygen solution powder used for calibration of oxygen sensors at zero oxygen. The data are presented as mean values ± SD of three independent experiments. * p < 0.05; ** p < 0.005, *** p < 0.001. Calculated original data are presented in Supplementary file S1.
Genes 12 00638 g002
Genes 12 00638 g003 550
Figure 3. The numbers of D. discoideum mitochondrial proteins assigned to the nine functional groups. The classification of identified proteins was performed at the base of the dictybase.org and proteins belonging to the same functional groups were counted in order to assess the size of the classes for details (see Methods Section 2.5 and Supplementary file S2).
Figure 3. The numbers of D. discoideum mitochondrial proteins assigned to the nine functional groups. The classification of identified proteins was performed at the base of the dictybase.org and proteins belonging to the same functional groups were counted in order to assess the size of the classes for details (see Methods Section 2.5 and Supplementary file S2).
Genes 12 00638 g003
Genes 12 00638 g004 550
Figure 4. The mitochondrial proteins upregulated and downregulated in stage S when compared to stage U during D. discoideum starvation. (A) Upregulated mitochondrial proteins; (B) downregulated mitochondrial proteins. The protein upregulation and downregulation were estimated using qualitative comparative analysis between stages U and S (see Methods Section 2.5). (C) Diagram of downregulated (blue) and upregulated (orange) mitochondrial proteins assigned to cellular processes. The cell metabolic adaptation under glucose depletion includes attenuation of the Krebs cycle, fatty acid β oxidation, and oxidative phosphorylation, and enhancement of the amino acids’ metabolism, being known as a fuel for gluconeogenesis. The early steps of gluconeogenesis is regulated by mitochondrial phosphoenolpyruvate carboxykinase (PCK2), which allows for alternative glucose/biosynthesis intermediates required for cell aggregation. TCA (tricarboxylic acid cycle, Krebs cycle); PEP (phosphoenolopyruvate); SOD (superoxide dismutase); caf17 (putative mitochondrial transferase); scsA (succinate-CoA ligase); mdhB (mitochondrial malate dehydrogenase); idhB (isocitrate dehydrogenase (NAD+)); α-KG (α-ketoglutarate); OM (mitochondrial outer membrane); IM (mitochondrial inner membrane).
Figure 4. The mitochondrial proteins upregulated and downregulated in stage S when compared to stage U during D. discoideum starvation. (A) Upregulated mitochondrial proteins; (B) downregulated mitochondrial proteins. The protein upregulation and downregulation were estimated using qualitative comparative analysis between stages U and S (see Methods Section 2.5). (C) Diagram of downregulated (blue) and upregulated (orange) mitochondrial proteins assigned to cellular processes. The cell metabolic adaptation under glucose depletion includes attenuation of the Krebs cycle, fatty acid β oxidation, and oxidative phosphorylation, and enhancement of the amino acids’ metabolism, being known as a fuel for gluconeogenesis. The early steps of gluconeogenesis is regulated by mitochondrial phosphoenolpyruvate carboxykinase (PCK2), which allows for alternative glucose/biosynthesis intermediates required for cell aggregation. TCA (tricarboxylic acid cycle, Krebs cycle); PEP (phosphoenolopyruvate); SOD (superoxide dismutase); caf17 (putative mitochondrial transferase); scsA (succinate-CoA ligase); mdhB (mitochondrial malate dehydrogenase); idhB (isocitrate dehydrogenase (NAD+)); α-KG (α-ketoglutarate); OM (mitochondrial outer membrane); IM (mitochondrial inner membrane).
Genes 12 00638 g004
Table
Table 1. Results of the comparative mitoproteomes of the U (unicellular), A (aggregation), and S (streams) stages. D. discoideum stages regarding the distinguished functional groups: A/U, S/U, and A/S.
Table 1. Results of the comparative mitoproteomes of the U (unicellular), A (aggregation), and S (streams) stages. D. discoideum stages regarding the distinguished functional groups: A/U, S/U, and A/S.
RefSeq
Database		Protein NamePhase to Phase Ratio Quotient
A/U		Phase to Phase Ratio Quotient
S/U		Phase to Phase Ratio Quotient
A/S
q ValueRatioFold ChangePeptidesq ValueRatioFold ChangePeptidesq ValueRatioFold ChangePeptides
Krebs
XP_643860.1citrate synthase, mitochondrial0.000730.561.7846 0.000230.631.645
XP_628920.1isocitrate dehydrogenase (NAD+) (IdhB)0.048420.661.52250.012530.611.6425
XP_636263.1succinate-CoA ligase (scsA)0.008580.482.116 0.005990.422.3716
XP_629516.1malate dehydrogenase (mdhB)0.000270.581.7145 0.000060.581.7345
-
XP_636911.13-oxoacid CoA-transferase0.005690.551.8229 0.040590.581.7329
XP_645396.1phosphoenolpyruvate carboxykinase (PCK2)0.040450.541.8630 0.000550.412.4332
OXPHOS
XP_640649.1cytochrome c oxidase subunit IV 0.032470.492.0619
Protein biosynthesis
XP_643155.1heat shock protein Hsp70 family protein0.000153.763.76320.000073.473.4733
XP_637813.1acetylglutamate kinase0.008860.611.6537 0.000610.541.8637
XP_638096.1branched-chain amino acid aminotransferase0.000150.452.2122 0.010410.571.7521
XP_647552.14-aminobutyrate transaminase0.005120.581.7343 0.003680.571.7443
Protein import
XP_629204.1Stress-70 protein, mitochondrial (mtHsp70)0.00040.61.6762 0.000060.531.8862
XP_645160.1mitochondrial substrate carrier family protein Q0.007393.113.1110 0.001912.252.2510
Other known processes
FA metabolism
XP_645587.1acetyl-CoA C-acetyltransferase0.000150.422.4137 0.005490.611.6535
XP_643323.1acyl-CoA oxidase A0.013252.162.1625 0.000065.545.5425
XP_635485.1electron transfer flavoprotein alpha subunit0.019320.511.9724 0.006150.422.3623
XP_642058.1electron transfer flavoprotein beta subunit0.035240.61.6521 0.0220.541.8621
XP_645815.1Redox
superoxide dismutase/SOD		0.01603 0.462.1810 0.004040.42.5212
XP_639996.1Uncharacterized proteins
putative mitochondrial transferase caf17		0.005420.412.4124 0.000060.283.5924
XP_639145.1putative delta-1-pyrroline-5-carboxylate dehydrogenase0.040230.631.5938 0.000060.492.0437
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mazur, M.; Wojciechowska, D.; Sitkiewicz, E.; Malinowska, A.; Świderska, B.; Kmita, H.; Wojtkowska, M. Mitochondrial Processes during Early Development of Dictyostelium discoideum: From Bioenergetic to Proteomic Studies. Genes 2021, 12, 638. https://doi.org/10.3390/genes12050638

AMA Style

Mazur M, Wojciechowska D, Sitkiewicz E, Malinowska A, Świderska B, Kmita H, Wojtkowska M. Mitochondrial Processes during Early Development of Dictyostelium discoideum: From Bioenergetic to Proteomic Studies. Genes. 2021; 12(5):638. https://doi.org/10.3390/genes12050638

Chicago/Turabian Style

Mazur, Monika, Daria Wojciechowska, Ewa Sitkiewicz, Agata Malinowska, Bianka Świderska, Hanna Kmita, and Małgorzata Wojtkowska. 2021. "Mitochondrial Processes during Early Development of Dictyostelium discoideum: From Bioenergetic to Proteomic Studies" Genes 12, no. 5: 638. https://doi.org/10.3390/genes12050638

APA Style

Mazur, M., Wojciechowska, D., Sitkiewicz, E., Malinowska, A., Świderska, B., Kmita, H., & Wojtkowska, M. (2021). Mitochondrial Processes during Early Development of Dictyostelium discoideum: From Bioenergetic to Proteomic Studies. Genes, 12(5), 638. https://doi.org/10.3390/genes12050638

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