Lipidomics profiling of goose granulosa cell model of stearoyl-CoA desaturase function identifies a pattern of lipid droplets associated with follicle development
Cell & Bioscience volume 11, Article number: 95 (2021)
Despite their important functions and nearly ubiquitous presence in cells, an understanding of the biology of intracellular lipid droplets (LDs) in goose follicle development remains limited. An integrated study of lipidomic and transcriptomic analyses was performed in a cellular model of stearoyl-CoA desaturase (SCD) function, to determine the effects of intracellular LDs on follicle development in geese.
Numerous internalized LDs, which were generally spherical in shape, were dispersed throughout the cytoplasm of granulosa cells (GCs), as determined using confocal microscopy analysis, with altered SCD expression affecting LD content. GC lipidomic profiling showed that the majority of the differentially abundant lipid classes were glycerophospholipids, including PA, PC, PE, PG, PI, and PS, and glycerolipids, including DG and TG, which enriched glycerophospholipid, sphingolipid, and glycerolipid metabolisms. Furthermore, transcriptomics identified differentially expressed genes (DEGs), some of which were assigned to lipid-related Gene Ontology slim terms. More DEGs were assigned in the SCD-knockdown group than in the SCD-overexpression group. Integration of the significant differentially expressed genes and lipids based on pathway enrichment analysis identified potentially targetable pathways related to glycerolipid/glycerophospholipid metabolism.
This study demonstrated the importance of lipids in understanding follicle development, thus providing a potential foundation to decipher the underlying mechanisms of lipid-mediated follicle development.
Lipids are a highly diverse class of molecules key to many cellular processes, including energy storage, structural moieties, signaling molecules, and protein modification [1, 2]. Lipids have been loosely defined as biological substances that are generally hydrophobic in nature and, in many cases, soluble in organic solvents. These chemical properties include a broad range of molecules, including fatty acids, phospholipids, sterols, sphingolipids, terpenes, with each having a role in living systems . Lipid metabolism involves numerous interwoven pathways; these networks are essential for regulating cell function, which is influenced by many external factors. Imbalanced changes in lipid composition can even affect various physiological functions of the cell . Lipid overload results in lipotoxicity, which initially manifests as organelle dysfunction, ultimately leading to cell death and tissue damage [5, 6]. Evidently, deregulation of this network contributes to the onset of pathology as lipids are tightly regulated, both spatially and temporally, in various parts of the organism [7, 8]. Many lipid molecular species have been identified as potential biomarkers involved in several diseases, including, epidemic diseases, cancers, inflammation, and genetic diseases .
Intracellular lipid droplets (LDs) occur through the deposition of lipid complex molecular pathways. They are found in a wide variety of cell types, and are viewed as highly dynamic organelles with unique spherical structures, of which primary neutral lipids are mainly sterol esters and triacylglycerols (TGs) [10, 11]. LDs affect cell and organ health and contribute to the dynamics of biological membranes, energy storage and expenditure, and conveying messages at all cellular levels, including in the nucleus [12, 13]. Despite their important functions and nearly ubiquitous presence in cells, many aspects of LD biology are unknown. We have previously confirmed that de novo lipogenesis occurs in goose GCs  and that LD accumulation capacity depends on the stage of granulosa cells (GCs) during goose follicle development . It has been reported that lipid metabolism in bovine, sheep, and human GCs is pivotal for oocyte maturation, fertilization, and early embryonic development [16,17,18,19]. Lipid content in follicles differs among species, as reported in studies performed in porcine , bovine , and canine oocytes , wherein LD number, size, and distribution were analyzed. However, an understanding of the biology of LDs in GCs remains limited. Recently, the role of stearoyl-CoA desaturase (SCD) as a regulator of fatty acid desaturation in goose GCs was identified, whereby metabolite changes in lipid-related metabolism were found to be closely related to follicular development . Although this observation indicated that such links might exist, it is unclear whether LDs play a role in goose follicle development and which lipid species or classes might be affected. It is therefore hypothesized that lipidomic analysis of the goose GCs at the level of individual lipid molecular species might help to elucidate such a mechanism.
Lipidomics comprises the characterization of an organism’s differences in individual lipid molecular species and their biological roles with respect to the expression of proteins involved in lipid metabolism and function [24, 25]. Thus, to obtain a detailed view of lipid metabolism, it is crucial to analyze the full lipid profiles of the individual lipids using a lipidomics approach. This is the first time that a combination of lipidomics and transcriptomics analysis has been applied to the GC model of SCD function to identify alterations in the regulation of specific gene, lipid classes, and lipid-mediated signaling processes that are involved in LDs associated with follicle development in geese. The detailed characterization of lipid classes and species provided in this study might be useful to understand the mechanism of intramyocellular LDs, and an analysis based on integrating two-omics datasets will broaden our understanding of lipid alterations with respect to lipid-gene networks in GCs associated with follicle development.
LD distribution in goose GCs
Previous studies have demonstrated that de novo lipogenesis occurs in goose ovarian GCs, which have lipid accumulation capacity . To investigate the localization of LDs in goose GCs, confocal microscopy was used to analyze LD location in goose GCs stained with BODIPY (500/510, green). GCs contained numerous internalized LDs that were generally spherical in shape and dispersed throughout the cytoplasm (Fig. 1). The dispersed LDs clumped to form clusters. The mobility of the dispersed LDs was from the perinuclear regions to the extracellular domain.
Altered SCD expression affects LD content
To further investigate the effect of SCD expression on LDs in GCs, the LDs were stained with Oil Red O after SCD-overexpression and knockdown (Fig. 2a, b). The accumulation of LD content was quantified using appropriate equations (see “Materials and methods” section for details). Compared with that in the control group, the overexpression of SCD significantly stimulated the accumulation of LD content. However, SCD knockdown significantly inhibited the accumulation of LDs (Fig. 2c).
Analysis of lipid classes
In total, 662 species of lipids were identified in the GC samples using LC–MS/MS analysis (Additional file 1: Fig. S1). There were 183 different lipids identified in the control group (referred to as LN) vs. SCD-overexpression group (referred to as LS) comparison (165 upregulated and 18 downregulated) based on fold-change analysis (variable importance in projection [VIP] > 1, fold-change [FC] > 1.2), 85 lipids in the GFP vector group (referred to as LG) vs. LS comparison (36 upregulated and 49 downregulated), 90 lipids in the scrambled siRNA group (referred to as LC) vs. siRNA-210 (referred to as LT) comparison (45 upregulated and 45 downregulated), and 106 lipids in the LC vs. siRNA-470 (referred to as LF) comparison (86 upregulated and 20 downregulated). Glycerophospholipids and glycerolipids were the predominant lipid species (Fig. 3; Additional file 6: Table S1). Venn diagram analysis revealed 26 overlapping lipids in the SCD-overexpression group, 35 overlapping lipids in the SCD-knockdown group, and nine overlapping lipids in the SCD-overexpression and SCD-knockdown groups (Additional file 2: Fig. S2). As shown in Table 1 and Fig. 4, 19 lipids followed the same trend, whereas seven lipids showed the opposite pattern, as compared to that of the two overexpressing groups (LN vs. LS comparison and LG vs. LS comparison). Thirty-three lipids followed the same trend, whereas two lipids showed the opposite pattern, as compared to that of the two knockdown groups (LC vs. LT comparison and LC vs. LF comparison). Seven lipids followed the same trend, whereas two lipids showed the opposite pattern between the overexpressed and knockdown groups.
Cross comparisons between and within each group
With the aid of the KEGG database, metabolic pathways are summarized in Additional file 7: Table S2 based on the identified lipids. In total, 183 lipids were associated with 57 metabolic pathways in the LN vs. LS comparison, 85 lipids with 40 metabolic pathways in the LG vs. LS comparison, 90 lipids with 46 metabolic pathways in the LC vs. LT comparison, and 106 lipids with 42 metabolic pathways in the LC vs. LF comparison. Figures 5a, b represent the perturbed SCD-associated pathways, the results of which are as follows: for LN vs. LS, glycerophospholipid, arachidonic, sphingolipid, and glycerolipid metabolisms altered significantly; for LG vs. LS, glycerophospholipid, sphingolipid, glycerolipid, and inositol phosphate metabolisms altered significantly; for LC vs. LT, glycerophospholipid, sphingolipid, glycerolipid, and inositol phosphate metabolisms altered significantly; and for LC vs. LF, glycerophospholipid, sphingolipid, glycerolipid, and inositol phosphate metabolisms altered significantly.
RNA-seq investigation of overexpression and knockdown of SCD expression in GCs
To obtain a comprehensive view of the role of the SCD gene in GCs, the six groups of cells (SCD-overexpression group [referred to as OS], GFP vector group [referred to as OG], control group [referred to as LN], siRNA-210 group [referred to as ST], siRNA-470 group [referred to as SF], scrambled siRNA group [referred to as LC]) were subjected to RNA-seq analysis. According to the RNA-seq data, 51.56, 53.78, 55.70, 57.70, 51.65, and 56.19 million clean reads were obtained from the NC, OS, OG, SC, ST, and SF samples, respectively. For each sample, 37.69, 38.67, 40.04, 42.60, 37.88, and 41.30 million unique reads could be mapped to the current version of the goose genome (Anser cygnoides) (Table 2). The basic characteristics of the tags in each sample library are shown in Additional file 8: Table S3.
Identification and functional enrichment of DEGs
In total, 679 DEGs were identified in the GCs after SCD overexpression and knockdown (Additional file 3: Fig. S3), including 125 DEGs (41 upregulated and 84 downregulated) in the NC vs. OS comparison, 121 DEGs (66 upregulated and 55 downregulated) in the OG vs. OS comparison, 199 DEGs (57 upregulated and 142 downregulated) in the SC vs. ST comparison, and 306 DEGs (92 upregulated and 214 downregulated) in the SC vs. SF comparison (Fig. 6a, b).
To characterize differences in lipid-related genes between the transcriptomes of the SCD-overexpression and SCD-knockdown groups, the GO slim terms of the complex Gene Ontology annotations were used. The results showed that DEGs were assigned to 20 lipid-related GO terms in the NC vs. OS comparison, 21 lipid-related GO terms in the OG vs. OS comparison, 63 lipid-related GO terms in the SC vs. ST comparison, and 63 lipid-related GO terms in the SC vs. SF comparison (Fig. 7; Additional file 9: Table S4). The top 20 pathways of KEGG functional analysis showed that SCD-overexpression had a significant effect on one carbon pool by folate, glycosphingolipid biosynthesis, ferroptosis, and influenza A. Knockdown of SCD had a significant effect on sphingolipid metabolism, proteasome, mTOR signaling pathway, hedgehog signaling pathway, glycosphingolipid biosynthesis (metabolism), ether lipid metabolism, amino sugar and nucleotide sugar metabolism, and protein processing in the endoplasmic reticulum (Additional file 4: Fig. S4a, b).
Level of coordination between transcript and lipidomic data
Monitored responses in parallel on the transcriptome and the lipid level were carried out, thus allowing one to compare the level of coordination between both molecular readouts. To perform this analysis, a holistic co-clustering approach was performed, as shown in Fig. 8. Perturbed SCD in GCs influenced mRNAs and lipids in glycerolipid/glycerophospholipid metabolism, glycine, serine, and threonine metabolism, and sphingolipid metabolism, which are associated with lipid-related metabolism pathways. However, SCD knockdown in GCs influenced more complicated lipid-related metabolism pathways, including fatty acid degradation/elongation, inositol phosphate metabolism, phosphatidylinositol signaling system, alpha-linolenic acid metabolism, arachidonic acid metabolism, and linoleic acid metabolism. Additionally, a targeted approach using prior biological knowledge in conjunction with canonical correlation analysis (CCA) revealed that lipids from glycerolipid/glycerophospholipid metabolism were subjected to a CCA together with transcript data of all those pathways as derived from SCD knockdown in GCs (Fig. 9, Additional file 10: Table S5). We were also interested in other lipid-related pathways subjected to a CCA together with transcript data (Additional file 4: Fig. S4). However, no similar CCA results were observed with SCD overexpression in GCs.
In the present study, we used a GC model of SCD function and performed a detailed study on cellular lipid metabolism using lipidomic profiling to evaluate the lipid molecular changes in a detailed manner. The method is based on ultra-high performance liquid chromatography combined with quadrupole-time-of-flight mass spectrometry and gas chromatography mass spectrometry, which allow for the analysis of hundreds of molecular lipids from a single sample run . More than 660 lipids were identified, including phospholipids, ceramides, diacylglycerol and triacylglycerols, cholesterol esters, and fatty acid compositions (Fig. 3). Among these lipids, glycerophospholipids and glycerolipids were the predominant species. Additionally, altered SCD expression had the greatest impact on glycerophospholipid metabolism, this lipid-related pathway contains the most key lipids, with more than six representative metabolites as follows: PC (16:0/20:3); PC (18:1/18:4); PC (18:2/20:3); LPC (18:3/0:0); PE (22:2/12:0); and PS (20:4/20:0). FFA (22:2) and TG (14:0/22:0/22:2) were representative metabolites of glycerolipid metabolism, and sphingolipid metabolism (Figs. 4, 5).
There is great diversity in the structures and functions of the vast array of lipid species, with intracellular levels of lipids tightly regulated by a network of metabolic pathways to sustain normal cellular functions . Glycerophospholipids make up phospholipids, which play both structural and metabolic roles in cells. Therefore, the level and species of phospholipids might reflect the biochemical and physiological conditions of cells . Phospholipids are important for determining LD formation and are attributed to the formation of large LD fusions . Glycerolipid metabolism is altered by fatty acids and, induces LD formation . Blocking PC synthesis causes elevated SREBP-1-dependent gene transcription and LD accumulation in mouse liver and human cells, and the ratio of PE to PC in the LD surface tends to affect the LD number and size [31, 32]. The present study provided the lipid species composition and lipid metabolism of the GC model of SCD function. Given the important biological activities of many lipids, such information could provide potential insights into the mechanisms of internalized lipid-mediated follicle development.
We further investigated whether altered SCD expression of GCs affects gene expression involved in lipid metabolism. We observed that altered SCD expression resulted in many DEGs assigned to lipid-related GO terms, including lipid binding, lipid transport, lipid localization, and LDs. Compared with those in the SCD-overexpression group, more DEGs were identified and assigned to lipid-related GO terms in the SCD-knockdown group (Fig. 7; Additional file 9: Table S4). These results indicated that SCD regulates the lipid metabolic processes during goose follicle development. Particularly, SCD knockdown might mobilize more genes associated with lipid metabolism. Combined analysis of lipid and gene expression profiles showed that more complicated lipid-related metabolism pathways were significantly clustered, especially in the SCD-knockdown group.
Numerous internalized LDs were dispersed throughout the cytoplasm of GCs (Figs. 1, 2). However, in steroidogenic tissues such as the ovary, the ovarian follicles of almost all mammals contain LDs, which are prominently and generally non-pathogenic, except in cases of genetic disorders of steroidogenesis . We speculate that there is a relationship between LDs and energy supply or the storage of cholesteryl esters during goose follicle development. Cholesterol esters, particularly CE (20:4), were significantly upregulated in SCD-overexpression and knockdown groups (Fig. 4c). Previous results demonstrated intracellular metabolite changes in which cholesterol was mostly upregulated in the SCD-overexpression group, whereas androsterone was upregulated in the SCD-knockdown group .
LDs make dynamic contact with nearly every other organelle. Some contact activity might exist among the ER, endosomes, mitochondria, and peroxisomes in the cell. This contact is believed to mediate the transfer of lipids (including phospholipids, sterols, and fatty acids) between the different cellular compartments [34, 35]. Pharmacological inhibition of the enzymatic activity of fatty acid synthase significantly reduces the production of progesterone in bovine GCs, suggesting that ovarian steroidogenesis relies on lipid metabolism in GCs [16, 36]. Additionally, LD accumulation capacity of goose GCs depends on the different stages of follicle development, with the stage of the highest accumulation capacity being F1 GCs, followed by F4–F2 GCs and pre-hierarchical GCs . As ovarian follicles grow with changes in both the size and number of LDs in bovines , the total intracellular LD content increases gradually with follicular growth. This exponential increase was demonstrated in large antral follicles in rats and goats [38, 39], with an increase in LD numbers and size in GCs within one day after ovulation in rabbits . These results confirm that LDs are considered a source of energy for oocyte maturation [41, 42].
SCD catalyzes the rate-limiting step in the production of monounsaturated fatty acids, which are major substrates for the synthesis of triglycerides, wax esters, cholesterol esters, and phospholipids . Recent research has revealed that alterations in SCD expression change the fatty acid profile of these lipids and have a significant role in lipid metabolism . The integrated optical density of the LDs is related to the amount of triglycerides in the LDs. Our results showed that altered SCD expression affected LDs in GCs. The overexpression of SCD stimulated the accumulation of LD content, whereas SCD knockdown inhibited the accumulation of LD content (Fig. 2c). Previous research has demonstrated that SCD is an important control point in lipid homeostasis in goose GCs . Similar results suggest that increased hepatic SCD activity accelerates the last stage of TG synthesis and stimulates lipid accumulation . These results provide direct evidence that SCD expression plays a major role in the triacylglycerol and phospholipid secretion processes. Much evidence indicates that the direct anti-steatotic effect of SCD deficiency stems from increased fatty acid oxidation and reduced lipid synthesis  and that the inhibition of SCD decreases the proliferation and apoptosis of cells . Additionally, a link between SCD activity and LD size was previously observed in the nematode Caenorhabditis elegans ; and experiments on cell lines cultured from patients with Berardinelli-Seip congenital lipodystrophy  demonstrated that SCD activity is necessary for large-sized LDs. Recent studies have indicated that the size of LDs might influence the balance of lipogenesis and lipolysis in adipocytes and hepatic cells , with decreased SCD expression resulting in reduced accumulation of oleic and palmitoleic acids and consequently fat cells with small LDs . However, the functional implications of the diverse sizes of LDs in goose GCs are poorly understood.
In summary, these results provide new evidence that altered SCD expression changes GC lipid profiles, genes, and regulated relevant metabolic pathways, which might implicate SCD as one of the major genes in the regulation of lipid secretion by goose GCs.
Animals and granulosa layer isolation
Geese (from a maternal line of Tianfu goose) were raised under natural temperature and light conditions at the experimental station of waterfowl breeding at Sichuan Agricultural University. Six geese with regular laying sequences were randomly selected as experimental samples and sacrificed by cervical bleeding under anesthesia, 2 h after oviposition. A pool of ovarian follicles was immediately collected from the goose abdominal cavities and divided into prehierarchal follicles and hierarchal preovulatory follicles according to previously reported nomenclature . For granulosa layer collection, the outer connective tissue was removed from the preovulatory follicles, and the preovulatory follicles were dissected to allow the yolk and adhered granulosa layer to flow out .
Cell culture and transfection
The granulosa layer was dispersed by incubation in 0.1% type II collagenase (Sigma, USA) in Dulbecco’s modified Eagle’s medium (DMEM, HyClone, USA) for 10 min in a 37 °C water bath. After incubation, cells were dispersed with a pipette and pelleted by centrifugation at 1000×g for 10 min (20 °C). The supernatant was discarded and the cells were resuspended in 3 mL of fresh basic medium without collagenase and centrifuged. The washing procedure was repeated twice. The GCs were dispersed in DMEM supplemented with 1% antibiotic/antimycotic solution (Solarbio, China) and 3% fetal bovine serum (Gibco, USA). Transient transfections based on the GC cellular model of SCD function (SCD-specific overexpression and knockdown) were performed using Lipofectamine® 3000 and Lipofectamine™ RNAiMAX Transfection Reagent (Invitrogen), according to manufacturer’s protocol. SCD siRNA-210 and siRNA-405 were used to achieve SCD mRNA knockdown, with scrambled siRNA as a negative control. The extent of SCD mRNA knockdown was measured as a percentage compared to that the scrambled siRNA. SCD-specific overexpression was used to achieve SCD mRNA overexpression, termed pEGFP-N1/SCD, and the pEGFP-N1/empty vector (GFP vector) served as the negative control, and another control with no transfection components was also included. The extent of SCD mRNA overexpression was measured as a percentage compared to that with the GFP vector and no transfection components conditions. Transfection efficiency was determined as previously described .
Oil red O staining
To determine LD accumulation in GCs after transfection, the cells were stained with Oil Red O solution (Sigma Chemical Co., St. Louis, MO, USA). After 48 h of transfection, the cultured GCs were washed with phosphate-buffered saline (PBS) and fixed with 10% formaldehyde at room temperature. Then, the cells were washed twice with PBS, and were subsequently stained with 0.5 μg/mL Oil Red O solution in the dark for 1 h at room temperature. Then, cells were washed with PBS and photographed using an optical microscope system (Nikon ECLIPSE 90i) at 200 × magnification. Finally, the LDs were dissolved in isopropanol and the optical density (OD) value was determined at a wavelength of 540 nm using a microplate reader (Thermo Varioskan, USA). The relative lipid content was calculated using the following equation: relative LD content (%) = sample OD/mg × mL−1 protein.
Immunofluorescence staining and confocal microscopy
For the photographic documentation of LD localization in goose GCs, live cells were stained with 4,4-difluoro-3a,4a-diaza-s-indacene (Bodipy 500/510, Thermo Fisher), wheat germ agglutinin (WGA, Alexa Fluor 594 conjugate, Invitrogen), and Hoechst dyes (33258, Invitrogen) for labeling the LDs, membrane, and nucleus, respectively. Briefly, cells on culture plates were rinsed with PBS and stained for 5 min with 2 μg/mL Bodipy 500/510. After staining, cells were washed three times with PBS for 5 min. The cells were then stained for 5 min with 5 μg/mL WGA and 5 μg/mL Hoechst 33,342. After staining and incubation, the cells were washed thrice with PBS for 5 min each. Finally, cells were fixed in formaldehyde for 10 min. The fixed cells were imaged using a confocal laser scanning microscope (FV1200, Olympus) mounted on an IX83 inverted microscope (Olympus). For treble-color images, a 550-nm laser line was used to image cells stained with Bodipy 500/510, a 635-nm laser line was used to image cells stained with WGA, and a 543-nm laser line was used to image cells stained with Hoechst 33258.
Sample preparation and extraction
For lipidomic analysis, each group of cells (LS, LG, LN, and two independent siRNA groups [siRNA-210 and siRNA-470, referred to as LT and LF, respectively], as well as the LC) were placed in liquid nitrogen for 2 min, then thawed on ice for 5 min and vortexed. The cells were centrifuged at 12,000 rpm at 4 °C for 10 min. The supernatants (300 µL) were homogenized with a 1 mL mixture containing methanol, MTBE, and an internal standard. This mixture was stirred for 2 min, followed by the addition of 500 µL water, stirring for 1 min, and centrifugation at 12,000 rpm at 4 °C for 10 min. A total of 500 µL of the supernatant obtained was extracted and concentrated. The powder was dissolved with 100 µL mobile phase B (0.1% acetic acid), and the supernatant (200 μL) was transferred to an LC–MS sampling vial with an inner liner, for LC–MS analysis. Quality control (QC) samples were produced by pooling equal aliquots taken from each individual sample in the analytical run together. These QC samples were used to monitor the stability and reproducibility of the analytes in the samples during the analysis.
For RNA-sequencing analysis, each group of cells was extracted from three independent biological replicates using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instructions. RNA quality was confirmed on 1% agarose gels. RNA purity and concentration were determined using a NanoPhotometer® spectrophotometer (IMPLEN, CA, USA) and Qubit® RNA Assay Kit with a Qubit® 2.0 Fluorometer (Life Technologies, CA, USA) according to the manufacturer’s instructions. RNA integrity was assessed using an RNA 6000 Nano Assay Kit (Bioanalyzer 2100 system; Agilent Technologies, CA, USA).
Lipidomics profiling was performed by MetWare (Wuhan, China) using an LC–ESI–MS/MS system (UPLC, Shim-pack UFLC SHIMADZU CBM A system, https://www.shimadzu.com/; MS, QTRAP® System, https://sciex.com/). LIT and triple quadrupole (QQQ) scans were acquired on a triple quadrupole-linear ion trap mass spectrometer (QTRAP®) LC–MS/MS System, equipped with an ESI Turbo Ion-Spray interface, operating in positive and negative ion mode, and controlled by Analyst 1.6.3 software (Sciex). Qualitative analysis of primary and secondary MS data was carried out by a comparison of the accurate precursor ions (Q1), product ions (Q3), retention time (RT), and fragmentation patterns with those obtained by injecting standards at the same conditions if the standards were available (Sigma-Aldrich, USA), or was conducted using the self-compiled database MWDB (MetWare, China). The quantitative analysis of metabolites was based on the MRM mode. The characteristic ions of each metabolite were screened using a QQQ mass spectrometer to obtain the signal strengths. Integration and correction of chromatographic peaks was performed using Progenesis QI software (Waters Co., Milford, MA, USA). The corresponding relative metabolite contents were represented as chromatographic peak area integrals. In addition, accurate masses of features representing significant differences were searched against the METLIN, Kyoto Encyclopedia of Genes and Genomes (KEGG), HMDB, and LIPIDMAPS databases.
Sequencing libraries were generated by Novogene (Beijing, China) using a NEBNext® Ultra™ RNA Library Prep Kit (Illumina®, NEB, USA) following the manufacturer’s instructions. Clustering of the index-coded samples was performed on a cBot Cluster Generation System by using a TruSeq PE Cluster Kit v3-cBot-HS (Illumina) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina HiSeq (Illumina, USA) platform and 125 bp/150 bp paired-end reads were generated.
Bioinformatic and statistical analyses of data
Raw LC–MS data were filtered, identified, integrated, corrected, aligned, and normalized using Progenesis QI software (Waters Co., Milford, MA, USA). A data matrix of RT, mass-to-charge ratio, and peak intensity was obtained. The processed data were analyzed using principal component analysis and orthogonal correction partial least squares discriminant analysis (PC) using SIMCA-P14.0 software (Umetrics, Umeå, Sweden). Differentially abundant metabolites between dietary treatments were identified from variable importance in projection (VIP) from OPLS-DA and Student’s t tests (VIP > 1 and P < 0.05). Metabolites were identified from public databases, including MassBank (http://www.massbank.jp/), KNApSAcK (http://kanaya.naist.jp/KNApSAcK/), Human Metabolome Database (http://www.hmdb.ca/), and METLIN (https://metlin.scripps.edu/). The KEGG database (http://www.genome.jp/kegg/) was used to view the enriched pathways of the different metabolites. Hierarchical clustering analysis and heat map analysis were conducted using the R package, version 3.3.1. The raw read counts were normalized considering both the different depths of sequencing among the samples and the gene GC content. Normalization was performed using the EDASeq package. We considered all DEGs with a P value < 0.05 after false discovery rate correction. The putative functions of DEGs were investigated by Gene Ontology (GO) enrichment analysis using the R package GOseq , in which gene length bias was corrected. GO terms with a corrected P value < 0.05 were considered significantly enriched. KOBAS software was used to test the statistical enrichment of DEGs in KEGG pathways . Statistical plots were prepared using Origin version 6.1.
Availability of data and materials
RNA-seq data used in the current study are available at NCBI Bioproject database under the accession number PRJNA720652.
Intracellular lipid droplets
Canonical correlation analysis
Differentially expressed genes
Dulbecco’s modified Eagle’s medium
Wheat germ agglutinin
Kyoto Encyclopedia of Genes and Genomes
Importance in projection
Storck EM, Ozbalci C, Eggert US. Lipid cell biology: a focus on lipids in cell division. Annu Rev Biochem. 2018;87(1):839–69.
Santos AL, Preta G. Lipids in the cell: organisation regulates function. Cell Mol Life Sci CMLS. 2018;75(11):1909–27.
Brown HA, Marnett LJ. Introduction to lipid biochemistry, metabolism, and signaling. Chem Rev. 2011;111(10):5817–20.
Bochkov VN, Oskolkova OV, Birukov KG, Levonen AL, Binder CJ, Stöckl J. Generation and biological activities of oxidized phospholipids. Antioxidants Redox Signaling. 2010;12(8):1009–59.
Ozcan U, Yilmaz E, Ozcan L, Furuhashi M, Vaillancourt E, Smith RO, Görgün CZ, Hotamisligil GS. Chemical chaperones reduce ER stress and restore glucose homeostasis in a mouse model of type 2 diabetes. Science (New York, NY). 2006;313(5790):1137–40.
Feldstein AE, Werneburg NW, Canbay A, Guicciardi ME, Bronk SF, Rydzewski R, Burgart LJ, Gores GJ. Free fatty acids promote hepatic lipotoxicity by stimulating TNF-alpha expression via a lysosomal pathway. Hepatology (Baltimore, MD). 2004;40(1):185–94.
Goldberg IJ, Karen R, Abumrad NA, Bickel PE, Sarah C, Fisher EA, Galis ZS, Granneman JG, Douglas LE, Robert M. Deciphering the role of lipid droplets in cardiovascular disease: a report from the 2017 National Heart, Lung, and Blood Institute Workshop. Circulation. 2018;138(3):305–15.
Liu Z, Xun H. Lipid metabolism in Drosophila: development and disease. Acta Biochim Biophys Sin. 2013;45:44–50.
Lydic TA, Goo YH. Lipidomics unveils the complexity of the lipidome in metabolic diseases. Clin Transl Med. 2018;7(1):4.
Olzmann JA, Carvalho P. Dynamics and functions of lipid droplets. Nat Rev Mol Cell Biol. 2019;20(3):137–55.
Den Brok MH, Raaijmakers TK, Colladocamps E, Adema GJ. Lipid droplets as immune modulators in myeloid cells. Trends Immunol. 2018;39(5):380–92.
Onal G, Kutlu O, Gozuacik D, Dokmeci ES. Lipid droplets in health and disease. Lipids Health Dis. 2017;16(1):128.
Mike WH, Reese ML, Goodman JM. The assembly of lipid droplets and their roles in challenged cells. EMBO J. 2019.
Wen R, Gan X, Hu S, Gao S, Deng Y, Qiu J, Sun W, Li L, Han C, Hu J. Evidence for the existence of de novo lipogenesis in goose granulosa cells. Poult Sci. 2019;98(2):1023–30.
Gao S, Gan X, He H, Hu S, Deng Y, Chen X, Li L, Hu J, Li L, Wang J. Dynamic characteristics of lipid metabolism in cultured granulosa cells from geese follicles at different developmental stages. Biosci Rep. 2019;39(12).
Elis S, Desmarchais A, Maillard V, Uzbekova S, Monget P, Dupont J. Cell proliferation and progesterone synthesis depend on lipid metabolism in bovine granulosa cells. Theriogenology. 2015;83(5):840–53.
Campbell BK, Onions V, Kendall NR, Guo L, Scaramuzzi RJ. The effect of monosaccharide sugars and pyruvate on the differentiation and metabolism of sheep granulosa cells in vitro. Reproduction (Cambridge, England). 2010;140(4):541–50.
Hu W, Qiao J. Expression and regulation of adipocyte fatty acid binding protein in granulosa cells and its relation with clinical characteristics of polycystic ovary syndrome. Endocrine. 2011;40(2):196–202.
Mcevoy TG, Coull GD, Broadbent PJ, Hutchinson JS, Speake BK. Fatty acid composition of lipids in immature cattle, pig and sheep oocytes with intact zona pellucida. J Reprod Fertil. 2000;118(1):163–70.
Campagna C, Bailey JL, Sirard M-A, Ayotte P, Maddox-Hyttel P. An environmentally-relevant mixture of organochlorines and its vehicle control, dimethylsulfoxide, induce ultrastructural alterations in porcine oocytes. Mol Reprod Dev. 2010;73(1):83–91.
Abe H, Yamashita S, Satoh T, Hoshi H. Accumulation of cytoplasmic lipid droplets in bovine embryos and cryotolerance of embryos developed in different culture systems using serum-free or serum-containing media. Mol Reprod Dev. 2002;61:57–66.
De Lesegno CV, Reynaud K, Pechoux C, Thoumire S, Chastant-Maillard S. Ultrastructure of canine oocytes during in vivo maturation. Mol Reprod Dev. 2008;75(1):115–25.
Yuan X, Hu S, Li L, Liu H, He H, Wang J. Metabolomic analysis of SCD during goose follicular development: implications for lipid metabolism. Genes (Basel). 2020;11(9):1001.
Liu H. Special topic: metabolomics and lipidomics. J Anal Testing. 2017;1(3):185–6.
Lange M, Ni Z, Criscuolo A, Fedorova M. Liquid chromatography techniques in lipidomics research. Chromatographia. 2018;82(2):77–100.
Chen W, Gong L, Guo Z, Wang W, Zhang H, Liu X, Yu S, Xiong L, Luo J. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: application in the study of rice metabolomics. Mol Plant. 2013;6(6):1769–80.
Tang QQ. Lipid metabolism and diseases 科学通报: 英文版. 2016;61(19):1471.
Dowhan W. The role of phospholipids in cell function. In: Gross RW, editor. Advances in lipobiology (JAI), vol. 2. Elsevier: Amsterdam; 1997. p. 79–107.
Han L, Pang K, Li XI, Loor JJ, Yang GY, Gao T. Lipidomic profiling analysis of the phospholipid molecules in SCAP-induced lipid droplet formation in bovine mammary epithelial cells. Prostaglandins Other Lipid Media. 2020;149:106420.
Pathmaperuma AN, Maña P, Cheung SN, Kugathas K, Josiah A, Koina ME, Broomfield A, Delghingaro-Augusto V, Ellwood DA, Dahlstrom JE, et al. Fatty acids alter glycerolipid metabolism and induce lipid droplet formation, syncytialisation and cytokine production in human trophoblasts with minimal glucose effect or interaction. Placenta. 2010;31(3):230–9.
Walker AK, Jacobs RL, Watts JL, Rottiers V, Jiang K, Finnegan DM, Shioda T, Hansen M, Yang F, Niebergall L. A conserved SREBP-1/phosphatidylcholine feedback circuit regulates lipogenesis in metazoans. Cell. 2011;147(4):840–52.
Listenberger L, Townsend E, Rickertsen C, Hains A, Brown E, Inwards EG, Stoeckman AK, Matis MP, Sampathkumar RS, Osna NA, et al. Decreasing phosphatidylcholine on the surface of the lipid droplet correlates with altered protein binding and steatosis. Cells. 2018;7(12):230.
Paulini F, Silva RC, Rolo JLJdP, Lucci CM. Ultrastructural changes in oocytes during folliculogenesis in domestic mammals. J Ovarian Res. 2014;7:102.
Gao G, Chen FJ, Zhou L, Su L, Xu D, Xu L, Li P. Control of lipid droplet fusion and growth by CIDE family proteins. Biochim Biophys Acta Mol Cell Biol Lipids. 2017;1862:1197–204.
Gallardo-Montejano VI, Saxena G, Kusminski CM, Yang C, Mcafee JL, Hahner L, Hoch K, Dubinsky W, Narkar VA, Bickel PE. Nuclear Perilipin 5 integrates lipid droplet lipolysis with PGC-1α/SIRT1-dependent transcriptional regulation of mitochondrial function. Nat Commun. 2016;7:12723.
Yenuganti VR, Viergutz T, Vanselow J. Oleic acid induces specific alterations in the morphology, gene expression and steroid hormone production of cultured bovine granulosa cells. Gen Comp Endocrinol. 2016;232:134–44.
Dadarwal D, Adams GP, Hyttel P, Brogliatti GM, Caldwell S, Singh J. Organelle reorganization in bovine oocytes during dominant follicle growth and regression. Reprod Biol Endocrinol. 2015;13:124–124.
Sangha GK, Guraya SS. Biochemical changes in lipids during follicular growth and corpora lutea formation and regression in rat ovary. Indian J Exp Biol. 1989;27(11):998–1000.
Sharma RK, Vats R, Sawhney AK. Changes in follicular lipids during follicular growth in the goat (Capra hircus) ovary. Small Ruminant Res. 1996;20(2):177–80.
Guraya SS. A histochemical study of pre-ovulatory and post-ovulatory follicles in the rabbit ovary. Reprod Fertil. 1968;15(3):381–7.
Silva RC, Báo SN, Jivago JLPR, Lucci CM. Ultrastructural characterization of porcine oocytes and adjacent follicular cells during follicle development: lipid component evolution. Theriogenology. 2011;76(9):1647–57.
Sturmey RG, O’Toole PJ, Leese HJ. Fluorescence resonance energy transfer analysis of mitochondrial:lipid association in the porcine oocyte. Reproduction (Cambridge, England). 2006;132(6):829–37.
Khor VK, Shen W, Kraemer FB. Lipid droplet metabolism. Curr Opin Clin Nutr Metab Care. 2013;16(6):632–7.
Matthew T, Flowers, James M, Ntambi. Role of stearoyl-coenzyme A desaturase in regulating lipid metabolism. Curr Opin Lipidol. 2008.
Dobrzyn A, Dobrzyn P, Lee SH, Miyazaki M, Cohen P, Asilmaz E, Hardie DG, Friedman JM, Ntambi JM. Stearoyl-CoA desaturase-1 deficiency reduces ceramide synthesis by downregulating serine palmitoyltransferase and increasing beta-oxidation in skeletal muscle. Am J Physiol Endocrinol Metab. 2005;288(3):E599.
Iwai T, Kume S, Chin-Kanasaki M, Kuwagata S, Araki H, Takeda N, Sugaya T, Uzu T, Maegawa H, Araki SI. Stearoyl-CoA desaturase-1 protects cells against lipotoxicity-mediated apoptosis in proximal tubular cells. Int J Mol Sci. 2016;17(11):1868.
Xun S, Juan L, Xiaoju Z, Joel G, Rodkaer SV, Faergeman NJ, Bin L, Watts JL. Regulation of lipid droplet size and phospholipid composition by stearoyl-CoA desaturase. J Lipid Res. 2013;54(9):2504–14.
Boutet E, El Mourabit H, Prot M, Nemani M, Khallouf E, Colard O, Maurice M, Durand-Schneider A-M, Chrétien Y, Grès S, et al. Seipin deficiency alters fatty acid Delta9 desaturation and lipid droplet formation in Berardinelli-Seip congenital lipodystrophy. Biochimie. 2009;91:796–803.
Welte MA, Gould AP. Lipid droplet functions beyond energy storage. Biochim Biophys Acta. 2017;1862(10 Part B):1260–72.
Choi Y, Kim YC, Han YB, Park Y, Pariza MW, Ntambi JM. The trans-10, cis-12 isomer of conjugated linoleic acid downregulates stearoyl-CoA desaturase 1 gene expression in 3T3-L1 adipocytes. J Nutr. 2000;130(8):1920–4.
Gilbert A. The Ovary. In: Bell, D. J, and Freeman, BM eds. Physiology and Biochemistry of the Domestic Fowl. Vol. 3. Academic Press, London, New York; 1971
Gilbert A, Evans A, Perry M, Davidson M. A method for separating the granulosa cells, the basal lamina and the theca of the preovulatory ovarian follicle of the domestic fowl (Gallus domesticus). J Reprod Fertil. 1977;50(1):179–81.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):1–12.
Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787.
This study was supported by a grant from the National Natural Science Foundation of China (No. 31672424 and 035Z2425), the China Agricultural Research System (No. CARS-42-4), the Open Fund of Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province (SNDK-KF-201802), the Project of National Science and Technology Plan for the Rural Development in China (No. 2015BAD03B06), and the Key Technology Support Program of Sichuan Province (No. 2016NYZ0044).
Ethics approval and consent to participate
The animal experiments were reviewed and approved by the Committee of the School of Farm Animal Genetic Resources Exploration and Innovation Key Laboratory, College of Animal Science and Technology, Sichuan Agricultural University, China (Permit No. DKY20170913). The experiments were performed in accordance with the Regulations for the Administration of Affairs Concerning Experimental Animals (China 1988). This study was conducted in accordance with the requirements of the Beijing Animal Welfare Committee. All efforts were made to minimize animal suffering.
Consent for publication
All co-authors agreed to publish the final version of the present manuscript.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Figure S1. Unsupervised hierarchical clustering of all lipids in each group; each column denotes one group. Increased lipids concentrations are shown in red, whereas decreased lipids concentrations are shown in blue.
: Figure S2. Venn diagram of overlapping and unique of lipids altered in each group.
: Figure S3. Unsupervised hierarchical clustering of all DEGs in each group; each column denotes one group. Increased DEGs concentrations are shown in red, whereas decreased DEGs concentrations are shown in blue.
: Figure S4. Topology analysis of top 20 pathways identified of transcriptomics. (a) Overexpressed group, NC vs. OS comparison and OG vs. OS comparison; (b) Knockdown group, SC vs. ST comparison and SC vs. SF comparison. Advanced bubble chart shows the enrichment of differentially abundant DEGs in pathways. The x-axis represents the rich factor (rich factor = number of different DEGs enriched in the pathway/number of all DEGs in the background DEGs set). The y-axis represents represent the enriched pathways. Size of the bubble represent the number of different abundant DEGs enriched in the pathway, and the color represents enrichment significance.
: Figure S5. Visualization of the CCA results of lipids and DEGs involved in (a) inositol phosphate metabolism, (b) arachidonic acid metabolism, (c) sphingolipid metabolism, (d) phosphatidylinositol signaling system and (e) glycine, serine and threonine metabolism under SCD knockdown.
. Significant differential lipids involved in LN vs. LS comparison, LG vs. LS comparison, LC vs. LT comparison, and LC vs. LF comparison.
. A functional analysis of pathways related to the differentially abundant lipids.
. Basic characteristics of tags in each sample libraries and data of sequencing reads mapping to the reference genome.
. The detail of DEGs were assigned with lipid-related GO slim terms.
. The detail of lipids from glycerolipid/glycerophospholipid metabolism were subjected to a CCA together with transcript data of all those pathways as derived from SCD knockdown in GCs.
About this article
Cite this article
Yuan, X., Hu, S., Li, L. et al. Lipidomics profiling of goose granulosa cell model of stearoyl-CoA desaturase function identifies a pattern of lipid droplets associated with follicle development. Cell Biosci 11, 95 (2021). https://doi.org/10.1186/s13578-021-00604-6
- Gene expression
- Goose granulosa cells
- Ovarian follicles
- Lipid droplets