Next Article in Journal
Could Targeted Pharmacotherapies Exert a “Disease Modification Effect” in Patients with Chronic Plaque Psoriasis?
Next Article in Special Issue
Basic Leucine Zipper Transcription Factors as Important Regulators of Leydig Cells’ Functions
Previous Article in Journal
New Insights into Psoriasis
Previous Article in Special Issue
MnHR4 Functions during Molting of Macrobrachium nipponense by Regulating 20E Synthesis and Mediating 20E Signaling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 21
10.3390/ijms232112846
ijms-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

Comprehensive and Quantitative Analysis of the Changes in Proteomic and Phosphoproteomic Profiles during Stimulation and Repression of Steroidogenesis in MA-10 Leydig Cells †

by
Zoheir B. Demmouche
1 and
Jacques J. Tremblay
1,2,*
1
Reproduction, Mother and Child Health, Room T3-67, CHU de Québec—Université Laval Research Centre, Québec, QC G1V 4G2, Canada
2
Centre for Research in Reproduction, Development and Intergenerational Health, Department of Obstetrics, Gynecology and Reproduction, Faculty of Medicine, Université Laval, Québec, QC G1V 0A6, Canada
*
Author to whom correspondence should be addressed.
†
This research was conducted as part of the requirements for a PhD degree.
Int. J. Mol. Sci. 2022, 23(21), 12846; https://doi.org/10.3390/ijms232112846
Submission received: 29 September 2022 / Revised: 21 October 2022 / Accepted: 21 October 2022 / Published: 25 October 2022
(This article belongs to the Special Issue Molecular Mechanisms of Steroid Hormone Biosynthesis and Action)
Browse Figures
Versions Notes

Abstract

:
Leydig cells produce testosterone, a hormone essential for male sex differentiation and spermatogenesis. The pituitary hormone, LH, stimulates testosterone production in Leydig cells by increasing the intracellular cAMP levels, which leads to the activation of various kinases and transcription factors, ultimately stimulating the expression of the genes involved in steroidogenesis. The second messenger, cAMP, is subsequently degraded to AMP, and the increase in the intracellular AMP levels activates AMP-dependent protein kinase (AMPK). Activated AMPK potently represses steroidogenesis. Despite the key roles played by the various stimulatory and inhibitory kinases, the proteins phosphorylated by these kinases during steroidogenesis remain poorly characterized. In the present study, we have used a quantitative LC-MS/MS approach, using total and phosphopeptide-enriched proteins to identify the global changes that occur in the proteome and phosphoproteome of MA-10 Leydig cells during both the stimulatory phase (Fsk/cAMP treatment) and inhibitory phase (AICAR-mediated activation of AMPK) of steroidogenesis. The phosphorylation levels of several proteins, including some never before described in Leydig cells, were significantly altered during the stimulation and inhibition of steroidogenesis. Our data also provide new key insights into the finely tuned and dynamic processes that ensure adequate steroid hormone production.

1. Introduction

Leydig cells are located in the interstitial space between the seminiferous tubules of the mammalian testis [1,2]. These cells are the source of androgens, the main one being testosterone [3]. Steroidogenesis is the biological process of converting cholesterol into steroid hormones, which involves cholesterol transport into the mitochondria, where steroidogenesis is initiated. In males, androgens are essential for male sex differentiation during fetal life and for initiating and maintaining spermatogenesis from puberty onwards. Androgens are also needed to acquire male secondary sex characteristics during puberty. Inadequate androgen production is associated with some cases of differences/disorders of sex development (DSD) in males [4].
Since androgens have pleiotropic roles in male reproductive function and overall health, their synthesis is regulated tightly. Steroidogenesis in Leydig cells is stimulated mainly by the pituitary luteinizing hormone (LH) [5]. The binding of LH to its G protein-coupled receptor (LHCGR) on Leydig cells activates adenylate cyclase, which leads to an increase in intracellular cAMP levels [6]. LH/cAMP-induced steroidogenesis then triggers the activation of several kinases, including protein kinase A (PKA), mitogen-activated protein kinase (MAPK), and calcium/calmodulin-dependent protein kinase I (CAMKI), which in turn phosphorylates several proteins required for increased steroid hormone synthesis (reviewed in [7]). This includes several transcription factors, such as the cAMP response element-binding protein (CREB)/CRE modulator (CREM), GATA4, and SF1 (reviewed in [8]). Furthermore, the de novo synthesis of the NR4A1/NUR77 (nuclear receptor subfamily 4, group A, member 1) transcription factor is required for maximal hormone-induced steroidogenesis [9]. The induction of Nr4a1 expression and the stimulation of steroidogenesis in Leydig cells also requires the release of Ca2+ from internal stores through the ryanodine receptors, leading to the activation of CAMKI [10]. Once adequate testosterone levels are reached, stimulation of steroidogenesis is blunted by two mechanisms. The first is the classic negative feedback loop, where testosterone acts at the level of the hypothalamus and pituitary to inhibit LH production [11]. The second is at the level of the Leydig cell itself. In Leydig cells, cAMP is degraded into AMP by phosphodiesterase (PDE) 8A, PDE8B, and PDE4 [12]. The ensuing increase in the intracellular AMP levels activates AMP-activated protein kinase (AMPK), a ubiquitous serine/threonine kinase best known as an energy balance sensor. AMPK is a heterotrimeric complex containing one catalytic subunit (alpha) and two regulatory subunits (beta and gamma) [13]. Once activated, AMPK potently blunts the LH/cAMP-induced steroidogenesis in Leydig cells [14].
We previously compared the transcriptome of Leydig cells that were either untreated, stimulated with Forskolin (Fsk, an agonist of adenylate cyclase), or co-treated with Fsk+AICAR (an agonist of AMPK). This led to the identification and characterization of several genes that were upregulated by the LH/cAMP stimulatory pathway and subsequently downregulated by the AMPK repressive pathway [14]. In addition to the changes in the transcriptome, the protein levels are also affected by the stimulation/repression of steroidogenesis in Leydig cells. However, global changes in protein and phosphoprotein levels have never been reported in this context. In the present work, we have used a quantitative mass spectrometry approach to elucidate the global and dynamic changes in the phosphoproteome of Leydig cells in response to stimulatory (Fsk/cAMP) and inhibitory (AICAR/AMPK) treatments.

2. Results

2.1. Validation of MA-10 Leydig Cell Responsiveness

Before analyzing the samples by quantitative LC-MS/MS, we first validated the responsiveness of the MA-10 Leydig cells to the different treatments. MA-10 Leydig cells were treated for 1 h with the vehicle (DMSO), Forskolin alone (Fsk), or Fsk+AICAR. Fsk is an agonist of adenylate cyclase, leading to increased intracellular cAMP levels and stimulation of steroidogenesis. AICAR is an agonist of the AMPK kinase, which we have identified as a potent repressor of hormone-activated steroidogenesis [14]. Although 1 h is sufficient to detect changes in the protein phosphorylation levels, it is usually too short to detect changes in the total protein levels by Western blot. We, therefore, isolated total RNA, which was used in a qPCR to quantify the mRNA levels for the steroidogenic acute regulatory (Star) protein. The Star gene codes for the steroidogenic acute regulatory (STAR) protein, a protein essential for hormone-induced cholesterol transport into mitochondria and, consequently, steroidogenesis [15]. The Star gene is an excellent marker of the dynamic steroidogenic process in Leydig cells, as its expression is strongly induced by LH/Fsk/cAMP and repressed by AICAR/AMPK ([14] and reviewed in [8]). As shown in Figure 1, the Star mRNA levels from the three samples used in LC-MS/MS (described below) were increased six- to ten-fold in the presence of Fsk. As expected [14], this increase was potently repressed when the MA-10 Leydig cells were co-treated with AICAR in addition to Fsk (Figure 1). These results confirm the responsiveness of the MA-10 Leydig cells to Fsk and AICAR and validate that the protein samples isolated from these cells are suitable for LC-MS/MS analysis.

2.2. Treatment of MA-10 Leydig Cells with Fsk or Fsk+AICAR Significantly Affects the Levels of 20 Proteins

Total proteins from the MA-10 Leydig cells treated for 1 h with either the vehicle (DMSO), Fsk alone, or Fsk+AICAR were extracted, digested, and quantitively analyzed by LC-MS/MS. A total of 5887 proteins were identified, and 4819 were quantified (data not shown). The level of the majority of these proteins was not significantly changed by the treatments due to the short treatment time. Only 20 proteins (listed in Table 1) were significantly affected between treatments (DMSO vs. Fsk, Fsk vs. Fsk+AICAR, and DMSO vs. Fsk+AICAR). The protein levels either increased (shown in red in Table 1) or decreased (shown in blue in Table 1). In addition, the same protein was sometimes found in different comparison groups. For instance, the protein levels of the orphan nuclear receptor NR4A1 (NUR77), a known regulator of the hormone-induced steroidogenic gene expression in Leydig cells [10,16], was increased by 2.4-fold by the Fsk treatment (DMSO vs. Fsk group), and this increase was blunted after the AICAR treatment (down by 2.34-fold; Fsk vs. Fsk+AICAR group), as previously reported [14]. Since the levels of the NR4A1 protein are back down to the control levels after treatment with Fsk+AICAR, the NR4A1 protein levels were not significantly changed in the DMSO vs. Fsk+AICAR group. These changes in NR4A1 protein levels (increased by Fsk/cAMP and reduced upon activation of AMPK by AICAR) serve as a positive control and validate the quantitative LC-MS/MS approach used.
The level of the protein Dymeclin, a Golgi-associated protein believed to be involved in intracellular trafficking [17,18,19], was the most increased after 1 h of Fsk treatment (23-fold; Table 1). Similar to NR4A1, the Fsk-dependent increase of Dymeclin was prevented when AMPK was activated by AICAR (decreased by 23-fold and back to control levels). Another protein affected by both Fsk and AICAR is the Mpv17-like protein 2 (MPV17L2). MPV17L2 is required for the assembly and stability of the mitochondrial ribosome, and in its absence, protein synthesis in mitochondria is impaired, and mitochondrial DNA aggregates [20]. Contrary to NR4A1 and Dymeclin, the MPV17L2 levels were reduced by 5.35-fold in the presence of Fsk and increased by 5.23-fold by AICAR-mediated AMPK activation (Table 1).
As summarized in Table 1, in addition to Dymeclin (up by 23-fold) and NR4A1 (up by 2.4-fold), the top Fsk-induced proteins included the homologous-pairing protein 2 homolog (PSMC3IP; up by 5-fold), acyl-coenzyme A dehydrogenase family member 12 (ACAD12; up by 4-fold), and N-acetyltransferase domain containing 1 (NATD1; up by 3.6-fold). In addition to MPV17L2, the most induced proteins in the presence of AICAR (activation of AMPK) included the cyclic AMP-responsive element-binding protein 1 (CREB1; up by 18-fold), while decreased proteins, in addition to Dymeclin and NR4A1, included the Lemur serine/threonine-protein kinase 2 (LMTK2; down by 2.8-fold). These results showed that after a short 1-h treatment, the level of some proteins involved in mitochondrial function and in steroidogenesis was already altered.

2.3. Identification of Proteins Differentially Phosphorylated in MA-10 Leydig Cells in Response to Fsk and AICAR

Our quantitative analysis of the phosphoproteome of MA-10 Leydig cells revealed significant changes in the phosphorylation levels of several phosphopeptides corresponding to various proteins in response to the treatments. The treatments included DMSO (the vehicle used as the control), Fsk to increase the cAMP levels mimicking LH-stimulated steroidogenesis, and Fsk+AICAR simultaneously, where AICAR activates AMPK, a kinase previously reported to potently inhibit hormone-induced steroidogenesis [14]. The top 10 differentially phosphorylated (up and down) phosphopeptides for each pairwise comparison (Fsk vs. DMSO, Fsk+AICAR vs. Fsk, Fsk+AICAR vs. DMSO) are presented in Table 2.

2.3.1. Proteins Differentially Phosphorylated in Response to Fsk

We first compared the phosphopeptides between the control and Fsk-treated MA-10 Leydig cells to identify changes in protein phosphorylation during the stimulatory phase of steroidogenesis. This led to the identification of 12,125 phosphopeptides, of which 8471 were quantified. As shown in Figure 2A, the quantified phosphopeptides were mapped on a volcano plot based on the significance and the ratio between the Fsk-treated and control (DMSO) MA-10 Leydig cells. A total of 61 phosphopeptides had significantly different phosphorylation levels between these two conditions; the phosphorylation of 35 phosphopeptides was increased, while the phosphorylation of 26 phosphopeptides was decreased following the Fsk treatment. These 61 phosphopeptides correspond to 60 unique phosphoproteins and are shown along with their relative phosphorylation levels between the two conditions on the heatmap in Figure 2B.
The proteins exhibiting increased phosphorylation after the Fsk treatment included the nuclear receptor NR4A1, the signal recognition particle SRP14, the protein kinase A anchoring protein AKAP1, and the serine/threonine kinase RAF1. In comparison, the proteins showing a reduction in their phosphorylation level in response to Fsk included the microtubule-associated protein tau (MAPT) and the Par-3 family cell polarity regulator protein PARD3 (Figure 2A,B). In response to Fsk in MA-10 Leydig cells, the differentially phosphorylated proteins are mostly implicated in the metabolic processes and the regulation of translation (Figure 2C). A pathway analysis revealed that these proteins are implicated mainly in the mTOR and Hippo signaling pathways, two pathways involved in metabolism (Figure 2D). As shown in Figure 2E, the differentially phosphorylated proteins are associated with various cellular components. Finally, the 60 differentially phosphorylated proteins are associated with various molecular functions, including binding to protein phosphatase, phosphatidyl inositol phosphate, and nucleic acid (Figure 2F).

2.3.2. Proteins Differentially Phosphorylated in Response to AMPK Activation

The activation of AMPK with AICAR represses LH/cAMP-induced steroidogenesis but does not significantly affect the unstimulated steroidogenic cells [14]. Therefore, to identify potential target proteins of AMPK, we compared all phosphopeptides between the Fsk- and Fsk+AICAR-treated MA-10 Leydig cells. The differences between those two conditions should reflect mainly the proteins phosphorylated by AMPK. This comparison identified 12,125 phosphopeptides, of which 8581 phosphopeptides were quantified. It also revealed that the phosphorylated amino acid in cells exposed to AICAR was, by far, a serine, as shown by the sequence motif presented in Figure 3A. Quantified phosphopeptides were then mapped on a volcano plot based on the significance and the ratio between the Fsk+AICAR- and Fsk-treated cells (Figure 3B). A total of 46 phosphopeptides had significantly different phosphorylation levels between these two conditions; in the Fsk+AICAR group, the phosphorylation of 27 phosphopeptides was increased, while the phosphorylation of 19 phosphopeptides was decreased compared to the Fsk group (Figure 3B,C). These 46 phosphopeptides correspond to 44 unique phosphoproteins, as some proteins have more than one phosphorylated residue (phosphopeptide sequence). These 44 phosphoproteins are presented along with their relative phosphorylation levels between the two conditions on the heatmap presented in Figure 3C.
The proteins with increased phosphorylation in the AICAR group included NFATC3 (transcription factor), PARD3, and MAPT; an example of a protein with reduced phosphorylation in the AICAR group is the nuclear mitotic apparatus protein NUMA1 (Figure 3A,B). The proteins differentially phosphorylated between the Fsk+AICAR and Fsk groups are mostly implicated in the microtubule polymerization or its organization (Figure 3D). This is supported by the localization of these proteins, mostly with microtubules (Figure 3F), with their main molecular functions being microtubule binding, histone binding, and deubiquitinase activity (Figure 3G). The pathway analysis revealed that the 44 proteins are mainly associated with the MAPK signaling pathway and mRNA surveillance pathway (Figure 3E).

2.3.3. Proteins Differentially Phosphorylated in Response to Combined Fsk+AICAR

Comparisons were made between the Fsk+AICAR and the control (DMSO) groups to identify all the changes in the protein phosphorylation triggered by the simultaneous treatment of the MA-10 Leydig cells with Fsk (stimulatory) and AICAR (inhibitory). A total of 12,125 phosphopeptides were identified, of which 8711 were quantified. As shown in Figure 4A, the quantified phosphopeptides were mapped on a volcano plot based on the significance and the ratio between the Fsk+AICAR and control (DMSO) cells. A total of 73 phosphopeptides exhibited significantly different phosphorylation levels between these two conditions; the phosphorylation of 38 phosphopeptides was increased in the MA-10 Leydig cells treated with Fsk+AICAR, while the phosphorylation of 35 phosphopeptides was decreased. These 73 phosphopeptides correspond to 68 unique phosphoproteins, as some proteins harbor more than one phosphorylation site, such as MAP1A. The 68 phosphoproteins identified, along with their relative phosphorylation levels between the two conditions, are shown on the heatmap presented in Figure 4B.
Some of the proteins with increased phosphorylation in the MA-10 Leydig cells in response to Fsk+AICAR included NR4A1, NUMA1, the protein GOLGA5 involved in maintaining the Golgi structure, and the microtubule-associated proteins MAP1A and MAP1B. MAP1A was also present amongst the proteins with reduced phosphorylation. This is explained by the fact that MAP1A contains two distinct phosphopeptides: one showed increased phosphorylation, while the other showed reduced phosphorylation in the presence of Fsk+AICAR (Figure 4A). In terms of the biological processes, differentially phosphorylated proteins were mainly associated with microtubule polymerization, depolymerization, and organization (Figure 4C). Consistent with this, these proteins are mostly localized with microtubules, and their main functions are to bind to the microtubules, actin, and tubulin (Figure 4F). The pathway analysis associated the differentially phosphorylated proteins with the axon guidance pathway (Figure 4D).
Taken together, these results provide the first description of the global changes in the phosphoproteome of MA-10 steroidogenic Leydig cells after a stimulatory treatment (Fsk) and an inhibitory treatment (activation of AMPK by AICAR).

2.4. Phosphorylation Trajectory of Representative Proteins

Several proteins exhibited an altered phosphorylation status in response to Fsk and to Fsk+AICAR in the MA-10 Leydig cells (Figure 2, Figure 3 and Figure 4). Although the protein levels remained unchanged for the majority of proteins (Table 1), changes in the phosphorylation level in response to Fsk and to Fsk+AICAR could still be the result of a concomitant change in the protein level. Therefore, the phosphorylation status of the representative proteins was plotted and compared to their protein levels, allowing for better visualization of the phosphorylation trajectory of a given protein.
The phosphorylation level of NR4A1, an orphan nuclear receptor known to regulate the expression of several steroidogenic genes, was increased by more than 20-fold after the Fsk treatment, while the total NR4A1 protein level was only increased by 2.4-fold. This increase tended to be attenuated when cells were treated with AICAR in addition to Fsk (Figure 5A). Interestingly, the phosphorylation levels of the progesterone membrane receptor components 2 (PGRMC2), a protein belonging to the membrane-associated progesterone receptor (MAPR) protein family, was significantly increased when cells were treated with Fsk and with Fsk+AICAR (Figure 5B). No change in the total PGRMC2 protein level was noted. The NUMA1 exhibited the same phosphorylation pattern as PGRMC2 (Figure 5C), with an increase in phosphorylation of the phosphopeptide QAASSQEPSELEELR. On the other hand, a different phosphopeptide (IASSSSENNFLSGSPSSPMGDILQTPQFQMR) of the NUMA1 protein remained unchanged in response to Fsk, whereas AICAR caused a significant reduction in its phosphorylation level (Figure 5D). Finally, the phosphorylation level of two proteins, known to be downstream of AMPK, MAPT, and PARD3, was significantly reduced when the MA-10 Leydig cells were treated with Fsk. However, this decrease in phosphorylation was no longer apparent when the cells were co-treated with Fsk+AICAR, suggesting that AMPK phosphorylates these proteins (Figure 5E,F). The total protein level of PGRMC2, NUMA1, MAPT, and PARD3 was not affected by the treatments. These results validate AMPK activation in the MA-10 Leydig cells treated with AICAR.
Altogether, these results reveal that protein phosphorylation, in response to Fsk and Fsk+AICAR in MA-10 Leydig cells, is a very dynamic process and is mainly independent of the changes in protein levels.

3. Discussion

Several studies have reported changes in the expression of a single gene or its protein levels in response to hormone-induced steroidogenesis. However, Leydig cell steroidogenesis is a strictly regulated process that requires changes in gene expression, protein levels, and protein phosphorylation of many genes and proteins. Therefore, the main objective of this work was to elucidate the global phosphoproteomic changes that occur in MA-10 Leydig cells in response to Forskolin-induced steroidogenesis, followed by AICAR-activated AMPK, leading to repression of steroidogenesis.

3.1. Global Variation in Protein Levels in MA-10 Leydig Cells 1 h after Fsk or Fsk+AICAR Treatment

We first determined the changes in the global protein levels and found that the level of only 20 proteins was significantly modified in the Leydig cells treated for 1 h with either Forskolin alone (stimulatory) or Forskolin+AICAR (inhibitory). One of those proteins was the nuclear receptor NR4A1. Our global proteomic analysis revealed that the NR4A1 protein levels were increased by 2.4-fold after the Forskolin treatment, and this increase was blunted when cells were co-treated with AICAR. This is consistent with previous studies from our group that reported the pattern of Nr4a1 mRNA levels [14], as well as a time course of NR4A1 protein levels [10], in response to (Bu)2-cAMP treatment of mouse MA-10 Leydig cells and primary Leydig cells from rats. Therefore, the observed changes in the NR4A1 protein levels in our current study validate our experimental approach.
Another protein affected by the treatments was Dymeclin, which was induced 23-fold in response to Fsk and then lost upon AMPK activation. Although it has never been reported in Leydig cells, Dymeclin is believed to be implicated in vesicle trafficking, as it shuttles between the cytosol and the Golgi apparatus in a highly dynamic manner [18,19]. Vesicle trafficking is an important process for steroidogenesis. For instance, the life cycle of lipid droplets, which are a source of cholesterol substrate for steroid hormones [21], involves intracellular trafficking (reviewed in [22]). The significant and rapid changes in Dymeclin protein levels in response to Fsk and Fsk+AICAR suggest that Dymeclin might play an active role in hormone-regulated steroidogenesis in Leydig cells.
The transcription factor, CREB, is known to activate steroidogenic gene expression in Leydig cells [23,24,25,26,27,28]. Despite this, the CREB protein levels remained unchanged in response to the Fsk/cAMP treatment, which is identical to what we observed in our current global proteomic analysis. Unexpectedly, the CREB protein levels were increased by 18-fold in the MA-10 Leydig cells treated with AICAR to turn on the repressive kinase, AMPK (Fsk+AICAR-treated MA-10 Leydig cells compared to the Fsk-treated cells). This would appear to suggest a repressive role for CREB in steroidogenesis. Although mainly known as an activator, CREB has also been reported to act as a repressor. For instance, CREB represses the expression of the cFos gene [29]. Furthermore, alternate exon usage is also known to switch CREB from an activator to a repressor [30]. It is, therefore, possible that activated AMPK influences CREB exon usage and that CREB would, therefore, repress steroidogenic gene expression in that context.
The MPV17L2 protein levels were also modulated by the treatments. MPV17L2 is an integral mitochondrial inner membrane protein [20], and as such, it is implicated in mitochondrial ribosome assembly [20] and in reactive oxygen species (ROS) production [31], two processes known to influence steroidogenesis in Leydig cells [32]. We found that the MPV17L2 protein levels increased 5-fold in the presence of AICAR, while it was reduced by 5-fold in the Fsk-treated cells. Our data are consistent with transcriptomic data from the pancreas, where the Mpv17l2 mRNA levels are decreased by 1.6-fold in AMPK-deficient mice [33], indicating that activated AMPK increases its expression.
Overall, our data show that globally, the levels of 20 proteins, including some never before reported in Leydig cells, are rapidly and differentially modulated in response to Fsk (stimulatory) and to Fsk+AICAR (inhibitory).

3.2. Comprehensive Analysis of the Phosphoproteome during Stimulation and Inhibition of Steroidogenesis in MA-10 Leydig Cells

In Leydig cells, the response to hormone stimulation involves the activation of various signaling pathways, leading to increased activity of various kinases, including PKA, MAPK, CAMKI, and AMPK. These kinases then phosphorylate various proteins, thus ensuring a proper cellular response to the stimulus. Although some phosphorylated proteins have been identified, the majority of studies have focused on a single protein at the time. One study reported the changes in the phosphoproteome of Leydig cells in which the cAMP levels were maintained artificially high by inhibiting all phosphodiesterase activity [34]. This study identified alterations in several phosphosites, some of which were also detected in our current study. However, constitutively high cAMP levels might trigger non-physiological responses. In our present work, we used Fsk, an agonist of adenylate cyclase that leads to a typical increase in cAMP levels, triggering a stimulatory response, and AICAR, an agonist of AMPK, leading to an inhibitory response. In both cases, the levels of cAMP, the main second messenger in Leydig cells, followed a normal increase–decrease rate typical of a more classic response [14]. Using this treatment scheme and a quantitative phosphoproteomics approach, we now report the global changes in the phosphoproteome that occur in Leydig cells in response to both a stimulatory and an inhibitory signal. We found that the phosphorylation level of a large number of proteins was modified in response to either Fsk (stimulatory response) or Fsk+AICAR (inhibitory response).
The proteins displaying increased phosphorylation levels in the MA-10 Leydig cells treated with Fsk include AKAP1 (also known as AKAP121), RAF1, NR4A1, and SRP14. The phosphorylation of both AKAP1 and RAF1 was previously reported as being altered in the MA-10 Leydig cells depleted of PDE4 and PDE8 activity [34]. PDE activity is required for normal Leydig cell function. Indeed, PDEs regulate steroidogenesis by their ability to degrade cAMP into AMP, thus inactivating several cAMP-dependent pathways and kinases, such as PKA, while simultaneously activating the repressive kinase AMPK [12]. AKAP1 is a member of the AKAP family of scaffold proteins. AKAP1 is known to interact with PKA and therefore regulates its intracellular localization, especially at the outer mitochondrial membrane (OMM) (reviewed in [35]), where PKA can phosphorylate its target proteins [36]. A key target of PKA at the OMM for the stimulation of steroidogenesis is the STAR protein [36,37]. PKA-dependent phosphorylation of STAR increases its cholesterol shuttling activity, leading to enhanced steroidogenesis [36,37]. AKAP1 is, therefore, considered an important mitochondrial signaling hub, and the mitochondrion is an essential organelle for steroidogenesis in Leydig cells. To date, AKAP1 is not known to be phosphorylated. Our current work demonstrates that AKAP1 is phosphorylated in response to Fsk. The physiological implication of AKAP1 phosphorylation remains to be established, although a previous study suggested that phosphorylation of AKAP1 can modulate the AKAP1-PKA association [38].
Another protein exhibiting increased phosphorylation in response to Fsk is SRP14. SRP14 is involved in protein translation and, more specifically, in the elongation arrest and efficient translocation of proteins into the endoplasmic reticulum [39]. This is consistent with the requirement of a de novo protein synthesis in hormone-induced steroidogenic cells [10]. Interestingly, the level of SRP14 protein in glioblastoma cells is increased by CYP17A1 [40], an enzyme implicated in steroidogenesis. As we observed in the Fsk-treated MA-10 Leydig cells, SRP14 was also found to be phosphorylated during a large screen of phosphorylated proteins in Ras-transformed cells [41]. The significance of SRP14 phosphorylation, however, remains unknown.
The proteins displaying reduced phosphorylation levels after Fsk treatment include MAPT and PARD3. Phosphorylation of both proteins was restored when AMPK was activated (Fsk+AICAR-treated cells). MAPT, a protein involved in microtubule assembly and stabilization, is a known target of AMPK [42]. Furthermore, phosphorylation of MAPT is reduced in the MA-10 Leydig cells deficient in PDE4 and PDE8 activity to maintain high cAMP levels [34]. These data are consistent with our present data showing a reduction in MAPT phosphorylation in the context of high cAMP levels (Fsk-treated cells) and increased phosphorylation in AMPK-activated cells. Phosphorylation of MAPT is known to influence protein stability and degradation (reviewed in [43]).
The protein PARD3 is part of the PAR3/PAR6/aPKC complex implicated in tight-junction assembly [44]. Activation of AMPK facilitates tight-junction assembly [45]. Whether activated AMPK interacts with and directly phosphorylates the PAR3/PAR6/aPKC complex remains uncertain [45]. However, our current data now indicate that PARD3 is phosphorylated upon AMPK activation. Furthermore, PARD3 is known to be phosphorylated by the PAR1 kinase, leading to PARD3 translocation from the tight junction to the cytosol [46]. Once in the cytosol, PARD3 promotes the interaction between PP1A phosphatase and LATS1, a mediator of the Hippo signaling pathway [46]. This results in the dephosphorylation and inactivation of LATS1 and, consequently, dephosphorylation, activation, and nuclear translocation of the TAZ transcriptional co-activator [46]. Once in the nucleus, TAZ represses Leydig cell steroidogenesis by suppressing the expression of several steroidogenic genes via an interaction with the nuclear receptor NR4A1, preventing it from binding to DNA and activating its target genes [47]. It is tempting to speculate that the phosphorylation of PARD3, which we observed in the AMPK-activated Leydig cells, might lead to the same events, therefore contributing to a repression of steroidogenesis.
The transcription factors belonging to the nuclear factor of activated T cells (NFAT) family are well-known transcriptional activators (reviewed in [48,49]). In our phosphoproteomic screen, the phosphorylation of the transcription factor NFATC3 remained unchanged in the presence of Fsk alone. However, it was significantly increased in the MA-10 Leydig cells treated with Fsk+AICAR, indicating that NFATC3 becomes phosphorylated upon activation of AMPK. Interestingly, phosphorylated NFAT transcription factors are known to be sequestered in the cytoplasm since phosphorylation masks the nuclear localization sequence [50,51]. The phosphorylation of NFATC3 that occurs upon activation of AMPK in MA-10 Leydig cells would be consistent with the exclusion of NFATC3 from the nucleus and a reduction in its ability to activate gene expression, therefore contributing to reduced steroidogenesis.
Another group of differentially phosphorylated proteins includes those where the phosphorylation level was significantly different only in the Fsk+AICAR-treated MA-10 Leydig cells vs. control cells (vehicle-treated). This included the NUMA1 protein, which was differentially phosphorylated on two different peptides; one exhibited increased phosphorylation, while the phosphorylation of another residue of the protein was significantly decreased. NUMA1 is a component of the nuclear matrix (reviewed in [52]) and is implicated in ciliogenesis and autophagy, as is AMPK [53]. NUMA1 was indeed proposed as a potential target of AMPK, although this remains to be validated [54]. Our current results support direct phosphorylation of NUMA1 by AMPK.
The PGRMC2 protein also exhibited increased phosphorylation in the MA-10 Leydig cells treated with Fsk+AICAR vs. controls. PGRMC2 is a membrane-associated progesterone receptor (MAPR) present in Leydig cells and is believed to play a role in the progesterone auto/paracrine action on Leydig cells [55,56] in addition to potentially regulating the activity of some cytochrome P450 enzymes [57]. We showed that PGRMC2 is present in MA-10 Leydig cells and that it is phosphorylated after AMPK activation.
The MA-10 Leydig cell line is the gold standard for studying Leydig cells, although they are not identical to the cultured primary Leydig cells. The MA-10 Leydig cell line was originally established from a Leydig cell tumor in a mouse (M5480P) [58] and has since been meticulously characterized. Treatment of MA-10 cells with the luteinizing hormone, Forskolin, or cAMP analogs, increases steroid hormone production in the same way as in normal Leydig cells ([58] and reviewed in [59]). The rate-limiting step in hormone-induced steroidogenesis is the transport of cholesterol into the mitochondria, where steroidogenesis is initiated (reviewed in [60]); this is also true for MA-10 Leydig cells [10,15,61]. MA-10 cells mainly produce progesterone because of a mutation in the Cyp17a1 coding sequence [58]. MA-10 Leydig cells are, therefore, suitable for the study of hormone-induced steroidogenesis and, more specifically, the early steps of steroidogenesis [62]. Since our current work focuses on the modulation of protein phosphorylation that occurs early following hormone treatment, the MA-10 Leydig cell line is a suitable and convenient model to use.
In conclusion, using a quantitative approach, our present study determined the global phosphoproteomic profile of MA-10 Leydig cells in response to Fsk (stimulation of steroidogenesis) and AICAR-mediated AMPK activation (repression of steroidogenesis). Our results indicate that steroidogenesis is a very dynamic process that involves differential protein phosphorylation by various kinases, altogether contributing to the fine-tuned regulation that is needed to achieve the appropriate production of steroid hormones.

4. Materials and Methods

4.1. Cell Culture

The MA-10 cell line was obtained from ATCC (Cat# CRL-3050, RRID:CVCL_D789. ATCC, Manassas, VA, USA). The MA-10 cells were grown in a DMEM/F12 medium supplemented with 2.438 g/L sodium bicarbonate, 3.57 g/L HEPES, and 15% horse serum on gelatin-coated plates. Penicillin and streptomycin sulphate were added to the cell culture media to a final concentration of 50 mg/L, and the cells were kept at 37 °C, 5% CO2, in a humidified incubator. The MA-10 Leydig cell line was validated by morphology and by quantifying the progesterone output, as previously described.

4.2. Chemicals

The AMPK agonist, AICAR, was obtained from Tocris Bioscience (Minneapolis, MN, USA). Forskolin (Fsk) was purchased from Sigma-Aldrich, Canada (Oakville, ON, Canada).

4.3. RNA Isolation, Reverse Transcription, and Quantitative PCR

The MA-10 Leydig cells were cultured in the presence of either DMSO (vehicle), Fsk (10 µM), or Fsk+AICAR (10 µM + 1 mM) for 60 min. Isolation of RNA, cDNA synthesis, and reverse transcription-quantitative PCR (RT-qPCR) were performed as previously described [63,64]. Briefly, the total RNA from the MA-10 Leydig cells grown and treated, as described above, was isolated using TRIZOL (Life Technologies, Burlington, ON, Canada) and reverse-transcribed using the iScript Advanced cDNA Synthesis Kit (Bio-Rad Laboratories, ON, Canada). A quantitative real-time PCR was performed using a CFX96™ Real-Time PCR Detection System (Bio-Rad Laboratories, ON, Canada) along with the SsoAdvanced Universal SYBR Green Supermix kit (Bio-Rad Laboratories, ON, Canada) according to the manufacturer’s protocols. Relative expression was normalized to the expression of Rpl19, which was used as an internal control. The mouse Star primers were: forward 5′-GTT CCT CGC TAC GTT CAA GC-3′ and reverse 5′-GAA ACA CCT TGC CCA CAT CT-3′. The mouse Rpl19 primers were: forward 5′-CTG AAG GTC AAA GGG AAT GTG-3′ and reverse 5′-GGA CAG AGT CTT GAT GAT CTC-3′. The melting temperature of both primer couples was 62.6 °C.

4.4. Sample Preparation for LC-MS/MS

The MA-10 Leydig cells (5.8 million) were plated in 15-cm plates and then treated at 80% confluence with either the vehicle (DMSO), Fsk (10 μM) or Fsk+AICAR (10 µM + 1 mM) for 60 min at 37 °C and 5% CO2. This experiment was repeated three times in triplicate.
Protein extraction, digestion, phosphopeptide enrichment, and mass spectrometry analyses were performed by the proteomics platform of the CHU de Quebec Research Centre (Quebec City, QC, Canada). Cell pellets were resuspended in an extraction buffer (50 mM ammonium bicarbonate, 0.5% sodium deoxycholate (SDC), 50 mM DTT, protease inhibitor cocktail, phosphatase inhibitor PhosSTOP, and 1 μM pepstatin) and vortexed for 1 min. They were then sonicated with a microprobe (Sonic Dismembrator 550, Fisher Scientific) 20 times for 1 s (on/off). The extract was centrifuged at 16,000× g for 15 min, the supernatant was collected, and then acetone-precipitated. The pellet was resuspended in ammonium bicarbonate 50 mM/SDC 1% and treated with bioruptor 15 times for 30 s (on/off). The protein concentration was determined by performing a Bradford assay.
Protein aliquots of 10 μg for the global protein analysis or 2 mg for the phosphopeptide analysis for each condition (vehicle, Fsk, Fsk+AICAR) were used. The proteins were first heated at 95 °C for 5 min and reduced with 0.2 mM DTT at 37 °C 30 min. They were then alkylated with 0.8 mM IAA (iodoacetamide) for 30 min at 37 °C in the dark. This was followed by digestion with trypsin (1:50) at 37 °C overnight. Digestion was stopped with the addition of formic acid; the sodium deoxycholate was removed after centrifugation. The peptides were desalted using a stage-tip column (for the global protein analysis) or HLB column (for the phosphopeptide analysis).
Phosphopeptide enrichment was performed using the High Select Phosphopeptide Enrichment Kits & Reagents kit (Fisher Scientific Canada, Catalog number A32993) according to the manufacturer’s protocol. Two (2) mg of dry protein from each condition (vehicle, Fsk, Fsk+AICAR) were resuspended in 150 μL of binding/equilibration buffer (pH < 3). The TiO2 spin-tip columns were washed and then equilibrated before adding 150 μL of the resuspended samples. The TiO2 spin-tip columns loaded with the samples were centrifuged twice at 1000× g for 5 min to bind the phosphopeptides. The TiO2 spin-tip columns were then washed twice with the binding/equilibration buffer, once with the wash buffer and then with HPLC water. The phosphopeptides were finally eluted with 50 μL of phosphopeptide elution buffer twice.

4.5. Quantitative Sample Analysis by LC-MS/MS

Aliquots of 1 μg (for global protein analysis) and 1 mg of enriched phosphopeptide were analyzed by nanoLC-MS/MS using a Dionex UltiMate 3000 nanoRSLC chromatography system (Thermo Fisher Scientific, Ontario, Canada), connected to an Orbitrap Fusion mass spectrometer (Thermo Fisher Scientific Ontario, Canada). Peptides were trapped at 20 μL/min in a loading solvent (2% acetonitrile, 0.05% TFA) on a 5 mm × 300 μm C18 Pepmap cartridge pre-column (Thermo Fisher Scientific, Ontario, Canada) for 5 min. The pre-column was then switched online with a Pepmap Acclaim column (Thermo Fisher Scientific, Ontario, Canada) 50 cm × 75 μm internal diameter separation column. The peptides were eluted with a linear gradient of 5–40% solvent B (A: 0.1% formic acid, B: 80% acetonitrile, 0.1% formic acid) over 270 min for a total of a 300-min run at 300 nL/min for the global protein analysis. For the phosphopeptide analysis, elution with the solvent gradient was completed in 90 min for a total run of 120 min. Mass spectra were acquired using a data-dependent acquisition mode using Thermo Xcalibur software, version 4.1.50. Full-scan mass spectra (350 to 1800 m/z) were acquired in the orbitrap using an AGC target of 4e5, a maximum injection time of 50 ms, and a resolution of 120,000. An internal calibration using lock mass on the m/z 445.12003 siloxane ion was used. Each MS scan was followed by the acquisition of fragmentation involving MS/MS spectra of the most intense ions for a total cycle time of 3 s (top speed mode). The selected ions were isolated using the quadrupole analyzer in a window of 1.6 m/z and fragmented by higher energy collision-induced dissociation (HCD) with 35% of collision energy. The resulting fragments were detected by the linear ion trap at a rapid scan rate with an AGC target of 1 × 104 and a maximum injection time of 50 ms. The dynamic exclusion of previously fragmented peptides was set for a period of 30 s and a tolerance of 10 ppm.

4.6. Quantitative Data Analyses

Spectra were searched against the Uniprot Reference mus musculus (61,295 entries) using the Andromeda module of the MaxQuant software, v. 1.6.10.43. The trypsin/P enzyme parameter was selected with two possible missed cleavages. The carbamidomethylation of cysteines was set as a fixed modification, while methionine oxidation and phosphorylation (Ser, Thr, Tyr) were set as variable modifications. The mass search tolerances were 5 ppm and 0.5 Da for MS and MS/MS, respectively. For the protein validation, a maximum false discovery rate of 1% at the peptide and protein levels was used based on a target/decoy search. MaxQuant was also used for label-free quantifications. The ‘match between runs’ option was used with a 20 min value as the alignment time window and 0.7 min as the match time window. For the global protein analysis, the protein group files were used, while for the phosphopeptide analysis, the modificationSpecifPeptides file was used. RStudio 1.2.5019 was used for data processing. A normalization step was performed using the median of the median intensities of each condition. The missing peptide intensity values were replaced by a noise value corresponding to the 1% percentile of the normalized value for each condition. A peptide was considered quantifiable only if at least three intensity values in one of the two conditions were present and with a minimum of two peptides (for the global protein analysis).

4.7. Statistical Analysis

For the RT-qPCR, statistical analyses were carried out using a nonparametric one-way ANOVA on ranks via the Kruskal–Wallis test, followed by a post-hoc Mann–Whitney U test. A p-value of < 0.05 was considered significant. For the LC-MS/MS analysis, a p-value limma test, a limma q-value (Benjamin–Hochberg correction), and a z-score were calculated. Phosphopeptides or proteins were considered variants if the q-value was < 0.05 and the z-score was +/− 1.96. All statistical analyses were performed using the OriginPro Version 2021 software (www.originlab.com, accessed on 10 August 2022) (OriginLab Corporation, Northampton, MA, USA).

Author Contributions

Z.B.D. performed all the experiments. Z.B.D. and J.J.T. analyzed and interpreted the data. J.J.T. conceived the study and coordinated and supervised the project. Z.B.D. drafted the manuscript with the assistance of J.J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a grant from the Canadian Institutes of Health Research (CIHR) (funding reference number PJT-148738) to JJT.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All mass spectrometry data (raw files and MaxQuant search result files) are publicly available on the ProteomeXchange repository (www.proteomexchange.org) with the identifier PXD037514.

Acknowledgments

We would like to thank Sylvie Bourassa (proteomics platform of the CHU de Quebec Research Centre, Quebec City) for her assistance with the quantitative protein and phosphoprotein mass spectrometry approach and Lisa Lahens for her assistance with the OriginPro software.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leydig, F. Zur anatomie der mannlichen geschlechtsorgane und analdrusen der saugethiere. Z. Wiss. Zool. 1850, 2, 1–57. [Google Scholar]
  2. Bouin, P.; Ancel, P. Recherches sur les cellules interstitielles du testicule des mammiferes. Arch. Zool. Exp. Gen. 1903, 1, 437–523. [Google Scholar]
  3. Hall, P.F.; Irby, D.C.; De Kretser, D.M. Conversion of cholesterol to androgens by rat testes: Comparison of interstitial cells and seminiferous tubules. Endocrinology 1969, 84, 488–496. [Google Scholar] [CrossRef]
  4. Hughes, I.A. Consensus statement on management of intersex disorders. Arch. Dis. Child. 2005, 91, 554–563. [Google Scholar] [CrossRef]
  5. Ewing, L.L.; Eik-Nes, K.B. On the formation of testosterone by the perfused rabbit testis. Can. J. Biochem. 1966, 44, 1327–1344. [Google Scholar] [CrossRef] [PubMed]
  6. Dufau, M.L.; Tsuruhara, T.; Horner, K.A.; Podesta, E.; Catt, K.J. Intermediate role of adenosine 3′:5′-cyclic monophosphate and protein kinase during gonadotropin-induced steroidogenesis in testicular interstitial cells. Proc. Natl. Acad. Sci. USA 1977, 74, 3419–3423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Tremblay, J.J. Molecular regulation of steroidogenesis in endocrine Leydig cells. Steroids 2015, 103, 3–10. [Google Scholar] [CrossRef]
  8. de Mattos, K.; Viger, R.S.; Tremblay, J.J. Transcription factors in the regulation of Leydig cell gene expression and function. Front. Endocrinol. 2022, 13, 881309. [Google Scholar] [CrossRef]
  9. Martin, L.J.; Tremblay, J.J. The human 3β-hydroxysteroid dehydrogenase/Δ5-Δ4 isomerase type 2 promoter is a novel target for the immediate early orphan nuclear receptor Nur77 in steroidogenic cells. Endocrinology 2005, 146, 861–869. [Google Scholar] [CrossRef]
  10. Martin, L.J.; Boucher, N.; Brousseau, C.; Tremblay, J.J. The orphan nuclear receptor NUR77 regulates hormone-induced StAR transcription in Leydig cells through cooperation with Ca2+/calmodulin-dependent protein kinase I. Mol. Endocrinol. 2008, 22, 2021–2037. [Google Scholar] [CrossRef] [Green Version]
  11. Darney, K.J.; Ewing, L. Autoregulation of testosterone secretion in perfused rat testes. Endocrinology 1981, 109, 993–995. [Google Scholar] [CrossRef] [PubMed]
  12. Shimizu-Albergine, M.; Tsai, L.-C.L.; Patrucco, E.; Beavo, J.A. cAMP-specific phosphodiesterases 8A and 8B, essential regulators of Leydig cell steroidogenesis. Mol. Pharmacol. 2012, 81, 556–566. [Google Scholar] [CrossRef] [PubMed]
  13. Xiao, B.; Sanders, M.J.; Underwood, E.; Heath, R.; Mayer, F.V.; Carmena, D.; Jing, C.; Walker, P.A.; Eccleston, J.F.; Haire, L.F.; et al. Structure of mammalian AMPK and its regulation by ADP. Nature 2011, 472, 230–233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Abdou, H.S.; Bergeron, F.; Tremblay, J.J. A cell-autonomous molecular cascade initiated by AMP-activated protein kinase represses steroidogenesis. Mol. Cell. Biol. 2014, 34, 4257–4271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Clark, B.J.; Wells, J.; King, S.R.; Stocco, D.M. The purification, cloning, and expression of a novel luteinizing hormone-induced mitochondrial protein in MA-10 mouse Leydig tumor cells. Characterization of the steroidogenic acute regulatory protein (StAR). J. Biol. Chem. 1994, 269, 28314–28322. [Google Scholar] [CrossRef]
  16. Abdou, H.S.; Villeneuve, G.; Tremblay, J.J. The calcium signaling pathway regulates Leydig cell steroidogenesis through a transcriptional cascade involving the nuclear receptor NR4A1 and the steroidogenic acute regulatory protein. Endocrinology 2013, 154, 511–520. [Google Scholar] [CrossRef] [Green Version]
  17. Dupuis, N.; Fafouri, A.; Bayot, A.; Kumar, M.; Lecharpentier, T.; Ball, G.; Edwards, D.; Bernard, V.; Dournaud, P.; Drunat, S.; et al. Dymeclin deficiency causes postnatal microcephaly, hypomyelination and reticulum-to-Golgi trafficking defects in mice and humans. Hum. Mol. Genet. 2015, 24, 2771–2783. [Google Scholar] [CrossRef]
  18. Dimitrov, A.; Paupe, V.; Gueudry, C.; Sibarita, J.-B.; Raposo, G.; Vielemeyer, O.; Gilbert, T.; Csaba, Z.; Attie-Bitach, T.; Cormier-Daire, V.; et al. The gene responsible for Dyggve-Melchior-Clausen syndrome encodes a novel peripheral membrane protein dynamically associated with the Golgi apparatus. Hum. Mol. Genet. 2009, 18, 440–453. [Google Scholar] [CrossRef] [Green Version]
  19. Osipovich, A.B.; Jennings, J.L.; Lin, Q.; Link, A.J.; Ruley, H.E. Dyggve–Melchior–Clausen syndrome: Chondrodysplasia resulting from defects in intracellular vesicle traffic. Proc. Natl. Acad. Sci. USA 2008, 105, 16171–16176. [Google Scholar] [CrossRef] [Green Version]
  20. Dalla Rosa, I.; Durigon, R.; Pearce, S.F.; Rorbach, J.; Hirst, E.M.A.; Vidoni, S.; Reyes, A.; Brea-Calvo, G.; Minczuk, M.; Woellhaf, M.W.; et al. MPV17L2 is required for ribosome assembly in mitochondria. Nucleic Acids Res. 2014, 42, 8500–8515. [Google Scholar] [CrossRef] [Green Version]
  21. Miller, W.L.; Bose, H.S. Early steps in steroidogenesis: Intracellular cholesterol trafficking. J. Lipid Res. 2011, 52, 2111–2135. [Google Scholar] [CrossRef] [Green Version]
  22. Dejgaard, S.; Presley, J. Interactions of lipid droplets with the intracellular transport machinery. Int. J. Mol. Sci. 2021, 22, 2776. [Google Scholar] [CrossRef] [PubMed]
  23. Kumar, S.; Kang, H.; Park, E.; Park, H.-S.; Lee, K. The expression of CKLFSF2B is regulated by GATA1 and CREB in the Leydig cells, which modulates testicular steroidogenesis. Biochim. Biophys. Acta Gene Regul. Mech. 2018, 1861, 1063–1075. [Google Scholar] [CrossRef] [PubMed]
  24. Orlando, U.; Cooke, M.; Maciel, F.C.; Papadopoulos, V.; Podestá, E.J.; Maloberti, P. Characterization of the mouse promoter region of the acyl-CoA synthetase 4 gene: Role of Sp1 and CREB. Mol. Cell. Endocrinol. 2013, 369, 15–26. [Google Scholar] [CrossRef]
  25. Manna, P.R.; Stocco, D.M. Crosstalk of CREB and Fos/Jun on a single cis-element: Transcriptional repression of the steroidogenic acute regulatory protein gene. J. Mol. Endocrinol. 2007, 39, 261–277. [Google Scholar] [CrossRef] [PubMed]
  26. Clem, B.F.; Hudson, E.A.; Clark, B.J. Cyclic adenosine 3′,5′-monophosphate (cAMP) enhances cAMP-responsive element binding (CREB) protein phosphorylation and phospho-CREB interaction with the mouse steroidogenic acute regulatory protein gene promoter. Endocrinology 2005, 146, 1348–1356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Manna, P.; Eubank, D.; Lalli, E.; Sassone-Corsi, P.; Stocco, D. Transcriptional regulation of the mouse steroidogenic acute regulatory protein gene by the cAMP response-element binding protein and steroidogenic factor 1. J. Mol. Endocrinol. 2003, 30, 381–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Manna, P.R.; Dyson, M.T.; Eubank, D.W.; Clark, B.J.; Lalli, E.; Sassone-Corsi, P.; Zeleznik, A.J.; Stocco, D.M. Regulation of Steroidogenesis and the Steroidogenic Acute Regulatory Protein by a Member of the cAMP Response-Element Binding Protein Family. Mol. Endocrinol. 2002, 16, 184–199. [Google Scholar] [CrossRef]
  29. Ofir, R.; Dwarki, V.J.; Rashid, D.; Verma, I.M. CREB represses transcription of fos promoter: Role of phosphorylation. Gene Expr. 1991, 1, 55–60. [Google Scholar]
  30. Walker, W.H.; Girardet, C.; Habener, J.F. Alternative exon splicing controls a translational switch from activator to repressor isoforms of transcription factor CREB during spermatogenesis. J. Biol. Chem. 1996, 271, 20145–21050. [Google Scholar] [CrossRef]
  31. Yi, W.-R.; Tu, M.-J.; Yu, A.-X.; Lin, J.; Yu, A.-M. Bioengineered miR-34a modulates mitochondrial inner membrane protein 17 like 2 (MPV17L2) expression toward the control of cancer cell mitochondrial functions. Bioengineered 2022, 13, 12489–12503. [Google Scholar] [CrossRef] [PubMed]
  32. Li, J.; Gao, L.; Chen, J.; Zhang, W.-W.; Zhang, X.-Y.; Wang, B.; Zhang, C.; Wang, Y.; Huang, Y.-C.; Wang, H.; et al. Mitochondrial ROS-mediated ribosome stalling and GCN2 activation are partially involved in 1-nitropyrene-induced steroidogenic inhibition in testes. Environ. Int. 2022, 167, 107393. [Google Scholar] [CrossRef] [PubMed]
  33. Kone, M.; Pullen, T.J.; Sun, G.; Ibberson, M.; Martinez-Sanchez, A.; Sayers, S.; Nguyen-Tu, M.S.; Kantor, C.; Swisa, A.; Dor, Y.; et al. LKB1 and AMPK differentially regulate pancreatic β-cell identity. FASEB J. 2014, 28, 4972–4985. [Google Scholar] [CrossRef] [PubMed]
  34. Golkowski, M.; Shimizu-Albergine, M.; Suh, H.W.; Beavo, J.A.; Ong, S.-E. Studying mechanisms of cAMP and cyclic nucleotide phosphodiesterase signaling in Leydig cell function with phosphoproteomics. Cell. Signal. 2016, 28, 764. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, Y.; Merrill, R.A.; Strack, S. A-kinase anchoring protein 1: Emerging roles in regulating mitochondrial form and function in health and disease. Cells 2020, 9, 298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Dyson, M.T.; Jones, J.K.; Kowalewski, M.P.; Manna, P.R.; Alonso, M.; Gottesman, M.E.; Stocco, U.M. Mitochondrial A-kinase anchoring protein 121 binds type II protein kinase A and enhances steroidogenic acute regulatory protein-mediated steroidogenesis in MA-10 mouse Leydig tumor cells. Biol. Reprod. 2008, 78, 267–277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Arakane, F.; King, S.R.; Du, Y.; Kallen, C.B.; Walsh, L.P.; Watari, H.; Stocco, D.M.; Strauss, J.F. Phosphorylation of steroidogenic acute regulatory protein (StAR) modulates its steroidogenic activity. J. Biol. Chem. 1997, 272, 32656–32662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Pryde, K.R.; Smith, H.L.; Chau, K.-Y.; Schapira, A.H. PINK1 disables the anti-fission machinery to segregate damaged mitochondria for mitophagy. J. Cell Biol. 2016, 213, 163–171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Lakkaraju, A.K.; Mary, C.; Scherrer, A.; Johnson, A.E.; Strub, K. SRP keeps polypeptides translocation-competent by slowing translation to match limiting ER-targeting sites. Cell 2008, 133, 440–451. [Google Scholar] [CrossRef] [Green Version]
  40. Lin, H.-Y.; Ko, C.-Y.; Kao, T.-J.; Yang, W.-B.; Tsai, Y.-T.; Chuang, J.-Y.; Hu, S.-L.; Yang, P.-Y.; Lo, W.-L.; Hsu, T.-I. CYP17A1 maintains the survival of glioblastomas by regulating SAR1-mediated endoplasmic reticulum health and redox homeostasis. Cancers 2019, 11, 1378. [Google Scholar] [CrossRef] [Green Version]
  41. Anton, K.; Sinclair, J.; Ohoka, A.; Kajita, M.; Ishikawa, S.; Benz, P.M.; Renne, T.; Balda, M.; Jorgensen, C.; Matter, K.; et al. PKA-regulated VASP phosphorylation promotes extrusion of transformed cells from the epithelium. J. Cell Sci. 2014, 127, 3425–3433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Thornton, C.; Bright, N.J.; Sastre, M.; Muckett, P.J.; Carling, D. AMP-activated protein kinase (AMPK) is a tau kinase, activated in response to amyloid β-peptide exposure. Biochem. J. 2011, 434, 503–512. [Google Scholar] [CrossRef] [Green Version]
  43. Alquezar, C.; Arya, S.; Kao, A.W. Tau Post-translational modifications: Dynamic transformers of Tau function, degradation, and aggregation. Front. Neurol. 2020, 11, 595532. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, X.; Macara, I.G. Par-3 controls tight junction assembly through the Rac exchange factor Tiam1. Nat. Cell Biol. 2005, 7, 262–269. [Google Scholar] [CrossRef]
  45. Zhang, L.; Li, J.; Young, L.H.; Caplan, M.J. AMP-activated protein kinase regulates the assembly of epithelial tight junctions. Proc. Natl. Acad. Sci. USA 2006, 103, 17272–17277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Lv, X.; Liu, C.; Wang, Z.; Sun, Y.; Xiong, Y.; Lei, Q.; Guan, K. PARD3 induces TAZ activation and cell growth by promotingLATS1 andPP1 interaction. EMBO Rep. 2015, 16, 975–985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Shin, J.H.; Lee, G.; Jeong, M.G.; Kim, H.K.; Won, H.Y.; Choi, Y.; Lee, J.; Nam, M.; Choi, C.S.; Hwang, G.; et al. Transcriptional coactivator with PDZ-binding motif suppresses the expression of steroidogenic enzymes by nuclear receptor 4 A1 in Leydig cells. FASEB J. 2020, 34, 5332–5347. [Google Scholar] [CrossRef] [Green Version]
  48. Horsley, V.; Pavlath, G.K. Nfat: Ubiquitous regulator of cell differentiation and adaptation. J. Cell Biol. 2002, 156, 771–774. [Google Scholar] [CrossRef] [Green Version]
  49. Rao, A.; Luo, C.; Hogan, P.G. Transcription factors of the NFAT family: Regulation and Function. Annu. Rev. Immunol. 1997, 15, 707–747. [Google Scholar] [CrossRef]
  50. Okamura, H.; Aramburu, J.; Garcia-Rodriguez, C.; Viola, J.; Raghavan, A.; Tahiliani, M.; Zhang, X.; Qin, J.; Hogan, P.G.; Rao, A. concerted dephosphorylation of the transcription factor NFAT1 induces a conformational switch that regulates transcriptional activity. Mol. Cell 2000, 6, 539–550. [Google Scholar] [CrossRef]
  51. Okamura, H.; Garcia-Rodriguez, C.; Martinson, H.; Qin, J.; Virshup, D.M.; Rao, A. A conserved docking motif for CK1 binding controls the nuclear localization of NFAT1. Mol. Cell. Biol. 2004, 24, 4184–4195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Radulescu, A.E.; Cleveland, D.W. NuMA after 30 years: The matrix revisited. Trends Cell Biol. 2010, 20, 214–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Akhshi, T.; Trimble, W.S. A non-canonical Hedgehog pathway initiates ciliogenesis and autophagy. J. Cell Biol. 2021, 220, e202004179. [Google Scholar] [CrossRef]
  54. Schaffer, B.E.; Levin, R.S.; Hertz, N.T.; Maures, T.J.; Schoof, M.L.; Hollstein, P.E.; Benayoun, B.A.; Banko, M.R.; Shaw, R.J.; Shokat, K.M.; et al. Identification of AMPK phosphorylation sites reveals a network of proteins involved in cell invasion and facilitates large-scale substrate prediction. Cell Metab. 2015, 22, 907–921. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Braun, B.C.; Okuyama, M.W.; Müller, K.; Dehnhard, M.; Jewgenow, K. Steroidogenic enzymes, their products and sex steroid receptors during testis development and spermatogenesis in the domestic cat (Felis catus). J. Steroid Biochem. Mol. Biol. 2018, 178, 135–149. [Google Scholar] [CrossRef]
  56. Ponikwicka-Tyszko, D.; Chrusciel, M.; Pulawska, K.; Bernaczyk, P.; Sztachelska, M.; Guo, P.; Li, X.; Toppari, J.; Huhtaniemi, I.T.; Wolczynski, S.; et al. Mifepristone treatment promotes testicular Leydig cell tumor progression in transgenic mice. Cancers 2020, 12, E3263. [Google Scholar] [CrossRef]
  57. Wendler, A.; Wehling, M. PGRMC2, a yet uncharacterized protein with potential as tumor suppressor, migration inhibitor, and regulator of cytochrome P450 enzyme activity. Steroids 2013, 78, 555–558. [Google Scholar] [CrossRef]
  58. Ascoli, M. Characterization of several clonal lines of cultured Leydig tumor cells: Gonadotropin receptors and steroidogenic responses. Endocrinology 1981, 108, 88–95. [Google Scholar] [CrossRef]
  59. Rahman, N.A.; Huhtaniemi, I.T. Testicular cell lines. Mol. Cell. Endocrinol. 2004, 228, 53–65. [Google Scholar] [CrossRef]
  60. Selvaraj, V.; Stocco, D.M.; Clark, B.J. Current knowledge on the acute regulation of steroidogenesis. Biol. Reprod. 2018, 99, 13–26. [Google Scholar] [CrossRef] [Green Version]
  61. King, S.R.; Ronen-Fuhrmann, T.; Timberg, R.; Clark, B.J.; Orly, J.; Stocco, D.M. Steroid production after in vitro transcription, translation, and mitochondrial processing of protein products of complementary deoxyribonucleic acid for steroidogenic acute regulatory protein. Endocrinology 1995, 136, 5165–5176. [Google Scholar] [CrossRef] [PubMed]
  62. Engeli, R.T.; Fürstenberger, C.; Kratschmar, D.V.; Odermatt, A. Currently available murine Leydig cell lines can be applied to study early steps of steroidogenesis but not testosterone synthesis. Heliyon 2018, 4, e00527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Mehanovic, S.; E Mendoza-Villarroel, R.; de Mattos, K.; Talbot, P.; Viger, R.S.; Tremblay, J.J. Identification of novel genes and pathways regulated by the orphan nuclear receptor COUP-TFII in mouse MA-10 Leydig cells. Biol. Reprod. 2021, 105, 1283–1306. [Google Scholar] [CrossRef]
  64. Hébert-Mercier, P.-O.; Bergeron, F.; Robert, N.M.; Mehanovic, S.; Pierre, K.J.; Mendoza-Villarroel, R.E.; de Mattos, K.; Brousseau, C.; Tremblay, J.J. Growth hormone-induced STAT5B regulates Star gene expression through a cooperation with cJUN in mouse MA-10 Leydig cells. Endocrinology 2022, 163, bqab267. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 12846 g001 550
Figure 1. Validation of treatments on MA-10 Leydig cells by assessing Star mRNA levels. MA-10 Leydig cells were treated with either DMSO (control, grey bars), Forskolin alone (Fsk, 10 µM, blue bars) or Forskolin and AICAR (Fsk+AICAR, 10 µM, and 1 mM, respectively, red bars) for 1 h. Total RNA was extracted and reverse-transcribed, and qPCR was performed to quantify Star mRNA levels. Rpl19 was used to normalize the data. The numbers on the x-axis refer to the three different samples used in the LC-MS/MS analysis. Results are displayed as the mean of three individual experiments, each performed in duplicate. For a given experiment, different letters indicate a statistically significant difference between groups (p < 0.05).
Figure 1. Validation of treatments on MA-10 Leydig cells by assessing Star mRNA levels. MA-10 Leydig cells were treated with either DMSO (control, grey bars), Forskolin alone (Fsk, 10 µM, blue bars) or Forskolin and AICAR (Fsk+AICAR, 10 µM, and 1 mM, respectively, red bars) for 1 h. Total RNA was extracted and reverse-transcribed, and qPCR was performed to quantify Star mRNA levels. Rpl19 was used to normalize the data. The numbers on the x-axis refer to the three different samples used in the LC-MS/MS analysis. Results are displayed as the mean of three individual experiments, each performed in duplicate. For a given experiment, different letters indicate a statistically significant difference between groups (p < 0.05).
Ijms 23 12846 g001
Ijms 23 12846 g002 550
Figure 2. Analysis of differentially phosphorylated proteins in MA-10 Leydig cells in response Forskolin. LC-MS/MS was performed on enriched phosphopeptides extracted and digested from MA-10 Leydig cells treated with either DMSO (control) or Forskolin (Fsk, 10 µm) for 1 h. Phosphorylated proteins were then quantified and compared between treatments. Results are the mean of three individual experiments, each performed in triplicate. (A) Volcano plot comparing the level of phosphorylated proteins detected by LC-MS/MS in control and Fsk-treated MA-10 Leydig cells. Each dot represents a distinct phosphoprotein. A Log2 fold change of +/−1.17 (corresponds to 2.25-fold) with a p-value of 0.05 was used and is represented by the blue dotted lines. In response to Fsk, blue dots depict a decrease in phosphorylation level (two proteins are identified in that group: MAPT and PARD3), while red dots correspond to an increase in phosphorylation level (four proteins are identified in that group: SRP14, NR4A1, RAF1, AKAP1). Grey dots correspond to proteins with no significant changes in phosphorylation levels. (B) Heatmap of protein phosphorylation levels for differentially phosphorylated proteins between treatment groups (control DMSO vs. Fsk). The name of the protein is indicated. The scale on the right represents expression levels (Log10). Biological Process (C), Pathway Analysis (D), Cellular Component (E), and Molecular Function (F) of the differentially phosphorylated proteins are shown.
Figure 2. Analysis of differentially phosphorylated proteins in MA-10 Leydig cells in response Forskolin. LC-MS/MS was performed on enriched phosphopeptides extracted and digested from MA-10 Leydig cells treated with either DMSO (control) or Forskolin (Fsk, 10 µm) for 1 h. Phosphorylated proteins were then quantified and compared between treatments. Results are the mean of three individual experiments, each performed in triplicate. (A) Volcano plot comparing the level of phosphorylated proteins detected by LC-MS/MS in control and Fsk-treated MA-10 Leydig cells. Each dot represents a distinct phosphoprotein. A Log2 fold change of +/−1.17 (corresponds to 2.25-fold) with a p-value of 0.05 was used and is represented by the blue dotted lines. In response to Fsk, blue dots depict a decrease in phosphorylation level (two proteins are identified in that group: MAPT and PARD3), while red dots correspond to an increase in phosphorylation level (four proteins are identified in that group: SRP14, NR4A1, RAF1, AKAP1). Grey dots correspond to proteins with no significant changes in phosphorylation levels. (B) Heatmap of protein phosphorylation levels for differentially phosphorylated proteins between treatment groups (control DMSO vs. Fsk). The name of the protein is indicated. The scale on the right represents expression levels (Log10). Biological Process (C), Pathway Analysis (D), Cellular Component (E), and Molecular Function (F) of the differentially phosphorylated proteins are shown.
Ijms 23 12846 g002
Ijms 23 12846 g003 550
Figure 3. Analysis of differentially phosphorylated proteins in MA-10 Leydig cells treated with Forskolin or Forskolin+AICAR. LC-MS/MS was performed on enriched phosphopeptides extracted and digested from MA-10 Leydig cells treated with Forskolin (Fsk, 10 µm) or Fsk+AICAR (10 µM and 1 mM, respectively) for 1 h. Phosphorylated proteins were then quantified and compared between treatments. Results are the mean of three individual experiments, each performed in triplicate. (A) Consensus sequence motif of phosphorylated peptides. The amino acid chosen is the one with the highest probability of being phosphorylated. (B) Volcano plot comparing the level of phosphorylated proteins detected by LC-MS/MS in Fsk-treated vs. Fsk+AICAR-treated MA-10 Leydig cells (to observe the effects of AICAR). Each dot represents a distinct phosphoprotein. A Log2 fold change of +/−1.17 (corresponds to 2.25-fold) with a p-value of 0.05 was used and is represented by the blue dotted lines. In response to AICAR, blue dots depict a decrease in phosphorylation level (one protein is identified in that group: NUMA1), while red dots correspond to an increase in phosphorylation level (three proteins are identified in that group: MAPT, PARD3, NFATC3). Grey dots correspond to proteins with no significant changes in phosphorylation levels. (C) Heatmap of protein phosphorylation levels for differentially phosphorylated proteins between treatments (Fsk vs. Fsk+AICAR). The name of the protein in indicated. The scale on the right represents expression levels (Log10). Biological Process (D), Pathway Analysis (E), Cellular Component (F), and Molecular Function (G) of the differentially phosphorylated proteins are shown.
Figure 3. Analysis of differentially phosphorylated proteins in MA-10 Leydig cells treated with Forskolin or Forskolin+AICAR. LC-MS/MS was performed on enriched phosphopeptides extracted and digested from MA-10 Leydig cells treated with Forskolin (Fsk, 10 µm) or Fsk+AICAR (10 µM and 1 mM, respectively) for 1 h. Phosphorylated proteins were then quantified and compared between treatments. Results are the mean of three individual experiments, each performed in triplicate. (A) Consensus sequence motif of phosphorylated peptides. The amino acid chosen is the one with the highest probability of being phosphorylated. (B) Volcano plot comparing the level of phosphorylated proteins detected by LC-MS/MS in Fsk-treated vs. Fsk+AICAR-treated MA-10 Leydig cells (to observe the effects of AICAR). Each dot represents a distinct phosphoprotein. A Log2 fold change of +/−1.17 (corresponds to 2.25-fold) with a p-value of 0.05 was used and is represented by the blue dotted lines. In response to AICAR, blue dots depict a decrease in phosphorylation level (one protein is identified in that group: NUMA1), while red dots correspond to an increase in phosphorylation level (three proteins are identified in that group: MAPT, PARD3, NFATC3). Grey dots correspond to proteins with no significant changes in phosphorylation levels. (C) Heatmap of protein phosphorylation levels for differentially phosphorylated proteins between treatments (Fsk vs. Fsk+AICAR). The name of the protein in indicated. The scale on the right represents expression levels (Log10). Biological Process (D), Pathway Analysis (E), Cellular Component (F), and Molecular Function (G) of the differentially phosphorylated proteins are shown.
Ijms 23 12846 g003
Ijms 23 12846 g004 550
Figure 4. Analysis of differentially phosphorylated proteins in MA-10 Leydig cells in response to Forskolin+AICAR. LC-MS/MS was performed on enriched phosphopeptides extracted and digested from MA-10 Leydig cells treated with DMSO (control) or Fsk+AICAR (10 µM and 1 mM, respectively) for 1 h. Phosphorylated proteins were then quantified and compared between treatments. Results are the mean of three individual experiments, each performed in triplicate. (A) Volcano plot comparing the level of phosphorylated proteins detected by LC-MS/MS in control (DMSO) vs. Fsk+AICAR-treated MA-10 Leydig cells. Each dot represents a distinct phosphoprotein. A Log2 fold change of +/−1.17 (corresponds to 2.25-fold) with a p-value of 0.05 was used and is represented by the blue dotted lines. In response to Fsk+AICAR, blue dots depict a decrease in phosphorylation level (one protein is identified in that group: MAP1A), while red dots correspond to an increase in phosphorylation level (five proteins are identified in that group: MAP1A, MAP1B, NR4A1, NUMA1, GOLGA5). Grey dots correspond to proteins with no significant changes in phosphorylation levels. (B) Heatmap of protein phosphorylation levels for differentially phosphorylated proteins between treatments (control DMSO vs. Fsk+AICAR). The name of the protein is indicated. The scale on the right represents expression levels (Log10). Biological Process (C), Pathway Analysis (D), Cellular Component (E) and Molecular Function (F) of the differentially phosphorylated proteins are shown.
Figure 4. Analysis of differentially phosphorylated proteins in MA-10 Leydig cells in response to Forskolin+AICAR. LC-MS/MS was performed on enriched phosphopeptides extracted and digested from MA-10 Leydig cells treated with DMSO (control) or Fsk+AICAR (10 µM and 1 mM, respectively) for 1 h. Phosphorylated proteins were then quantified and compared between treatments. Results are the mean of three individual experiments, each performed in triplicate. (A) Volcano plot comparing the level of phosphorylated proteins detected by LC-MS/MS in control (DMSO) vs. Fsk+AICAR-treated MA-10 Leydig cells. Each dot represents a distinct phosphoprotein. A Log2 fold change of +/−1.17 (corresponds to 2.25-fold) with a p-value of 0.05 was used and is represented by the blue dotted lines. In response to Fsk+AICAR, blue dots depict a decrease in phosphorylation level (one protein is identified in that group: MAP1A), while red dots correspond to an increase in phosphorylation level (five proteins are identified in that group: MAP1A, MAP1B, NR4A1, NUMA1, GOLGA5). Grey dots correspond to proteins with no significant changes in phosphorylation levels. (B) Heatmap of protein phosphorylation levels for differentially phosphorylated proteins between treatments (control DMSO vs. Fsk+AICAR). The name of the protein is indicated. The scale on the right represents expression levels (Log10). Biological Process (C), Pathway Analysis (D), Cellular Component (E) and Molecular Function (F) of the differentially phosphorylated proteins are shown.
Ijms 23 12846 g004
Ijms 23 12846 g005 550
Figure 5. Representative proteins that are differentially phosphorylated in MA-10 Leydig cells in response to Forskolin and to Forskolin+AICAR. Intensity plots of representative phosphoproteins (blue lines) and total proteins (red lines). Ctl, DMSO-treated cells; Fsk, Forskolin-treated cells; Fsk+AICAR, cells treated with Forskolin and AICAR. Results are shown as Fold Change compared to the control. Results are the mean of three individual experiments +/−SEM. The phosphopeptide sequence is shown beneath the graph. For a given protein, different letters indicate a statistically significant difference between treatment groups (p < 0.05).
Figure 5. Representative proteins that are differentially phosphorylated in MA-10 Leydig cells in response to Forskolin and to Forskolin+AICAR. Intensity plots of representative phosphoproteins (blue lines) and total proteins (red lines). Ctl, DMSO-treated cells; Fsk, Forskolin-treated cells; Fsk+AICAR, cells treated with Forskolin and AICAR. Results are shown as Fold Change compared to the control. Results are the mean of three individual experiments +/−SEM. The phosphopeptide sequence is shown beneath the graph. For a given protein, different letters indicate a statistically significant difference between treatment groups (p < 0.05).
Ijms 23 12846 g005
Table
Table 1. Variation in protein levels from MA-10 Leydig cells treated for 1 h with vehicle (DMSO), Fsk, or Fsk+AICAR.
Table 1. Variation in protein levels from MA-10 Leydig cells treated for 1 h with vehicle (DMSO), Fsk, or Fsk+AICAR.
DMSO vs. Fsk
Upregulated by Forskolin
Protein nameProteinGene nameFskDMSOFsk/DMSOq-ValueFsk/DMSO
DymeclinQ8CHY3Dym3.64 × 1071.57 × 10623.25.08 × 10-6 Up
Homologous-pairing protein 2 homologO35047Psmc3ip7.93 × 1061.57 × 1065.061.07 × 10-4 Up
Acyl-Coenzyme A dehydrogenase family member 12D3Z7X0Acad128.76 × 1062.22 × 1063.942.56 × 10-2 Up
N-acetyltransferase domain containing 1Q9DBW3Natd15.61 × 1061.57 × 1063.581.58 × 10-4 Up
Nuclear receptor subfamily 4 group A member 1P12813Nr4a13.76 × 1061.57 × 1062.41.41 × 10-3 Up
Downregulated by Forskolin
Protein nameProteinGene nameFskDMSODMSO/Fskq-ValueFsk/DMSO
H/ACA ribonucleoprotein complex subunit 3Q9CQS2Nop101.21 × 1067.91 × 1066.551.44 × 10-2Down
Mpv17-like protein 2Q8VIK2Mpv17l21.21 × 1066.46 × 1065.358.13 × 10-3Down
Overexpressed in colon carcinoma 1 protein homologP0C913Occ11.21 × 1066.36 × 1065.278.32 × 10-3Down
Glutathione peroxidase 7Q99LJ6Gpx71.21 × 1064.79 × 1063.973.60 × 10-2Down
Fsk vs. Fsk+AICAR
Upregulated by AICAR
Protein nameProteinGene nameFsk+AICARFskFsk+AICAR/Fskq-ValueFsk+AICAR/Fsk
Cyclic AMP-responsive element-binding protein 1Q01147Creb12.19 × 1071.21 × 10618.153.83 × 10-4 Up
Mpv17-like protein 2Q8VIK2Mpv17l26.32 × 1061.21 × 1065.235.72 × 10-3 Up
Overexpressed in colon carcinoma 1 protein homologP0C913Occ15.95 × 1061.21 × 1064.937.06 × 10-3 Up
HAUS augmin-like complex subunit 3Q8QZX2Haus33.98 × 1061.21 × 1063.33.10 × 10-2 Up
Downregulated by AICAR
Protein nameProteinGene nameFsk+AICARFskFsk/Fsk+AICARq-ValueFsk+AICAR/Fsk
DymeclinQ8CHY3Dym1.61 × 1063.64 × 10722.648.11x 10-6Down
Ubinuclein-2Q80WC1Ubn21.61 × 1061.56 × 1079.71.71 × 10-5Down
FerrochelataseP22315Fech1.61 × 1066.69 × 1064.164.02 × 10-4Down
Serine/threonine-protein kinase LMTK2Q3TYD6Lmtk21.61 × 1064.56 × 1062.846.67 × 10-3Down
Armadillo-like helical domain-containing protein 3Q6PD19Armh31.61 × 1064.52 × 1062.818.76 × 10-4Down
RibokinaseQ8R1Q9Rbks1.61 × 1064.05 × 1062.521.13 × 10-3Down
Nuclear receptor subfamily 4 group A member 1P12813Nr4a11.61 × 1063.76 × 1062.341.48 × 10-3Down
DMSO vs. Fsk+AICAR
Upregulated by Forskolin + AICAR
Protein nameProteinGene nameFsk+AICARDMSOFsk+AICAR/DMSOq-ValueFsk+AICAR/DMSO
Glutamyl-tRNA(Gln) amidotransferase subunit B. mitochondrialQ99JT1Gatb3.56 × 1061.57 × 1062.278.86 × 10-4 Up
Decapping and exoribonuclease proteinO70348Dxo2.89 × 1061.57 × 1061.843.70 × 10-2 Up
Downregulated by Forskolin + AICAR
Protein nameProteinGene nameFsk+AICARDMSODMSO/Fsk+AICARq-ValueFsk+AICAR/DMSO
SAGA-associated factor 29Q9DA08Sgf291.61 × 1061.39 × 109864.223.00 × 10-5Down
Coiled-coil domain-containing protein 187Q8C5V8Ccdc1871.61 × 1061.40 × 1078.713.87 × 10-4Down
FerrochelataseP22315Fech1.61 × 1066.73 × 1064.191.34 × 10-4Down
Armadillo-like helical domain-containing protein 3Q6PD19Armh31.61 × 1065.59 × 1063.483.67 × 10-2Down
Table
Table 2. Top 10 phosphopeptides differentially phosphorylated in MA-10 Leydig cells treated with vehicle (DMS), Fsk, or Fsk+AICAR (pair comparisons).
Table 2. Top 10 phosphopeptides differentially phosphorylated in MA-10 Leydig cells treated with vehicle (DMS), Fsk, or Fsk+AICAR (pair comparisons).
Fsk vs. DMSO
Phosphorylation decreased by Forskolin
Phosphopeptide SequenceProteinProtein NameGene
Name		FskDMSODMSO/
Fsk		q-ValueFsk/
DMSO
IQTSNVTNKNDPKQ64012RNA-binding protein RalyRaly2.12 × 1063.47 × 108163.64.54 × 10-4Down
DLPPFEDESEGLLGTEGPMEEEEDGEELIGDGMERP97310DNA replication licensing factor MCM2Mcm22.12 × 1065.46 × 10725.71.63 × 10-3Down
IDEPNTPYHNMIGDDEDAYSDSEGNEVMTPDILAKQ9DCL8Protein phosphatase inhibitor 2Ppp1r22.12 × 1065.14 × 10724.191.53 × 10-3Down
HLSNVSSTGSIDMVDSPQLATLADEVSASLAKP10637Microtubule-associated protein tauMapt2.12 × 1064.06 × 10719.141.53 × 10-3Down
SSSPSASLTEHEVSDSPGDEPSESPYESADETQTEASVSSKQ8BVE8Histone-lysine N-methyltransferase NSD2Whsc12.12 × 1063.51 × 10716.548.63 × 10-3Down
ETGGTYPPSPPPHSSPTPAATVAATVSTAVPGEPLLPRQ80U72Protein scribble homologScrib2.12 × 1062.36 × 10711.131.55 × 10-2Down
RISHSLYSGIEGLDESPTRQ99NH2Partitioning defective 3 homologPard32.12 × 1062.19 × 10710.332.37 × 10-2Down
KSMYSRVPECQVTTYYYVGFAYLMMRQ8QZY1Eukaryotic translation initiation factor 3 subunit LEif3l2.12 × 1062.12 × 1079.993.82 × 10-3Down
ALEETPPDSPAAEQENSVNCVDPLRQ8K1K3TERF1-interacting nuclear factor 2Tinf22.12 × 1062.03 × 1079.543.36 × 10-3Down
LAAQESSEAEDVTVDRQ6PGL7WASH complex subunit FAM21Fam212.12 × 1061.99 × 1079.364.73 × 10-2Down
Phosphorylation increased by Forskolin
Phosphopeptide SequenceProteinProtein NameGene
Name		FskDMSOFsk/
DMSO		q-ValueFsk/
DMSO
KSSVEGLEPAENKP16254Signal recognition particle 14 kDa proteinSrp142.63 × 1082.90 × 10690.64.54 × 10-4 Up
STSQGSINSPVYSRQ8K4G5Actin-binding LIM protein 1Ablim17.08 × 1072.12 × 10633.481.63 × 10-3 Up
QNPEQSADEDAEKNEEDSEGSSDEDEDEDGVGNTTFLKQ8R1B4Eukaryotic translation initiation factor 3 subunit CEif3c7.06 × 1072.12 × 10633.391.53 × 10-3 Up
PCIQAQYGTPATSPGPRP12813Nuclear receptor subfamily 4 group A member 1Nr4a14.52 × 1072.12 × 10621.361.53 × 10-3 Up
NTFTAWSEEDSDYEIDDRQ6ZQ58La-related protein 1Larp11.66 × 1081.02 × 10716.289.61 × 10-3 Up
EKVESAGPGGDSEPTGSTGALAHTPRO08550Histone-lysine N-methyltransferase 2BKmt2b2.42 × 1072.12 × 10611.443.17 × 10-3 Up
SPDLSNQNSDQANEEWETASESSDFASERQ7TSC1Protein PRRC2APrrc2a1.77 × 1072.12 × 1068.358.52 × 10-3 Up
KVSVEPQDSHQDAQPRQ8BXB6Solute carrier organic anion transporter family member 2B1Slco2b11.97 × 1082.47 × 1077.981.30 × 10-2 Up
KSSAAAAAAAAAEGALLPQTPPSPRQ9ERD6Ras-specific guanine nucleotide-releasing factor RalGPS2Ralgps21.64 × 1072.12 × 1067.779.61 × 10-3 Up
SHSESASPSALSSSPNNLSPTGWSQPKQ99N57RAF proto-oncogene serine/threonine-protein kinaseRaf11.61 × 1072.12 × 1067.631.30 × 10-2 Up
Fsk+AICAR vs. Fsk
Phosphorylation decreased by AICAR
Phosphopeptide SequenceProteinProtein NameGene
Name		Fsk + AICARFskFsk/
Fsk+AICAR		q-ValueFsk+AICAR/Fsk
YDEKTGLALQTEEFIPNYYCTDERRK7N6T2Vomeronasal 2 receptor 67Vmn2r672.15 × 1061.56 × 109727.231.81 × 10-4Down
GASAATGIPLESDEDSNDNDNDLENENCMHTNQ8R5H1Ubiquitin carboxyl-terminal hydrolase 15Usp152.15 × 1066.62 × 10730.832.50 × 10-3Down
KGSDDDGGDSPVQDIDTPEVDLYQLQVNTLRO88574Histone deacetylase complex subunit SAP30Sap302.15 × 1063.75 × 10717.431.12 × 10-2Down
WGQPPSPTPVPRPPDADPNTPSPKQ6PDQ2Chromodomain-helicase-DNA-binding protein 4Chd42.15 × 1063.33 × 10715.528.64 × 10-3Down
TNSMGSATGPLPGTKQ8BZ47Zinc finger protein 609Znf6092.15 × 1063.21 × 10714.944.99 × 10-3Down
NLATSADTPPSTIPGTGKQ8BKX6Serine/threonine-protein kinase SMG1Smg12.15 × 1061.63 × 1077.574.56 × 10-3Down
SPHDSKSPLDHRSPLERE9PZM4Chromodomain-helicase-DNA-binding protein 2Chd22.15 × 1061.58 × 1077.354.63 × 10-3Down
IASSSSENNFLSGSPSSPMGDILQTPQFQMRE9Q7G0Nuclear mitotic apparatus protein 1Numa12.15 × 1061.40 × 1076.521.14 × 10-2Down
NQENVSHLSVSSASPTSSVASAAGSVTSSSLQKO54826Protein AF-10Mllt102.15 × 1061.39 × 1076.452.00 × 10-2Down
GSSGEGLPFAEEGNLTIKQ8VHK1Caskin-2Caskin22.15 × 1061.23 × 1075.751.02 × 10-2Down
Phosphorylation increased by AICAR
Phosphopeptide SequenceProteinProtein NameGene
Name		Fsk + AICARFskFsk+AICAR/
Fsk		q-ValueFsk+AICAR/Fsk
ARPAQAPVSEELPPSPKPGKQ6PGL7WASH complex subunit FAM21Fam216.22 × 1072.12 × 10629.291.54 × 10-2 Up
HLSNVSSTGSIDMVDSPQLATLADEVSASLAKP10637Microtubule-associated protein tauMapt4.88 × 1072.12 × 10622.992.50 × 10-3 Up
APSPEPPTEEVAAETNSTPDDLEAQDALSPETTEEKQ9ESJ4NCK-interacting protein with SH3 domainNckipsd4.73 × 1072.12 × 10622.262.50 × 10-3 Up
SSSPSASLTEHEVSDSPGDEPSESPYESADETQTEASVSSKQ8BVE8Histone-lysine N-methyltransferase NSD2Whsc13.40 × 1072.12 × 10616.004.63 × 10-3 Up
VLLTHEVMCSRO08792Transcription factor COE2Ebf23.04 × 1072.12 × 10614.332.50 × 10-3 Up
VTETEDDSDSDSDDDEDDVHVTIGDIKQ9D824Pre-mRNA 3’-end-processing factor FIP1Fip1l12.40 × 1072.12 × 10611.312.50 × 10-3 Up
SRSVEDDEEGHLICQSGDVLSARP22518Dual specificity protein kinase CLK1Clk12.04 × 1072.12 × 1069.638.25 × 10-3 Up
RISHSLYSGIEGLDESPTRQ99NH2Partitioning defective 3 homologPard31.94 × 1072.12 × 1069.163.87 × 10-3 Up
DAEDLSPCLPSSSQEDTAVPSSPGPSDEVSNTEAEARG5E8P0Gamma-tubulin complex component 6Tubgcp61.92 × 1072.12 × 1069.063.47 × 10-2 Up
SVDLKTASPESGRSGFQDEESFRE9Q4Y4Centrosomal protein 192Cep1921.89 × 1072.12 × 1068.914.99 × 10-3 Up
Fsk+AICAR vs. DMSO
Phosphorylation decreased by Forskolin + AICAR
Phosphopeptide SequenceProteinProtein NameGene
Name		Fsk + AICARDMSODMSO/Fsk+AICARq-ValueFsk+AICAR/DMSO
AADPPAENSSAPEAEQGGAEP62960Nuclease-sensitive element-binding protein 1Ybx12.15 × 1061.60 × 10874.332.97 × 10-3Down
DLPPFEDESEGLLGTEGPMEEEEDGEELIGDGMERP97310DNA replication licensing factor MCM2Mcm22.15 × 1065.46 × 10725.42.27 × 10-3Down
SLTISVDSASTSRQ99MZ6Unconventional myosin-VIIbMyo7b2.15 × 1064.54 × 10721.123.65 × 10-2Down
LDGESDKEQFDDDQKQ8C1B1Calmodulin-regulated spectrin-associated protein 2Camsap22.15 × 1063.27 × 10715.236.41 × 10-3Down
ALEETPPDSPAAEQENSVNCVDPLRQ8K1K3TERF1-interacting nuclear factor 2Tinf22.15 × 1062.03 × 1079.434.79 × 10-3Down
QATESPAYGIPLKDGSEQTDEEAEGPFSDDEMVTHKQ99KK1Receptor expression-enhancing protein 3Reep32.15 × 1061.89 × 1078.82.18 × 10-2Down
HSPIKEEPCGSISETVCKQ6PDM1Male-specific lethal 1 homologMsl12.15 × 1061.88 × 1078.752.63 × 10-2Down
QDVETCRPSSPFGRQ6PG16Holliday junction recognition proteinHjurp2.15 × 1061.78 × 1078.271.96 × 10-2Down
SSPDPAVNPVPKQ924W7Suppression of tumorigenicity 5 proteinSt52.15 × 1061.73 × 1078.036.41 × 10-3Down
SEVNSEDSDIQEVLPVPKQ7TQI8Testis-specific Y-encoded-like protein 2Tspyl22.83 × 1062.20 × 1077.81.83 × 10-2Down
Phosphorylation increased by Forskolin + AICAR
SequenceProteinProtein NameGene
Name		Fsk + AICARDMSOFsk+AICAR/
DMSO		q-ValueFsk+AICAR/DMSO
KSSVEGLEPAENKP16254Signal recognition particle 14 kDa proteinSrp142.76 × 1082.90 × 10695.061.47 × 10-3 Up
TPLGASLDEQSSGTPKQ922B9Sperm-specific antigen 2 homologSsfa21.37 × 1082.35 × 10658.371.47 × 10-3 Up
TVSTQHSTESQDNDQPDYDSVASDEDTDVETRQ9JLQ2ARF GTPase-activating protein GIT2Git21.84 × 1083.32 × 10655.336.87 × 10-3 Up
QNPEQSADEDAEKNEEDSEGSSDEDEDEDGVGNTTFLKQ8R1B4Eukaryotic translation initiation factor 3 subunit CEif3c8.66 × 1072.12 × 10640.931.74 × 10-3 Up
STSQGSINSPVYSRQ8K4G5Actin-binding LIM protein 1Ablim18.20 × 1072.12 × 10638.773.41 × 10-3 Up
LKNDSDLFGLGLEEMGPKESSDEDRQ5U3K5Rab-like protein 6Rabl64.16 × 1072.12 × 10619.671.76 × 10-3 Up
NTFTAWSEEDSDYEIDDRQ6ZQ58La-related protein 1Larp11.88 × 1081.02 × 10718.396.41 × 10-3 Up
AELEEMEEVHPSDEEEEETKAESFYQKQ9QYR6Microtubule-associated protein 1AMap1a3.86 × 1072.12 × 10618.261.76 × 10-3 Up
QSNASSDVEVEEKETNVSKEDTDQEEKQ99JF8PC4 and SFRS1-interacting proteinPsip13.81 × 1072.12 × 10618.016.41 × 10-3 Up
PCIQAQYGTPATSPGPRP12813Nuclear receptor subfamily 4 group A member 1Nr4a12.46 × 1072.12 × 10611.622.39 × 10-2 Up
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Demmouche, Z.B.; Tremblay, J.J. Comprehensive and Quantitative Analysis of the Changes in Proteomic and Phosphoproteomic Profiles during Stimulation and Repression of Steroidogenesis in MA-10 Leydig Cells. Int. J. Mol. Sci. 2022, 23, 12846. https://doi.org/10.3390/ijms232112846

AMA Style

Demmouche ZB, Tremblay JJ. Comprehensive and Quantitative Analysis of the Changes in Proteomic and Phosphoproteomic Profiles during Stimulation and Repression of Steroidogenesis in MA-10 Leydig Cells. International Journal of Molecular Sciences. 2022; 23(21):12846. https://doi.org/10.3390/ijms232112846

Chicago/Turabian Style

Demmouche, Zoheir B., and Jacques J. Tremblay. 2022. "Comprehensive and Quantitative Analysis of the Changes in Proteomic and Phosphoproteomic Profiles during Stimulation and Repression of Steroidogenesis in MA-10 Leydig Cells" International Journal of Molecular Sciences 23, no. 21: 12846. https://doi.org/10.3390/ijms232112846

APA Style

Demmouche, Z. B., & Tremblay, J. J. (2022). Comprehensive and Quantitative Analysis of the Changes in Proteomic and Phosphoproteomic Profiles during Stimulation and Repression of Steroidogenesis in MA-10 Leydig Cells. International Journal of Molecular Sciences, 23(21), 12846. https://doi.org/10.3390/ijms232112846

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