Cocultured porcine granulosa cells respond to excess non-esterified fatty acids during in vitro maturation
Journal of Ovarian Research volume 14, Article number: 142 (2021)
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 . 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 . High levels of NEFAs during oocyte maturation and early embryo development have hazardous effects on development competence by affecting metabolism and cell fate . 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 .
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 . The study of Pawlak et al.  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 . 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.
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).
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 , 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).
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.
The KEGG analysis showed that differentially expressed genes were markedly enriched for  metabolism, including carbohydrate, amino acid, and antibiotics;  inflammatory responses, including “Toll-like receptor signaling pathway” and “Rheumatoid arthritis”;  signal transduction, including “TGF-beta signaling pathway”, “mTOR signaling pathway” and “Thyroid hormone signaling pathway”; and  catabolism, including “Peroxisome” and “Endocytosis” (Fig. 5).
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.
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  (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 . 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 . 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 . 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 .
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 . Besides, increased OA concentrations decreased the production of estradiol in bovine GCs . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . Instead of a lipotoxic impact, unsaturated fatty acids generally have a moderate impact on somatic cells . Lipid storage and the β-oxidation of NEFAs in GCs  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.
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 . 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  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.
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.
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 . GO and pathway enrichment analyses were performed using the DAVID Database (V6.8) with Sus scrofa background . Data were also analyzed using Ingenuity Pathway Analysis (IPA, Ingenuity® Systems, www.ingenuity.com/) 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), , 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.
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
In Vitro Maturation
Reactive Oxygen Species
North Carolina State University
Tyrode’s Lactate-buffered HEPES supplemented with 0.1 % Polyvinyl Alcohol
Dulbecco’s Minimal Essential Medium/Ham’s F12
Epidermal Growth Factor
Differentially Expressed Genes
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
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
Hepatocyte Growth Factor
Transforming Growth Factor-beta
Specificity Protein 1
Jagged Canonical Notch Ligand 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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Chaput C, Sirard M. Embryonic response to high beta-hydroxybutyrate (BHB) levels in postpartum dairy cows. Domest Anim Endocrinol. 2020;72:106431.
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.
Baddela VS, Sharma A, Vanselow J. Non-esterified fatty acids in the ovary: friends or foes? Reproductive Biology Endocrinology. 2020;18:1–14.
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.
Silvestris E, de Pergola G, Rosania R, Loverro G. Obesity as disruptor of the female fertility. Reproductive Biology Endocrinology. 2018;16(1):1–13.
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.
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.
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.
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.
Vanselow J, Conley AJ, Berger T. Aromatase and the three little paralogs. Biol Reprod. 2021;105(1):5–6.
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.
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.
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.
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.
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-.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009;4(1):44.
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.
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.
Ethics approval and consent to participate
Consent for publication
We declare that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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). https://doi.org/10.1186/s13048-021-00904-y
- Granulosa cells