Skip to main content

Cocultured porcine granulosa cells respond to excess non-esterified fatty acids during in vitro maturation



Non-esterified fatty acids (NEFAs) are one of the main lipid components of follicular fluid at concentrations that depend on circulating levels. Elevated levels of NEFAs impair oocyte quality, development potential, and may subsequently influence the metabolism and reproductive fitness of offspring. Granulosa cells (GCs) are the follicular cells that are closely communicating with the oocyte. However, the responses of GCs exposed to high levels of NEFAs when cocultured with cumulus-oocyte complexes (COCs), and how they attenuate the negative effects of NEFAs on oocytes, are unclear.


To better understand this protective effect, monolayers of porcine GCs were cocultured with COCs during in vitro maturation (IVM) in the presence of elevated levels of NEFAs. Genomic expression analysis was conducted to explore the responses of the GCs to the elevated levels of NEFAs. After limma algorithm analysis, 1,013 genes were differentially expressed between GCs cultured with and without elevated NEFAs. Among them, 438 genes were upregulated and 575 were downregulated. The differentially expressed genes were enriched in pathways related to metabolism, inflammation, and epithelial-mesenchymal transition.


The pathways and upstream regulators suggested that the cocultured GCs responded to the elevated NEFAs with (1) inhibition of the transition from granulosa to luteal cell, (2) interactions of metabolism change, anti-inflammation, mitochondrial function, and cell transition, (3) intercommunication with cocultured COCs of anti-inflammatory factors.


In mammals, especially in bovine and humans, the contexts of negative energy balance and obesity result in elevated levels of circulating non-esterified fatty acids (NEFAs) and may affect fertility [1, 2]. The elevated NEFAs in circulation originate mainly from lipolysis in adipocytes and from the diet [3]. In vitro studies showed that increased NEFA levels can affect follicular growth through impaired viability and steroidogenesis of granulosa cells (GCs) [4, 5]. Moreover, the early presence of high NEFA levels in the follicular fluid can result in subnormal corpus luteum function [6]. High levels of NEFAs during oocyte maturation and early embryo development have hazardous effects on development competence by affecting metabolism and cell fate [3]. In our previous study, exposure to elevated NEFAs level during IVM impaired the development of porcine blastocysts by increasing endoplasmic reticulum (ER) stress, reactive oxygen species (ROS), mitochondrial dysfunction, and by inducing inflammation [7].

The coculture of follicle shell-like GCs encapsulated in an agarose matrix with cumulus-oocyte complexes (COCs) indicated that GCs could improve oocyte developmental competence and blastocyst formation, reduce the intra-oocyte ROS content, and reduce apoptosis [8]. The study of Pawlak et al. [9] indicated that there is bidirectional cooperation of lipid metabolism between oocytes and surrounding cumulus cells. Besides, the coculture of porcine COCs with a monolayer of GCs reduced the stress induced by the presence of high level of NEFAs by activation of anti-inflammatory factors in GCs (unpublished data). However, the mechanisms of the positive effects of GCs on blastocysts obtained from cocultured oocytes exposed to high level of NEFAs are unclear.

In vitro studies on human follicles and embryos are limited because of ethical concerns and limitations. The pig is considered one of the optimal experimental models for human biomedical and reproductive research due to physiological, digestive, and reproductive similarities including the duration of oocyte in vitro maturation, the kinetics of cell division, and the timing of the embryo genome activation [10]. In this study, gene expression in porcine GCs cocultured with COCs and exposed to elevated level of NEFAs during IVM was analyzed to explore the pathways involved in protecting oocytes and improving blastocyst quality after in vitro fertilization.


Differentially expressed genes

Following loess and quantile normalization and the limma algorithm with Flexaray, 1,381 out of the 43,794 probes were differentially expressed in cocultured GCs with NEFAs relative to the control (P<0.05). Of these, 1,295 probes were differentially expressed with the cut off criterion of 1.2-fold change (Fig. 1), and only 39 probes were differentially expressed with the cut off criterion of 1.5-fold change. With the verification of the database, 1,013 probes were annotated as differentially expressed genes and 438 of them were significantly upregulated in NEFA exposed GCs compared to the non-NEFA group, and 575 were downregulated.

Fig. 1
figure 1

Volcano plot visualizing the statistical difference of the DEGs. The horizontal axis represents the p-value of detected transcripts, and the horizontal blue line at –log10 (P-value) of 1.3 is the threshold of 0.05. The vertical axis represents the fold-change (FC) of detected transcripts and the vertical blue lines are the cut-off criterion of 1.2. Each data point represents a gene or a variant of a gene. The upregulated and downregulated genes in the NEFAs treated group are illustrated as the dark red spots upper right and upper left, respectively


The other four granulosa cell replicates collected were used for SYBR green-based qRT-PCR to validate the expression level of six DEGs identified by microarrays. We observed a significant correlation between the qRT-PCR results and the microarray data for the six DEGs. The expression of Mitogen-Activated Protein Kinase Kinase Kinase 8 (MAP3K8), Toll Like Receptor 6 (TLR6), Interleukin 33 (IL33), and Glyceraldehyde-3-Phosphate Dehydrogenase (GAPDH) had similar trends, and TLR6 and GAPDH were significantly expressed by both methods (P<0.05), confirming the reliability of the microarray analysis (Fig. 2).

Fig. 2
figure 2

The abundance of different transcripts from genes affected by the elevated level of NEFAs (upregulated: AGPS, MAP3K8, and TLR6 and downregulated: AGO1, IL33, and GAPDH) was analyzed by real-time RT-qPCR. The RT-qPCR results showed that MAP3K8, TLR6, AGO1, IL33, and GAPDH had the same direction as in the microarray data, and two of them (TLR6 and GAPDH) showed significant differences relative to the respective controls. * represents p-value < 0.05, ** represents p-value < 0.01. GC samples were analyzed after IVM under the high level of NEFAs throughout the IVM culture or control treatment. NC is the control group, and FA is the group exposed to the elevated level of NEFAs. AGPS: Alkylglycerone Phosphate Synthase; MAP3K8: Mitogen-Activated Protein Kinase Kinase Kinase 8; TLR6: Toll Like Receptor 6; AGO1: Argonaute RISC Component 1; IL33: Interleukin 33; GAPDH: Glyceraldehyde-3-Phosphate Dehydrogenase

The expression of genes related to progesterone and estradiol production was also investigated. There was no significant change in the expression of Steroidogenic Acute Regulatory Protein [11], Hydroxy-Delta-5-Steroid Dehydrogenase, 3 Beta- And Steroid Delta-Isomerase 1 (HSD3B1), and Cytochrome P450 Family 19 Subfamily A Member 3 (CYP19A3) in response to exposure to excessive NEFAs, and only the abundance of Cytochrome P450 Family 11 Subfamily A Member 1 (CYP11A1) was increased (Fig. 3) (P<0.05).

Fig. 3
figure 3

Expression of genes related to steroidogenesis. GCs were cultured for 7 days to form a monolayer. On day 8, a monolayer of GCs was cocultured with COCs either with or without the supplementation with the high level of NEFAs during in vitro maturation. The transcript abundance of different genes related to estradiol and progesterone was determined after IVM. The expressions of STAR, HSD3B1, and CYP19A3 was not affected by the high level of NEFAs. NC is the control group, and FA is the group exposed to the elevated level of NEFA. STAR: Steroidogenic Acute Regulatory Protein; CYP11A1: Cytochrome P450 Family 11 Subfamily A Member 1; HSD3B1: Hydroxy-Delta-5-Steroid Dehydrogenase, 3 Beta- And Steroid Delta-Isomerase 1; CYP19A3: Cytochrome P450 Family 19 Subfamily A 3

Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways

To further explore the biological functions associated with the differentially expressed genes, GO annotations were applied. With the annotation of DAVID, 827 IDs were recognized and categorized into 35 biological processes, 18 cellular components, and 15 molecular functions (Fig. 4)(Supplemental Table S1). In terms of percentage of associated genes (one of the enrichment factors), the top significantly enriched biological processes were transcription-related regulation, cellular processes including cell proliferation, growth and cell cycle, as well as localization, catabolism, and protein modification. Interestingly, cell redox hemostasis and apoptosis were also biological processes affected by the high level of NEFAs. In the category cellular components, genes were clustered into terms related to integrin complex, ubiquitin ligase complex, and intracellular structures. The dominant molecular functions were binding of zinc, chromatin and lipid, and transcription regulation.

Fig. 4
figure 4

Gene ontology (GO) analysis was performed using DEGs between cocultured granulosa cells supplemented or not with NEFAs. The 2D bar graph contains biological process (green), cellular component (orange) and molecular function (purple) categories. The vertical axis shows terms of GO categories, and the horizontal axis shows the number of DEGs

The KEGG analysis showed that differentially expressed genes were markedly enriched for [1] metabolism, including carbohydrate, amino acid, and antibiotics; [2] inflammatory responses, including “Toll-like receptor signaling pathway” and “Rheumatoid arthritis”; [3] signal transduction, including “TGF-beta signaling pathway”, “mTOR signaling pathway” and “Thyroid hormone signaling pathway”; and [4] catabolism, including “Peroxisome” and “Endocytosis” (Fig. 5).

Fig. 5
figure 5

Bubble chart of potential signaling pathways generated by KEGG. Pathway analysis was performed to associate the unique DEGs with pathways using the KEGG database. The size and color of each bubble represent the number of DEGs in the pathway and the P-value, respectively. Fold enrichment was calculated as the ratio of genes (the number of genes from the dataset enriched in the pathway in question divided by all the genes from the dataset) divided by the ratio of background (the number of annotated genes from the background enriched in the pathway in question divided by all annotated genes from background)

IPA canonical pathway, tox list, and upstream regulators

Ingenuity Pathway Analysis with human, mouse, and rat data only recognized 876 probes, including 374 up-regulated and 502 down-regulated (Supplemental Table S2). The pathways most affected by the NEFA treatment (pathways which had low p-values) were “Clathrin-mediated endocytosis signaling”, “phospholipase C signaling”, “glucocorticoid receptor signaling”, “semaphoring signaling in neurons”, “estrogen-dependent breast cancer signaling”, “HOTAIR regulatory pathway”, “regulation of eIF4 and p70S6K signaling”, “EGF signaling”, “Huntington’s disease signaling”, and “Insulin secretion signaling pathway” (Fig. 6) (Supplemental Table S3). These highly significant pathways are related to endocytosis, signaling transduction, and metabolism, in agreement with the results from the KEGG database analysis. Besides, differentially expressed genes affected by excessive NEFA exposure were enriched in the tox lists related to “cell death” and “mitochondrial dysfunction” including “decreases respiration of mitochondria” and “decreases transmembrane potential of mitochondria and mitochondrial membrane” (Fig. 7) (Supplemental Table S4). The most significant upstream regulators identified were beta-estradiol, erb-B2 receptor tyrosine kinase 2 (ERBB2), trichostatin A, tumor protein P53 (TP53), and dexamethasone (Supplemental Table S5) Fig. 8.

Fig. 6
figure 6

Bubble chart of the top 20 canonical pathways generated by IPA analysis. The size and color of each bubble represent the number of DEGs in each pathway and the P-value, respectively. The ratio was calculated as the number of genes from our dataset that overlap with the canonical pathway in question divided by the total number of genes that are represented in that canonical pathway

Fig. 7
figure 7

Bubble chart of tox lists affected by the high level of NEFAs with a p-value less than 0.05, which was generated by IPA analysis. The size and color of each bubble represent the number of DEGs in each pathway and the P-value, respectively. The ratio was calculated as the number of genes of our dataset that overlap with the tox list divided by the total number of genes that are represented in that tox list

Fig. 8
figure 8

Upstream regulators with predicted activation or inhibition (with an absolute value of Z-sore greater than 2) affected by the high level of NEFAs from IPA analysis. These upstream regulators were categorized into 9 types including chemicals, cytokine, enzyme, group, growth factor, kinase, other, transcription factor, and transporter. The vertical axis shows categories and names of regulators, and the horizontal axis shows the Z-score (red represents activated, green represents inhibited)


Elevated levels of NEFAs are caused by either lipolysis or fat intake; therefore, the concentration of NEFAs varies importantly according to body condition and nutrition. Elevated NEFAs activated ROS-mediated inflammation, apoptosis, mitochondrial dysfunction, and insulin resistance in hepatocytes [12,13,14]. Moreover, high levels of NEFAs reduced reproductive performance in humans and bovine [1, 15]. High levels of NEFAs also induced apoptosis and impaired cell viability, and steroidogenic capacity in cultured ovarian cells [3,4,5]. Moreover, cumulus cells were shown to protect oocytes from the stress of excess NEFAs and to allow better blastocysts to be generated [16, 17]. In accordance with this, our previous study found that blastocysts which originated from oocytes exposed to elevated NEFA levels during IVM had impaired mitochondrial functions and increased inflammatory response [7] (unpublished data). Besides, some of the affected pathways in blastocysts from oocytes cocultured with GCs during IVM were associated with reduction of inflammation, indicating that the GCs were protecting the cocultured COCs and especially the oocytes (unpublished data). In this study, the original results provided some explanation on how fatty acids may rapidly impact the functions of GCs and how the cells adapt to their presence by triggering a defense mechanism.

Inhibited luteinization of cocultured GCs under elevated NEFAs levels

Granulosa cells are one of the epithelial cell types that display polarity, cellular adhesion, and cell-cell junctions. Following the LH surge, GCs lose their epithelial nature and undergo an epithelial-mesenchymal transition (EMT) known as “luteinization” [18,19,20]. Luteinization of GCs can be induced by physiological doses of LH, FSH, or 3’,5’-cAMP in vitro [21]. In this study, hepatocyte growth factor (HGF), EGF, and transforming growth factor-beta (TGFβ) signaling pathways were significantly affected by the elevated NEFA level (showed in IPA canonical pathways and KEGG pathways), and these factors are involved in the initiation and regulation of the EMT [20, 22]. Moreover, IPA analysis results suggested that the EMT process was downregulated in GCs exposed to high NEFAs concentration. The MYC proto-oncogene (MYC), specificity protein 1 (SP1), and jagged canonical notch ligand 1 (JAG1) are involved in inducing the EMT [23,24,25]. Prediction of their down-regulation in this study implies that the EMT from GCs to luteal cells was possibly inhibited in vitro in the presence of COCs and NEFAs. A recent study suggested that increased concentrations of oleic acid (OA) may promote the EMT process in bovine GCs [26]. Oleic acid was demonstrated to have protective effects against lipotoxicity in oocytes and early embryos, contrarily to palmitic acid (PA), stearic acid (SA), and probably NEFAs (a combination of PA, SA, and OA)[9, 27, 28]. The promotion of luteinization by OA partly supports our result that the high level of NEFAs inhibited luteinization to some extent. The inhibition of EMT is somewhat similar to what was observed in post-partum dairy cows where high NEFAs were associated with slow follicular growth and decreased differentiation leading to delayed ovulation and infertility [29, 30]. Besides, elevated levels of NEFAs affected follicle development and inhibited ovulation in bovine [31]. In humans, extreme obesity is associated with amenorrhea and polycystic ovary syndrome, and both conditions can lead to impairment of the final differentiation of the follicle prior to ovulation [32, 33]. Moreover, the coculture of COCs may also interfere with the EMT transition of GCs because the maintenance of efficient bidirectional interaction between granulosa cells and COCs can reduce apoptosis in GCs [34].

Moreover, elevated NEFA levels were also reported to decrease progesterone production in bovine thecal cells, and inhibit testosterone production in mouse Leydig cells [35, 36]. These results can be explained by the inhibitory effect of NEFAs on the hydroxylation of cholesterol [36]. Besides, increased OA concentrations decreased the production of estradiol in bovine GCs [26]. In the present study, the expression of HSD3B1, which is involved in the conversion from pregnenolone to progesterone, was not affected by the elevated level of NEFAs, suggesting that progesterone production was not increased. The protein STAR is responsible for the transportation of cholesterol, the substrate for steroidogenesis; and CYP11A1 encodes P450scc, which catalyzes the formation of pregnenolone, the precursor for both progesterone and estradiol [37]. In the pig genome, there are three paralogs (CYP19A1, CYP19A2, and CYP19A3) of the aromatase-encoding gene CYP19A, the key gene of estradiol production, and CYP19A3 is highly expressed in the ovary [38, 39]. The maintained and increased levels of these two genes indicate that cocultured GCs can maintain estradiol production which keeps the follicle in development, and that luteinization of cocultured GCs was not stimulated under the physiologically elevated NEFA levels used here. The contradictory results may be explained by the counteracting effects of OA against PA, SA, and NEFAs. Whether higher NEFA levels (either one saturated fatty acids or mixed fatty acids) would inhibit progesterone production or accelerate luteinization in pig GCs needs to be explored.

Apoptosis and Inflammation

When cells are in the presence of an excess of NEFAs, these NEFAs are stored as lipid droplets, or they enter mitochondria to be β-oxidated to produce ATP. However, elevated NEFAs remaining in the cytosol may undergo lipid peroxidation, generating ROS and depleting protective glutathione levels [40]. The regulation of peroxisome, glycolysis/gluconeogenesis, and pyruvate metabolism highlighted in this study indicates that the original balance between β-oxidation and glycolysis in GCs was disrupted by the elevated NEFAs. Besides, the affected cell redox hemostasis suggests an increase in ROS production. In hepatocytes, NEFAs can activate oxidative stress-mediated NF-κB signaling to induce an inflammatory response [13]. In our previous study, exposure to elevated NEFAs level during IVM resulted in blastocysts with increased inflammation, apoptosis, and mitochondrial dysfunction associated with ROS production [7]. In line with this fact, endocytosis, rheumatoid arthritis, and toll-like receptor signaling pathways were significantly affected in GCs in the present study. Besides, NEFAs induced apoptosis and inhibited the proliferation of bovine and porcine GCs [5, 41], and directly caused ER stress and increased the intracellular Ca2+ level in oocytes [40]. The most significantly affected upstream regulators in this study, ERBB2 and TP53, indicated that elevated NEFAs induced inflammation and apoptosis in GCs similar to what we observed in blastocysts obtained from NEFA-exposed oocytes (unpublished data). Cells can adapt to the increased lipid stress to survive, but prolonged stress initiates inflammation and apoptosis [42]. This study showed that NEFAs induced inflammation and apoptosis in cocultured GCs.

Impaired mitochondrial function

The mitochondrion is the principal organelle that produces ATP, and its well-organized functioning guarantees the efficiency of cellular processes. The increased ROS content due to the β-oxidation of excess fatty acids may cause mitochondrial dysfunction. A balance of pyruvate and fatty acid oxidation is required in healthy mammalian oocyte mitochondria to maintain a low level of ROS production [42]. Besides, the disrupted balance mentioned above was induced by the high level of NEFAs suggesting increased production of ROS. Excessive palmitic acid induced ROS formation, caspase 3 activated apoptosis, and deterioration of mitochondria, which contributed to the degeneration of the cumulus cell layers [43]. In oocytes, NEFAs can induce oocyte mitochondrial dysfunction with abnormal ATP production and high ROS content [41, 44, 45]. Moreover, our previous studies also demonstrated NEFA induced mitochondrial dysfunction in blastocysts from oocytes exposed to NEFAs [7]. Consistent with the NEFA induced effects, mitochondrial functions and cell redox hemostasis were also affected in this study, suggesting that the oxidation of excessive NEFAs created more ROS and impaired the mitochondrial function in cocultured GCs.

GCs protect COCs against exposure to excessive NEFA levels

There is a possible protection mechanism in follicular cells that can desaturate the potentially toxic saturated fatty acids into mono-unsaturated fatty acids [17]. Instead of a lipotoxic impact, unsaturated fatty acids generally have a moderate impact on somatic cells [46]. Lipid storage and the β-oxidation of NEFAs in GCs [40] may be able to partly explain how the GCs protected the cocultured COCs and especially the oocytes. The genes STAT6, MYC, SP1 were affected potential upstream regulators which are involved in inflammation, apoptosis, and ER stress [47, 48]. Downregulation of these regulators implies that the inflammatory and apoptotic responses, which should have been induced by elevated NEFAs, were attenuated to a certain degree in cocultured GCs. When we compared blastocysts obtained from oocytes exposed to NEFAs alone to blastocysts from oocytes exposed to NEFAs in the presence of GCs, we observed that the co-culture reduced the inflammation reaction. However, the precise factors involved remain unknown and more research is needed.

Effects of the high levels of NEFAs on different ovarian cell types

When simultaneously exposed to elevated levels of NEFAs either in vivo (pathological nutrition status including NEB and obesity) or in vitro, the different ovarian cell types respond differently. Commonly, apoptosis is induced by elevated levels of NEFAs in pig, human, bovine, and mouse granulosa cells [4, 49, 50], and apoptosis and inflammation are also increased in blastocysts originated from oocytes exposed to NEFAs [2, 7, 51]. Since granulosa cells and blastocysts go through mitosis, the proliferation of granulosa cells and the cell number of blastocysts are also affected by high levels of NEFAs [35, 41, 51, 52]. Lipid content is increased in granulosa cells, cumulus cells, and oocytes [16, 41, 44, 53, 54], and mitochondrial dysfunctions are induced in granulosa cells, cumulus cells, oocytes, and blastocysts [42,43,44, 55, 56]. The cellular metabolism of granulosa cells and blastocysts originated from oocytes exposed to NEFAs is also changed because fatty acids can be oxidized and produce ATP [2, 7, 51]. The more specific responses to elevated levels of NEFAs are more dependent on the type of cells, including the development ability of oocytes and blastocysts (maturation, fertilization, and development), the steroidogenesis activity of granulosa cells and theca cells, and the epithelial-mesenchymal transition of granulosa cells to luteal cells [2, 5, 7, 9, 27, 28, 35, 52, 57, 58].


The model for the mechanisms of the effects of high levels of NEFAs on GCs and how the GCs protect COCs is shown in Fig. 9. High levels of NEFAs in the IVM medium may inhibit the luteinization of porcine GCs cocultured with COCs, potentially explaining some delayed ovulation problems in other species. Besides, elevated NEFA levels increased the ROS-mediated inflammatory response in GCs, but this reaction may protect oocytes and explain the protective effect of GCs during IVM as described before.

Fig. 9
figure 9

The response of GCs to high levels of NEFAs and how GCs protect COCs exposed to high levels of NEFAs. The cell survival, proliferation, and steroidogenesis of granulosa cells are affected by exposure to high levels of NEFAs. Apoptosis, ER stress, inflammation, and lipid accumulation in granulosa cells are increased. The epithelial-mesenchymal transition of granulosa cells to luteal cells is also negatively affected. Because the bidirectional communication between granulosa cells, cumulus cells, and oocytes ensures follicle development and oocyte maturation, granulosa cells and cumulus cells protect the oocyte when they are exposed to hazardous environments. The quantity of NEFAs that the oocyte is exposed to is reduced because some of the NEFAs are either stored as lipid droplets or oxidized in granulosa cells and cumulus cells. Besides, granulosa cells can produce and deliver some anti-inflammatory factors to cumulus cells and the oocyte. Therefore, blastocysts from oocytes exposed to the high level of NEFAs still exhibit the anti-inflammatory response, and the blastocyst rate from oocytes cocultured with GCs while exposed to NEFAs was decreased to a lesser extent than the blastocyst rate from exposed oocytes cultured without granulosa cells

Materials and methods


All supplements and reagents used in this study were purchased from sigma-Aldrich (St. Louis, MO, United State) unless otherwise specified.

Preparation of NEFAs

PA, SA, and OA were dissolved in absolute ethanol at 468 mM, 194mM, and 534mM, respectively to make the stock solution as described previously [7]. The NEFA mixture was added to North Carolina State University (NCSU)-23 medium with 3 % fatty-acid-free bovine serum albumin and was sonicated for 30 min with an Arrayit ultrasonic bath sonicator (Sunnyvale, CA). The medium containing the NEFAs was then incubated at 38.5 °C in a humidified atmosphere of 5 % CO2 in air for 24 h for equilibration before use.

Collection of GCs and COCs

Porcine ovaries were collected from a local slaughterhouse and transported within 2 h to the laboratory in physiological saline supplemented with 100 IU/L of penicillin G and 100 mg/mL of streptomycin sulfate at approximately 35 °C. The ovaries were washed three times with saline warmed to 37 °C and COCs, GCs, and follicular fluid were aspirated from follicles 3-6 mm in diameter with an 18-gauge needle attached to a 10-ml syringe, avoiding blood contamination. The cell mixture was washed twice in warmed Tyrode’s lactate-buffered HEPES supplemented with 0.1 % polyvinyl alcohol (PVA-TL-HEPES).

Preparation of granulosa cell monolayers

After the COCs had fallen to the bottom of the container due to their larger size, the granulosa cell-containing superstratum of follicular fluid was passed through a Falcon 40-µm mesh cell strainer (BD Biosciences, Stockholm, Sweden), centrifuged at 800 × g for 8 min and washed three times with warm phosphate-buffered saline (PBS, Corning, New York, USA) containing penicillin-streptomycin (100 U per mL, ThermoFisher Scientific, Waltham, USA) and 50 µg/mL gentamycin. The cells were re-suspended at 5 × 105 per mL in DMEM/F12 (ThermoFisher Scientific, Waltham, USA) supplemented with 10 % fetal calf serum (Corning, New York, USA), penicillin-streptomycin (100 U per mL) and 50 µg/mL gentamycin and cultured in 4-well plates (ThermoFisher Scientific, Waltham, USA) at 38.5 °C under a humidified atmosphere (5 % CO2). At the 24-h point, cells were washed with PBS twice and re-supplemented with fresh Dulbecco’s minimal essential medium/Ham’s F12 (DMEM/F12) with 10 % fetal calf serum for 6 more days of culture during which the medium (70 %) was changed every 48 h.

Co-culture of COCs and granulosa cell monolayers

After 7 days of culture, the DMEM/F12 culture medium was replaced with NEFA-containing NCSU-23 (or NCSU-23 without NEFA but with same amount of BSA and ethanol as control) supplemented with 10 ng/mL epidermal growth factor (EGF), 0.5 g/mL follicle-stimulating hormone (FSH), 0.5 g/mL luteinizing hormone (LH), 1 mM dibutyryl-cyclic adenosine monophosphate, 1 % (v/v) purified porcine follicular fluid, and 50 µg/mL gentamycin. COCs were collected from the follicular fluid of a different batch of ovaries then washed with PVA-TL-HEPES. Those with cytoplasmic uniformity and at least three layers of cumulus cells fully surrounding the oocyte [59] were selected for maturation in 4-well plates containing granulosa cell monolayers. After 22 h, COCs were washed, put back in the same plate with fresh NEFA-containing NCSU-23 (or control NCSU-23) supplemented with only EGF and follicular fluid and incubated for another 22 h. After maturation, COCs were removed from the culture medium and the GC monolayers were washed three times with PBS. Then, Trizol (ThermoFisher Scientific) was used to capture the adherent cells and the mixture was stored at – 80 °C until used.

RNA extraction

The Direct-zolTM RNA Miniprep Plus (Zymo Research, CA, USA) was used to extract the RNA of GCs in Trizol. Cells from four wells of one 4-well plate were pooled as one replicate, and eight biological replicates were collected on different days (four for microarray analysis and four for quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) validation). An equal volume of ethanol was added to cells lysed in Trizol and mixed thoroughly. The mixture was passed through a Zymo-SpinTM IIIGG column with a DNase I treatment, followed by the washing buffer. The quality and concentration of RNA were analyzed with a Bioanalyzer (Agilent, Mississauga, Ontario, Canada). All the extracted RNA samples were of good quality, with integrity numbers > 8.5. Samples were stored at −80 °C.

Manipulation of RNA microarrays

Five ng of RNA from four replicates from two conditions were amplified linearly using a RiboAmp® HS PLUS RNA Kit (Applied Biosystems). Antisense RNA (aRNA) was produced using T7 RNA polymerase and labeled with Cy3 and Cy5 using a ULS Fluorescent Labeling Kit (Leica, Germany). Two corresponding treatments (900 ng aRNA per replicate) were mixed in a two-color dye-swap design and hybridized on an Agilent-manufactured EmbryoGENE 4 × 44 K porcine transcriptome microarray slide (Agilent Technologies, CA, USA) for 17 h at 65 °C. The slide was washed for 1 min in Expression Wash Buffer 1 (RT), 3 min in Expression Wash Buffer 2 (42 °C), 10 s in 100 % acetonitrile (RT), and 30 s in Stabilization and Drying Solution. Microarray slide was then scanned with the PowerScanner (Tecan, Switzerland) and features were extracted with Array-pro 6.3 (MediaCybernetics, MD). Intensity files were normalized and analyzed with FlexArray. The background was subtracted from raw fluorescence intensity prior to normalization within Cy-3/Cy-5 and between each array with Loess and quantile, respectively.

Data analysis

Treatment effects were analyzed with the limma algorithm, which attributed to each probe a probability of differing significantly from the control. The differentially expressed genes (DEGs) were then annotated by the latest version (Revision 6) of the Affymetrix porcine annotation database combined with the current annotation data of pig (Sus scrofa 11.1) in NCBI [37]. GO and pathway enrichment analyses were performed using the DAVID Database (V6.8) with Sus scrofa background [60]. Data were also analyzed using Ingenuity Pathway Analysis (IPA, Ingenuity® Systems, with human, mouse, and rat backgrounds. IPA analysis considers the general tendency in gene expression variations, including the direction in which each tends to change and the concordance between the variations of the significantly affected genes and the signaling pathways. It then determines the most plausible upstream regulators (explaining the observed gene expression changes) and tox list (linking to pharmacological responses) by taking into consideration the results and the usual context in which the same genes are affected. The extremely high similarity between the porcine and human genomes and the availability of more annotated functions in human studies made the IPA analysis a necessary complement.

Quantitative Real-Time RT-qPCR

Total RNA extracted from the other four biological replicates was reverse transcribed using oligo-dT and random primers with the qScriptTM Flex cDNA Synthesis Kit (Quanta Biosciences, MD). Real-time PCR reactions were conducted on a LightCycler 480 System (Roche, Germany) with LightCycler 480® SYBR Green I Master. The efficiency of the primers and the standard curve for each gene were deduced from serial dilutions (1 pg to 0.1 fg) of the corresponding cDNA fragment used as template. The geometric mean of two housekeeping genes, succinate dehydrogenase complex flavoprotein subunit A (SDHA) and peptidylprolyl isomerase A (PPIA), [61], was used to normalize gene expression with the GeNorm VBA applet software. Primer information is provided in Table 1. Statistical analysis was conducted in Prism (GraphPad Software, version 5.00 for Windows, San Diego California USA) with the independent-samples t-test.

Table 1 Information on RT-qPCR primers

Availability of data and materials

The data sets of transcriptomic microarray results were deposited in National Center for Biotechnology Information’s Gene Expression Omnibus and are accessible through the accession number GSE168957.



Non-Esterified Fatty Acids


Granulosa Cells


Cumulus-Oocyte Complexes


In Vitro Maturation


Reactive Oxygen Species


Endoplasmic Reticulum


Palmitic Acid


Stearic Acid


Oleic Acid


North Carolina State University


Tyrode’s Lactate-buffered HEPES supplemented with 0.1 % Polyvinyl Alcohol


Phosphate-Buffered Saline


Dulbecco’s Minimal Essential Medium/Ham’s F12


Epidermal Growth Factor


Follicle-Stimulating Hormone


Luteinizing Hormone


Antisense RNA


Room Temperature


Differentially Expressed Genes


Gene Ontology


Ingenuity Pathway Analysis


Succinate Dehydrogenase Complex Flavoprotein Subunit A


Peptidylprolyl Isomerase A


Quantitative Real-time Reverse Transcription Polymerase Chain Reaction


Mitogen-Activated Protein Kinase Kinase Kinase 8


Toll Like Receptor 6


Interleukin 33


Glyceraldehyde-3-Phosphate Dehydrogenase


Steroidogenic Acute Regulatory Protein


Hydroxy-Delta-5-Steroid Dehydrogenase, 3 Beta- And Steroid Delta-Isomerase 1


Cytochrome P450 Family 11 Subfamily A Member 1


Kyoto Encyclopedia of Genes and Genomes


Erb-B2 Receptor Tyrosine Kinase 2


Tumor Protein P53


Epithelial-Mesenchymal Transition


Hepatocyte Growth Factor


Transforming Growth Factor-beta


MYC Proto-oncogene


Specificity Protein 1


Jagged Canonical Notch Ligand 1


  1. Valck × SD, Arias-Alvarez M, De Pauw I, Fievez V, Vlaeminck B, Fransen E, et al. Fatty acid composition of the follicular fluid of normal weight, overweight and obese women undergoing assisted reproductive treatment: a descriptive cross-sectional study. Reproductive Biology Endocrinology. 2014;12(1):13.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Desmet KL, Van Hoeck V, Gagne D, Fournier E, Thakur A, O’Doherty AM, et al. 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(1):1004.

    Article  CAS  Google Scholar 

  3. Van Hoeck V, Bols PE, Binelli M, Leroy JL. Reduced oocyte and embryo quality in response to elevated non-esterified fatty acid concentrations: a possible pathway to subfertility? Anim Reprod Sci. 2014;149(1-2):19–29.

    Article  PubMed  CAS  Google Scholar 

  4. Mu Y-M, Yanase T, Nishi Y, Tanaka A, Saito M, Jin C-H, et al. Saturated FFAs, palmitic acid and stearic acid, induce apoptosis in human granulosa cells. Endocrinology. 2001;142(8):3590–7.

    Article  CAS  PubMed  Google Scholar 

  5. Vanholder T, Leroy J, Van Soom A, Opsomer G, Maes D, Coryn M, et al. Effect of non-esterified fatty acids on bovine granulosa cell steroidogenesis and proliferation in vitro. Anim Reprod Sci. 2005;87(1):33–44.

    Article  CAS  PubMed  Google Scholar 

  6. Yung M, VandeHaar M, Fogwell R, Sharma B. Effect of energy balance and somatotropin on insulin-like growth factor I in serum and on weight and progesterone of corpus luteum in heifers. Journal of animal science. 1996;74(9):2239–44.

    Article  CAS  PubMed  Google Scholar 

  7. Shi M, Sirard M-A. Transcriptome and epigenome analysis of porcine embryos from non-esterified fatty acids exposed oocyte. Domest Anim Endocrinol. 2021;76:106605.

  8. Park JE, Lee J, Lee ST, Lee E. In vitro maturation on ovarian granulosa cells encapsulated in agarose matrix improves developmental competence of porcine oocytes. Theriogenology. 2021;164:42–50.

  9. Pawlak P, Malyszka N, Szczerbal I, Kolodziejski P. Fatty acid induced lipolysis influence embryo development, gene expression and lipid droplet formation in the porcine cumulus cells. Biol Reprod. 2020;103(1):36–48.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Madeja ZE, Pawlak P, Piliszek A. Beyond the mouse: Non-rodent animal models for study of early mammalian development and biomedical research. Int J Dev Biol. 2019;63(3-4-5):187–201.

    Article  CAS  PubMed  Google Scholar 

  11. Testart J, Gautier E, Brami C, Rolet F, Sedbon E, Thebault A. Genetics. Intracytoplasmic sperm injection in infertile patients with structural chromosome abnormalities. Hum Reprod. 1996;11(12):2609–12.

    Article  CAS  PubMed  Google Scholar 

  12. Song Y, Li X, Li Y, Li N, Shi X, Ding H, et al. Non-esterified fatty acids activate the ROS–p38–p53/Nrf2 signaling pathway to induce bovine hepatocyte apoptosis in vitro. Apoptosis. 2014;19(6):984–97.

    Article  CAS  PubMed  Google Scholar 

  13. Shi X, Li D, Deng Q, Li Y, Sun G, Yuan X, et al. NEFAs activate the oxidative stress-mediated NF-κB signaling pathway to induce inflammatory response in calf hepatocytes. J Steroid Biochem Mol Biol. 2015;145:103–12.

    Article  CAS  PubMed  Google Scholar 

  14. Gao W, Du X, Lei L, Wang H, Zhang M, Wang Z, et al. NEFA-induced ROS impaired insulin signalling through the JNK and p38MAPK pathways in non‐alcoholic steatohepatitis. J Cell Mol Med. 2018;22(7):3408–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Leroy J, Vanholder T, Mateusen B, Christophe A, Opsomer G, de Kruif A, et al. Non-esterified fatty acids in follicular fluid of dairy cows and their effect on developmental capacity of bovine oocytes in vitro. Reproduction. 2005;130(4):485–95.

    Article  CAS  PubMed  Google Scholar 

  16. Aardema H, Lolicato F, van de Lest CH, Brouwers JF, Vaandrager AB, van Tol HT, et al. Bovine cumulus cells protect maturing oocytes from increased fatty acid levels by massive intracellular lipid storage. Biol Reprod. 2013;88(6):164, 1–15.

    Article  PubMed  Google Scholar 

  17. Aardema H, van Tol HTA, Wubbolts RW, Brouwers JFHM, Gadella BM, Roelen BAJ. Stearoyl-CoA desaturase activity in bovine cumulus cells protects the oocyte against saturated fatty acid stress. Biol Reprod. 2017;96(5):982–92.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Childs AJ, McNeilly AS. Epithelial-to-mesenchymal transition in granulosa cells: a key to activation of follicle growth? Biol Reprod. 2012;86(5):152, 1–2.

    Article  PubMed  CAS  Google Scholar 

  19. Irving-Rodgers HF, Harland ML, Rodgers RJ. A novel basal lamina matrix of the stratified epithelium of the ovarian follicle. Matrix biology. 2004;23(4):207–17.

    Article  CAS  PubMed  Google Scholar 

  20. Nivet A-L, Dufort I, Gilbert I, Sirard M-A. Short-term effect of FSH on gene expression in bovine granulosa cells in vitro. Reproduction Fertility Development. 2018;30(8):1154–60.

    Article  CAS  PubMed  Google Scholar 

  21. Channing CP, Tsai V, Sachs D. Role of insulin, thyroxin and cortisol in luteinization of porcine granulosa cells grown in chemically defined media. Biol Reprod. 1976;15(2):235–47.

    Article  CAS  PubMed  Google Scholar 

  22. Stone RC, Pastar I, Ojeh N, Chen V, Liu S, Garzon KI, et al. Epithelial-mesenchymal transition in tissue repair and fibrosis. Cell tissue research. 2016;365(3):495–506.

    Article  CAS  PubMed  Google Scholar 

  23. Zavadil J, Cermak L, Soto-Nieves N, Böttinger EP. Integration of TGF‐β/Smad and Jagged1/Notch signalling in epithelial‐to‐mesenchymal transition. EMBO J. 2004;23(5):1155–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Cho KB, Cho MK, Lee WY, Kang KW. Overexpression of c-myc induces epithelial mesenchymal transition in mammary epithelial cells. Cancer letters. 2010;293(2):230–9.

    Article  CAS  PubMed  Google Scholar 

  25. Kolesnikoff N, Attema JL, Roslan S, Bert AG, Schwarz QP, Gregory PA, et al. Specificity protein 1 (Sp1) maintains basal epithelial expression of the miR-200 family: implications for epithelial-mesenchymal transition. J Biol Chem. 2014;289(16):11194–205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Yenuganti VR, Koczan D, Vanselow J. Genome wide effects of oleic acid on cultured bovine granulosa cells: evidence for the activation of pathways favoring folliculo-luteal transition. BMC Genomics. 2021;22(1):1–10.

    Google Scholar 

  27. Aardema H, Vos PL, Lolicato F, Roelen BA, Knijn HM, Vaandrager AB, et al. Oleic acid prevents detrimental effects of saturated fatty acids on bovine oocyte developmental competence. Biol Reprod. 2011;85(1):62–9.

    Article  CAS  PubMed  Google Scholar 

  28. Sharma A, Baddela VS, Becker F, Dannenberger D, Viergutz T, Vanselow J. Elevated free fatty acids affect bovine granulosa cell function: a molecular cue for compromised reproduction during negative energy balance. Endocrine connections. 2019;8(5):493–505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Chaput C, Sirard M. Embryonic response to high beta-hydroxybutyrate (BHB) levels in postpartum dairy cows. Domest Anim Endocrinol. 2020;72:106431.

  30. Girard A, Dufort I, Sirard M-A. The effect of energy balance on the transcriptome of bovine granulosa cells at 60 days postpartum. Theriogenology. 2015;84(8):1350–61.e6.

    Article  CAS  PubMed  Google Scholar 

  31. Baddela VS, Sharma A, Vanselow J. Non-esterified fatty acids in the ovary: friends or foes? Reproductive Biology Endocrinology. 2020;18:1–14.

    Article  CAS  Google Scholar 

  32. Foong SC, Abbott DH, Zschunke MA, Lesnick TG, Phy JL, Dumesic DA. Follicle luteinization in hyperandrogenic follicles of polycystic ovary syndrome patients undergoing gonadotropin therapy for in vitro fertilization. The Journal of Clinical Endocrinology Metabolism. 2006;91(6):2327–33.

    Article  CAS  PubMed  Google Scholar 

  33. Silvestris E, de Pergola G, Rosania R, Loverro G. Obesity as disruptor of the female fertility. Reproductive Biology Endocrinology. 2018;16(1):1–13.

    Article  CAS  Google Scholar 

  34. Palmerini MG, Nottola SA, Tunjung WAS, Kadowaki A, Bianchi S, Cecconi S, et al. EGF-FSH supplementation reduces apoptosis of pig granulosa cells in co-culture with cumulus-oocyte complexes. Biochem Biophys Res Commun. 2016;481(1-2):159–64.

    Article  CAS  PubMed  Google Scholar 

  35. Vanholder T, Leroy JL, Van Soom A, Maes D, Coryn M, Fiers T, et al. Effect of non-esterified fatty acids on bovine theca cell steroidogenesis and proliferation in vitro. Anim Reprod Sci. 2006;92(1-2):51–63.

    Article  CAS  PubMed  Google Scholar 

  36. Meikle AW, Benson SJ, Liu X, Boam WD, Stringham JD. Nonesterified fatty acids modulate steroidogenesis in mouse Leydig cells. American Journal of Physiology-Endocrinology Metabolism. 1989;257(6):E937-E42.

    Article  Google Scholar 

  37. Zhang J, Liu Y, Yao W, Li Q, Liu H, Pan Z. Initiation of follicular atresia: gene networks during early atresia in pig ovaries. Reproduction. 2018;156(1):23–33.

    Article  CAS  PubMed  Google Scholar 

  38. Vanselow J, Conley AJ, Berger T. Aromatase and the three little paralogs. Biol Reprod. 2021;105(1):5–6.

  39. Vanselow J, Conley AJ, Corbin CJ, Berger T. Genomic Structure of the Porcine CYP19 Locus and Expression of the CYP19A3 Paralog. Genes. 2021;12(4):533.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Aardema H, van Tol HT, Vos PL. An overview on how cumulus cells interact with the oocyte in a condition with elevated NEFA levels in dairy cows. Anim Reprod Sci. 2019;207:131–7.

    Article  CAS  PubMed  Google Scholar 

  41. Shibahara H, Ishiguro A, Inoue Y, Koumei S, Kuwayama T, Iwata H. Mechanism of palmitic acid-induced deterioration of in vitro development of porcine oocytes and granulosa cells. Theriogenology. 2020;141:54–61.

    Article  CAS  PubMed  Google Scholar 

  42. Bradley J, Swann K. Mitochondria and lipid metabolism in mammalian oocytes and early embryos. Int J Dev Biol. 2019;63(3-4-5):93–103.

    Article  CAS  PubMed  Google Scholar 

  43. Lolicato F, Brouwers JF, de Lest CHv, Wubbolts R, Aardema H, Priore P, et al. The cumulus cell layer protects the bovine maturing oocyte against fatty acid-induced lipotoxicity. Biol Reprod. 2015;92(1):16, 1-.

    Article  PubMed  CAS  Google Scholar 

  44. Itami N, Shirasuna K, Kuwayama T, Iwata H. Palmitic acid induces ceramide accumulation, mitochondrial protein hyperacetylation, and mitochondrial dysfunction in porcine oocytes. Biol Reprod. 2018;98(5):644–53.

    Article  PubMed  Google Scholar 

  45. Ogawa K, Itami N, Ueda M, Kansaku K, Shirasuna K, Kuwayama T, et al. Non-esterified fatty acid‐associated ability of follicular fluid to support porcine oocyte maturation and development. Reproductive medicine biology. 2018;17(2):155–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Coll T, Eyre E, Rodríguez-Calvo R, Palomer X, Sánchez RM, Merlos M, et al. Oleate reverses palmitate-induced insulin resistance and inflammation in skeletal muscle cells. Journal of biological chemistry. 2008;283(17):11107–16.

    Article  CAS  PubMed  Google Scholar 

  47. Leon-Cabrera SA, Molina-Guzman E, Delgado-Ramirez YG, Vázquez-Sandoval A, Ledesma-Soto Y, Pérez-Plasencia CG, et al. Lack of STAT6 attenuates inflammation and drives protection against early steps of colitis-associated colon cancer. Cancer immunology research. 2017;5(5):385–96.

    Article  CAS  PubMed  Google Scholar 

  48. Sarosiek KA, Fraser C, Muthalagu N, Bhola PD, Chang W, McBrayer SK, et al. Developmental regulation of mitochondrial apoptosis by c-Myc governs age-and tissue-specific sensitivity to cancer therapeutics. Cancer cell. 2017;31(1):142–56.

    Article  CAS  PubMed  Google Scholar 

  49. Chen Z, Lei L, Wen D, Yang L. Melatonin attenuates palmitic acid-induced mouse granulosa cells apoptosis via endoplasmic reticulum stress. Journal of ovarian research. 2019;12(1):43.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Wang Y, Li C, Li J, Wang G, Li L. Non-esterified Fatty Acid-Induced Reactive Oxygen Species Mediated Granulosa Cells Apoptosis Is Regulated by Nrf2/p53 Signaling Pathway. Antioxidants. 2020;9(6):523.

    Article  CAS  PubMed Central  Google Scholar 

  51. Van Hoeck V, Sturmey RG, Bermejo-Alvarez P, Rizos D, Gutierrez-Adan A, Leese HJ, et al. Elevated non-esterified fatty acid concentrations during bovine oocyte maturation compromise early embryo physiology. PloS one. 2011;6(8):e23183.

  52. Jorritsma R, Cesar M, Hermans J, Kruitwagen C, Vos P, Kruip T. Effects of non-esterified fatty acids on bovine granulosa cells and developmental potential of oocytes in vitro. Anim Reprod Sci. 2004;81(3-4):225–35.

    Article  CAS  PubMed  Google Scholar 

  53. Carro M, Buschiazzo J, Rios GL, Oresti GM, Alberio R. Linoleic acid stimulates neutral lipid accumulation in lipid droplets of maturing bovine oocytes. Theriogenology. 2013;79(4):687–94.

    Article  CAS  PubMed  Google Scholar 

  54. Yang X, Wu LL, Chura LR, Liang X, Lane M, Norman RJ, et al. Exposure to lipid-rich follicular fluid is associated with endoplasmic reticulum stress and impaired oocyte maturation in cumulus-oocyte complexes. Fertility sterility. 2012;97(6):1438–43.

    Article  CAS  PubMed  Google Scholar 

  55. Igosheva N, Abramov AY, Poston L, Eckert JJ, Fleming TP, Duchen MR, et al. Maternal diet-induced obesity alters mitochondrial activity and redox status in mouse oocytes and zygotes. PloS one. 2010;5(4):e10074.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Marei WF, De Bie J, Mohey-Elsaeed O, Wydooghe E, Bols PE, Leroy JL. Alpha-linolenic acid protects the developmental capacity of bovine cumulus–oocyte complexes matured under lipotoxic conditions in vitro. Biol Reprod. 2017;96(6):1181–96.

    Article  PubMed  Google Scholar 

  57. Valck × SD, Van Hoeck V, Arias-Alvarez M, Maillo V, Lopez-Cardona AP, Gutierrez-Adan A, et al. Elevated non-esterified fatty acid concentrations during in vitro murine follicle growth alter follicular physiology and reduce oocyte developmental competence. Fertility sterility. 2014;102(6):1769–76. e1.

    Article  CAS  PubMed  Google Scholar 

  58. Van Hoeck V, Rizos D, Gutierrez-Adan A, Pintelon I, Jorssen E, Dufort I, et al. Interaction between differential gene expression profile and phenotype in bovine blastocysts originating from oocytes exposed to elevated non-esterified fatty acid concentrations. Reprod Fertil Dev. 2015;27(2):372–84.

    Article  PubMed  CAS  Google Scholar 

  59. Funahashi H, Cantley TC, Day BN. Synchronization of meiosis in porcine oocytes by exposure to dibutyryl cyclic adenosine monophosphate improves developmental competence following in vitro fertilization. Biol Reprod. 1997;57(1):49–53.

    Article  CAS  PubMed  Google Scholar 

  60. Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009;4(1):44.

    Article  PubMed  CAS  Google Scholar 

  61. Lee W-J, Jang S-J, Lee S-C, Park J-S, Jeon R-H, Subbarao RB, et al. Selection of reference genes for quantitative real-time polymerase chain reaction in porcine embryos. Reprod Fertility Dev. 2017;29(2):357–67.

    Article  CAS  Google Scholar 

Download references


Thanks to Dominic Gagné for technical support and Sylvie Bilodeau for language editing.


This study was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC grant number 445230-12). Fonds de recherche du Québec – Nature et technologies (FRQNT) provided a fellowship to Meihong Shi.

Author information

Authors and Affiliations



Marc-André Sirard designed the study, reviewed, and edited the paper. Meihong Shi designed the study, performed the experiments, analyzed the data, and wrote the paper. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Marc-André Sirard.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

We declare that there are no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information


Additional file 1.


Additional file 2.


Additional file 3.


Additional file 4.


Additional file 5.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shi, M., Sirard, MA. Cocultured porcine granulosa cells respond to excess non-esterified fatty acids during in vitro maturation. J Ovarian Res 14, 142 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • NEFAs
  • Granulosa cells
  • IVM
  • Transcriptome