3.1 Bioinformatics
From the GSE193812 dataset, a total of 472 differentially expressed genes in adipose tissue were identified, as shown in Figure 1A (Table 1). In the GSE10946 dataset, a total of 61 DEGs were found in the granule cell cluster, as shown in Figure 1B (Table 2). The common differentially expressed gene was SLC16A6 (Figure 1C). The GSE10946 dataset was selected to verify the expression of SLC16A6 and perform a ROC curve analysis. SLC16A6 was significantly increased in the granulosa cells of the PCOS group (P < 0.05) (Figure 1D). The AUC was 0.866, indicating that SLC16A6 expression in the granulosa cells of the PCOS group was strongly reliable (Figure 1E). A total of 356 DEGs were detected between the SLC16A6 high- and low-expression groups (Figure 1F). GO and KEGG analyses were performed on 356 DEGs, and a total of 380 significant GO terms were identified (Table 3). The top 30 GO terms are shown in Figure 1G. These pathways included "insulin-like growth factor receptor binding", "insulin-like growth factor II binding", insulin-like growth factor I binding, estradiol 17-beta-dehydrogenase activity, steroid dehydrogenase activity, etc., and "sex hormone-related pathways" (Figure 1G). Twenty significantly different pathways were identified according to the KEGG enrichment results (Table 4): "TGF-beta signaling pathway", "ovarian steroidogenesis", "proline signaling pathway", "regulation of lipolysis in" "adipocytes", "type II diabetes mellitus" and other inflammation-, intraovarian sterol-, and diabetes-related pathways (Figure 1H).
3.2, Network pharmacology
A total of 15 corresponding components associated with SLC16A6 (Table 5), which contained three coffee-related substances, namely, "coffee acid", "caffeic", and "trans-caffeic acid", were obtained from the herb database. Seven effective substances in coffee were collected, namely, "coffee acid", "caffeic", "trans-caffeic acid", "cafestol", "trigonelline", "caffeine", and "chlorogenic acid". A total of 147 targets corresponding to 7 substances were identified (Table 6) (Figure 2A). The database revealed 3240 PCOS genes and 73 intersecting genes between active coffee ingredients and PCOS (Table 7) (Figure 2B). The PPI network diagram of the intersecting genes was constructed (Figure 2C), and the core targets of the PPI network diagram were analyzed. The top three targets were "INS", "TP53", and "TNF" (Figure 2D). A network diagram of coffee components and therapeutic targets for PCOS treatment was constructed (Figure 2E), and the corresponding targets of each component were statistically analyzed. Forty-six targets corresponding to caffeine were identified (Figure 2F), which was far more than the number of other components, indicating its core role. KEGG and GO analyses were performed on the 73 intersecting genes. Figure 3A shows the involvement of genes in each pathway, and the INS genes were involved in multiple pathways. Figure 3B shows the top 30 GO terms, including "regulation of insulin secretion", "glucose transmembrane transporter activity", "D-glucose transmembrane transporter" activity" and "other insulin or glucose-related pathways". After the KEGG analysis, the intersecting gene enrichment pathways were closely related to fat metabolism, insulin, and oocyte development, as shown in Figure 3C.
3.3. Experimental verification
3.3.1 Reproductive and metabolic phenotypes of obese PCOS rats
In this study, the obese PCOS group was confirmed by changes in body weight. By the end of the induction (4 weeks), the body weight of the model group increased significantly (339.5±29.8 g) compared with that of the control group (234.7±11.8 g) (Fig. 4A). Compared with the control group (35.04%), the weight change in the model group was significantly greater (95.35%) (Figure 4B).The vaginal pictures of the control group suggested that the control group had a regular, complete estrous cycle of 4 to 5 days, while the vaginal smear of the model group showed aperiorism, mainly in the interestrus period (Figure 4C.4D). The HE-stained sections from the control group had a normal ovarian morphology. Microscopy revealed that the ovaries of the control group had multiple corpora lutea and follicles at various stages of development, and multilayer granulosa cells were neatly arranged. The HE-stained sections of the ovaries in the model group showed an obvious polycystic morphology with typical polycystic pathological changes under the microscope, with a large number of dilated cystic follicles, no or little luteal tissue, corona radiata and oocytes, and a significant reduction in the number of granulosa cell layers (Fig. 4E).
In conclusion, rats subjected to long-term intragastric administration of letrozole and a high-fat diet (model group) exhibited obvious obese PCOS-like changes.
3.3.2 Caffeine treatment restored the estrous cycle in obese PCOS rats
Vaginal exfoliative cytology was performed on the four groups of rats, and continuous observations were recorded until the end of the experiment. The results of vaginal smears showed that all rats in the control group had a regular and complete estrous cycle, which lasted for 4-5 days. In the model group, the estrous cycle was disrupted in 6 rats, and the normal rate of the estrous cycle was only 25%, which was mainly manifested as the prolongation of interestrus. The CFYN group and GLP-1 group were similar to the control group. Five rats in the CFYN group had a normal estrous cycle after the experiment, and the normal rate of the estrous cycle was 83%. Six rats in the GLP-1 group had normal estrous cycles at the end of the experiment, and the percentage of normal estrous cycles was 75% (Fig. 6A). Based on histological HE staining, the study revealed that the model group still exhibited significant PCOS-like changes in the ovaries, and the CFYN group and the GLP-1 group exhibited a reduced number of cystic follicles, an increased number of corpus lutea, and improvements in polycystic changes (Figure 6B).
3.3.3 Effects of caffeine on body weight and endocrine metabolism in obese PCOS rats
In this study, the successful establishment of an obese PCOS model was confirmed by changes in body weight. After the establishment of the model (after 4 weeks), the obese PCOS rats were randomly divided into three groups: the untreated obese PCOS group (model group), caffeine gavage treatment group (CFYN group), and oral semaglutide gavage treatment group (GLP-1 group). After 3 weeks of treatment, no significant difference in body weight was observed between the CFYN group (373.3±66.2 g) and the model group (401.5±31.7 g), while the GLP-1 group (281.8±36.3 g) and the model group (401.5±31.7 g) showed a significant reduction in body weight (Fig. 5A).
Serum samples were collected after a 12 h fast. The GLP-1, fasting plasma glucose and fasting insulin levels were compared among the model, CFYN, GLP-1 and control groups, and the corresponding insulin resistance index was calculated. A significant difference was observed between the model group and the control group (P<0.05). The FINS level of the CFYN group was significantly improved compared with the model group (P<0.05), and a significant difference was not observed between the CFYN group and the GLP-1 treatment group (Fig. 5B-E).
In this study, the serum testosterone, FSH, and LH levels were analyzed in each group (Figure 5F-H). The data of this study showed that the testosterone level of the model group (11.74±1.2 nmol/L) was significantly higher than that of the control group (7.23±1.35 nmol/L) (P<0.05), which showed the typical hyperandrogenemia characteristics of PCOS. A significant difference in testosterone levels was not observed between the CFYN group and the GLP-1 group. (P<0.05) (Fig. 5F). The FSH level in the model group (5.9±0.84 U/l) was significantly higher than that in the control group (2.79±0.38 U/l) (P<0.05), and a significant difference was detected between the CFYN group and the model group (P<0.05); however, no significant difference was observed between the CFYN group and the GLP-1 group (Fig. 5G). Moreover, the LH level in the model group (8.69±1.49 mIU/ml) was significantly higher than that in the control group (3.69±1.27 mIU/ml) (P<0.05), which is a typical characteristic of PCOS. A significant difference was observed between the CFYN group and the model group (P<0.05), but no significant difference was detected between the CFYN group and the GLP-1 group (Fig. 5H), which indicated that CFYN could significantly restore the ovarian function of PCOS rats and improve their symptoms.
3.3.4 Effect of caffeine on the expression of the SLC16A6 transporter gene in obese PCOS rats
RT‒PCR was used to detect the gene expression of the SLC16A6 transporter in different groups of rats. Compared with the normal control group, the expression of the SLC16A6 gene in the model group was significantly decreased, and significant differences were not observed among the control group, CFYN group and GLP-1 group. However, significant differences were noted between the CON group, CFYN group and GLP-1 group and the model group (P<0.05) (Fig. 7).