Next Article in Journal
Near-Infrared Markers based on Bacterial Phytochromes with Phycocyanobilin as a Chromophore
Next Article in Special Issue
Proteomics Analysis of Early Developmental Stages of Zebrafish Embryos
Previous Article in Journal
An Ancient Mutation in the TPH1 Gene is Consistent with the Changes in Mammalian Reproductive Rhythm
Previous Article in Special Issue
Expression Analysis of mRNA Decay of Maternal Genes during Bombyx mori Maternal-to-Zygotic Transition
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 23
10.3390/ijms20236066
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

The PPARδ Agonist GW501516 Improves Lipolytic/Lipogenic Balance through CPT1 and PEPCK during the Development of Pre-Implantation Bovine Embryos

by
Muhammad Idrees
1,
Lianguang Xu
1,
Marwa El Sheikh
1,
Tabinda Sidrat
1,
Seok-Hwan Song
1,
Myeong-Don Joo
1,
Kyeong-Lim Lee
2 and
Il-Keun Kong
1,2,3,*
1
Division of Applied Life Science (BK21 Plus), Gyeongsang National University, Jinju 52828, Gyeongnam Province, Korea
2
The King Kong Ltd., Jinju 52828, Korea
3
Institute of Agriculture and Life Science, Gyeongsang National University, Jinju 52828, Gyeongnam Province, Korea
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(23), 6066; https://doi.org/10.3390/ijms20236066
Submission received: 31 October 2019 / Revised: 20 November 2019 / Accepted: 29 November 2019 / Published: 2 December 2019
(This article belongs to the Special Issue Regulation of Gene Expression During Embryonic Development)
Browse Figures
Versions Notes

Abstract

:
The PPARs (peroxisome proliferator-activated receptors) play critical roles in the regulation of lipid and glucose metabolism. PPARδ, a member of the PPARs family, is associated with decreased susceptibility to ectopic lipid deposition and is implicated in the regulation of mitochondrial processes. The current study aimed to determine the role of PPARδ in fatty acid β-oxidation and its influence on PEPCK for the lipogenic/lipolytic balance during in vitro bovine oocyte maturation and embryo development. Activation of PPARδ by GW501516, but not 2-BP, was indicated by intact embryonic PEPCK (cytosolic) and CPT1 expression and the balance between free fatty acids and mitochondrial β-oxidation that reduced ROS and inhibited p-NF-κB nuclear localization. Genes involved in lipolysis, fatty acid oxidation, and apoptosis showed significant differences after the GW501516 treatment relative to the control- and 2-BP-treated embryos. GSK3787 reversed the PPARδ-induced effects by reducing PEPCK and CPT1 expression and the mitochondrial membrane potential, revealing the importance of PPARδ/PEPCK and PPARδ/CPT1 for controlling lipolysis during embryo development. In conclusion, GW501516-activated PPARδ maintained the correlation between lipolysis and lipogenesis by enhancing PEPCK and CPT1 to improve bovine embryo quality.

Graphical Abstract

1. Introduction

Lipids are a potent source of energy, so considerable attention has been paid to their metabolism during in vitro oocyte maturation and early embryonic development [1]. Lipid contents and lipid-derived free fatty acid regulation have proven to be of great importance to oocyte developmental competence and early embryo physiology [2]. In bovine oocytes, lipids are primary triglycerides of specific fatty acids stored in distinct droplet organelles that re-localize during oocyte maturation [1,3]. In vitro cultured embryos are frequently associated with mitochondrial dysfunction and high ROS levels due to elevated lipid contents [4]. Accumulation of cytoplasmic lipids is one of the major drawbacks of in vitro bovine embryo production, and low cryotolerance is associated with enhanced cytoplasmic lipids compared with that of in vivo derived counterparts [5,6].
The potential of PPARs (peroxisome proliferator-activated receptors) in the management of metabolic disorders has been recognized by many researchers owing to their positive roles in both lipid and glucose utilization [7]. The PPARs, ligand-activated transcription factors belonging to the nuclear receptor superfamily, mainly exist in three subtypes, PPARα, PPARγ, and PPARβ/δ, all of which are regulated by fatty acids and their derivatives [8]. PPARs heterodimerize with RXR (retinoid-X-receptor), bind to PPREs (peroxisome proliferator response elements), and form PPAR-RXR heterodimers, which stimulate the expression of genes involved in glucose and fatty acid metabolism, insulin signaling, cell differentiation, and cell proliferation [8,9]. PPARδ, a ubiquitously expressed member of the PPARs family, plays a vital role in the transcription of many proteins related to lipid homeostasis [10]. Its activation significantly upregulates lipolysis by activating AMPK and the transcription of fatty acid β-oxidation-related genes [11,12]. PPARδ also improves glyceroneogenesis by transcribing PEPCK (phosphoenolpyruvate carboxykinase), an enzyme well known for the regulation of TCA cycle flux [8,13,14]. On the other hand, PEPCK is highly associated with energy production and plays an important role in the homeostasis of lipids by re-esterification of free fatty acids to generate triglycerides [15,16,17]. Among PPARs, PPARγ and δ are strongly expressed in the granulosa and theca cells of rat oocytes, whereas their knockout or deletion exerts reproductive defects [18,19]. Deletion of PPARδ from mice oocytes reduced embryo development, cell proliferation, and implantation potential [20,21]. Prostacyclin plays a critical role in embryo development and implantation, and PPARδ activation is essential for prostacyclin-mediated blastocyst formation and hatching [21]. Activated PPARδ represents a novel therapeutic target to improve in vitro fertilization and enhanced implantation of cultured embryos [21]. Regulation/termination of PPARs in vivo or in vitro indicated that each receptor has physiological roles in gamete maturation and embryo development [22].
In this study, we sought to determine the effects of the PPARδ synthetic agonist GW501516 on lipid metabolism during in vitro bovine oocyte maturation and its consequences on ROS levels, mitochondrial functioning, and the rate of developed bovine embryos. We analyzed PPARδ/PEPCK and PPARδ/CPT1 in control-, 2-bromopalmitate (2-BP)-, and GW501516-treated embryos. The PPARδ-specific antagonist GSK3787 (4-chloro-N-(2-ethyl) benz amide) [23] was used to determine the mitochondrial membrane potential, the expression of PEPCK and CPT1, and blastocyst development and hatching.

2. Results

2.1. Dynamic Changes in PPARs and PEPCK Expression during In Vitro Oocyte Maturation and Embryo Development

To investigate the expression pattern of PPARs, we first analyzed the mRNA levels of PPAR genes (PPARα, β/δ, and γ) at different developmental stages of bovine embryos. The results revealed that all PPAR members were maternally expressed and showed their expression in COCs (cumulus oocyte complexes), MII oocytes, two-cell embryos, 3.5-day embryos, and day-8 blastocysts (Figure 1A). The PPAR genes were minimally expressed during the MII stage oocyte, but the expression began to be enhanced as the developmental stages proceeded from the two-cell stage to day-8 blastocysts [24]. Meanwhile, PPARδ showed reduced expression relative to other PPAR members at the 3.5-day embryonic stage (Figure 1B).
To confirm the PCR results and examine PPAR localization, we co-stained the development stage-dependent samples with PPARα and PPARγ antibodies for immunofluorescence (Figure 1C). The results showed the protein expressions of PPARα and PPARγ in GV oocytes as well as in surrounding cumulus cells. At the MII stage, the protein expressions of both the PPARs were cytoplasmic. From two-cell to day-8 blastocysts, the expression of PPARα and γ was highly enhanced and found all over the embryo as well as in the inner cell mass. Thereafter, we analyzed PPARδ and PEPCK stage-dependent protein expressions through immunofluorescence. To do this, we co-stained COCs, MII, and PN (Pronuclear) with PPARδ and PEPCK antibodies (Figure 1D). PPARδ and PEPCK were expressed in cumulus cells as well as in GV oocytes, but at the MII stage, the expression of both the proteins were reduced. Surprisingly, PPARδ showed nuclear localization (as oocytes are transcriptionally inactive), while PEPCK was cytoplasmic. At the pronuclear stage (zygote), both proteins showed cytoplasmic localization. The dynamic changes in the expression of PPARδ and PEPCK during maturation and fertilization compelled us to analyze these proteins in detail (Figure 1E). From the two-cell to the eight-cell embryonic stage, PEPCK and PPARδ were expressed in the nucleus as well as in the cytoplasm, but in 3.5-day embryos and day-8 blastocysts, the expression of PEPCK was completely cytoplasmic. PPARδ was expressed in both the cytoplasm and the nucleus of bovine embryonic cells.

2.2. Effect of 2-BP and GW501516 on PPARs, Lipid Metabolism, and ROS Level during GVBD Induction

To evaluate the effects of 2-BP and GW501516 on oocytes, we first examined the concentration- and stage-dependent responses of both the compounds during embryonic development (Figure 2A,B). After obtaining the effective concentrations, we analyzed the influence of these compounds on fatty acid β-oxidation in MII oocytes. To do this, we checked CPT1 protein expression through Western blotting, and significant enhancement in GW501516-treated oocytes was observed relative to the control and 2-BP groups (Figure 2C) [25,26,27]. Next we analyzed the effects of 2-BP and GW501516 on lipid contents in MII oocytes (matured). For this we performed Nile red staining in control-, 2-BP-, and GW501516-treated in vitro maturation (IVM) media (Figure 2D). The Nile red staining revealed high lipid contents in control and 2-BP samples relative to GW501516. High lipid contents enhance ROS levels, and this phenomenon was observed via H2DCFDA staining in MII oocytes (Figure 2E) [28]. The level of ROS was markedly reduced in GW501516-treated oocytes relative to 2-BP and control groups.

2.3. PPARδ/PEPCK Expression and Mitochondrial β-Oxidation in Bovine Embryos

To determine the roles of PPARδ in preimplantation embryonic development, we investigated the effect of 2-BP and GW501516 on the development of bovine embryos (Table S2). The obtained blastocysts were immunostained with PPARδ antibody, and the results showed significantly high nuclear localization in the 2-BP and GW501516 groups relative to the control (Figure 3A). After that, we analyzed the PEPCK expression, which was markedly reduced in the 2-BP group relative to the control and GW501516 groups (Figure 3B). Immunofluorescent expression of the mitochondrial β-oxidation surrogate marker CPT1 also showed downregulation in the 2-BP-treated group (Figure 3C). 2-BP significantly translocated PPARδ to the nucleus, but repressed PEPCK and CPT1 expression by inhibiting fatty acid conversion to triglycerides and mitochondrial β-oxidation [29,30].
To further explore the effects of GW501516 and 2-BP on lipolysis, we assessed the mRNA expression of AGTL and PLIN2 through RT-qPCR in bovine day-8 blastocysts (Figure 3D). The mRNA expressions of all the genes showed upregulation in the GW501516 group relative to the 2-BP and control groups. The expression of ATGL and PLIN2 was considerably reduced in the 2-BP group, which showed its effect on lipolysis, as perilipins interact with ATGL for the hydrolysis of triglycerides [31,32]. To examine the effect of the GW501516-activated PPARδ and PEPCK pathways on mitochondria and lipid contents in bovine blastocysts, we combined Mitotracker green and Nile red staining (Figure 3E). The mitochondrial distribution was markedly higher in the GW501516 and control groups compared to 2-BP-treated blastocysts, whereas the lipid contents were significantly enhanced in the 2-BP and control groups. Reduced mitochondrial activity enhanced ROS production, and this phenomenon was observed in 2-BP-treated blastocysts (Figure 3F).

2.4. Disturbance in Lipolytic/Lipogenic Balance Enhanced Apoptosis in Bovine Blastocysts

The subtle balance between lipolysis and lipogenesis is essential for metabolic homeostasis during oocyte maturation and embryo development [33]. Disturbance of lipolytic/lipogenic balance enhances ROS levels and cellular apoptosis by inducing p-NF-κB DNA binding [34]. To investigate 2-BP-induced disturbance in lipolytic/lipogenic balance, we examined p-NF-κB localization in control-, 2-BP- and GW501516-treated blastocysts through immunofluorescence. As shown in Figure 4A, the p-NF-κB nuclear localization was much higher in the 2-BP group compared to the other two groups. The nuclear localized p-NF-κB signals were observed all over the blastocysts and initiated apoptosis as shown in a TUNEL assay (Figure 4B). To further confirm that 2-BP disturbance enhanced apoptosis, we analyzed caspase-3, a major marker of apoptosis initiation. To do this, we combined Nile red staining with caspase-3 immunofluorescent staining, and the results showed that 2-BP significantly enhanced lipids and apoptosis in the bovine day-8 blastocysts (Figure 4C). The above findings indicate that PPARδ activation by GW501516 sustained a balance between lipolysis and lipogenesis due to enhanced PEPCK and CPT1 expression, while 2-BP reduced PEPCK and CPT1 expression and interrupted lipid metabolism and also deregulated mitochondrial activity.

2.5. PPARδ Inhibition Reduced Embryo Development and Hatching

We then sought to examine the effect of specific inhibitor (GSK3787) of PPARδ on PEPCK and CPT1 expression during bovine embryo development. GSK3787 antagonizes ligand-induced changes in PPARδ-dependent gene expression [23]. The dose-dependent responses of GSK3787, in terms of bovine embryo cleavage and blastocyst development, are shown in Figure 5A. To investigate the effects of PPARδ inhibition on oocyte maturation, we used aceto-orcein staining. The results showed the non-significant effect of GSK3787 on oocyte maturation relative to GW501516- and control-treated oocytes (Figure 5B). The difference in the rate of in vitro fertilization of MII oocytes was also nonsignificant in inhibitor-treated oocytes as compared to GW501516 and control groups (Figure 5C) [21]. After that, GSK3787 was added to the IVC media, and during the first 84 h of in vitro culture, the two-cell and 3.5-day embryo rates were significantly affected by inhibitor treatment (Figure 5D and Table 1). Furthermore, we analyzed the blastocyst hatching rate, which showed a marked reduction in the GSK3787 group (19.4%) compared to GW501516- (47.0%) and control-treated (35.9%) blastocysts (Figure 5E).

2.6. PPARδ Reversibility Affected Lipid Metabolism and Embryo Survival

We speculated that PPARδ inhibition would reduce lipid metabolism and mitochondrial β-oxidation. For this, we investigated PEPCK and CPT1 through immunofluorescence, and both had significantly lower protein expression with GSK3787 treatment relative to control and GW501516-treated blastocysts (Figure 6A). After that to further explore the link between PPARδ and lipolysis in bovine embryos, we assessed the mRNA expression of ATGL, LMF1, LMF2, and LPL genes in day-8 blastocysts through RT-qPCR. The results showed that all lipolysis-related genes were significantly reduced with PPARδ inhibition (Figure 6B). Previously, it was identified that PPARδ increased the 5’ promoter activity of the SIRT1 gene and enhanced mTOR through PEPCK regulation of central carbon metabolism [13,35,36]. Through Western blot analysis, we found that GW501516 significantly enhanced SIRT1 and p-mTOR expression relative to the control and inhibitor groups (Figure 6C). These results identified that PPARδ activation/inhibition has significant effects on lipid metabolism in bovine embryos.

2.7. PPARδ Effects on Mitochondria and Implantation Potential of Bovine Day-8 Blastocysts

Oocytes and early cleaved embryos are critically dependent on mitochondria for proper development from ovulation to the compacted morula stage [37]. Mitochondria are not only important for free fatty acid β-oxidation, but also an essential biomarker for blastocyst implantation potential [31,38]. Through JC-1 staining, we analyzed the active mitochondria (red), which were significantly lower in the GSK3787 group relative to the control and GW501516 groups in bovine day-8 blastocysts (Figure 7A). Furthermore, to confirm mitochondrial functioning, the BCL-2 and cytochrome c proteins were analyzed through immunofluorescence. The results showed that BCL-2 was significantly enhanced in GW501516 relative to the inhibitor-treated group, in which the expression of cytochrome c was markedly upregulated (Figure 7B). Thereafter, the blastocyst implantation potential was analyzed through an invasion assay, and the invasion area plus proliferation rate were significantly enhanced in GW501516 relative to the inhibitor and control groups (Figure 7C).

3. Discussion

Management of energy stores is critical during oocyte maturation and early embryo development, and PPARs function as critical regulators of lipid and fatty acid homeostasis [1,10]. Previously, it was reported that PPARδ activation considerably enhanced mouse embryo development, hatching, and live birth rates, but the effects of activated PPARδ on embryo lipid metabolism and mitochondrial functioning were not explored [25]. Here, we demonstrated that PPARδ activation by a synthetic ligand (GW501516) enhanced PEPCK and CPT1 protein expression in bovine embryos, and, as a result, the balance between lipolysis and fatty acid β-oxidation improved embryo quality. Our work collectively suggests that PPARδ activation facilitates bovine blastocyst development, hatching, and cryosurvival via enhanced CPT1 and PEPCK expression.
Lipid metabolism is recognized as crucial to oocyte maturation and embryo development. Many studies using a variety of techniques have identified that triglycerides are the major constituent of lipid contents in bovine oocytes and early embryos [32]. According to Ferguson and Leese, triglyceride contents in bovine oocytes decreased from 59 ng to 46 ng during maturation and 34 ng post fertilization [32,39]. This decrease in triglycerides is due to the breakdown (lipolysis) and endogenous formation of free fatty acids for mitochondrial β-oxidation to act as a source of energy [1,40]. Previously, it has been reported that the inhibition of lipolysis reduced oocyte maturation and embryo development [41]. Similarly, the inhibition of fatty acid oxidation, by inhibiting CPT1 with methyl palmoxirate, also markedly reduces bovine oocyte maturation and fertilization [40]. Many studies reported that fatty acids are ligands for PPARs to transcribe genes for mitochondrial fatty acid oxidation, and PPAR inhibition reduces fatty acid oxidation as well as lipolysis [8,42]. Furthermore, elevated concentrations of non-esterified fatty acids and glucose enhance lipid contents and reduce lipolysis and fatty acid β-oxidation [43,44]. Increased lipid contents were correlated with apoptosis in blastocysts [38]. High lipid content is one of the causes of reduced cryotolerance in in vitro developed embryos [45]. Despite the known deleterious effects of lipid contents on cryopreservation, these are an important energy source required by all cells [5,6]. Specifically, the metabolism of fatty acids through mitochondrial β-oxidation is one of the leading pathways of oocyte nuclear and cytoplasmic maturation and embryo development [46].
Until recently, the physiological functions of PPARδ remained elusive in bovine oocytes and embryos. The utilization of subtype-selective agonists and specific antagonists revealed that PPARδ played an important role in the metabolic adaptation of bovine embryos to in vitro conditions. In the current study, for the first time, all PPAR family proteins were detected in bovine GV oocytes, MII oocytes, and from two-cell to day-8 blastocysts (Figure 1). PPARδ, a member of the PPAR nuclear receptor family, plays an important role in lipolysis and free fatty acid β-oxidation, generated by lipolysis for the production of ATP [8]. Prostacyclin (PGI2) plays a vital role in embryo development and hatching, and PPARδ is essential for PGI2 stimulation of mouse embryo development [21]. Bovine embryos exposed to PGI2 also showed enhanced blastocyst development and hatching [47]. Previously, it has been stated that free fatty acids activate PPARδ and translocate it to the nucleus; as a result, lipolysis-related genes are transcribed [8]. We found that 2-bromo palmitate (nonmetabolized fatty acid) significantly enhanced PPARδ nuclear localization but did not reduce lipid contents, while GW501516-based PPARδ activation significantly reduced oocyte lipid contents and ROS levels (Figure 2) [42].
PPARδ transcribed PEPCK in the liver, the encoding protein of PCK1, and played an important role in free fatty acid reduction [48]. PEPCK is highly associated with energy production and is a key factor in the regulation of the TCA cycle [13,16]. PEPCK is also required for steroidogenesis in Leydig cells, and its expression in COCs is significantly correlated with successful embryo potential and pregnancy outcome [17,49]. An important role of PEPCK is the maintenance of triglyceride/free fatty acid flux and also glucose and lipid homeostasis [15]. Additionally, bovine embryos exposed to NEFA (non-esterified fatty acids) downregulate PEPCK expression [50]. Furthermore, inhibition of lipolysis with a lipase inhibitor reduces mRNA expression of PEPCK [51]. Previous work and our studies suggest that 2-BP significantly reduced PEPCK and CPT expression relative to GW501516 (Figure 3) [52]. High lipid contents and dysfunctional mitochondrial also enhance apoptosis, and this phenomenon was found in 2-BP-treated bovine blastocysts (Figure 4) [53]. In accordance with previous studies, PPARδ enhances and improves glycolysis, but through an unknown mechanism; here we found that PPARδ, through PEPCK, promotes cellular growth by activating mTOR via coordination of the regulation of central carbon metabolism (Figure 6) [13,36]. Our results are quite similar in that PPARδ inhibition (GSK3787) had no effect on oocyte maturation; however, in contrast to previous studies, we found reduced fertilization with PPARδ inhibition [25]. The reductions in embryo cleavage and blastocyst hatching are quite similar to previous work (Figure 5) [25]. Likewise, PPARδ inhibition significantly reduced mitochondrial activity, and we found that PPARδ inhibition not only reduced the mitochondrial membrane potential, but also the implantation potential of bovine embryos (Figure 7) [30,31].

4. Materials and Methods

All experiments were conducted with slaughterhouse-derived materials. The Gyeongsang National University Institute of Animal Care Committee approved all experiments, including surgical procedures (GNU-130902-A0059). All of the chemicals and reagents were obtained from Sigma-Aldrich (St. Louis, MO, USA), unless otherwise noted.

4.1. Experimental Design

4.1.1. Experiment 1 (PPARδ Activation)

GW501516 (cat. #SML1491), and 2-BP (cat. #238422) were separately added to the IVM and IVC media. To optimize the concentration of GW501516, the COCs were divided into 12 experimental groups and treated with 0 (control group), 0.1, 0.2, 0.3, 0.5, 1, 1.5, 2, 3, 4, 6, 8, and 10 µM GW501516 in IVM and IVC media. Based on the blastocyst development rate, the 1 µM concentration was selected for GW501516. Similarly, 2-BP in various concentrations (0, 1, 2, 4, 5, 6, 8, 10, and 20 µM) was applied to COCs and embryos in IVM and IVC media, and 5 µM was selected as the minimal effective concentration for further studies. Four biological replicates were performed to obtain the optimal working concentration. After selecting the effective concentrations, the protein expression levels of PPARα, PPARγ, PPARδ, CPT1, PEPCK, and p-NF-κB were analyzed in MII oocytes and day-8 blastocysts using Western blot and immunofluorescence. Lipid concentration and apoptotic signals were analyzed using Nile red staining and a TUNEL assay. Real-time RT-qPCR was performed to confirm the transcription levels of genes (ATGL and PLIN2) related to lipid metabolism.

4.1.2. Experiment 2 (PPARδ Inhibition)

The PPARδ-specific inhibitor GSK3787 (cat. #G7423) in different concentrations was added in IVM media to COCs and in IVC media to embryos (0 (control group), 1, 2, 5, 10, 15, and 20 µM). Ten micro-molar was selected as the minimal effective concentration after comparison with the control group. The parameters were oocyte maturation rate, fertilization rate, percentage of cleavage, developed blastocysts, and hatching percentages. Five biological replicates were performed to obtain the optimal working concentrations. After selecting the effective inhibitor concentration, it was compared with the GW501516 and control groups. The protein expression levels of PPARδ, CPT1, PEPCK, Sirt1, mTOR, BCL2, and cytochrome c were analyzed through Western blot and immunofluorescence. Real-time RT-qPCR was performed for genes (ATGL, LMF 1, LMF 2, and LPL) related to lipolysis. An in vitro invasion assay was also used to analyze the embryo quality.

4.2. Oocyte Collection

Bovine ovaries were collected at a local abattoir as previously described [54]. Three- to six-millimeter diameter COCs were aspirated using an 18-gauge disposable needle. TL-HEPES (10 mM HEPES (H-6147), 10 mM sodium lactate, 114 mM sodium chloride (S-5886), 2 mM sodium bicarbonate (S-5761), 0.5 mM magnesium chloride (M-2393), 3.2 mM potassium chloride (P-5405), 0.34 mM sodium phosphate monobasic (S-5011), 2 mM calcium chloride (C-7902), 1 µL/mL phenol red, 100 IU/mL penicillin, and 0.1 mg/mL streptomycin solution was used for the collection of COCs. After aspiration, the COCs were allowed to settle down as sediment in 15-mL conical tubes at 37 °C for 5 min. The COCs of ≥3 uniform layers were collected under a stereomicroscope, while expended or denuded oocytes were discarded and washed three times with TLH-PVA (P-8136).

4.3. In Vitro Maturation (IVM)

A four-well Nunc dish (Nunc, Roskilde, Denmark) was used for oocyte maturation. Approximately 50 COCs were placed into each well containing 700 μL of IVM medium (TCM199; Invitrogen Corp., Carlsbad, CA, USA) with 10% (v/v) FBS (fetal bovine serum) (Gibco BRL, Life Technologies, Grand Island, NY, USA; cat. #16000-044), 1 µg/mL oestradiol-17β, 10 µg/mL follicle-stimulating hormone, 0.6 mM cysteine, 10 ng/mL epidermal growth factor, and 0.2 mM sodium pyruvate (Gibco BRL, Life Technologies, Grand Island, NY, USA, cat. #11360-070). The oocytes selected were allowed to undergo in vitro maturation at 38.5 °C in a humidified atmosphere of 5% CO2 in air for 22–24 h [55].

4.4. In Vitro Fertilization and Culture

The IVM matured oocytes were fertilized in vitro with frozen–thawed bovine sperm, as described previously [55]. In brief, the sperm were thawed at 38.0 °C for 1 min and diluted in D-PBS. After dilution, the sperm were centrifuged at 750× g for 5 min at room temperature. Following centrifugation, sperm pellets were resuspended in 500 µL of heparin (20 μg/mL) in IVF media (Tyrode’s lactate solution supplemented with 6 mg/mL bovine serum albumin, 22 µg/mL sodium pyruvate, 100 IU/mL penicillin, and 0.1 mg/mL streptomycin). To facilitate capacitation, the sperm were incubated at 38.5 °C in a humidified atmosphere of 5% CO2 air for 15 min. Thereafter, heparin-treated sperm were diluted in an IVF medium (final density of 1 × 106 sperms/mL) and co-cultured with mature oocytes for 20 h. The presumed zygotes were thoroughly pipetted to remove cumulus cells and were washed and cultured in four-well dishes for eight days of embryonic development. The culturing conditions were 38.5 °C and 5% CO2 air, while the medium composition was SOF-BE1 [28] medium supplemented with 5 μg/mL insulin, 5 μg/mL transferrin, 5 ng/mL sodium selenite (Sigma cat. #11074547001), 4 mg/mL fatty-acid-free BSA, and 100 ng/mL epidermal growth factor (EGF).

4.5. Immunofluorescence

Immunofluorescence staining was performed as previously described [56]. In brief, bovine oocytes or blastocysts were fixed in 4% (v/v) paraformaldehyde prepared in 1 M phosphate-buffered saline (PBS) and preserved at 4 °C for a minimum of 30 min. On the staining day, oocytes or blastocysts were taken in four-well dishes and washed twice in PVA-PBS (0.3% PVA 1× PBS) for 10 min. Proteinase K solution was then added for 5 min to retrieve the antigen. Subsequently, the blastocysts were incubated for 30 min in blocking solution containing normal bovine or donkey serum and 0.1% Triton X-100 in PBS. Primary antibodies were applied, and the four-well dishes were kept at 4 °C overnight. The next day, the blastocysts or oocytes were washed twice with PVA-PBS for 10 min. After washing, secondary antibodies (FITC and TRITC conjugated, Santa Cruz Biotechnology, Dallas, TX, USA) were applied at room temperature for an additional 90 min. Blastocysts and oocytes were again washed thrice with PBS for 5 min. After that, oocytes or blastocysts were treated with 4′, 6′-diamidino-2-phenylindole (DAPI) at 10 µg/mL for 5 min to stain the nuclei, and fixed on slides. Images were captured with a confocal laser-scanning microscope (Fluoview FV 1000, Olympus, Tokyo, Japan). The relative integrated density of the signal and area were measured with the ImageJ analysis program (National Institutes of Health, Bethesda, MD, USA; https://imagej.nih.gov/ij).

4.6. TUNEL Assay

An in situ cell death detection kit, TMR red (St. Louis, MO, USA cat. #12156792910) was used to perform the TUNEL assay according to the manufacturer’s protocols. Briefly, paraformaldehyde-fixed embryos were washed twice with PVA-PBS [54]. After washing, the embryos were permeabilized (0.5% (v/v) Triton X-100 and 0.1% (w/v) sodium citrate) at room temperature for 30 min. Then, permeabilized embryos were washed with PVA-PBS and incubated in the dark with fluorescent-conjugated terminal deoxynucleotidyl transferase dUTP at 37 °C for 1 h. The stained embryos were again washed with PVP-PBS and incubated in DAPI at 10 µg/mL for 5 min. After washing with PVP-PBS for 5 min, the blastocysts were mounted onto a glass slide, and their nuclear configuration was analyzed. The number of cells per blastocyst was determined by counting DAPI-stained cells (blue), and TUNEL-positive cells were labeled bright red under an epifluorescence microscope (Olympus IX71, Olympus, Tokyo, Japan) equipped with a mercury lamp.

4.7. H2DCFDA Assay for ROS Detection

ROS was measured with 2, 7, dichlorodihydrofluorescein diacetate (H2DCFDA) (cat. #D6883) as previously described [57]. In brief, fresh blastocysts were incubated in PBS containing 10 µM H2DCFDA for 30 min in a humidified atmosphere of 5% (v/v) CO2 in air at 38.5 °C. When incubation was finished, the blastocysts were washed three times with PBS, mounted onto glass slides, and examined under an epifluorescence microscope (Olympus IX71) under 490 nm excitation and 525 nm emission filters.

4.8. Extraction of mRNA and cDNA Synthesis

Total mRNA was extracted as previously described [55]. In brief, mRNA was extracted from different biological replicates, with 10 COCs, 20 mature oocytes/two cells stage embryos, 10 embryos (day 3.5) and five blastocysts (day 8) per replicate, using a Dynabeads mRNA direct kit (Dynal AS, Oslo, Norway). In 100 μL of lysis buffer, oocytes or blastocysts were suspended and vortexed at room temperature for 2 min. The lysate was mixed with prewashed Dynabeads oligo (dT) (20 μL) and annealed by rotation at room temperature for 3 min. To remove the supernatant, a Dynal MPC magnetic particle concentrator (Dynal AS, Oslo, Norway) was used. Three hundred microliters of washing buffer A were used to wash the magnetic beads harboring the hybridized mRNA and oligo (dT). After that, 150 μL of washing buffer B was applied. To denature the secondary structures, bound mRNAs were resuspended in 8 μL of 10 mM Tris-HCl and heated at 65 °C for 5 min, followed by rapid quenching of the reaction on ice for 3 min. Superscript III reverse transcriptase was used for mRNA to reverse-transcribe into first-strand cDNA (Bio-Rad Laboratories Hercules, CA, USA cat. #1708891). The primers and PCR conditions for each gene are given below (Table S1).

4.9. Real-Time Polymerase Chain Reaction

Quantitative RT-PCR was performed using a CFX98 instrument (Bio-Rad Laboratories) as previously described [56]. A total of 10 μL of reaction mixture containing 3 μL of diluted cDNA, 0.2 mM bovine-specific primer, and 1× iQ SYBR Green Super mix (iQ SYBR Green Super mix kit, Bio-Rad Laboratories cat. #170-8882) was used. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) primer was used for normalization and to detect variation in the expression of this internal control gene in all cDNA samples. After confirming the nonsignificant difference in GAPDH among the samples, all transcripts were quantified using independent real-time PCR analyses. Quantitative PCR programing comprised pre denaturation at 95 °C for 3 min, followed by 44 cycles of 95 °C for 15 s, 57 °C for 20 s, and 72 °C for 30 s, and a final extension at 72 °C for 5 min. Amplification was followed by melting curve analysis using progressive denaturation, during which the temperature was increased from 65 °C to 95 °C at a rate of 0.2 °C per second. Quantitative analysis of the targeted genes was performed using the ΔΔC (t) method, and the results were reported as a relative expression compared with the calibrator after standardization of the transcript to the normal estimation of the endogenous control GAPDH. The coefficients of variation of the intra- and interassay variance were calculated according to the formula s.d./mean × 100, for all genes profiled by qPCR.

4.10. Protein Extraction and Western Blot Analysis

Mature oocytes (100 per extract) or day-8 blastocysts (20 per extract) were washed with PBS, dissolved thoroughly in pro-prep™ (iNtRON Biotechnology, Inc., Burlington, NJ, USA cat. #17081), sonicated to make cell lysates, and then centrifuged at 13,200 rpm at 4 °C for 25 min. The supernatant proteins were collected and quantified with a Bradford assay as previously described [54], with some minor modifications. Sodium dodecyl sulphate polyacrylamide gel (SDS-PAGE) fractionation was used to quantify equal amounts of protein (10 µg) using a Bio-Rad protein assay kit (cat. #5000002). SDS-PAGE-containing proteins were transferred to a PVDF (cat. #GE 10600023) membrane and blocked in 5% skim milk or 5% BSA for 1 h. Thereafter, the proteins were incubated with primary antibodies overnight at 4 °C. The next day, after washing, the membrane was incubated with secondary antibody at room temperature for 90 min. ECL (Pierce TM ECL Western Blotting Substrate, Thermo Fisher Scientific, Waltham, MA USA) detection reagent was used according to the manufacturer’s instructions to detect proteins. To detect the molecular weights of the proteins, prestained protein ladders (Abcam, Cambridge, Cambs, UK cat. #ab116029) were used. X-ray films (iNtRON Biotechnology, Inc., Burlington, NJ, USA) were scanned, and the ImageJ program was used to detect the optical densities of the bands.

4.11. Mitochondrial–Lipid Dual Staining

Dual staining of lipid and mitochondria was performed as previously described [58]. In brief, the live blastocyst mitochondria were stained with 2 µM Mitotracker Green FM (Invitrogen, Eugene, OR, USA) (diluted in PBS and 0.4% BSA) for 30 min at 38.7 °C. After that, the blastocysts were washed thrice in PVA-PBS and fixed in 4% formaldehyde for 24 h at 4 °C. After fixation, the oocytes or blastocysts were washed with PVA-PBS and stained with Nile red dye (Invitrogen, Molecular Probes) at 10 µM (diluted in PBS) applied at room temperature for 20 min. The samples were washed thrice in PVA-PBS and stained with DAPI (1:100 (v/v) in D-PBS) for 10 min. The oocytes or blastocysts were mounted on a glass slides without compression in a mounting medium and covered with coverslips. A confocal laser-scanning Olympus Fluoview FV1000 microscope (Olympus, Tokyo, Japan) was used to excite the green fluorescence excitation wavelength at 594 nm, and the emission was read at 608 nm for Mitotracker and lipophilic fluorescent dye NR at 485 nm. For image analysis, the intensities of green fluorescence (mitochondria) and red fluorescence (lipids) were measured using ImageJ software (version 1.50, National Institutes of Health, Bethesda, MD, USA; https://imagej.nih.gov/ij) after normalization through subtraction of the background intensity from each image of the experimental groups.

4.12. Antibodies

The following antibodies from Santa Cruz Biotechnology (Dallas, TX, USA) were used in this study: CPT1 (cat. #sc-393070), p-NF-κB (cat. #sc-271908), PEPCK (cat. #sc-166778), PPARα (cat. #sc-9000), BCL-2 (cat. #sc-783), Caspase-3 (cat. #sc-1225), Sirt 1 (cat. #sc-15404) mouse β-actin (cat. #sc-47778). While the other antibodies, PPARγ Abcam (cat. #ab45036), cytochrome C (cat. #ab110325), p-mTOR Abcam (cat. #ab84400) and PPARδ (LS bioscience cat. #LS-C437498, Seattle, WA, USA).

4.13. Statistical Analysis

A computer-based Sigma Gel System (SPSS Software, Inc., Chicago, IL, USA) was used for embryo development analysis. To analyze the density and integral optical density (IOD) of scanned X-ray films of Western blot and immunofluorescence images, GraphPad Prism 6 (GraphPad Software, San Diego, CA, USA) and ImageJ software were used. To determine the statistical significance (p-value), one-way ANOVA followed by Student’s t-test was used to analyze the data. The density values of the data are expressed as the mean ± SEM of three independent experiments. Significance: * p < 0.05, ** p < 0.01, and *** p < 0.001.

5. Conclusions

In conclusion, we noted that bovine embryonic PEPCK expression was regulated by PPARδ and played a significant role in the lipolytic/lipogenic balance. Bovine embryos treated with GW501516 had low ROS levels and apoptosis relative to control- and 2-BP-treated embryos. Moreover, we found an increased hatching rate of bovine embryos treated with GW501516, which was most likely due to PPARδ-regulated lipolysis. This approach may be of significant value for the development of more efficacious protocols for in vitro-generated bovine embryos.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/20/23/6066/s1.

Author Contributions

Designed research, M.I.; performed research, M.I. and L.X.; analyzed data, M.I., S.-H.S. and L.X.; provided reagents and helped with experiments, S.-H.S., M.-D.J., K.-L.L., M.E.S., and T.S.; wrote the paper, M.I. and L.X.; review and supervision, I.-K.K.

Funding

This work was supported by IPET through the Agri-Bio Industry Technology Development Program, funded by MAFRA (Grant number: 117029-3 and 315017-5). Muhammad Idrees, Lianguang Xu, Marwa El Sheikh, Tabinda Sidrat, Seok-Hwan Song, Myeong-Don Joo, and Kyeong-Lim Lee were supported by scholarships from the BK21 plus the Ministry of Education, Republic of Korea.

Acknowledgments

The authors are grateful to Ahmed Mesalam for revising the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PPARδPeroxisome Proliferator-activated receptor δ
PEPCKPhosphoenolpyruvate carboxykinase
CPT1Carnitine palmitoyltransferase I
p-NF-κBPhosphorylated-nuclear factor kappa light chain enhancer of activated B cells
2-BP2-Bromo Palmitate
GSK37874-chloro-N-(2-ethyl) Benz amide
GV oocyteGerminal Vesicle oocyte
ROSReactive Oxygen Species

References

  1. Dunning, K.R.; Russell, D.L.; Robker, R.L. Lipids and oocyte developmental competence: The role of fatty acids and beta-oxidation. Reproduction 2014, 148, 15–27. [Google Scholar] [CrossRef] [PubMed]
  2. Aardema, H.; Vos, P.L.; Lolicato, F.; Roelen, B.A.; Knijn, H.M.; Vaandrager, A.B.; Helms, J.B.; Gadella, B.M. Oleic acid prevents detrimental effects of saturated fatty acids on bovine oocyte developmental competence. Biol. Reprod. 2011, 85, 62–69. [Google Scholar] [CrossRef] [PubMed]
  3. Cetica, P.; Pintos, L.; Dalvit, G.; Beconi, M. Activity of key enzymes involved in glucose and triglyceride catabolism during bovine oocyte maturation in vitro. Reproduction 2002, 124, 675–681. [Google Scholar] [CrossRef] [PubMed]
  4. Prastowo, S.; Amin, A.; Rings, F.; Held, E.; Wondim, D.S.; Gad, A.; Neuhoff, C.; Tholen, E.; Looft, C.; Schellander, K.; et al. Fateful triad of reactive oxygen species, mitochondrial dysfunction and lipid accumulation is associated with expression outline of the AMP-activated protein kinase pathway in bovine blastocysts. Reprod. Fertil. Dev. 2016, 29, 890–905. [Google Scholar] [CrossRef] [PubMed]
  5. Abe, H.; Yamashita, S.; Satoh, T.; Hoshi, H. Accumulation of cytoplasmic lipid droplets in bovine embryos and cryotolerance of embryos developed in different culture systems using serum-free or serum-containing media. Mol. Reprod. Dev. 2002, 61, 57–66. [Google Scholar] [CrossRef] [PubMed]
  6. Schwarz, K.R.L.; Botigelli, R.C.; Del Collado, M.; de Castro, F.C.; Fernandes, H.; Paschoal, D.M.; Leal, C.L.V. Effects of fetal calf serum on cGMP pathway and oocyte lipid metabolism in vitro. Reprod. Fertil. Dev. 2017, 29, 1593–1601. [Google Scholar] [CrossRef]
  7. Luquet, S.; Gaudel, C.; Holst, D.; Lopez-Soriano, J.; Jehl-Pietri, C.; Fredenrich, A.; Grimaldi, P.A. Roles of PPAR delta in lipid absorption and metabolism: A new target for the treatment of type 2 diabetes. Biochim. Biophys. Acta 2005, 1740, 313–317. [Google Scholar] [CrossRef]
  8. Ravnskjaer, K.; Frigerio, F.; Boergesen, M.; Nielsen, T.; Maechler, P.; Mandrup, S. PPARdelta is a fatty acid sensor that enhances mitochondrial oxidation in insulin-secreting cells and protects against fatty acid-induced dysfunction. J. Lipid Res. 2010, 51, 1370–1379. [Google Scholar] [CrossRef]
  9. Liu, Y.; Colby, J.K.; Zuo, X.; Jaoude, J.; Wei, D.; Shureiqi, I. The Role of PPAR-delta in Metabolism, Inflammation, and Cancer: Many Characters of a Critical Transcription Factor. Int. J. Mol. Sci. 2018, 19, 3339. [Google Scholar] [CrossRef]
  10. De Lange, P.; Lombardi, A.; Silvestri, E.; Goglia, F.; Lanni, A.; Moreno, M. Peroxisome Proliferator-Activated Receptor Delta: A Conserved Director of Lipid Homeostasis through Regulation of the Oxidative Capacity of Muscle. PPAR Res. 2008, 2008, 172676. [Google Scholar] [CrossRef]
  11. Kramer, D.K.; Al-Khalili, L.; Guigas, B.; Leng, Y.; Garcia-Roves, P.M.; Krook, A. Role of AMP kinase and PPARdelta in the regulation of lipid and glucose metabolism in human skeletal muscle. J. Biol. Chem. 2007, 282, 19313–19320. [Google Scholar] [CrossRef] [PubMed]
  12. Rogers, N.H.; Perfield II, J.W.; Strissel, K.J.; Obin, M.S.; Greenberg, A.S. Loss of ovarian function in mice results in abrogated skeletal muscle PPARdelta and FoxO1-mediated gene expression. Biochem. Biophys. Res. Commun. 2010, 392, 1–3. [Google Scholar] [CrossRef] [PubMed]
  13. Montal, E.D.; Dewi, R.; Bhalla, K.; Ou, L.; Hwang, B.J.; Ropell, A.E.; Gordon, C.; Liu, W.J.; DeBerardinis, R.J.; Sudderth, J.; et al. PEPCK Coordinates the Regulation of Central Carbon Metabolism to Promote Cancer Cell Growth. Mol. Cell 2015, 60, 571–583. [Google Scholar] [CrossRef] [PubMed]
  14. Fan, W.; Waizenegger, W.; Lin, C.S.; Sorrentino, V.; He, M.X.; Wall, C.E.; Li, H.; Liddle, C.; Ruth, T.Y.; Atkins, A.R.; et al. PPARdelta Promotes Running Endurance by Preserving Glucose. Cell Metab. 2017, 25, 1186–1193. [Google Scholar] [CrossRef] [PubMed]
  15. Millward, C.A.; DeSantis, D.; Hsieh, C.W.; Heaney, J.D.; Pisano, S.; Olswang, Y.; Reshef, L.; Beidelschies, M.; Puchowicz, M.; Croniger, C.M. Phosphoenolpyruvate carboxykinase (Pck1) helps regulate the triglyceride/fatty acid cycle and development of insulin resistance in mice. J. Lipid Res. 2010, 51, 1452–1463. [Google Scholar] [CrossRef] [PubMed]
  16. Assou, S.; Boumela, I.; Haouzi, D.; Anahory, T.; Dechaud, H.; De Vos, J.; Hamamah, S. Dynamic changes in gene expression during human early embryo development: From fundamental aspects to clinical applications. Hum. Reprod. Update 2011, 17, 272–290. [Google Scholar] [CrossRef]
  17. Assou, S.; Haouzi, D.; Mahmoud, K.; Aouacheria, A.; Guillemin, Y.; Pantesco, V.; Reme, T.; Dechaud, H.; De Vos, J.; Hamamah, S. A non-invasive test for assessing embryo potential by gene expression profiles of human cumulus cells: A proof of concept study. Mol. Hum. Reprod. 2008, 14, 711–719. [Google Scholar] [CrossRef]
  18. Barak, Y.; Liao, D.; He, W.; Ong, E.S.; Nelson, M.C.; Olefsky, J.M.; Boland, R.; Evans, R.M. Effects of peroxisome proliferator-activated receptor delta on placentation, adiposity, and colorectal cancer. Proc. Natl. Acad. Sci. USA 2002, 99, 303–308. [Google Scholar] [CrossRef]
  19. Komar, C.M.; Braissant, O.; Wahli, W.; Curry, T.E., Jr. Expression and localization of PPARs in the rat ovary during follicular development and the periovulatory period. Endocrinology 2001, 142, 4831–4838. [Google Scholar] [CrossRef]
  20. MacLaren, L.A.; Guzeloglu, A.; Michel, F.; Thatcher, W.W. Peroxisome proliferator-activated receptor (PPAR) expression in cultured bovine endometrial cells and response to omega-3 fatty acid, growth hormone and agonist stimulation in relation to series 2 prostaglandin production. Domest. Anim. Endocrinol. 2006, 30, 155–169. [Google Scholar] [CrossRef]
  21. Huang, J.C.; Wun, W.S.; Goldsby, J.S.; Wun, I.C.; Noorhasan, D.; Wu, K.K. Stimulation of embryo hatching and implantation by prostacyclin and peroxisome proliferator-activated receptor delta activation: Implication in IVF. Hum. Reprod. 2007, 22, 807–814. [Google Scholar] [CrossRef]
  22. Huang, J.C. The role of peroxisome proliferator-activated receptors in the development and physiology of gametes and preimplantation embryos. PPAR Res. 2008, 2008, 732303. [Google Scholar] [CrossRef] [PubMed]
  23. Palkar, P.S.; Borland, M.G.; Naruhn, S.; Ferry, C.H.; Lee, C.; Sk, U.H.; Sharma, A.K.; Amin, S.; Murray, I.A.; Anderson, C.R.; et al. Cellular and pharmacological selectivity of the peroxisome proliferator-activated receptor-beta/delta antagonist GSK3787. Mol. Pharmacol. 2010, 78, 419–430. [Google Scholar] [CrossRef] [PubMed]
  24. Memili, E.; First, N.L. Control of gene expression at the onset of bovine embryonic development. Biol. Reprod. 1999, 61, 1198–1207. [Google Scholar] [CrossRef] [PubMed]
  25. Kang, H.J.; Kang, H.J.; Hwang, S.J.; Lim, H.; Song, H. Activation of peroxisome proliferators-activated receptor delta (PPARdelta) promotes blastocyst hatching in mice. Mol. Hum. Reprod. 2011, 17, 653–660. [Google Scholar] [CrossRef] [PubMed]
  26. Sudano, M.J.; Paschoal, D.M.; da Silva Rascado, T.; Crocomo, L.F.; Magalhães, L.C.O.; Junior, A.M.; Machado, R.; da Cruz Landim-Alvarenga, F. Crucial surviving aspects for vitrified in vitro-produced bovine embryos. Zygote 2014, 22, 124–131. [Google Scholar] [CrossRef]
  27. Wan, J.; Jiang, L.; Lü, Q.; Ke, L.; Li, X.; Tong, N. Activation of PPARdelta up-regulates fatty acid oxidation and energy uncoupling genes of mitochondria and reduces palmitate-induced apoptosis in pancreatic beta-cells. Biochem. Biophys. Res. Commun. 2010, 391, 1567–1572. [Google Scholar] [CrossRef]
  28. Somfai, T.; Kaneda, M.; Akagi, S.; Watanabe, S.; Haraguchi, S.; Mizutani, E.; Dang-Nguyen, T.Q.; Geshi, M.; Kikuchi, K.; Nagai, T. Enhancement of lipid metabolism with L-carnitine during in vitro maturation improves nuclear maturation and cleavage ability of follicular porcine oocytes. Reprod. Fertil. Dev. 2011, 23, 912–920. [Google Scholar] [CrossRef]
  29. Bickel, P.E.; Tansey, J.T.; Welte, M.A. PAT proteins, an ancient family of lipid droplet proteins that regulate cellular lipid stores. Biochim. Biophys. Acta 2009, 1791, 419–440. [Google Scholar] [CrossRef]
  30. Szkudelski, T.; Szkudelska, K. Short-term effects of palmitate and 2-bromopalmitate on the lipolytic activity of rat adipocytes. Life Sci. 2011, 89, 450–455. [Google Scholar] [CrossRef]
  31. Benkhalifa, M.; Ferreira, Y.J.; Chahine, H.; Louanjli, N.; Miron, P.; Merviel, P.; Copin, H. Mitochondria: Participation to infertility as source of energy and cause of senescence. Int. J. Biochem. Cell Biol. 2014, 55, 60–64. [Google Scholar] [CrossRef] [PubMed]
  32. Ferguson, E.M.; Leese, H.J. Triglyceride content of bovine oocytes and early embryos. J. Reprod. Fertil. 1999, 116, 373–378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jin, J.X.; Lee, S.; Taweechaipaisankul, A.; Kim, G.A.; Lee, B.C. Melatonin regulates lipid metabolism in porcine oocytes. J. Pineal. Res. 2017, 62, e12388. [Google Scholar] [CrossRef] [PubMed]
  34. Sadeghi, A.; Seyyed Ebrahimi, S.S.; Golestani, A.; Meshkani, R. Resveratrol Ameliorates Palmitate-Induced Inflammation in Skeletal Muscle Cells by Attenuating Oxidative Stress and JNK/NF-κB Pathway in a SIRT1-Independent Mechanism. J. Cell. Biochem. 2017, 118, 2654–2663. [Google Scholar] [CrossRef] [PubMed]
  35. Okazaki, M.; Iwasaki, Y.; Nishiyama, M.; Taguchi, T.; Tsugita, M.; Nakayama, S.; Kambayashi, M.; Hashimoto, K.; Terada, Y. PPARbeta/delta regulates the human SIRT1 gene transcription via Sp1. Endocr. J. 2010, 57, 403–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lee, C.H.; Olson, P.; Hevener, A.; Mehl, I.; Chong, L.W.; Olefsky, J.M.; Gonzalez, F.J.; Ham, J.; Kang, H.; Peters, J.M.; et al. PPARdelta regulates glucose metabolism and insulin sensitivity. Proc. Natl. Acad. Sci. USA 2006, 103, 3444–3449. [Google Scholar] [CrossRef] [Green Version]
  37. Bradley, J.; Swann, K. Mitochondria and lipid metabolism in mammalian oocytes and early embryos. Int. J. Dev. Biol. 2019, 63, 93–103. [Google Scholar] [CrossRef] [Green Version]
  38. Lolicato, F.; Brouwers, J.F.; de Lest, C.H.V.; Wubbolts, R.; Aardema, H.; Priore, P.; Roelen, B.A.; Helms, J.B.; Gadella, B.M. The cumulus cell layer protects the bovine maturing oocyte against fatty acid-induced lipotoxicity. Biol. Reprod. 2015, 92, 16. [Google Scholar] [CrossRef]
  39. Kim, J.Y.; Kinoshita, M.; Ohnishi, M.; Fukui, Y. Lipid and fatty acid analysis of fresh and frozen-thawed immature and in vitro matured bovine oocytes. Reproduction 2001, 122, 131–138. [Google Scholar] [CrossRef]
  40. Ferguson, E.M.; Leese, H.J. A potential role for triglyceride as an energy source during bovine oocyte maturation and early embryo development. Mol. Reprod. Dev. 2006, 73, 1195–1201. [Google Scholar] [CrossRef]
  41. Schwarz, K.R.L.; de Castro, F.C.; Schefer, L.; Botigelli, R.C.; Paschoal, D.M.; Fernandes, H.; Leal, C.L. The role of cGMP as a mediator of lipolysis in bovine oocytes and its effects on embryo development and cryopreservation. PLoS ONE 2018, 13, e0191023. [Google Scholar]
  42. Forman, B.M.; Chen, J.; Evans, R.M. Hypolipidemic drugs, polyunsaturated fatty acids, and eicosanoids are ligands for peroxisome proliferator-activated receptors alpha and delta. Proc. Natl. Acad. Sci. USA 1997, 94, 4312–4317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Ogawa, K.; Itami, N.; Ueda, M.; Kansaku, K.; Shirasuna, K.; Kuwayama, T.; Iwata, H. Non-esterified fatty acid-associated ability of follicular fluid to support porcine oocyte maturation and development. Reprod. Med. Biol. 2018, 17, 155–163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Uhde, K.; van Tol, H.T.; Stout, T.A.; Roelen, B.A. Exposure to elevated glucose concentrations alters the metabolomic profile of bovine blastocysts. PLoS ONE 2018, 13, e0199310. [Google Scholar] [CrossRef] [PubMed]
  45. Abe, T.; Kawahara-Miki, R.; Hara, T.; Noguchi, T.; Hayashi, T.; Shirasuna, K.; Kuwayama, T.; Iwata, H. Modification of mitochondrial function, cytoplasmic lipid content and cryosensitivity of bovine embryos by resveratrol. J. Reprod. Dev. 2017, 63, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Paczkowski, M.; Silva, E.; Schoolcraft, W.B.; Krisher, R.L. Comparative importance of fatty acid beta-oxidation to nuclear maturation, gene expression, and glucose metabolism in mouse, bovine, and porcine cumulus oocyte complexes. Biol. Reprod. 2013, 88, 111. [Google Scholar] [CrossRef]
  47. Song, B.S.; Kim, J.S.; Koo, D.B.; Park, J.S.; Lee, K.K.; Han, Y.M. 178 Improved Developmental Competence of Bovine Embryos by Pgi2 Analog Treatment. Reprod. Fertil. Dev. 2006, 18, 197. [Google Scholar] [CrossRef]
  48. Elo, B.; Villano, C.M.; Govorko, D.; White, L.A. Larval zebrafish as a model for glucose metabolism: Expression of phosphoenolpyruvate carboxykinase as a marker for exposure to anti-diabetic compounds. J. Mol. Endocrinol. 2007, 38, 433–440. [Google Scholar] [CrossRef] [Green Version]
  49. Ahn, S.W.; Gang, G.T.; Tadi, S.; Nedumaran, B.; Kim, Y.D.; Park, J.H.; Kweon, G.R.; Koo, S.H.; Lee, K.; Ahn, R.S.; et al. Phosphoenolpyruvate carboxykinase and glucose-6-phosphatase are required for steroidogenesis in testicular Leydig cells. J. Biol. Chem. 2012, 287, 41875–41887. [Google Scholar] [CrossRef] [Green Version]
  50. Desmet, K.L.; Van Hoeck, V.; Gagné, D.; Fournier, E.; Thakur, A.; O’doherty, A.M.; Walsh, C.P.; Sirard, M.A.; Bols, P.E.J.; Leroy, J.L.M.R. Exposure of bovine oocytes and embryos to elevated non-esterified fatty acid concentrations: Integration of epigenetic and transcriptomic signatures in resultant blastocysts. BMC Genom. 2016, 17, 1004. [Google Scholar] [CrossRef]
  51. Wan, Z.; Matravadia, S.; Holloway, G.P.; Wright, D.C. FAT/CD36 regulates PEPCK expression in adipose tissue. Am. J. Physiol. Cell Physiol. 2013, 304, 478–484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Yuan, Y.G.; Xu, L.; Zhang, S.; Mesalam, A.; Lee, K.L.; Liu, H.; Joo, M.D.; Idrees, M.; Kong, I.K. Polydatin and I-CBP112 protects early bovine embryo against nicotinamide-induced mitochondrial dysfunction. Theriogenology 2019, 134, 1–10. [Google Scholar] [CrossRef] [PubMed]
  53. El Sheikh, M.; Mesalam, A.; Mesalam, A.A.; Idrees, M.; Lee, K.L.; Kong, I.K. Melatonin Abrogates the Anti-Developmental Effect of the AKT Inhibitor SH6 in Bovine Oocytes and Embryos. Int. J. Mol. Sci. 2019, 20, 2956. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Idrees, M.; Xu, L.; Song, S.H.; Joo, M.D.; Lee, K.L.; Muhammad, T.; El Sheikh, M.; Sidrat, T.; Kong, I.K. PTPN11 (SHP2) Is Indispensable for Growth Factors and Cytokine Signal Transduction During Bovine Oocyte Maturation and Blastocyst Development. Cells 2019, 8, 1272. [Google Scholar] [CrossRef] [Green Version]
  55. Chowdhury, M.M.R.; Mesalam, A.; Khan, I.; Joo, M.D.; Lee, K.L.; Xu, L.; Afrin, F.; Kong, I.K. Improved developmental competence in embryos treated with lycopene during in vitro culture system. Mol. Reprod. Dev. 2018, 85, 46–61. [Google Scholar] [CrossRef]
  56. Gao, M.; Li, X.; He, Y.; Han, L.; Qiu, D.; Ling, L.; Liu, H.; Liu, J.; Gu, L. SIRT7 functions in redox homeostasis and cytoskeletal organization during oocyte maturation. FASEB J. 2018, 32, 6228–6238. [Google Scholar] [CrossRef]
  57. Milakovic, I.; Jeseta, M.; Hanulakova, S.; Knitlova, D.; Hanzalova, K.; Hulinska, P.; Machal, L.; Kempisty, B.; Antosik, P.; Machatkova, M. Energy Status Characteristics of Porcine Oocytes During In Vitro Maturation is Influenced by Their Meiotic Competence. Reprod. Domest. Anim. 2015, 50, 812–819. [Google Scholar] [CrossRef]
  58. Sommerfeld, V.; Niemann, H. Cryopreservation of bovine in vitro produced embryos using ethylene glycol in controlled freezing or vitrification. Cryobiology 1999, 38, 95–105. [Google Scholar] [CrossRef]
Ijms 20 06066 g001 550
Figure 1. PPARs’ mRNA and protein expression during maturation and embryo development. (A) RT-PCR-based PPARα, PPARγ, and PPARδ expression at the COCs (cumulus oocyte complexes) (10 per sample), MII oocyte (20 per sample), two-cell embryo (20 per sample), 3.5-day embryo (10 per sample) and day-8 blastocyst stages (five per sample). (B) Relative mRNA expressions of PPARα, PPARγ and PPARδ from the COCs stage to day-8 blastocysts. (C) Development stage-dependent immunofluorescent expression of PPARα and PPARγ (n = 15 per sample). (D) Immunofluorescent co-localization of PPARδ and PEPCK from COCs to day-8 blastocysts (n = 15 per sample). The experiments were repeated three times, and the data are shown here as mean ± S.E.M. N.S., not significant. Significance: * p < 0.05, and ** p < 0.01. The original magnification is ×200.
Figure 1. PPARs’ mRNA and protein expression during maturation and embryo development. (A) RT-PCR-based PPARα, PPARγ, and PPARδ expression at the COCs (cumulus oocyte complexes) (10 per sample), MII oocyte (20 per sample), two-cell embryo (20 per sample), 3.5-day embryo (10 per sample) and day-8 blastocyst stages (five per sample). (B) Relative mRNA expressions of PPARα, PPARγ and PPARδ from the COCs stage to day-8 blastocysts. (C) Development stage-dependent immunofluorescent expression of PPARα and PPARγ (n = 15 per sample). (D) Immunofluorescent co-localization of PPARδ and PEPCK from COCs to day-8 blastocysts (n = 15 per sample). The experiments were repeated three times, and the data are shown here as mean ± S.E.M. N.S., not significant. Significance: * p < 0.05, and ** p < 0.01. The original magnification is ×200.
Ijms 20 06066 g001
Ijms 20 06066 g002 550
Figure 2. Effects of GW501516 and 2-BP on oocyte maturation and embryo development. (A) and (B) Development stage-dependent responses of GW501516 and 2-bromopalmitate on embryo cleavage percentage and blastocyst developmental percentage. (C) Western blot expression of CPT1 in MII oocytes treated with control, 2-BP, and GW501516 in IVM media (n = 100 per group); the experiment was repeated three times. (D) Immunofluorescent expression of lipid droplets in COCs and after treatment with control, 2-BP, and GW501516 in IVM media (n = 20 per group) The original magnification is ×200. (E) H2DCFDA staining in MII oocytes (20 per group) matured in IVM media treated with control, 2BP, and GW501516. The original magnification is ×40. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 2. Effects of GW501516 and 2-BP on oocyte maturation and embryo development. (A) and (B) Development stage-dependent responses of GW501516 and 2-bromopalmitate on embryo cleavage percentage and blastocyst developmental percentage. (C) Western blot expression of CPT1 in MII oocytes treated with control, 2-BP, and GW501516 in IVM media (n = 100 per group); the experiment was repeated three times. (D) Immunofluorescent expression of lipid droplets in COCs and after treatment with control, 2-BP, and GW501516 in IVM media (n = 20 per group) The original magnification is ×200. (E) H2DCFDA staining in MII oocytes (20 per group) matured in IVM media treated with control, 2BP, and GW501516. The original magnification is ×40. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Ijms 20 06066 g002
Ijms 20 06066 g003 550
Figure 3. Activated PPARδ enhanced PEPCK and mitochondrial β-oxidation. (A) Immunofluorescent expression of PPARδ, showing nuclear localized PPARδ in day-8 blastocysts (20 per group). (B) Immunoreactivity of PEPCK in day-8 blastocysts treated with control, 2-BP, and GW501516 in IVC media (15 per group). (C) Immunofluorescent expression of CPT1 in day-8 blastocysts in control, 2-BP, and GW501516 groups (15 per group). (D) ATGL and PLIN2 were examined through RT-qPCR in day-8 blastocysts (n = 5). (E) Mitotracker (Green) and Nile red dyes (Red) were examined in control-, 2-BP-, and GW501516-treated day-8 blastocysts. The original magnification is ×200. (F) H2DCFDA staining was performed to measure ROS in day-8 control-, 2-BP- and GW501516-treated blastocysts (20 per group). ImageJ software was used to quantify the signal intensity of the immunofluorescence images. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. N.S., not significant. * p < 0.05, ** p < 0.01, and *** p < 0.001. The original magnification is ×100.
Figure 3. Activated PPARδ enhanced PEPCK and mitochondrial β-oxidation. (A) Immunofluorescent expression of PPARδ, showing nuclear localized PPARδ in day-8 blastocysts (20 per group). (B) Immunoreactivity of PEPCK in day-8 blastocysts treated with control, 2-BP, and GW501516 in IVC media (15 per group). (C) Immunofluorescent expression of CPT1 in day-8 blastocysts in control, 2-BP, and GW501516 groups (15 per group). (D) ATGL and PLIN2 were examined through RT-qPCR in day-8 blastocysts (n = 5). (E) Mitotracker (Green) and Nile red dyes (Red) were examined in control-, 2-BP-, and GW501516-treated day-8 blastocysts. The original magnification is ×200. (F) H2DCFDA staining was performed to measure ROS in day-8 control-, 2-BP- and GW501516-treated blastocysts (20 per group). ImageJ software was used to quantify the signal intensity of the immunofluorescence images. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. N.S., not significant. * p < 0.05, ** p < 0.01, and *** p < 0.001. The original magnification is ×100.
Ijms 20 06066 g003
Ijms 20 06066 g004 550
Figure 4. Lipolytic/lipogenic disturbance enhanced apoptosis and reduced embryo quality. (A) Immunofluorescent analysis of p-NF-κB expression in control-, 2-BP-, and GW501516-treated bovine day-8 blastocysts. (B) A TUNEL assay was performed to examine apoptosis in blastocysts, and red signals show the apoptosis-positive cells (20 per group). The original magnification is ×100. (C) Immunofluorescent expression of Caspase 3 (Green) and Nile red (Red) in day-8 blastocysts (15 per group). ImageJ software was used to quantify the signal intensity of immunofluorescent images. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. * p < 0.05, ** p < 0.01, and *** p < 0.001. The original magnification is ×200.
Figure 4. Lipolytic/lipogenic disturbance enhanced apoptosis and reduced embryo quality. (A) Immunofluorescent analysis of p-NF-κB expression in control-, 2-BP-, and GW501516-treated bovine day-8 blastocysts. (B) A TUNEL assay was performed to examine apoptosis in blastocysts, and red signals show the apoptosis-positive cells (20 per group). The original magnification is ×100. (C) Immunofluorescent expression of Caspase 3 (Green) and Nile red (Red) in day-8 blastocysts (15 per group). ImageJ software was used to quantify the signal intensity of immunofluorescent images. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. * p < 0.05, ** p < 0.01, and *** p < 0.001. The original magnification is ×200.
Ijms 20 06066 g004
Ijms 20 06066 g005 550
Figure 5. PPARδ inhibition reduced embryo development and hatching. (A) Dose-dependent effect of GSK3787 (PPARδ inhibitor) on bovine embryo cleavage percentage and blastocyst development percentage. (B) Aceto-orcein staining was performed in MII oocytes to analyze the percent of maturation in control-, GSK3787-, and GW501516 IVM-treated oocytes. The original magnification is ×40. (C) DAPI staining was performed for fertilized oocytes (pronuclear PN), and the percent of fertilized zygotes is presented as a histogram. The original magnification is ×200. (D) Percent of two-cell embryos from five independent experiments; the data presented here are the mean ± S.E.M. The experiments were repeated three times, and the data are shown here as the mean ± S.E.M. N.S., not significant. * p < 0.05.
Figure 5. PPARδ inhibition reduced embryo development and hatching. (A) Dose-dependent effect of GSK3787 (PPARδ inhibitor) on bovine embryo cleavage percentage and blastocyst development percentage. (B) Aceto-orcein staining was performed in MII oocytes to analyze the percent of maturation in control-, GSK3787-, and GW501516 IVM-treated oocytes. The original magnification is ×40. (C) DAPI staining was performed for fertilized oocytes (pronuclear PN), and the percent of fertilized zygotes is presented as a histogram. The original magnification is ×200. (D) Percent of two-cell embryos from five independent experiments; the data presented here are the mean ± S.E.M. The experiments were repeated three times, and the data are shown here as the mean ± S.E.M. N.S., not significant. * p < 0.05.
Ijms 20 06066 g005
Ijms 20 06066 g006 550
Figure 6. PPARδ reversibility and lipid metabolism in bovine embryos. (A) Immunofluorescent co-localization of PEPCK and CPT1 in control-, GSK3787-, and GW501516 IVC-treated day-8 blastocysts (15 per group). The original magnification is ×200. (B) RT-qPCR-based mRNA quantification of ATGL, LMF1, LMF2, and LPL in day-8 blastocysts (five per sample). (C) Western blot analysis of SIRT1 and p-mTOR in control-, GSK3787-, and GW501516 IVC-treated blastocysts (20 per group). The experiments were repeated three times, and the data are shown here as mean ± S.E.M. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Figure 6. PPARδ reversibility and lipid metabolism in bovine embryos. (A) Immunofluorescent co-localization of PEPCK and CPT1 in control-, GSK3787-, and GW501516 IVC-treated day-8 blastocysts (15 per group). The original magnification is ×200. (B) RT-qPCR-based mRNA quantification of ATGL, LMF1, LMF2, and LPL in day-8 blastocysts (five per sample). (C) Western blot analysis of SIRT1 and p-mTOR in control-, GSK3787-, and GW501516 IVC-treated blastocysts (20 per group). The experiments were repeated three times, and the data are shown here as mean ± S.E.M. * p < 0.05, ** p < 0.01, and *** p < 0.001.
Ijms 20 06066 g006
Ijms 20 06066 g007 550
Figure 7. PPARδ effect on mitochondria and implantation potential of bovine blastocysts. (A) JC-1 staining in day-8 blastocysts from the control-, GSK3787-, and GW501516-treated groups (20 per group). (B) Blastocysts were co-stained with BCL-2 (red) and Cytochrome c (green) for immunofluorescence to analyze mitochondrial activity (15 per group). The original magnification is ×200. (C) The effects of PPARδ inhibition on blastocyst implantation potential were determined. Bright-field image showing the area of invasion and DAPI for migrant cells in day-8 blastocysts. The original magnification is ×40. ImageJ software was used to quantify the signal intensity of immunofluorescent images. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. N.S., not significant. * p < 0.05, and ** p < 0.01.
Figure 7. PPARδ effect on mitochondria and implantation potential of bovine blastocysts. (A) JC-1 staining in day-8 blastocysts from the control-, GSK3787-, and GW501516-treated groups (20 per group). (B) Blastocysts were co-stained with BCL-2 (red) and Cytochrome c (green) for immunofluorescence to analyze mitochondrial activity (15 per group). The original magnification is ×200. (C) The effects of PPARδ inhibition on blastocyst implantation potential were determined. Bright-field image showing the area of invasion and DAPI for migrant cells in day-8 blastocysts. The original magnification is ×40. ImageJ software was used to quantify the signal intensity of immunofluorescent images. The experiments were repeated three times, and the data are shown here as mean ± S.E.M. N.S., not significant. * p < 0.05, and ** p < 0.01.
Ijms 20 06066 g007
Table
Table 1. Cleavage and development percentage of bovine embryos, control, GW501516 (PPARδ activator), and GSK3787 (PPARδ inhibitor).
Table 1. Cleavage and development percentage of bovine embryos, control, GW501516 (PPARδ activator), and GSK3787 (PPARδ inhibitor).
GroupsNo. of Fertilized ZygotesNo. of Cleavage Embryo
(% ± SEM)		No. of Blastocysts
(% ± SEM)		No. of Hatched Blastocysts
(% ± SEM)
Control399306 (76.6 ± 1.2) b133 (33.2 ± 0.9) b47 (35.9 ± 2.8) b
GW501516399323 (81.1 ± 1.1) b154 (38.8 ± 1.6) c72 (47.0 ± 2.5) c
GSK3787375266 (70.8 ± 1.6) a102 (27.3 ± 1.1) a20 (19.4 ± 2.3) a
a,b,c Values with different superscripts in the same column are significantly different (p < 0.05). This experiment was completed in eight replicates.

Share and Cite

MDPI and ACS Style

Idrees, M.; Xu, L.; El Sheikh, M.; Sidrat, T.; Song, S.-H.; Joo, M.-D.; Lee, K.-L.; Kong, I.-K. The PPARδ Agonist GW501516 Improves Lipolytic/Lipogenic Balance through CPT1 and PEPCK during the Development of Pre-Implantation Bovine Embryos. Int. J. Mol. Sci. 2019, 20, 6066. https://doi.org/10.3390/ijms20236066

AMA Style

Idrees M, Xu L, El Sheikh M, Sidrat T, Song S-H, Joo M-D, Lee K-L, Kong I-K. The PPARδ Agonist GW501516 Improves Lipolytic/Lipogenic Balance through CPT1 and PEPCK during the Development of Pre-Implantation Bovine Embryos. International Journal of Molecular Sciences. 2019; 20(23):6066. https://doi.org/10.3390/ijms20236066

Chicago/Turabian Style

Idrees, Muhammad, Lianguang Xu, Marwa El Sheikh, Tabinda Sidrat, Seok-Hwan Song, Myeong-Don Joo, Kyeong-Lim Lee, and Il-Keun Kong. 2019. "The PPARδ Agonist GW501516 Improves Lipolytic/Lipogenic Balance through CPT1 and PEPCK during the Development of Pre-Implantation Bovine Embryos" International Journal of Molecular Sciences 20, no. 23: 6066. https://doi.org/10.3390/ijms20236066

APA Style

Idrees, M., Xu, L., El Sheikh, M., Sidrat, T., Song, S. -H., Joo, M. -D., Lee, K. -L., & Kong, I. -K. (2019). The PPARδ Agonist GW501516 Improves Lipolytic/Lipogenic Balance through CPT1 and PEPCK during the Development of Pre-Implantation Bovine Embryos. International Journal of Molecular Sciences, 20(23), 6066. https://doi.org/10.3390/ijms20236066

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