- 1Centre for Reproductive Medicine, Shanghai 10th People Hospital of Tongji University, Shanghai, China
- 2Department of Basic Medical Science, Faculty of Medicine Vajira Hospital, Navamindradhiraj University, Bangkok, Thailand
- 3Shanghai Clinical College, Anhui Medical University, Hefei, China
We evaluated metabolic profiles between cumulus cells (CCs) and mural granulosa cells (MGCs) derived from women with endometriosis to identify their correlations with oocyte quality. CCs and MGCs were collected from women with and without endometriosis undergoing in vitro fertilization/intracytoplasmic sperm injection treatment. The metabolomics of CCs and MGCs were measured by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) followed by a quantitative polymerase chain reaction to further confirm the genes involved in the metabolic results. LC-MS/MS analysis revealed differences in 24 metabolites of CCs and 71 metabolites of MGCs between groups. Among them, five metabolites were upregulated and 19 metabolites were downregulated in CCs with endometriosis, whereas three metabolites were upregulated and 68 metabolites were downregulated in MGCs with endometriosis. Metabolites related to sphingolipid metabolism, which included palmitic acid (PA) and docosahexaenoic acid, increased significantly only in CCs with endometriosis, whereas sphingosine and PA were significantly downregulated in MGCs with endometriosis compared with CCs and MGCs without endometriosis. Gene expression involved in ceramide synthesis (CERS1, SPTL1, and SMPD1) and autophagy (BECN1, LAMP, and PC3) were significantly higher in CCs with endometriosis according to FASN, BECN1, and LAMP protein expressions. However, gene expression involved in ceramide synthesis (SPHK1, ASAH1, and SGPP1) and autophagy (BECN1, LAMP, and PC3) were significantly lower in MGCs with endometriosis, whereas CERS1 and UGCG expression increased. There are differences in sphingolipid metabolites in CCs and MGCs with endometriosis compared with women without endometriosis. These differences seem to be involved in the regulation of autophagic cell death in preovulatory follicles.
1 Introduction
Endometriosis is an estrogen-dependent disorder affecting women of reproductive age and is characterized by the presence of endometrial tissue (glands and stroma) outside the uterine cavity and involved with chronic pelvic pain and infertility (1). Endometriosis affects fertility and several processes in the female reproductive system including folliculogenesis, ovulation, oocyte quality, and implantation (2). However, the mechanism of infertility is still unclear and insufficient to summarize the relation between endometriosis and its poor quality of oocytes (3).
Interestingly, recent studies demonstrated that endometriosis affects granulosa cells (GCs) in many aspects including an increased level of intracellular reactive oxygen species (ROS) generation, apoptosis, and dysregulation of the pathway involved in cell growth and development (4). During folliculogenesis, GCs play a critical role in the maturation of the oocyte. GCs differentiated into cumulus cells (CCs) and mural granulosa cells (MGCs) during the development of follicles. CCs are the cells in immediate contact with the oocyte, and MGCs line the follicular wall around the antrum. These cells play a critical role in follicular development and the maturation of oocytes (5).
Metabolomics is the new powerful tool for approaching disease progression via the new biomarker identification (6). Many reports demonstrated that metabolomic analysis in women with endometriosis in endometrial fluid (7), follicular fluid (8), urine (9), serum (10), and plasma (11) have represented the alteration of metabolite profiles. Palmitic acid (PA; 16:0) is the most common saturated fatty acid found in the human body (12) and can be provided the role of a precursor in sphingolipid metabolism-regulated several biological processes, including cell survival, migration, and apoptosis (13). Previous studies demonstrated that alteration of follicular free fatty acid (FFA) levels is involved with oocyte quality in both animals (14) and humans (15).
Recent studies also showed the excess palmitic acid-induced apoptosis of GCs in both animals (16) and humans (17), leading to impairment of ovarian follicular development and oocyte maturation. Moreover, it has been shown that sphingolipid metabolism changed in serum, peritoneal fluid, and endometrial tissue of women with endometriosis (13). In contrast, a report demonstrated that there is no specific metabolomic profile in the follicular fluid, which revealed no impairment of the cumulus–oocyte complex (COC) microenvironment of women with endometriosis (18).
In the present study, we aim to investigate the metabolic profiles in CCs and MGCs of women with and without endometriosis in order to find the correlation with the quality of the oocyte. We integrated the results of palmitic acid and sphingolipids’ balance to cell death at mRNA and protein levels via an autophagic pathway to evaluate the impact of endometriosis in the follicle related to the quality of the oocyte.
2 Materials and methods
2.1 Collection and isolation of cumulus cells and mural granulosa cells
After ovum pick up (OPU) by transvaginal ultrasound-guided needle aspiration, MGCs and CCs were collected from the follicular fluid under a stereomicroscope. The density gradient technique was used to purify MGCs and CCs by 40% and 80% gradient (SAGE), respectively. Then they were centrifuged for 15 min at 300 g, and the middle layer was collected. The remaining red blood cells were lysed by blood cell lysing buffer (Invitrogen by Thermo Fisher Scientific, Waltham, MA, USA). Then samples were washed by centrifugation 3 times with 1× phosphate-buffered saline (PBS). After that, cells were cryopreserved in liquid nitrogen until further metabolic analysis (19, 20).
2.2 Liquid chromatography coupled to tandem mass spectrometry
2.2.1 Chromatography–mass spectrometry analysis
2.2.1.1 Chromatographic conditions
CCs and MGCs were extracted by methanol/acetonitrile/aqueous solution (2:2:1, v/v) (21) and separated using Agilent 1290 Infinity LC ultra-high-performance liquid chromatography (UHPLC) hydrophilic interaction liquid chromatography (HILIC) column: column temperature was 25°C; flow rate was 0.3 ml/min. Mobile phase composition A = water +25 mM ammonium acetate +25 mM ammonia and B = acetonitrile. The gradient elution procedure was 0–0.5 min, 95% B; 0.5–7 min, linear change of B from 95% to 65%; 7–8 min, linear change of B from 65% to 40%; 8–9 min, B maintained at 40%; 9–9.1 min, linear change of B from 40% to 95%; and 9.1–12 min, B maintained at 95%. Samples during the entire analysis were placed in a 4°C autosampler (21, 22).
2.2.1.2 Quadrupole time-of-flight mass spectrometry conditions
After the sample was tested, the first- and second-level spectra were collected using the AB Triple TOF 6600 mass spectrometer. The electrospray ionization (ESI) source conditions were carried out according to the instructions after HILIC separation. Ion Source Gas1 (Gas1) was set to 60, Gas2 was also 60, and Curtain gas (23) was 30. IonSapary Voltage Floating (ISVF) was ±5,500 V (positive and negative modes); source temperature was 600°C; the m/z range of time-of-flight (TOF) MS scan and product ion scan was 60–1,000 and 25–1,000 Da, respectively; the accumulation time of TOF MS scan and product ion scan were 0.20 s/spectra and 0.05 s/spectra, respectively; the secondary mass spectrum was obtained by information-dependent acquisition (IDA) and adopts high sensitivity mode; declustering potential (DP) was set to ±60 V (positive and negative modes); collision energy was 35 ± 15 eV. IDA setting excludes isotopes within 4 Da, with candidate ions to monitor per cycle of 6 (24–26).
2.3 Quantitative PCR
Total RNA was extracted by TIANGEN RNA simple Total RNA Kit (TIANGEN), and RNA concentration was measured by NanoDrop ND 2000 spectrophotometer (Thermo Scientific). qPCR was carried according to the manufacturer’s instructions: cDNA (1 μl), forward and reverse primers (0.2 μl, 10 μM), RNase-free ddH2O (3.6 μl), and ChamQ Universal SYBR qPCR Master Mix (5 μl). Three-step PCR amplification protocol was followed: 95°C 30 s; 95°C 10 s → 60°C 30 s, a total of 40 cycles; 95°C 15 s → 60°C 60 s → 95°C 15 s. The housekeeping gene GAPDH was used as an internal reference, and the specific primer sequences are shown in Table 1 (27, 28).
2.4 Western blotting examination
Proteins were extracted by a nuclear protein and cytoplasmic protein extraction kit (Servicebio, Ghent, Belgium). Bicinchoninic acid (BCA) protein concentration assay kit (Servicebio) was used to measure protein concentration. Western blotting analyses were performed as usual procedures. Rabbit anti-human FASN polyclonal antibody (FASN) (abmart, Shanghai, China; 1:1,000), rabbit anti-human Beclin1 polyclonal antibody (Beclin1) (abmart, 1:1,000), and rabbit anti-human LAMP polyclonal antibody (Beclin1) (abmart, 1:1,000) were used as the primary antibodies for CCs. Rabbit anti-human ASAH polyclonal antibody (29) (abcam, Cambridge, UK; 1:1,000), rabbit anti-human Beclin1 polyclonal antibody (Beclin1) (abmart, 1:1,000), and rabbit anti-human LAMP polyclonal antibody (Beclin1) (abmart, 1:1,000) were used as the primary antibodies for MGCs. Goat anti-mouse IgG conjugated with horseradish peroxidase (HRP) was used as a secondary antibody (abclonal, 1:3,000). The density of the target bands was analyzed by the Alpha software processing system (30).
2.5 Data analysis
The original data of LC-MS/MS were converted into.mzXML format by ProteoWizard, and then the XCMS program was used for peak alignment, retention time correction, and peak area extraction. Data were preprocessed by Pareto scaling, and multi-dimensional and one-dimensional statistical analyses were performed. Statistical analysis includes Student’s t-test and multiple variation analysis, and volcano maps were drawn by R software (23, 31). The data of qPCR were analyzed using GraphPad PRISM 9 software (GraphPad Software Version 9, La Jolla, CA, USA) and free software for data manipulation, calculation, and graphical display, and the mRNA expression level of the target gene was expressed as the mean ± standard deviation. The independent-samples t-test was used to assess statistical significance, and p < 0.05 was considered statistically significant.
3 Result
3.1 Demographic data
Demographic data from women with and without endometriosis are summarized in Table 2. There were no significant differences in age, body mass index (BMI), and the other baseline clinical and endocrine profiles in women in the endometriosis group compared with the control group except that the level of estradiol was lower in the endometriosis group (609.3 pmol/L) when compared with the control group (172 pmol/L). The fertilization rate was lower in the endometriosis group (98.28%) when compared with the normal group (88.81%).
3.2 Metabolomic profiling of cumulus cells and mural granulosa cells in endometriosis and control women
Partial least squares discriminant analysis (PLS-DA) score plots showed that supervised multivariate analysis revealed more obvious differences between the metabolites between endometriosis and control, indicating a significant difference in metabolic components in positive ion mode and negative ion mode both in CCs (Figure 1A) and in MGCs (Figure 1C). A volcano plot showed the difference in metabolites between the two kinds of cells, calculated from data for all tested substances. Cluster analysis of all the samples further revealed that the repeat results for endometriosis and controls clustered together, indicating a significant difference in metabolism both in CCs (Figure 1B) and in MGCs (Figure 1D).
Figure 1 Metabolomic profiling of CCs and MGCs in endometriosis and control women by PLS-DA score plots and volcano plots. Score plots of principal component analysis of metabolome in the positive ion modes and the negative ion modes from CCs (A) and MGCs (C). These plots display a clear separation between normal and endometriosis metabolomes. The circle (blue and green) around each sample group represents the 95% confidence intervals. Volcano plot of CCs (B) and MGCs (D) showing the statistical significance (y-axis) and fold change (x-axis) for the difference between the metabolome of normal and endometriosis. p < 0.05, fold change ≥ 1.5 in positive ion mode and negative ion mode. CCs, cumulus cells; MGCs, mural granulosa cells; PLS-DA, partial least squares discriminant analysis.
3.3 Metabolomic analysis of cumulus cells in endometriosis and control women
3.3.1 Heat map and metabolite–metabolite correlation analysis
The results of the metabolomic analysis showed that there were differences in the metabolic profile of CCs and MGCs in endometriosis patients compared to controls. There were a total of 20,556 metabolite ion peaks: 11,784 in the positive ion mode and 8,772 in the negative ion mode. In CCs, 24 metabolites showed a difference in metabolite level between endometriosis and control, 14 different metabolites in positive ion mode, and 10 different metabolites in negative ion mode (Figure 2A). In MGCs, 71 metabolites showed a difference in metabolite level between endometriosis and control, 49 different metabolites in positive ion mode, and 22 different metabolites in negative ion mode (Figure 3A). The metabolites of variable importance in projection (VIP) >1 with multi-dimensional statistical analysis and p-value <0.05 with univariate statistical analysis as the significant differential metabolites were selected.
Figure 2 Metabolomic analysis of CCs in endometriosis and control women. Heat map of the changes in metabolites related to cumulus cells in endometriosis and control women. A heat map from cumulus cells in the positive ion modes and negative ion modes (A). These plots display a clear separation between normal and endometriosis metabolomes. Data of heap map plots reveal unit-variance scaling. The blue color represents the trend of reduction, and red color represents an increasing trend. Metabolite–metabolite correlation analysis of cumulus cells in endometriosis and control women in the positive ion modes and negative ion modes (B). Positive correlations are shown in blue, and negative correlations are shown in red. KEGG-enriched pathways changed between endometriosis patients and control (C). CCs, cumulus cells; KEGG, Kyoto Encyclopedia of Genes and Genomes.
Figure 3 Metabolomic analysis of MGCs in endometriosis and control women. Heat map of the changes in metabolites related to MGCs in endometriosis and control women. A heat map from MGCs in the positive ion modes and negative ion modes (A). These plots display a clear separation between normal and endometriosis metabolomes. Data of heap map plots reveal unit-variance scaling. The blue color represents the trend of reduction, and red color represents an increasing trend. Metabolite–metabolite correlation analysis of MGCs in endometriosis and control women in the positive ion modes and negative ion modes (B). Positive correlations are shown in blue, and negative correlations are shown in red. KEGG-enriched pathways changed between endometriosis patients and control (C). MGCs, mural granulosa cells; KEGG, Kyoto Encyclopedia of Genes and Genomes.
For preciseness, the role of cumulus cells on induced metabolites in endometriosis and multiple metabolites was screened. These metabolites were randomly divided into four different types including carbohydrate, amino acid, lipid, and other metabolites. There were differences in the metabolic profiles of CCs and MGCs between endometriosis and control. There were 3 metabolites upregulated in MGCs of endometriosis including, 1 lipid, 1 amino acid, and 1 other, whereas 68 metabolites were downregulated including 27 lipids, 8 amino acids, 4 carbohydrates, 6 energy, 9 nucleotides, and 9 nucleic acid pathways. All the detected metabolites were analyzed by heat mapping. In CCs, there were 5 metabolites upregulated in endometriosis, including 1 carbohydrate, 2 lipids, and 2 others, whereas 19 metabolites were downregulated including 7 lipids, 4 amino acids, 3 carbohydrates, 1 nucleic acid, and 4 other pathways. All the detected metabolites were analyzed by heat mapping. As shown in Figure 2A, the heat map indicated significantly increased levels of palmitic acid and docosahexaenoic acid (DHA) at 1.7859 and 1.8403, respectively, while a significantly decreased level of sphingomyelin at 0.445-fold change in CCs of endometriosis patients; in MGCs, there was a downregulation of sphingosine as 0.369-fold change (Figure 3A). Metabolite–metabolite correlations between the tissue of CCs and MGCs between endometriosis and control showed unique profiles. Metabolite–metabolite correlations showed significant correlation coefficients (p < 0.01) in CCs (Figure 2B) and in MGCs (Figure 3B). Moreover, palmitic acid, DHA, and sphingomyelin also played important roles in metabolite correlations in cumulus cells of endometriosis patients, while sphingosine showed played important roles in metabolite correlation of MGCs.
3.3.2 Kyoto encyclopedia of genes and genomes enrichment
Our results showed that CCs and MGCs differentially expressed metabolites in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis between endometriosis patients and control. We used the qualitative difference in metabolite expression quantity of each sample. In hierarchical clustering, the same group of samples can appear in the same cluster through clustering. Meanwhile, metabolites clustered in the same cluster have similar expression patterns and may be in relatively close reaction steps in the metabolic process. Figures 2C, 3C shows the results of the metabolite hierarchy clustering with significant differences in the sample group. The results showed that the choline metabolism in cancer and biosynthesis of unsaturated fatty acids were significantly changed between cumulus cells of endometriosis patients and control (Figure 2C), while the choline metabolism in cancer, glycerophospholipid metabolism, purine metabolism, and ABC transporters were significantly changed between MGCs of endometriosis patients and control (Figure 3C).
3.4 Gene and protein expression in cumulus cells and mural granulosa cells
Concerning the alteration of palmitic acid and sphingomyelin in CCs of endometriosis patients compared to controls, we evaluated the level of gene involvement in ceramide synthesis and autophagic pathway. In detail, genes involved in the ceramide synthesis including Cers1, sptl1, and SMPD1 increased significantly (with p < 0.05) in cumulus cells of women with endometriosis when compared to normal (1.5-fold, 2.4-fold, and 1.6-fold, respectively) (Figure 4). Moreover, the autophagic gene expression including BECN1, LAMP, and P62 increased significantly (with p < 0.05) in CCs of endometriosis (2.8-fold, 2.4-fold, and 2.6-fold, respectively) but not in MGCs of endometriosis (Figure 4A). Our Western blotting results showed that the expression of FASN, Beclin1, and LAMP in CCs of the endometriosis group was higher than that of the control group (Figure 4C).
Figure 4 qRT-PCR expression analysis and protein expression of the ceramide synthesis and the autophagic pathway in CCs and MGCs compared between normal and endometriosis patients. (A) The expression levels of a gene involved in ceramide synthesis (CERS1, SPTL1, SMPD1, SMPD2, and SMPD3) and genes involved in autophagy (BECN1, LAMP, and P62) were analyzed in CCs compared between normal and endometriosis patients. (B) The expression levels of a gene involved in ceramide synthesis (SPHK1, ASAH1, CERS1, UGCG, and SGPP1) and genes involved in autophagy (BECN1, LAMP, and P62) were analyzed in MGCs compared between normal and endometriosis patients. Data (technical triplicates of three biological experiments) are reported as means ± standard error. *, **, *** indicate a statistically significant difference (p< 0.05, 0.01 and 0.001 respectively). Protein expression in CCs (C) and MGCs (D) was determined by Western blotting using antibodies against FASN, ASAH1, Beclin1, LAMP, and Actin. The actin band indicates equal loading of proteins. CCs, cumulus cells; MGCs, mural granulosa cells.
Concerning the alteration of sphingosine in MGCs of endometriosis patients compared to controls, we evaluated the level of gene involvement in sphingosine synthesis, S1P synthesis, ceramide synthesis, and autophagic pathway. In detail, genes involved in sphingosine synthesis including ASAH1 and genes involved in S1P synthesis including SPH1K decreased significantly (with p < 0.05) (0.1-fold and 0.48-fold, respectively), while genes involved in the ceramide synthesis including Cers1 increased significantly (with p < 0.05). In addition, UGCG mRNA levels in MGCs were significantly higher in endometriosis patients than in controls (1.88-fold). However, the autophagic gene expression including BECN1, LAMP, and PC3 decreased in MGCs of endometriosis (Figure 4B). Our Western blotting results showed that the expression of ASAH1, Beclin1, and LAMP in the CC of the endometriosis group was higher than that of the control group (Figure 4D).
4 Discussion
In endometriosis, changes in the follicular microenvironment of the COC lead to poor oocyte quality. As our results showed, the rate of fertilization and the level of estradiol in endometriosis had a tendency to be lower than in healthy women. However, the role of CCs and MGCs in endometriosis is still unclear. CCs around the oocyte play a critical role in supporting the oocyte maturation (32) and in steroidogenesis by promoting P450 aromatase activity for estrogen synthesis, which is important for follicular development and generating a competent oocyte to reach the mature metaphase II (MII) stage and fertilization (33). However, MGCs line the follicular wall and perform endocrine functions as well as assist follicular development (34). Therefore, it is possible that alteration of fatty acids in the COC, especially in CCs, influences the oocyte and fertilization. Adverse effects of endometriosis in the pathophysiology of the CCs have been extensively indicated including alteration in the cell cycle (29) and disruption of signaling pathways of CC growth and development (35, 36). CCs are in transcripts related to metabolism and cell proliferation in mice, whereas MGCs are rich in transcripts related to cell signaling and differentiation. In mice, CCs are enhanced in transcripts related to metabolism and cell proliferation, whereas MGCs are enhanced in transcripts related to cell signaling and differentiation (37). Moreover, CCs were shown to be high in metabolites linked to cholesterol transport and estradiol generation in human preovulatory follicles, while MGCs were enriched in metabolites associated with anti-apoptosis (38). Additionally, previous studies demonstrated that many types of endometriosis also may negatively affect CCs steroidogenesis by reducing the level of P450 aromatase, increasing ROS, and inducing apoptosis (4, 39). A previous report has shown that endometriosis impaired CC steroidogenesis, leading to an imbalance in estrogen synthesis (40). Interestingly, cell death of the follicular cells can be induced through oxidative stress, hyper-androgenemia, and disturbance of gonadotropin hormone. Endometriosis impairs the cell cycle in granulosa cells (29). Recently, many reports indicated that estrogen receptor (ER) stress induction results in alteration of cumulus–oocyte complex (41), activation of ovarian fibrosis (42), and apoptosis of follicular cells (43). Moreover, a previous report found that CCs from endometriosis patients reveal senescence through ER stress-associated endometriosis-associated infertility and senescence phenotype (44).
Sphingolipid metabolism is important for maintaining normal follicular and oocyte maturation through modulating energy metabolism, cell proliferation, apoptosis, and steroid hormone synthesis (13). Previous studies reported that sphingolipid metabolism in endometriosis showed alteration in the serum, plasma, peritoneal fluid, follicular fluid, and endometrial tissue (13, 23, 45–47). However, no information is reported about CCs and MGCs correlated between sphingolipid metabolism and autophagic pathway in endometriosis. Our study focuses on the role of these follicular cells in terms of metabolism and integrated biochemical studies. Interestingly, our results demonstrated that the metabolites related to palmitic acid, sphingolipid metabolism, and autophagic cell death were enriched in only CCs, indicating that sphingolipid metabolism plays a critical role in follicle and oocyte growth. In the literature, previous evidence showed that sphingolipid metabolism regulated cell death at the molecular level. In this study, we used a metabolomic assay to evaluate the sphingolipid metabolism of CCs and MGCs in endometriosis systemically. Our results represented the alteration of the biomolecules presented in CCs and MGCs in endometriosis patients. The metabolic profiling of CCs has evidenced that the level of palmitic acid was increased in CCs but not in MGCs according to upregulated FASN gene expression in CCs. Therefore, the metabolite results suggested that palmitic acid may play an essential role in CCs and MGCs of endometriosis. Several studies suggested that elevation of palmitic acid leads to inflammation in several tissues during pathophysiologic progression (48, 49). Additionally, recent studies demonstrated that the accumulation of palmitic acid and its downstream metabolism leads to inflammation and cell death (50, 51). Y. M. Mu’s team showed that saturated FFAs, palmitic acid, and stearic acid induce apoptosis in human granulosa cells in vitro (17). In the COC, oocytes seem to be more sensitive to fatty acid stress. Moreover, several reports suggested elevation of fatty acids in the blood, and follicular fluid affected oocyte maturation (14, 15, 52–54). Previous reports on obesity elevated palmitic acid can cause insulin resistance (IR) and impairment of glucose metabolism in ovarian GCs (55). This evidence supports that palmitic acid is one of the key molecules in lipid metabolism. In addition, the alteration of sphingolipidomic data in serum, peritoneal fluid, and endometrial tissue in endometriosis patients reveals the critical role of lipid metabolism in pathological progression (13). Interestingly, we suggested that sphingolipid metabolism can be used by cumulus cells, which mediate its effect on oocyte maturation.
The key molecules in sphingolipid metabolism are ceramide and sphingosine-1-phosphate (S1P), which play a critical role in survival and cell death. Ceramide is a bio-effector molecule that mediates cell death, whereas S1P induces cell proliferation (56). Ceramide is a central molecule of sphingolipid metabolism and serves as the precursor for several sphingolipids including sphingomyelin (SM) and Cer-1-phosphate and glucosylceramide (GlcCer).
Recent studies also demonstrated that elevation and accumulation of GlcCer are associated with many human diseases including Gaucher’s disease, polycystic kidney disease, diabetes, and endometriosis, leading to overproliferation (57). Both in vitro and in vivo studies indicated that GlcCers inhibited GlcCer synthase (GCS) through several processes leading to cell death (58–60). Moreover, endometriosis patients reveal an association between GCS level and serum and PF GlcCer accumulation, indicating abnormalities in endometrial proliferation. Interestingly, an increasing level of Cer would have induced apoptosis in endometrial cells (61, 62). According to our results, the elevation of genes involved with the ceramide synthesis (Cers1, sptl1, and SMPD1) in CCs associated with elevated palmitic acid in metabolomic analysis supports that sphingolipid metabolism may be associated with the pathophysiology of endometriosis (13). However, sphingolipid metabolites in MGCs found that low levels of palmitic acids, sphingomyelin, and sphingosine are associated with the low level of the ceramide synthesis gene. Our studies indicated that sphingolipid metabolism is mediated pathway through regulating cell survival and cell death.
Autophagy is one of the types of cell death that modulate cellular homeostasis by supplying the cell with energy and metabolites and preventing oxidative stress to maintain normal cellular function under stress conditions (63, 64). A previous report demonstrated that the key molecule of sphingolipid metabolism, ceramide (Cer), induced autophagic cell death (65). Our results also revealed that cluster genes and proteins in autophagy (BECN1 and LAMP2) significantly increased in association with an increased level of palmitic acid and the ceramide synthesis gene in CCs under the pathological condition of endometriosis. This evidence indicated that CCs may undergo cell death, which may directly affect the oocyte quality. This evidence indicates that sphingolipid metabolism in surrounding cumulus cells plays an important role in the maturation of oocytes and steroidogenesis. However, in MGCs, this cluster of genes decreased in endometriosis, which correlated with the role of MGCs in the angiogenesis compartment of follicular cells.
A recent study demonstrates the association between apoptosis and survival molecules in cumulus cells can be used as a marker for estimating oocyte quality (66). Interestingly, our study indicated that the poor oocyte quality in patients with endometriosis is related to a higher level of sphingolipid metabolism correlated with the increase of autophagic gene and protein levels in the CCs but not in the MGCs. A balance in sphingolipid metabolism in CCs is associated with oocyte maturation and further development of an embryo.
The number of cumulus cells closely associated with energy sufficiency cell proliferation directly affects oocyte development. The cell depends on apoptosis and autophagic signaling pathways; therefore, saturating FFA metabolites may potentially promote cell death (17). We can suggest that this procedure could be integrated through biochemical studies in the follicular cells. It is better known that oocyte quality depends on the follicular microenvironment, and adequate bidirectional signaling between COC is important for both oocyte and cumulus cell competence acquisition (67). The most relevant and significant data of metabolites we obtained on the CCs and MGCs in endometriosis patients (autophagic genes and proteins decrease linked to the increase of lipid levels) could then be used to help to identify oocytes with higher capacity for development and fertilization to enhance the possibility of success of in vitro fertilization (IVF) and improving the probability of pregnancy in endometriosis patients.
The metabolic characterization of CCs and MGCs to integrate the multiple approaches from biochemical analysis to gene analysis in association with pathogenesis and disease progression has gained more value in understanding endometriosis. Our work demonstrates how integrating the different approaches of the ‘omic’ sciences from metabolomics to gene expression can help in better understanding the mechanisms of endometriosis progression and can lead to improvement of diagnosis. The advantage of systems biology is the integration of proteomics, transcriptomics, and metabolomic information to obtain a better understanding mechanism of human pathophysiology. Moreover, insight into the differences between the metabolic profile of CCs and MGCs may lead to the discovery of a new biomarker for estimating the quality of oocytes in endometriosis patients. It is possible that CCs and MGCs associated with competent oocytes have fulfilled their role of sustaining oocyte development, and their response to metabolic alteration is under consideration. Taken together, these findings may optimize the development of novel methods to improve oocyte selection strategies in assisted reproductive technology (ART) protocols and shed light on the molecular mechanisms governing oocyte–CC–MGC cross-talk. The functional roles of those metabolites in follicles related to oocyte growth, maturation, and subsequent embryonic developmental need to be further explored and studied.
5 Conclusion
The results indicate that CCs and MGCs reveal differences in terms of metabolites. CCs show the alteration of lipid metabolism correlated with autophagic cell death, while MGCs show no alteration, as shown in Figure 5. Our finding provides new insight into follicular cells in endometriosis.
Figure 5 Schematic representation of the metabolomics of sphingolipids and autophagic gene expression between the cumulus cells and the mural granulosa cells. Direction of arrows represents the direction of possible activity of sphingolipid pathways. Red colors of arrows represent increased level of metabolite or gene, and blue colors of arrows represent decreased level of metabolite or gene. SPT, serine palmitoyltransferase; CerS, ceramide synthase; DES, dihydroceramide desaturase; SMS, sphingomyelin synthase; SMase, sphingomyelinase; SPP, sphingosine-1-phosphate phosphatase; SPHIK, sphingosine kinase; UGCG, UDP-glucose ceramide glucosyltransferase; BECN1, Beclin1; LAMP, lysosome-associated membrane protein.
Data availability statement
The original contributions presented in the study are included in the article/supplementary materials. Further inquiries can be directed to the corresponding author.
Ethics statement
This study was reviewed and approved by the Medical and Life Science Ethics Committee of Tongji University (approval No. 2017yxy001). The patients/participants provided their written informed consent to participate in this study.
Author contributions
R-CC: study design, manuscript writing, and critical discussion. BT and E-MG.: study design, experiment execution, data analysis, manuscript writing, and critical discussion. KG: experiment execution, manuscript writing, and discussion. Y-BL, LW, and XD: experiment execution and manuscript review. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key R&D Program of China (No. 2017YFC1002003) and the Ministry of Science and Technology of China (No. 2017YFC1001601).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Parasar P, Ozcan P, Terry KL. Endometriosis: Epidemiology, diagnosis and clinical management. Curr Obstet Gynecol Rep (2017) 6(1):34–41. doi: 10.1007/s13669-017-0187-1
2. Macer ML, Taylor HS. Endometriosis and infertility: a review of the pathogenesis and treatment of endometriosis-associated infertility. Obstet Gynecol Clin North Am (2012) 39(4):535–49. doi: 10.1016/j.ogc.2012.10.002
3. Sanchez AM, Vanni VS, Bartiromo L, Papaleo E, Zilberberg E, Candiani M, et al. Is the oocyte quality affected by endometriosis? a review of the literature. J Ovarian Res (2017) 10(1):43. doi: 10.1186/s13048-017-0341-4
4. Sanchez AM, Somigliana E, Vercellini P, Pagliardini L, Candiani M, Vigano P. Endometriosis as a detrimental condition for granulosa cell steroidogenesis and development: From molecular alterations to clinical impact. J Steroid Biochem Mol Biol (2016) 155(Pt A):35–46. doi: 10.1016/j.jsbmb.2015.07.023
5. Scott R 3rd, Zhang M, Seli E. Metabolism of the oocyte and the preimplantation embryo: implications for assisted reproduction. Curr Opin Obstet Gynecol (2018) 30(3):163–70. doi: 10.1097/GCO.0000000000000455
6. Dutta M, Singh B, Joshi M, Das D, Subramani E, Maan M, et al. Metabolomics reveals perturbations in endometrium and serum of minimal and mild endometriosis. Sci Rep (2018) 8(1):6466. doi: 10.1038/s41598-018-23954-7
7. Domínguez F, Ferrando M, Díaz-Gimeno P, Quintana F, Fernández G, Castells I, et al. Lipidomic profiling of endometrial fluid in women with ovarian endometriosis†. Biol Reprod (2017) 96(4):772–9. doi: 10.1093/biolre/iox014
8. Cordeiro FB, Cataldi TR, Perkel KJ, do Vale Teixeira da Costa L, Rochetti RC, Stevanato J, et al. Lipidomics analysis of follicular fluid by ESI-MS reveals potential biomarkers for ovarian endometriosis. J Assist Reprod Genet (2015) 32(12):1817–25. doi: 10.1007/s10815-015-0592-1
9. Vicente-Muñoz S, Morcillo I, Puchades-Carrasco L, Payá V, Pellicer A, Pineda-Lucena A. Nuclear magnetic resonance metabolomic profiling of urine provides a noninvasive alternative to the identification of biomarkers associated with endometriosis. Fertil Steril (2015) 104(5):1202–9. doi: 10.1016/j.fertnstert.2015.07.1149
10. Dutta M, Joshi M, Srivastava S, Lodh I, Chakravarty B, Chaudhury K. A metabonomics approach as a means for identification of potential biomarkers for early diagnosis of endometriosis. Mol Biosyst (2012) 8(12):3281–7. doi: 10.1039/c2mb25353d
11. Lee YH, Cui L, Fang J, Chern BS, Tan HH, Chan JK. Limited value of pro-inflammatory oxylipins and cytokines as circulating biomarkers in endometriosis - a targeted 'omics study. Sci Rep (2016) 6:26117. doi: 10.1038/srep26117
12. Carta G, Murru E, Banni S, Manca C. Palmitic acid: Physiological role, metabolism and nutritional implications. Front Physiol (2017) 8:902. doi: 10.3389/fphys.2017.00902
13. Lee YH, Tan CW, Venkatratnam A, Tan CS, Cui L, Loh SF, et al. Dysregulated sphingolipid metabolism in endometriosis. J Clin Endocrinol Metab (2014) 99(10):E1913–21. doi: 10.1210/jc.2014-1340
14. Leroy JL, 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. doi: 10.1530/rep.1.00735
15. Jungheim ES, Macones GA, Odem RR, Patterson BW, Lanzendorf SE, Ratts VS, et al. Associations between free fatty acids, cumulus oocyte complex morphology and ovarian function during in vitro fertilization. Fertil Steril (2011) 95(6):1970–4. doi: 10.1016/j.fertnstert.2011.01.154
16. Vanholder T, Leroy JL, Soom AV, 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-2):33–44. doi: 10.1016/j.anireprosci.2004.09.006
17. Mu YM, Yanase T, Nishi Y, Tanaka A, Saito M, Jin CH, et al. Saturated FFAs, palmitic acid and stearic acid, induce apoptosis in human granulosa cells. Endocrinology (2001) 142(8):3590–7. doi: 10.1210/endo.142.8.8293
18. Bouet PE, Chao de la Barca JM, El Hachem H, Descamps P, Legendre G, Reynier P, et al. Metabolomics shows no impairment of the microenvironment of the cumulus-oocyte complex in women with isolated endometriosis. Reprod BioMed Online (2019) 39(6):885–92. doi: 10.1016/j.rbmo.2019.08.001
19. Chilvers RA, Bodenburg YH, Denner LA, Urban RJ. Development of a novel protocol for isolation and purification of human granulosa cells. J Assist Reprod Genet (2012) 29(6):547–56. doi: 10.1007/s10815-012-9739-5
20. Aghadavod E, Zarghami N, Farzadi L, Zare M, Barzegari A, Movassaghpour AA, et al. Isolation of granulosa cells from follicular fluid; applications in biomedical and molecular biology experiments. Adv BioMed Res (2015) 4:250. doi: 10.4103/2277-9175.170675
21. Dołowy M. Influence of different chromatographic conditions on the level of detection and quantitation of spironolactone by means of TLC-densitometry. J Anal Methods Chem (2019) 2019:8792783. doi: 10.1155/2019/8792783
22. Sadok I, Gamian A, Staniszewska MM. Chromatographic analysis of tryptophan metabolites. J Sep Sci (2017) 40(15):3020–45. doi: 10.1002/jssc.201700184
23. Marianna S, Alessia P, Susan C, Francesca C, Angela S, Francesca C, et al. Metabolomic profiling and biochemical evaluation of the follicular fluid of endometriosis patients. Mol Biosyst (2017) 13(6):1213–22. doi: 10.1039/C7MB00181A
24. Gómez-Caudillo L, Ortega-Lozano AJ, Martínez-Batallar ÁG, Rosas-Vargas H, Minauro-Sanmiguel F, Encarnación-Guevara S. Principal component analysis on LC−MS/MS and 2DE−MALDI−TOF in glioblastoma cell lines reveals that mitochondria act as organelle sensors of the metabolic state in glioblastoma. Oncol Rep (2020) 44(2):661–73. doi: 10.3892/or.2020.7625
25. Gaston KW, Limbach PA. The identification and characterization of non-coding and coding RNAs and their modified nucleosides by mass spectrometry. RNA Biol (2014) 11(12):1568–85. doi: 10.4161/15476286.2014.992280
26. Donnelly DP, Rawlins CM, DeHart CJ, Fornelli L, Schachner LF, Lin Z, et al. Best practices and benchmarks for intact protein analysis for top-down mass spectrometry. Nat Methods (2019) 16(7):587–94. doi: 10.1038/s41592-019-0457-0
27. Luddi A, Gori M, Marrocco C, Capaldo A, Pavone V, Bianchi L, et al. Matrix metalloproteinases and their inhibitors in human cumulus and granulosa cells as biomarkers for oocyte quality estimation. Fertil Steril (2018) 109(5):930–9.e3. doi: 10.1016/j.fertnstert.2018.01.030
28. Zhang J, Liu J, Zhu K, Hong Y, Sun Y, Zhao X, et al. Effects of BMAL1-SIRT1-positive cycle on estrogen synthesis in human ovarian granulosa cells: an implicative role of BMAL1 in PCOS. Endocrine (2016) 53(2):574–84. doi: 10.1007/s12020-016-0961-2
29. Toya M, Saito H, Ohta N, Saito T, Kaneko T, Hiroi M. Moderate and severe endometriosis is associated with alterations in the cell cycle of granulosa cells in patients undergoing in vitro fertilization and embryo transfer. Fertil Steril (2000) 73(2):344–50. doi: 10.1016/S0015-0282(99)00507-5
30. Read CC, Dyce PW. All-trans retinoic acid exposure increases connexin 43 expression in cumulus cells and improves embryo development in bovine oocytes. Mol Reprod Dev (2019) 86(12):1865–73. doi: 10.1002/mrd.23274
31. Chang L, Ni J, Beretov J, Wasinger VC, Hao J, Bucci J, et al. Identification of protein biomarkers and signaling pathways associated with prostate cancer radioresistance using label-free LC-MS/MS proteomic approach. Sci Rep (2017) 7:41834. doi: 10.1038/srep41834
32. Uyar A, Torrealday S, Seli E. Cumulus and granulosa cell markers of oocyte and embryo quality. Fertil Steril (2013) 99(4):979–97. doi: 10.1016/j.fertnstert.2013.01.129
33. Dumesic DA, Meldrum DR, Katz-Jaffe MG, Krisher RL, Schoolcraft WB. Oocyte environment: follicular fluid and cumulus cells are critical for oocyte health. Fertil Steril (2015) 103(2):303–16. doi: 10.1016/j.fertnstert.2014.11.015
34. Turathum B, Gao E-M, Chian R-C. The function of cumulus cells in oocyte growth and maturation and in subsequent ovulation and fertilization. Cells (2021) 10(9):2292. doi: 10.3390/cells10092292
35. Sanchez AM, Viganò P, Quattrone F, Pagliardini L, Papaleo E, Candiani M, et al. The WNT/β-catenin signaling pathway and expression of survival promoting genes in luteinized granulosa cells: endometriosis as a paradigm for a dysregulated apoptosis pathway. Fertil Steril (2014) 101(6):1688–96. doi: 10.1016/j.fertnstert.2014.02.040
36. Sanchez AM, Viganò P, Somigliana E, Panina-Bordignon P, Vercellini P, Candiani M. The distinguishing cellular and molecular features of the endometriotic ovarian cyst: from pathophysiology to the potential endometrioma-mediated damage to the ovary. Hum Reprod Update (2014) 20(2):217–30. doi: 10.1093/humupd/dmt053
37. Wigglesworth K, Lee K-B, Emori C, Sugiura K, Eppig JJ. Transcriptomic diversification of developing cumulus and mural granulosa cells in mouse ovarian Follicles1. Biol Reproduc (2015) 92(1):23-1. doi: 10.1095/biolreprod.114.121756
38. Gao E-M, Turathum B, Wang L, Zhang D, Liu Y-B, Tang R-X, et al. The differential metabolomes in cumulus and mural granulosa cells from human preovulatory follicles. Reprod Sci (2022) 29(4):1343–56. doi: 10.1007/s43032-021-00691-3
39. Nakahara K, Saito H, Saito T, Ito M, Ohta N, Takahashi T, et al. Ovarian fecundity in patients with endometriosis can be estimated by the incidence of apoptotic bodies. Fertil Steril (1998) 69(5):931–5. doi: 10.1016/S0015-0282(98)00038-7
40. Tummon IS, Maclin VM, Radwanska E, Binor Z, Dmowski WP. Occult ovulatory dysfunction in women with minimal endometriosis or unexplained infertility. Fertil Steril (1988) 50(5):716–20. doi: 10.1016/S0015-0282(16)60304-7
41. 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. Fertil Steril (2012) 97(6):1438–43. doi: 10.1016/j.fertnstert.2012.02.034
42. Takahashi N, Harada M, Hirota Y, Nose E, Azhary JM, Koike H, et al. Activation of endoplasmic reticulum stress in granulosa cells from patients with polycystic ovary syndrome contributes to ovarian fibrosis. Sci Rep (2017) 7(1):10824. doi: 10.1038/s41598-017-11252-7
43. Zeng Y, Sun H, Li Y, Shao M, Han P, Yu X, et al. Exposure to triptolide affects follicle development in NIH mice: Role of endoplasmic reticulum stress in granulosa cell apoptosis. Hum Exp Toxicol (2017) 36(1):82–92. doi: 10.1177/0960327116638725
44. Lin X, Dai Y, Tong X, Xu W, Huang Q, Jin X, et al. Excessive oxidative stress in cumulus granulosa cells induced cell senescence contributes to endometriosis-associated infertility. Redox Biol (2020) 30:101431. doi: 10.1016/j.redox.2020.101431
45. Lee YH, Yang JX, Allen JC, Tan CS, Chern BSM, Tan TY, et al. Elevated peritoneal fluid ceramides in human endometriosis-associated infertility and their effects on mouse oocyte maturation. Fertil Steril (2018) 110(4):767–77.e5. doi: 10.1016/j.fertnstert.2018.05.003
46. Vouk K, Hevir N, Ribić-Pucelj M, Haarpaintner G, Scherb H, Osredkar J, et al. Discovery of phosphatidylcholines and sphingomyelins as biomarkers for ovarian endometriosis. Hum Reprod (2012) 27(10):2955–65. doi: 10.1093/humrep/des152
47. Vouk K, Ribič-Pucelj M, Adamski J, Rižner TL. Altered levels of acylcarnitines, phosphatidylcholines, and sphingomyelins in peritoneal fluid from ovarian endometriosis patients. J Steroid Biochem Mol Biol (2016) 159:60–9. doi: 10.1016/j.jsbmb.2016.02.023
48. Wang Y, Qian Y, Fang Q, Zhong P, Li W, Wang L, et al. Saturated palmitic acid induces myocardial inflammatory injuries through direct binding to TLR4 accessory protein MD2. Nat Commun (2017) 8:13997. doi: 10.1038/ncomms13997
49. Schaffer JE. Lipotoxicity: when tissues overeat. Curr Opin Lipidol (2003) 14(3):281–7. doi: 10.1097/00041433-200306000-00008
50. Korbecki J, Bajdak-Rusinek K. The effect of palmitic acid on inflammatory response in macrophages: an overview of molecular mechanisms. Inflammation Res (2019) 68(11):915–32. doi: 10.1007/s00011-019-01273-5
51. Yu G, Luo H, Zhang N, Wang Y, Li Y, Huang H, et al. Loss of p53 sensitizes cells to palmitic acid-induced apoptosis by reactive oxygen species accumulation. Int J Mol Sci (2019) 20(24):6268. doi: 10.3390/ijms20246268
52. 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. doi: 10.1371/journal.pone.0023183
53. 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. doi: 10.1095/biolreprod.110.088815
54. 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(4):562–73. doi: 10.1210/me.2011-1362
55. Xu L, Wang W, Zhang X, Ke H, Qin Y, You L, et al. Palmitic acid causes insulin resistance in granulosa cells via activation of JNK. J Mol Endocrinol (2019) 62(4):197–206. doi: 10.1530/JME-18-0214
56. Ogretmen B. Sphingolipid metabolism in cancer signalling and therapy. Nat Rev Cancer (2018) 18(1):33–50. doi: 10.1038/nrc.2017.96
57. Marsh NL, Elias PM, Holleran WM. Glucosylceramides stimulate murine epidermal hyperproliferation. J Clin Invest (1995) 95(6):2903–9. doi: 10.1172/JCI117997
58. Natoli TA, Smith LA, Rogers KA, Wang B, Komarnitsky S, Budman Y, et al. Inhibition of glucosylceramide accumulation results in effective blockade of polycystic kidney disease in mouse models. Nat Med (2010) 16(7):788–92. doi: 10.1038/nm.2171
59. Uchida Y, Murata S, Schmuth M, Behne MJ, Lee JD, Ichikawa S, et al. Glucosylceramide synthesis and synthase expression protect against ceramide-induced stress. J Lipid Res (2002) 43(8):1293–302. doi: 10.1194/jlr.M100442-JLR200
60. Deng W, Li R, Guerrera M, Liu Y, Ladisch S. Transfection of glucosylceramide synthase antisense inhibits mouse melanoma formation. Glycobiology (2002) 12(3):145–52. doi: 10.1093/glycob/12.3.145
61. Dmowski WP, Gebel H, Braun DP. Decreased apoptosis and sensitivity to macrophage mediated cytolysis of endometrial cells in endometriosis. Hum Reprod Update (1998) 4(5):696–701. doi: 10.1093/humupd/4.5.696
62. Chrobak A, Sieradzka U, Sozański R, Chełmońska-Soyta A, Gabryś M, Jerzak M. Ectopic and eutopic stromal endometriotic cells have a damaged ceramide signaling pathway to apoptosis. Fertil Steril (2009) 92(6):1834–43. doi: 10.1016/j.fertnstert.2008.09.035
63. Levine B, Klionsky DJ. Development by self-digestion: molecular mechanisms and biological functions of autophagy. Dev Cell (2004) 6(4):463–77. doi: 10.1016/S1534-5807(04)00099-1
64. Singh R, Cuervo AM. Autophagy in the cellular energetic balance. Cell Metab (2011) 13(5):495–504. doi: 10.1016/j.cmet.2011.04.004
65. Young MM, Kester M, Wang HG. Sphingolipids: regulators of crosstalk between apoptosis and autophagy. J Lipid Res (2013) 54(1):5–19. doi: 10.1194/jlr.R031278
66. Bosco L, Chiarelli R, Roccheri MC, Matranga D, Ruvolo G. Relationship between apoptosis and survival molecules in human cumulus cells as markers of oocyte competence. Zygote (2017) 25(5):583–91. doi: 10.1017/S0967199417000429
Keywords: mural granulosa cells, cumulus cells, metabolomes, sphingolipid metabolism, autophagy
Citation: Turathum B, Gao E-M, Grataitong K, Liu Y-B, Wang L, Dai X and Chian R-C (2022) Dysregulated sphingolipid metabolism and autophagy in granulosa cells of women with endometriosis. Front. Endocrinol. 13:906570. doi: 10.3389/fendo.2022.906570
Received: 28 March 2022; Accepted: 07 July 2022;
Published: 03 August 2022.
Edited by:
Aniel Sanchez Puente, Lund University, SwedenReviewed by:
Kui Liu, The University of Hong Kong, Hong Kong SAR, ChinaKangpu Xu, Weill Cornell Medicine, Cornell University, United States
Copyright © 2022 Turathum, Gao, Grataitong, Liu, Wang, Dai and Chian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ri-Cheng Chian, cmNoaWFuQDEyNi5jb20=
†These authors have contributed equally to this work and share first authorship