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

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.


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]. 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 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.

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 concentrationand 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. 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.

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 concentrationand 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. (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.

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].  The original magnification is ×200. (F) H2DCFDA staining was performed to measure ROS in day-8 control-, 2-BPand 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.
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, 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.
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).

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. 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).

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.  (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.

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).
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).

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. ×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 twocell 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.

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 GW501516treated 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.

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).

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

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].

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.

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.

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.

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% CO 2 air for 15 min. Thereafter, heparin-treated sperm were diluted in an IVF medium (final density of 1 × 10 6 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% CO 2 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).

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).

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.

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) CO 2 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.

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).

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.

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.

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.

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.

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.