TMT-based proteomic and bioinformatic analyses of human granulosa cells from obese and normal-weight female subjects
Reproductive Biology and Endocrinology volume 19, Article number: 75 (2021)
Increasing evidence supports a relationship between obesity and either infertility or subfertility in women. Most previous omics studies were focused on determining if the serum and follicular fluid expression profiles of subjects afflicted with both obesity-related infertility and polycystic ovary syndrome (PCOS) are different than those in normal healthy controls. As granulosa cells (GCs) are essential for oocyte development and fertility, we determined here if the protein expression profiles in the GCs from obese subjects are different than those in their normal-weight counterpart.
GC samples were collected from obese female subjects (n = 14) and normal-weight female subjects (n = 12) who were infertile and underwent in vitro fertilization (IVF) treatment due to tubal pathology. A quantitative approach including tandem mass tag labeling and liquid chromatography tandem mass spectrometry (TMT) was employed to identify differentially expressed proteins. Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were then conducted to interrogate the functions and pathways of identified proteins. Clinical, hormonal, and biochemical parameters were also analyzed in both groups.
A total of 228 differentially expressed proteins were noted, including 138 that were upregulated whereas 90 others were downregulated. Significant pathways and GO terms associated with protein expression changes were also identified, especially within the mitochondrial electron transport chain. The levels of free fatty acids in both the serum and follicular fluid of obese subjects were significantly higher than those in matched normal-weight subjects.
In GCs obtained from obese subjects, their mitochondria were damaged and the endoplasmic reticulum stress response was accompanied by dysregulated hormonal synthesis whereas none of these changes occurred in normal-weight subjects. These alterations may be related to the high FFA and TG levels detected in human follicular fluid.
The prevalence of obesity is increasing all over the world, in consistence with the increase in obese subjects suffering from infertility . Compared to age-matched normal-weight women, obese women have an increased risk of ovulatory subfertility and anovulatory infertility [2,3,4]. While ovarian stimulation has been used to overcome anovulation, obese subjects usually require higher doses of gonadotropins. Decreased oocyte retrieval, lower oocyte quality, and reduced rates of pre-implantation embryo development are generally observed in obese subjects, in addition to an increased risk of cycle cancellation or miscarriage compared with normal-weight subjects [5,6,7,8,9]. A considerable number of clinical reports have highlighted the negative effect of body mass index (BMI) on the success rate of fertility treatments. However, multiple studies reported no significant association between BMI and assisted reproductive outcomes [10,11,12]. This controversy indicates that the effect of obesity on oocyte quality and fertility outcome is somewhat complex, and underlying mechanisms remain to be elucidated .
Folliculogenesis requires a carefully orchestrated crosstalk between the oocyte and surrounding somatic cells . Two granulosa cell types can be distinguished in the follicle, namely mural granulosa cells (GCs), lying on the basal membrane of the follicular wall, and cumulus GCs, which surround the oocyte . These different GCs provide various molecules, such as sugars, amino acids, signaling factors, and nucleotides, to facilitate oocyte development , maturation, meiosis, and oocyte cytoskeleton alteration . Hence, a better understanding of obesity-associated changes in GCs could allow us to identify the mechanisms underlying obesity-induced female infertility.
In the research field of obesity-related female infertility, omics studies were mostly limited to the analysis of follicular fluid or plasma in polycystic ovary syndrome (PCOS) [18, 19]. However, PCOS is a complex multigenic disorder, which is strongly influenced by various epigenetic and environmental factors . It is therefore difficult to determine the effect of obesity on intrafollicular cell function. The aim of the current study was to investigate how protein expression changes in the GCs of obese women in which PCOS was ruled out as a possible confounder. The results strongly indicate that the GC proteomic expression patterns undergo significant alteration in the obese subjects.
This study enrolled subjects visiting the Reproductive Medical Center of Shanghai Ruijin Hospital in vitro fertilization (IVF) Clinics between September 2019 and January 2020. The study was approved by the local ethics committee of Shanghai Jiao Tong University School of Medicine. All subjects gave informed consent before participation. Twenty-six women with normal menstruation cycles were recruited. The Chinese-specific cut-offs for general adiposity were used, with normal weight as BMI 18.5-23.9 kg/m2 and general obesity as BMI values equal to or more than 28.0 kg/m2 . Therefore, subjects with a normal BMI (20.38‒21.02 kg/m2) (n = 12) were assigned to the control group, and obese subjects with a BMI of 29.25–32.33 kg/m2 (n = 14) were assigned to the obesity group. The ages of the control group and the obesity group were 31.00 ± 1.86 years and 31.79 ± 1.37 years respectively. All participants were Chinese and experienced infertility due to tubal pathology. Subjects with endometriosis, PCOS, or other medical disorders that could affect folliculogenesis were excluded.
Purification of GCs from follicular fluid
All subjects underwent controlled ovarian stimulation for IVF procedures using standard gonadotropin-releasing hormone antagonist protocols. After receiving human chorionic gonadotrophin injections for 34‒36 h, follicular fluid was obtained from follicles (16-20 cm in diameter) that contained granulosa cells (GCs). The purification procedure details are as follow: when the mature follicles were obtainable, the cumulus-oocyte complexes were separated from the follicles for IVF treatment, and the remaining follicular fluid was collected and centrifuged at 250×g for 10 min. Then the resulting pellet was resuspended with 1X PBS and subsequently layered on 40 % Percoll® (Sigma-Aldrich; Merck KGaA) followed by centrifugation at 450×g for 20 min at room temperature. After performing this density gradient centrifugation, the mixture was divided into four layers from the button to the top surface for separating red blood cells, Percoll, granulosa cells and follicular fluid with PBS. The granulosa cell layer was collected and washed with 1 × PBS at 250×g for 10 min. The pellet was resuspended with 1X PBS,and then the red blood cell lysate was added to a volumetric ratio of 1:3 and kept at 4℃ for 5 min. Finally, three wash steps were carried out using 1X PBS at 250×g for 1 min at room temperature, and purified GCs were collected.
Total protein extraction
In order to establish protein content equivalence of 120 µg in each sample, the GCs of a single subject whose protein content was less than 120 µg were pooled with other GCs belonging to one or two subjects in the same group. This procedure generated 5 control and 5 obese samples in each group. Samples were ground individually in liquid nitrogen and lysed with lysis buffer, followed by ultrasonication on ice for 5 min. The mixture was then centrifuged at 12,000 × g for 15 min at 4 ℃, and the supernatant was collected. Extracts were reduced with 10 mM DTT for 1 h at 56 ℃ individually and were alkylated with iodoacetamide in the absence of light for 1 h at room temperature. The samples were then mixed with precooled acetone, followed by incubation and centrifugation. Next, the precipitate was collected, washed with cold acetone, and then the pellet was dissolved in dissolution buffer containing 0.1 M triethylammonium bicarbonate (TEAB, pH 8.5) and 6 M urea. Finally, the Bradford Protein Assay Kit (Beyotime Biotechnology, China) was used to determine the protein concentrations.
Protein digestion and tandem mass tag (TMT) labeling
A total of 120 µg per protein sample was added to 100 µL of dissolution buffer. Thereafter, 1.5 µg trypsin and 500 µL of 100 mM TEAB buffer were added, followed by digestion at 37 °C for 4 h. Subsequently, 1.5 µg trypsin and CaCl2 were added, and the sample was digested overnight. Formic acid was mixed with the digested sample, and the pH was adjusted to < 3, and then the sample was centrifuged at 12,000 × g for 5 min. After the supernatant flowing through a C18 desalting column, the column was washed thrice (0.1 % formic acid, 4 % acetonitrile) and then eluted with elution buffer (0.1 % formic acid, 75 % acetonitrile). The eluents were collected and lyophilized. Next, 100 µL of 0.1 M TEAB buffer was added for reconstitution, followed by the addition of 41 µL of acetonitrile-dissolved TMT labeling reagent (ThermoFisher, Massachusetts, USA). After mixing the samples by shaking for 2 h, the reaction was stopped by the addition of 8 % ammonia. Then, labeling samples were desalted and lyophilized.
Separation of fractions
To develop a gradient elution, mobile phase A (2 % acetonitrile, adjusted pH to 10.0) and B (98 % acetonitrile) were used. The lyophilized powder was dissolved in solution A. The sample was fractionated using a C18 column, and the column oven was set at 50 °C. At 214 nm, the eluates were monitored, followed by collection into a tube every minute, and eventually combination into 10 fractions. Under vacuum, all fractions were dried and then reconstituted in 0.1 % (v/v) formic acid in water.
Liquid chromatography tandem mass spectrometry (MS) analysis
For transition library construction, shotgun proteomics analyses were conducted using an EASY-nLCTM 1200 UHPLC system (Thermo Fisher, Massachusetts, USA) coupled with a Q Exactive HF(X) mass spectrometer (Thermo Fisher, Massachusetts, USA). The Q Exactive HF(X) MS allowed operation in a data-dependent acquisition (DDA) mode.
Identification and quantification of differentially abundant proteins
The resulting spectra from each run were searched against the SwissProt human database using the Proteome Discoverer v2.2 (ThermoFisher, Massachusetts, USA) search engine. A maximum of two miscleavage sites were allowed. Peptide spectrum matches (PSMs) with a credibility of more than 99 % were identified as PSMs. The identified proteins contained at least one unique peptide. PSMs and proteins with a false discovery rate (FDR) of no more than 1.0 % were retained. Differentially expressed proteins (DEPs) were identified using the limma package  with the following criteria: (1) fold change value > 1.2- or < 1/1.2, (2) P value < 0.05. DEP analysis was carried out in R 4.1.0.
Gene Ontology (GO) annotation of identified proteins was performed using the UniProt-GOA database and InterProScan program version 5 . Pathway annotation was carried out using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Enrichment analysis of functional annotations was conducted using the R package clusterProfiler . A P value < 0.05 was considered significant, and enriched GO terms or pathways were sorted by P values. Principle component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) clustering analyses  were conducted using the princomp() function or Rt-SNE package, respectively, in R. Hierarchical clustering analysis of DEPs was performed using the R gplots package. Protein-protein interactions (PPI) were derived using the STRING database  (https://string-db.org/). Network visualization of PPI was carried out using Cytoscape version 3.8.0 , and hub scores for each identified DEP were obtained using the cytohubba plugin . The modules of hub genes were filtered using the MCODE plugin in Cytoscape with the following criteria: degree cut-off > 2, node score cut-off > 0.2, K-score > 2, and max depth = 10.
Determination of lipids in serum and follicular fluid
Venous blood or follicular fluid samples were centrifuged at 3,000 × g for 10 min to isolate the upper layer for analysis. Serum sex hormones, including follicle-stimulating hormone, luteinizing hormone, estradiol, progesterone and testosterone, were analyzed using commercially available kits from the Unicel DXI 800 Access immunoassay system (Beckman Coulter). Total cholesterol, triglycerides, high-density lipoprotein, and low-density lipoprotein levels were determined using commercially available kits purchased from Randox Laboratories Ltd., Northern Ireland, United Kingdom. The assays were carried out on a Hitachi 7600–210 Automatic Biochemical Analyzer (Hitachi Corporation, Japan). Commercially available kits for free fatty acid (FFA) determination were purchased from Nanjing Jiancheng Bioengineering Institute (Jiangsu, China).
All analyses were performed using R version 4.1.0. Normally distributed continuous data are presented as the mean ± SE and were compared using the Student’s t-test, while non-parametric data were presented as the median (interquartile range) and were compared using the Mann–Whitney U test, when appropriate. Correlations were evaluated using Spearman’s correlation coefficient. A two-sided P-value < 0.05 was considered statistically significant for all tests.
Subject clinical characteristics
Descriptive statistics and ovarian stimulation characteristics were recorded (Table 1). There was no statistical difference in age, basal hormone levels, or ovarian stimulation procedure data between the two groups. The BMI in the obesity group ranged from 29.25 to 32.33 kg/m2, much higher than that in the control group (20.38‒21.02 kg/m2).
Subject lipid parameters
Table 2 summarizes lipid levels in the serum and follicular fluid of all participants. FFA and triglycerides levels in follicular fluid and FFA levels in serum of obese women were significantly higher than those in matched normal-weight women. No differences were observed in total cholesterol, high-density lipoprotein, and low-density lipoprotein levels between the two groups.
To obtain a global view of the GC proteome in normal-weight and obese subjects, sample clustering based on all identified proteins was performed. As shown in Fig. 1 a, PCA analysis revealed two significant clusters, which indicated distinct protein expression levels in obese subjects. The above result was further confirmed using t-SNE clustering analysis, which also showed two isolated clusters (Fig. 1b). Subsequently, DEPs between the two groups were screened. A total 228 DEPs were identified, including 138 upregulated and 90 downregulated proteins (Fig. 1 c). As shown in Fig. 1 c, the top 10 DEPs were highlighted and labeled, including RNASEH2C, INSL3, TAPBPL, TSPAN8, TUBB8, and SUMO3. Next, hierarchical clustering based on DEPs was performed, and the results revealed clusters similar to those of PCA and t-SNE analysis (Fig. 1d). The annotation and quantification information for the top 10 DEPs are presented in Table 3.
Changes in the mitochondrial electron transport chain in obese subject GCs
Functional enrichment analyses, including GO and KEGG, were carried out to investigate the roles of DEPs in GCs and to explore the effect of obesity on oocyte quality and fertility outcome. DEPs were enriched in biological processes related to the mitochondrial respiratory chain, including electron transport, ATP synthesis, cellular respiration, energy derivation by oxidation of organic compounds, and ATP metabolic process (Fig. 2 a). In the cellular component analysis, DEPs were found to be enriched in respirasome, mitochondrial respirasome, endoplasmic reticulum lumen, and other energy metabolism-related subcellular locations (Fig. 2b). Similar results were observed in molecular function analysis, in which electron transfer activity, protease binding, and oxidoreductase activity, among other terms, were significantly enriched (Fig. 2 c).
KEGG pathway analysis revealed that the most significantly enriched pathways included energy metabolism (oxidative phosphorylation, thermogenesis, and carbon metabolism), extracellular matrix (ECM) (ECM-receptor interaction, focal adhesion, and regulation of actin cytoskeleton), steroid hormone synthesis (steroid biosynthesis and ovarian steroidogenesis), diabetes complications (AGE-RAGE signaling pathway in diabetic complications), apoptosis, the PI3K-AKT signaling pathway, and the relaxin (INSL3) signaling pathway, among other terms (Fig. 2d).
Key genes and modules based on the PPI network
We built the PPI network of DEPs using the STRING database. All DEPs were filtered into the DEP PPI network complex, which contained 189 nodes and 623 edges (Fig. 3 a). Among the 189 nodes, 11 central node genes were selected as hub genes with the criterium of filtering degree ≥ 20 (i.e., each node having more than 20 connections/interactions). The hub genes were FN1, MAPK3, CYCS, CYC1, NDUFAB1, SDHC, COL1A1, COX4I1, FGF2, COL3A1, and UQCR10. To identify the core gene modules that might play an important role in related pathways, core gene modules were also analyzed, and a total of three hub modules are shown in Fig. 3b-d. Functional enrichment analysis revealed that the genes in hub modules 1, 2, and 3 were mainly associated with “ATP synthesis,” “leukocyte migration”, and “extracellular matrix organization”, respectively (Fig. 3 a).
Obesity and infertility are inextricably linked, and their relationship has been increasingly recognized . As a key factor, GCs play an important role in oocyte and follicle development. They secrete nutrients and hormones essential for primary and secondary follicle growth, providing a suitable microenvironment for oocyte meiosis and maturation [30, 31]. Therefore, a precise understanding of the effect of obesity on GCs may help to determine the mechanism of obesity-related female infertility or subfertility.
The exposure of ovarian cells to high levels of fatty acids could result in an inflammatory response within ovarian follicles , excessive ovarian androgen production through ineffective FFA β-oxidation in mitochondria , insulin resistance , and oxidative DNA damage caused by reactive oxygen species (ROS) reacting with DNA within the cell, including mitochondrial DNA (mtDNA) . In the current study, the concentration of FFA and triglycerides in the follicular fluid of obese subjects was higher than that in control group. This was consistent with the results of our previous study, in which we found that FFA and triglycerides, as important biomedical indicators of abnormal lipid metabolism, were elevated in the serum and follicular fluid of obese subjects with PCOS and associated with lower embryo quality .
Mass spectrometry is widely used in proteomics for accurate and high-throughput analysis, allowing for the identification and quantification of proteins in complex samples with high sensitivity . In this study, we used TMT-based high-throughput MS to investigate the proteome of GC samples from normal-weight and obese subjects. DEPs between the two groups were screened, and a total of 228 DEPs were identified, including 138 upregulated and 90 downregulated proteins. Significant pathways and GO terms associated with protein expression changes were also identified, especially those associated with the mitochondrial electron transport chain. With regard to the biological process category, almost all of the significant GO terms were associated with the mitochondrial respiratory chain.
Caloric excess and high-fat diet is the most common cause of obesity. High-fat diet can lead to nutritional excess, increased electron flux, rising oxidative stress, an accumulation of oxidized substrates, and eventually, mitochondrial injury . Physiologically, the mitochondrial matrix flowing within the inner membrane, which hosts several copies mtDNA as well as components of the tricarboxylic acid pathway and electron transport chain, is the site of ATP production . Excess FFA can increase mitochondrial activity, contributing to higher oxidative stress, which adversely affects the replication of mtDNA and leads to the dysfunction or dysregulation of key mtDNA gene products. Subsequently, mtDNA gene products control oxidative phosphorylation, thereby impairing the maturation and developmental potential of oocytes . Damage to mitochondrial function caused by high concentrations of FFA has been confirmed by both in vivo and in vitro studies in multiple tissue and cell types [41,42,43]. In the current study, a significant correlation between high FFA in follicular fluid and mitochondrial dysfunction in GCs from obese subjects was observed. Among molecular function GO terms, electron transfer activity, which is essential for mitochondrial energy production, was listed as the most significantly perturbed function. In addition, cytochrome-c-oxidase activity, which is also key factor for the respiratory electron transport chain, was altered in GCs from obese subjects.
The ECM is a diverse mixture of molecules, influencing almost all aspects of the follicle development process, including migration, division, differentiation, cell death, and cell anchorage . Various ECM types in the ovarian follicles have been characterized, including those associated with the follicular basal lamina, the stromal matrix of the thecal layers, the matrix of thecal blood vessels, and that of the cumulus oocyte complexes [45,46,47]. Changes relevant to the structure of the ECM were related to a large part of the obtained molecular function GO terms, including changes in extracellular matrix structural constituent, extracellular matrix structural constituent conferring tensile strength, protease binding, collagen binding, and scaffold protein binding. The current study findings indicated that obesity and high FFA levels may affect the function of ovarian GCs via alteration of the ECM.
The regulation of endoplasmic reticulum homeostasis is indispensable during folliculogenesis and oocyte maturation . Some physiological and pathological conditions, such as glucose deprivation, inflammation, oxidative stress, elevated FFA levels, aberrant Ca2+ regulation, and hypoxia disrupt endoplasmic reticulum function, which results in the accumulation of unfolded or misfolded proteins within the endoplasmic reticulum, leading to the induction of endoplasmic reticulum stress . Fatty acid-induced endoplasmic reticulum stress in mouse cumulus-oocyte complexes impairs protein secretion and mitochondrial activity, resulting in abnormal embryo development . Another mouse study also demonstrated a relationship between obesity and endoplasmic reticulum stress in cumulus-oocyte complexes, which was associated with reduced mitochondrial membrane potential, high autophagy levels, and high intracellular lipid levels . The mechanisms underlying endoplasmic reticulum changes in the GCs of obese female subjects in this study may be similar to those reported in animal studies [50, 51].
Regarding the cellular component category, the major enriched terms were related to the mitochondria, ECM, and endoplasmic reticulum, which were also consistent with the results for the biological process and molecular function categories. The mitochondria-related terms included mitochondrial respirasome, respiratory chain complex, cytochrome complex, and oxidoreductase complex. ECM terms included collagen trimer, fibrillar collagen trimer, banded collagen fibril, and complex of collagen trimers. Endoplasmic reticulum terms included endoplasmic reticulum lumen.
Several studies have shown that oxidative stress caused by increased levels of ROS can induce GC apoptosis, resulting in decreased estradiol 17β levels, lower ovulation rates, and compromised oocyte quality . In the KEGG pathway analysis, energy metabolism and ECM-related pathways were among the most significantly enriched, which were consistent with GO results. Thus, oxidative stress caused by mitochondrial dysfunction and a high fatty acid environment may affect steroid synthesis, in turn affecting the quality of oocytes.
In patients with PCOS, advanced glycation end-product (AGE) levels were significantly increased compared with those in matched normal-weight controls [53, 54]. AGEs can alter insulin receptor substrate and AKT phosphorylation, thereby stimulating the downstream PI3K signaling cascade and compromising the effects of insulin signaling [55, 56]. Excessive deposition of AGEs can impair the ECM structure, altering its biochemical characteristics and associated metabolism . The binding of AGE to its receptor RAGE produces NAPDH, which enhances oxidative stress, activates the NF-κB signaling pathway and further stimulates the production of cytokines and other factors. It finally contributes to cellular damage and even apoptosis [58, 59]. Moreover, the PI3K-AKT signaling pathway was associated with human GC damage and apoptosis [60, 61]. An interesting result of the KEGG pathway analysis was the significant enrichment of the AGE-RAGE signaling pathway in the current study. This is the first to demonstrate the association between AGEs and impaired GC function in obese non-PCOS subjects. It suggests that obesity-related infertility or subfertility is closely associated with the metabolic, inflammatory, structural, and functional changes in GCs.
TSPAN8 is a cell surface glycoprotein, which forms complexes with integrins. Previous genome-wide association studies identified TSPAN8 loci, which was correlated with an increased risk of type 2 diabetes (T2D) in humans [62,63,64]. Close examination of the haploblock structure with the inclusion of additional SNPs suggested that TSPAN8 is likely to be a T2D causal gene . Genetic ablation of TSPAN8 resulted in a reduction of body weight in mice fed a normal diet as well as in a resistance to body weight gain upon high-fat and high-carbohydrate diet feeding . Small ubiquitin-like modifier 3 (SUMO3) is a member of the SUMO family of eukaryotic proteins. Subcutaneous adipose tissue from cows incubated with insulin exhibited an upregulation of SUMO3, in parallel to an increased expression of other key genes associated with insulin signaling, adipogenesis, and lipogenesis . In the field of reproduction, previous literature reported that excessive SUMOylation (SUMO1/2/3) was a marker of defective spermatozoa . INSL3 is a member of the insulin-like hormone superfamily and is mainly produced in gonadal tissues. Within the follicle, INSL3 acts via its G protein-coupled receptor, RXFP2, in an autocrine/paracrine manner to drive the production of the major steroid precursor androstenedione as well as its conversion into estrogen by GCs . MED25 was previously shown to play a critical role in the endoplasmic reticulum stress response  Here, a number of the DEPs identified in the current work, including TSPAN8, SUMO3, INSL3, and MED25, were implicated in obesity-related infertility or subfertility. In addition, INSL3 was one of the most differentially expressed proteins, which was consistent with the enrichment of the relaxin signaling pathway in KEGG analysis. Other major DEPs, including RNASEH2C, TAPBPL, TUBB8, NFRKB, MASTL, and MYLK, have never been reported to be associated with infertility, metabolic disorders, or hormone imbalances. Thus, their role in the occurrence and development of female infertility requires further attention.
In summary, TMT-based proteomic and bioinformatic analyses indicated that mitochondrial damage, the endoplasmic reticulum stress response as well as dysregulated hormonal synthesis were exhibited in GCs obtained from obese subjects whereas none of these changes occurred in normal-weight subjects. These alterations may be related to the high FFA and TG levels detected in human follicular fluid. Therefore, the current research provides valuable insight for the further etiological investigation of obesity-related infertility or subfertility and the identification of therapeutic targets.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Body mass index
Free fatty acid
Principle component analysis
t-distributed stochastic neighbor embedding
Peptide spectrum matches
In vitro fertilization
Differentially expressed protein
Kyoto Encyclopedia of Genes and Genomes
Polycystic ovary syndrome
Reactive oxygen species
Small ubiquitin-related modifier
Advanced glycation end-product
Tandem mass tag
Best D, Avenell A, Bhattacharya S. How effective are weight-loss interventions for improving fertility in women and men who are overweight or obese? A systematic review and meta-analysis of the evidence. Hum Reprod Update. 2017;23:681–705.
Pandey S, Bhattacharya S. Impact of obesity on gynecology. Womens Health (Lond). 2010;6:107–17.
Kumbak B, Oral E, Bukulmez O. Female obesity and assisted reproductive technologies. Semin Reprod Med. 2012;30:507–16.
Penzias AS. Recurrent IVF failure: other factors. Fertil Steril. 2012;97:1033–8.
Broughton DE, Moley KH. Obesity and female infertility: potential mediators of obesity’s impact. Fertil Steril. 2017;107:840–47.
Shah DK, Missmer SA, Berry KF, Racowsky C, Ginsburg ES. Effect of obesity on oocyte and embryo quality in women undergoing in vitro fertilization. Obstet Gynecol. 2011;118:63–70.
Provost MP, Acharya KS, Acharya CR, Yeh JS, Steward RG, Eaton JL, Goldfarb JM, Muasher SJ. Pregnancy outcomes decline with increasing body mass index: analysis of 239,127 fresh autologous in vitro fertilization cycles from the 2008–2010 Society for Assisted Reproductive Technology registry. Fertil Steril. 2016;105:663–69.
Wang JX, Davies M, Norman RJ. Body mass and probability of pregnancy during assisted reproduction treatment: retrospective study. BMJ. 2000;321:1320–1.
Sermondade N, Huberlant S, Bourhis-Lefebvre V, Arbo E, Gallot V, Colombani M, Fréour T. Female obesity is negatively associated with live birth rate following IVF: a systematic review and meta-analysis. Hum Reprod Update. 2019;25:439–51.
Lashen H, Ledger W, Bernal AL, Barlow D. Extremes of body mass do not adversely affect the outcome of superovulation and in-vitro fertilization. Human reproduction (Oxford England). 1999;14:712–15.
Dechaud H, Anahory T, Reyftmann L, Loup V, Hamamah S, Hedon B. Obesity does not adversely affect results in patients who are undergoing in vitro fertilization and embryo transfer. Eur J Obstet Gynecol Reprod Biol. 2006;127:88–93.
Martinuzzi K, Ryan S, Luna M, Copperman AB. Elevated body mass index (BMI) does not adversely affect in vitro fertilization outcome in young women. J Assist Reprod Genet. 2008;25:169–75.
Valckx SDM, De Pauw I, De Neubourg D, Inion I, Berth M, Fransen E, Bols PEJ, Leroy JLMR. BMI-related metabolic composition of the follicular fluid of women undergoing assisted reproductive treatment and the consequences for oocyte and embryo quality. Human reproduction (Oxford, England). 2012; 27:3531-39.
Uyar A, Torrealday S, Seli E. Cumulus and granulosa cell markers of oocyte and embryo quality. Fertil Steril. 2013;99:979–97.
Chronowska E. High-throughput analysis of ovarian granulosa cell transcriptome. Biomed Res Int. 2014; 2014:213570.
Russell DL, Gilchrist RB, Brown HM, Thompson JG. Bidirectional communication between cumulus cells and the oocyte: Old hands and new players? Theriogenology. 2016;86:62–8.
Coticchio G, Dal Canto M, Mignini Renzini M, Guglielmo MC, Brambillasca F, Turchi D, Novara PV, Fadini R. Oocyte maturation: gamete-somatic cells interactions, meiotic resumption, cytoskeletal dynamics and cytoplasmic reorganization. Hum Reprod Update. 2015;21:427–54.
Sun L, Hu W, Liu Q, Hao Q, Sun B, Zhang Q, Mao S, Qiao J, Yan X. Metabonomics reveals plasma metabolic changes and inflammatory marker in polycystic ovary syndrome patients. J Proteome Res. 2012;11:2937–46.
Zhang X, Xu X, Li P, Zhou F, Kong L, Qiu J, Yuan Z, Tan J. TMT Based Proteomic Analysis of Human Follicular Fluid From Overweight/Obese and Normal-Weight Patients With Polycystic Ovary Syndrome. Front Endocrinol. 2019;10:821.
Escobar-Morreale HF. Polycystic ovary syndrome: definition, aetiology, diagnosis and treatment. Nature reviews Endocrinology. 2018;14:270–84.
Consultation WHOE. Appropriate body-mass index for Asian populations and its implications for policy and intervention strategies. Lancet. 2004;363:157–63.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
van der Maaten L, Hinton G. Visualizing Data using t-SNE. Journal of Machine Learning Research. 2008;9:2579–605.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607-D13.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):11.
Talmor A, Dunphy B. Female obesity and infertility. Best Pract Res Clin Obstet Gynaecol. 2015;29:498–506.
Komatsu K, Masubuchi S. Mouse oocytes connect with granulosa cells by fusing with cell membranes and form a large complex during follicle development. Biol Reprod. 2018;99:527–35.
Ernst EH, Lykke-Hartmann K. Transcripts encoding free radical scavengers in human granulosa cells from primordial and primary ovarian follicles. J Assist Reprod Genet. 2018;35:1787–98.
Gonzalez MB, Lane M, Knight EJ, Robker RL. Inflammatory markers in human follicular fluid correlate with lipid levels and Body Mass Index. J Reprod Immunol. 2018;130:25–9.
Gervais A, Battista MC, Carranza-Mamane B, Lavoie HB, Baillargeon JP. Follicular fluid concentrations of lipids and their metabolites are associated with intraovarian gonadotropin-stimulated androgen production in women undergoing in vitro fertilization. J Clin Endocrinol Metab. 2015;100:1845–54.
Xu L, Wang W, Zhang X, Ke H, Qin Y, You L, Li W, Lu G, Chan WY, Leung PCK, et al. Palmitic acid causes insulin resistance in granulosa cells via activation of JNK. J Mol Endocrinol. 2019;62:197–206.
Bradley J, Swann K. Mitochondria and lipid metabolism in mammalian oocytes and early embryos. Int J Dev Biol. 2019;63:93–103.
Niu Z, Lin N, Gu R, Sun Y, Feng Y. Associations between insulin resistance, free fatty acids, and oocyte quality in polycystic ovary syndrome during in vitro fertilization. J Clin Endocrinol Metab. 2014;99:E2269-76.
Li S, Tang H. Computational Methods in Mass Spectrometry-Based Proteomics. Adv Exp Med Biol. 2016;939:63–89.
Bhatraju NK, Agrawal A. Mitochondrial Dysfunction Linking Obesity and Asthma. Ann Am Thorac Soc. 2017;14:368-S73.
Kimura I, Ichimura A, Ohue-Kitano R, Igarashi M. Free Fatty Acid Receptors in Health and Disease. Physiol Rev. 2020;100:171–210.
Grindler NM, Moley KH. Maternal obesity, infertility and mitochondrial dysfunction: potential mechanisms emerging from mouse model systems. Mol Hum Reprod. 2013;19:486–94.
Sacca SC, Cutolo CA, Ferrari D, Corazza P, Traverso CE. The Eye, Oxidative Damage and Polyunsaturated Fatty Acids. Nutrients. 2018; 10.
Harasymowicz NS, Dicks A, Wu CL, Guilak F. Physiologic and pathologic effects of dietary free fatty acids on cells of the joint. Ann N Y Acad Sci. 2019;1440:36–53.
Reis AH. Acidemia and blood free fatty acids: analysis of cardiovascular risk factors in a new context. Discov Med. 2017;23:183–88.
Irving-Rodgers HF, Rodgers RJ. Extracellular matrix in ovarian follicular development and disease. Cell Tissue Res. 2005;322:89–98.
Rodgers RJ, Irving Rodgers HF. Extracellular matrix of the bovine ovarian membrana granulosa. Mol Cell Endocrinol. 2002;191:57–64.
Rodgers RJ, Irving-Rodgers HF, Russell DL. Extracellular matrix of the developing ovarian follicle. Reproduction. 2003;126:415–24.
Rodgers RJ, Irving-Rodgers HF, van Wezel IL. Extracellular matrix in ovarian follicles. Mol Cell Endocrinol. 2000;163:73–9.
Guzel E, Arlier S, Guzeloglu-Kayisli O, Tabak MS, Ekiz T, Semerci N, Larsen K, Schatz F, Lockwood CJ, Kayisli UA. Endoplasmic Reticulum Stress and Homeostasis in Reproductive Physiology and Pathology. Int J Mol Sci. 2017; 18.
Huang N, Yu Y, Qiao J. Dual role for the unfolded protein response in the ovary: adaption and apoptosis. Protein Cell. 2017;8:14–24.
Wu LL, Russell DL, Norman RJ, Robker RL. Endoplasmic reticulum (ER) stress in cumulus-oocyte complexes impairs pentraxin-3 secretion, mitochondrial membrane potential (DeltaPsi m), and embryo development. Mol Endocrinol. 2012;26:562–73.
Wu LL, Russell DL, Wong SL, Chen M, Tsai TS, St John JC, Norman RJ, Febbraio MA, Carroll J, Robker RL. Mitochondrial dysfunction in oocytes of obese mothers: transmission to offspring and reversal by pharmacological endoplasmic reticulum stress inhibitors. Development. 2015;142:681–91.
Chaube SK, Shrivastav TG, Tiwari M, Prasad S, Tripathi A, Pandey AK. Neem (Azadirachta indica L.) leaf extract deteriorates oocyte quality by inducing ROS-mediated apoptosis in mammals. Springerplus. 2014;3:464.
Azhary JMK, Harada M, Kunitomi C, Kusamoto A, Takahashi N, Nose E, Oi N, Wada-Hiraike O, Urata Y, Hirata T, et al. Androgens Increase Accumulation of Advanced Glycation End Products in Granulosa Cells by Activating ER Stress in PCOS. Endocrinology. 2020; 161.
Diamanti-Kandarakis E, Piperi C, Patsouris E, Korkolopoulou P, Panidis D, Pawelczyk L, Papavassiliou AG, Duleba AJ. Immunohistochemical localization of advanced glycation end-products (AGEs) and their receptor (RAGE) in polycystic and normal ovaries. Histochem Cell Biol. 2007;127:581–9.
Wu CH, Huang HW, Huang SM, Lin JA, Yeh CT, Yen GC. AGE-induced interference of glucose uptake and transport as a possible cause of insulin resistance in adipocytes. J Agric Food Chem. 2011;59:7978–84.
Diamanti-Kandarakis E, Chatzigeorgiou A, Papageorgiou E, Koundouras D, Koutsilieris M. Advanced glycation end-products and insulin signaling in granulosa cells. Exp Biol Med (Maywood). 2016;241:1438–45.
Bansode SB, Gacche RN. Glycation-induced modification of tissue-specific ECM proteins: A pathophysiological mechanism in degenerative diseases. Biochim Biophys Acta Gen Subj. 2019;1863:129411.
Wautier MP, Chappey O, Corda S, Stern DM, Schmidt AM, Wautier JL. Activation of NADPH oxidase by AGE links oxidant stress to altered gene expression via RAGE. Am J Physiol Endocrinol Metab. 2001;280:E685-94.
Faria A, Persaud SJ. Cardiac oxidative stress in diabetes: Mechanisms and therapeutic potential. Pharmacol Ther. 2017;172:50–62.
Zhao Y, Zhang C, Huang Y, Yu Y, Li R, Li M, Liu N, Liu P, Qiao J. Up-regulated expression of WNT5a increases inflammation and oxidative stress via PI3K/AKT/NF-kappaB signaling in the granulosa cells of PCOS patients. J Clin Endocrinol Metab. 2015;100:201–11.
Nakahara T, Iwase A, Nakamura T, Kondo M, Bayasula, Kobayashi H, Takikawa S, Manabe S, Goto M, Kotani T, Kikkawa F. Sphingosine-1-phosphate inhibits H2O2-induced granulosa cell apoptosis via the PI3K/Akt signaling pathway. Fertil Steril. 2012;98:1001-8 e1.
Stuebe AM, Lyon H, Herring AH, Ghosh J, Wise A, North KE, Siega-Riz AM. Obesity and diabetes genetic variants associated with gestational weight gain. Am J Obstet Gynecol. 2010;203:283 e1–17.
Kong X, Zhang X, Xing X, Zhang B, Hong J, Yang W. The Association of Type 2 Diabetes Loci Identified in Genome-Wide Association Studies with Metabolic Syndrome and Its Components in a Chinese Population with Type 2 Diabetes. PLoS One. 2015;10:e0143607.
Zhou DZ, Liu Y, Zhang D, Liu SM, Yu L, Yang YF, Zhao T, Chen Z, Kan MY, Zhang ZF, et al. Variations in/nearby genes coding for JAZF1, TSPAN8/LGR5 and HHEX-IDE and risk of type 2 diabetes in Han Chinese. J Hum Genet. 2010;55:810–5.
Champy MF, Le Voci L, Selloum M, Peterson LB, Cumiskey AM, Blom D. Reduced body weight in male Tspan8-deficient mice. Int J Obes (Lond). 2011;35:605–17.
Lopreiato V, Hosseini A, Rosa F, Zhou Z, Alharthi A, Trevisi E, Loor JJ. Dietary energy level affects adipose depot mass but does not impair in vitro subcutaneous adipose tissue response to short-term insulin and tumor necrosis factor-alpha challenge in nonlactating, nonpregnant Holstein cows. J Dairy Sci. 2018;101:10206–19.
Vigodner M, Shrivastava V, Gutstein LE, Schneider J, Nieves E, Goldstein M, Feliciano M, Callaway M. Localization and identification of sumoylated proteins in human sperm: excessive sumoylation is a marker of defective spermatozoa. Hum Reprod. 2013;28:210–23.
Ivell R, Anand-Ivell R. Insulin-like peptide 3 (INSL3) is a major regulator of female reproductive physiology. Hum Reprod Update. 2018;24:639–51.
Sela D, Conkright JJ, Chen L, Gilmore J, Washburn MP, Florens L, Conaway RC, Conaway JW. Role for human mediator subunit MED25 in recruitment of mediator to promoters by endoplasmic reticulum stress-responsive transcription factor ATF6alpha. J Biol Chem. 2013;288:26179–87.
The authors thank Prof. Peter Reinach for his editorial assistance.
This research project was financially supported by grants from the National Natural Science Foundation of China (No. 82071694, 81671530) and Science and Technology Commission of Shanghai Municipality (No. 21140904000, 201409005800).
Ethics approval and consent to participate
All participants signed informed consent forms approved by the Shanghai Jiao Tong University Committee on the Use of Human Subjects in Medical Research. Ethical approval of follicular fluid collection was obtained from the institutional Ethics Committee of Ruijin Hospital, School of Medicine, Shanghai Jiao Tong University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Si, C., Wang, N., Wang, M. et al. TMT-based proteomic and bioinformatic analyses of human granulosa cells from obese and normal-weight female subjects. Reprod Biol Endocrinol 19, 75 (2021). https://doi.org/10.1186/s12958-021-00760-x
- Granulosa cells
- Proteomic analysis
- Free fatty acids
- Electron transport chain