Involvement of ferroptosis in the granulosa cells proliferation of PCOS through the circRHBG/miR-515/SLC7A11 axis.

BACKGROUND: Polycystic ovary syndrome (PCOS) is the most common endocrine disease among women, but its etiology remains unknown. In recent years, many circular RNAs (circRNAs) have been confirmed to be related to the development of PCOS. The role and mechanism of circRNA in the development of PCOS need to be further explored.
METHODS: In the present study, we used the circRNA chip to detect the difference in the expression of circRNA in the granulosa cells of PCOS patients and controls. Five upregulated circRNAs were then selected for quantitative reverse transcription polymerase chain reaction (qRT-PCR) verification, and circRHBG was found to be upregulated in PCOS. Subsequently, Cell Counting Kit-8 assay and EdU assay were used to observe the effect of circRHBG on the proliferation of KGN and SVOG cells. Furthermore, the pairwise binding relationship between circRHBG/miR515-5 and miR515-5p/SLC7A11 was verified by luciferase reporter assay. The interaction between circRHBG and SLC7A11 was detected with qRT-PCR and Western blot.
RESULTS: CircRNA high-throughput chips and qRT-PCR verified that circRHBG was significantly upregulated in granular cells of PCOS patients. Knockdown of circRHBG inhibits KGN and SVOG cell proliferation. Luciferase reporter assays and Ago2-RIP detection showed that circRHBG competes with SLC7A11 to bind to miR-515-5p. Subsequent experiments verified knockdown of circRHBG promotes ferroptosis in PCOS.
CONCLUSIONS: circRHBG inhibits ferroptosis in PCOS cells through the circRHBG/miR-515-5p/SLC7A11 axis in PCOS, which may provide new diagnostic molecular markers and therapeutic targets for PCOS. 2021 Annals of Translational Medicine. All rights reserved.

Introduction

Polycystic ovary syndrome (PCOS) is a heterogeneous disease, which is defined as excessive androgen, ovarian dysfunction while other specific diagnoses need to be ruled out (1). There's no test to definitively diagnose PCOS, three PCOS diagnostic standards currently exist (1-3). The incidence of PCOS in premenopausal women is reported to be 6–20%, and it is the most common endocrine disease in adult women (4,5). The cause of PCOS has not yet been clarified to a large extent, and increasing evidence shows that PCOS is a complex polygenic genetic disease that is greatly affected by environmental factors (6,7). PCOS is related to abdominal obesity, insulin resistance, obesity, metabolic disorders, and cardiovascular risk factors (3,8). With the development of next-generation sequencing technology, more and more genes involved in regulating the progression of PCOS have been discovered (9,10).

Circular RNA (circRNA) can be synthesized by reverse splicing of coding gene exons or classical mRNA splicing without degradation of the intron lasso (11,12). In addition to the effect of exon reverse splicing on its parental gene mRNA splicing, studies have reported that circRNA can regulate gene expression in the nucleus, competitively binding miRNA or protein, and acts as the backbone of the circRNA-protein complex (13-16). Many studies have confirmed the role of circRNA in various diseases; for example, circRNA HIPK3 causes endothelial proliferation and vascular dysfunction in diabetic retinopathy (17). Cancer-specific PML-RARA translocation-derived fusion circRNA can also promote transformation and cell survival (18). The role and mechanism of circRNA in PCOS has gradually become a research hotspot (17,19).

As a common type of non-coding RNA, many microRNAs (miRNAs) have also been found to play a role in the occurrence and development of PCOS. Studies have found that, compared with non-PCOS women, the expression of miRNA in the follicular fluid of PCOS women is very different, and the expression of potential target genes related to the PCOS phenotype has also changed (20,21). Some researchers have constructed a circRNA-microRNA interaction network. The results showed that the differentially expressed circRNA in PCOS is related to inflammation, bacterial infection, and oxidative stress (22). However, to date, these studies have only investigated discrepancy of circRNAs between PCOS patients and controls but have not further studied the specific mechanism. The circRHBG-microRNA pathway in PCOS has not been well studied.

In the present study, we screened out circRHBG with elevated expression in PCOS patients through circRNA microarray and speculated that it plays an important role in PCOS. Bioinformatics analysis was then applied to predict the downstream target of circRHBG, and we found that miR-515 is a potential target of it. In addition, we found that the upregulation of circRHBG expression inhibited ferroptosis in KGN and SVOG cells through the circRHBG/miR-515/SLC7A11 axis. Our results partially proved the role of the circRHBG/miR-515/SLC7A11 axis in PCOS, which may help in the clinical diagnosis and treatment of PCOS. We present the following article in accordance with the MDAR reporting checklist (available at https://dx.doi.org/10.21037/atm-21-4174).

Methods

Clinical sample data and ethics approval

We recruited 12 infertility patients (6 PCOS and 6 controls) aged 20 and 38 who underwent in vitro fertilization (IVF) at the Department of Reproductive Medicine Center, the First Affiliated Hospital of Sun Yat-Sen University. We adopted the Rotterdam definition for PCOS, which involves diagnosis with 2 of 3 clinical manifestations, including clinical and/or biochemical hyperandrogenism, ovulation dysfunction, and PCOM after excluding patients with Cushing's syndrome, congenital adrenal hyperplasia, and androgen secreting tumors. Six women without PCOS at the similar age as that of PCOS patients were gathered as the normal control group. The patients with body mass index (BMI) >28 kg/m2, basal FSH level >10 mIU/mL, endometriosis, diminished ovarian reserve, chromosome abnormality, hydrosalpinx, systemic illness or endocrine disorders were excluded from this study. All patients provided signed written informed consent, and all experiments were performed under the approval of the local ethics committee of the First Affiliated Hospital of Sun Yat-Sen University ([2021]399-1). All procedures performed in this study involving human participants were in accordance with the Declaration of Helsinki (as revised in 2013).

All patients were treated the antagonist stimulation protocol. Briefly, on cycle day 2, ovarian stimulation was started by daily injection of recombinant FSH (r-FSH) (Gonal-F; Merck Serono). Gonadotropin-releasing hormone (GnRH) antagonist (Cetrotide, 0.25 mg; Merck Serono) injection was started from the day 6 of stimulation. When at least three follicles had reached 17 mm or two follicles had reached 18 mm in diameter, an intramuscular injection of 5,000 IU-10000 IU of human chorionic gonadotropin (hCG) is used for triggering final oocyte maturation. Oocyte retrieval was scheduled approximately 34–36 h after hCG injection by transvaginal ultrasound-guided follicular aspiration. Follicle aspiration of mature oocytes was used to obtain ovarian granulosa cells from all 6 samples. Clear follicular fluid containing granulosa cells was collected and centrifuged at 626 ×g for 10 min. The cells were resuspended in hyaluronidase and digested at 37 °C for 20 min. Finally, the lymphocyte separation solution was added and centrifuged at 626 ×g for 10 min. The white turbid clusters of cells in the middle were granular cells.

circRNA microarray

CapitalBio Technology Human CircRNA Array v2 was designed with 4 identical arrays per slide, with each array containing probes analyzing about 170,340 human circRNAs.

Cell culture and transfection

KGN and SVOG cells (GC cell line) were purchased from American Type Culture Collection (ATCC, Manassas, VA, USA). All cells were cultured in Dulbecco’s modified Eagle’s medium containing 10% fetal bovine serum (Gibco, Grand Island, NY, USA, Cat.# 16140071) and 1% gentamicin at 37 °C in 5% CO2. Lipofectamine 3000 reagent (Invitrogen, Waltham, USA) was used to transfect all plasmids, siRNA, and negative controls purchased from GeneChem (Shanghai, China), according to the manufacturer's instructions.

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis

TRIzol reagent (Life Technologies, Carlsbad, CA, USA) was used to isolate and extract total RNA from cells. SYBR Premix Ex Taq II kit (Takara Bio, Dalian, China, Cat.# RR036A) was then used to synthesize cDNA. SYBR Green PCR Master Mix (Takara Bio, Dalian, China, Cat.# RR420A) was used to amplify the target RNA on the CFX96 Real-time PCR System (Bio-Rad, Hercules, CA, USA), and β-actin was selected as the internal control. The expression change of the target gene was calculated using the 2-ΔΔCt method.

Luciferase reporter assay

The WT and MUT luciferase reporter vectors of circRHBG and SLC7A11 3'-untranslated region (UTR) were constructed. After co-transfection of miR-155-5p mimics or miR-155-5p locked nucleic acid (LNA) and their negative control, the dual-luciferase reporter gene detection kit (Promega, Madison, WI, USA) was used to evaluate the luciferase reporter gene activity.

Cell Counting Kit-8 (CCK-8) assay

The CCK-8 assay was used to assess the proliferation ability of cells. CCK-8 solution (10 µL) was added to each well 48 h after transfection. After incubating for 2 h at 37 °C, a microtiter plate reader (Epoch2; BioTek, Winooski, VT, USA) was used to measure the absorbance at 450 nM.

EdU assay

The Cell-Light EdU DNA Cell Proliferation Kit (RiboBio,Guangzhou, China) was used to test cell proliferation. Two days after transfection of KGN and SVOG cells with si-NC and si-RHBG, EdU was added for 48 h. The cells were then fixed with 4% paraformaldehyde and stained with Apollo dye solution (Beyotime, Shanghai, China). Nucleic acids were stained by DAPI. Fluorescence microscope and Image-Pro Plus software were used for the images and to evaluate cellular EdU incorporation.

Western blot

The total protein extracted from the cells was separated by 10% SDS-PAGE gel and transferred to a PVDF membrane. The protein was then placed in 5% skim milk and blocked at room temperature for 1 h. Anti-SLC7A11 (1:1,000; Abcam, Cambridge, MA, USA) or anti-GAPDH (1:1,000; Abcam, Cambridge, MA, USA) primary antibody was then added and incubated at 4 °C overnight. After washing 3 times with TBST, the protein was incubated with secondary antibody for 2 h at room temperature and quantified by chemiluminescence detection.

Detection of glutathione (GSH)/oxidized glutathione (GSSG) ratio

The GSH/GSSG Quantitative Kit II (Dojindo, Shanghai, China) was used to measure the GSH/GSSG ratio. The sample and standard solution were prepared, and concentration detected, according to the kit instructions. Finally, a microplate reader was used to calculate the relative level.

Statistical analysis

All data were expressed as the mean of 3 experiments (biological replicates) ± standard deviation (SD). Student’s t-test and 1-way analysis of variance were used to compare quantitative data. All statistical analyses were performed by using SPSS version 25.0 (IBM, SPSS, Chicago, IL, USA). A 2-sided P value <0.05 was considered statistically significant.

Results

circRHBG upregulation in PCOS

We used circRNA high-throughput chips to screen the differential expression profiles of related circRNA in 6 pairs of granular cells in the PCOS patients and control patients. According to the multiple of the difference greater than 1-fold and the P<0.05 criterion, it was found that there were 575 PCOS-related circRNAs, among which 362 were upregulated circRNA and 213 were downregulated circRNA (). Five upregulated circRNAs were then selected for qRT-PCR verification in 9 pairs of PCOS patients and non-PCOS patients. The results showed that circRNA hsa_circ_0014734 (circRHBG) derived from RHBG was upregulated in 9 cases of granular cells from PCOS patients (9/9: P<0.05), which was significantly higher than that of control group. This was consistent with the microarray results ().

Figure 1

circRHBG upregulation in PCOS. (A) Circular RNA (circRNA) chip detects the differential expression of PCOS-related circRNA in PCOS patients and non-PCOS controls. (B) Volcano plot using fold change >1 and P<0.05 to screen for the differential expression of circRNA. (C) Five upregulated circRNAs were verified by qRT-PCR in 9 pairs of PCOS patients’ and non-PCOS women’s granular cells. PCOS, polycystic ovary syndrome.

circRHBG knockdown inhibits KGN and SVOG cell proliferation

In order to explore the biological function of circRHBG in PCOS, we used CCK-8 assay and the EdU method to compare the viability of KGN and SVOG cells with different circRHBG expression levels. After transfection of the 2 siRNAs, the circRHBG expression level in KGN and SVOG cells decreased significantly, and the knockdown efficiency of si-circRHGB1 was higher (). The CCK-8 assay results showed that the cell viability of si-circRHBG-transfected cells was significantly lower than that of si-NC transfected cells, and the cell viability of si-circRHGB1-transfected cells was lower than that of si-circRHGB2-transfected cells (). We used si-circRHGB1 with higher knockdown efficiency for the EdU experiments. The EdU test results showed that the EdU incorporation rate of KGN and SVOG cells transfected with si-circRHGB1 was significantly lower than that of cells transfected with si-NC (). Briefly, after knocking down the expression of circRHBG, the proliferation ability of KGN and SVOG cells was significantly weakened.

Figure 2

circRHBG knockdown inhibits KGN and SVOG cell proliferation. (A) Expression level of circRHBG in KGN and SVOG cells after transfection of 2 siRNAs. (B) Cell Counting Kit-8 assay was used to evaluate the proliferation ability of KGN and SVOG cells. (C) EdU images of KGN and SVOG cells transfected with si-circRHBG1 and si-NC. **P<0.01.

circRHBG functions as a sponge for miR-515-5p

circRNA often regulates gene expression through an endogenous competitive RNA mechanism. In order to explore the downstream target molecules of circRHBG, we used circRNA interactome (https://circinteractome.nia.nih.gov) to predict miRNAs that may bind to it. The results showed that there is a sequence that binds to miR-515-5p on circRHBG (). To further verify the binding of these 2 molecules, we performed the dual-luciferase reporter assay in KGN and SVOG cells. The results showed that the fluorescence intensity of cells transfected with circRHGB-wt + miR-NC mimics was higher than that of cells transfected with circRHGB-wt + miR-515-5p mimics. The fluorescence intensity of cells transfected with circRHGB-mut + miR-NC mimics was not different from that of cells transfected with circRHGB-mut + miR-515-5p mimics (). To further verify the binding of circRHBG to miR-515-5p, we performed MS2bp-MS2bs-based RIP assay in KGN and SVOG cells. The results showed that miR-515-5p was enriched in the MS2bs-circRHBG group (). qRT-PCR showed that, in KGN and SVOG cells, the expression level of miR-515-5p in cells transfected with si-circRHBG was significantly higher than that in cells transfected with si-NC (). These results indicate that circRHBG acts as a miR-515-5p sponge in PCOS.

Figure 3

circRHBG functions as a sponge for miR-515-5p. (A) Prediction of the binding site of miR-515-5p on circRHBG. (B) Luciferase activity of KGN and SVOG cells transfected with circRHBG-wt + miR-NC mimics, circRHBG-wt + miR-515-5p mimics, circRHBG-mut + miR-NC mimics, and circRHBG-mut + iR-515-5p mimics. (C) MS2-based RNA immunoprecipitation assay in KGN and SVOG cells transfected with MS2bs-circRHBG, MS2bs-circRHBGmt, or MS2bs-Rluc. (D) Expression of mir-515-5p in KGN and SVOG cells transfected with si-circRHBG and si-NC, respectively. **P<0.01.

circRHBG competes with SLC7A11 to bind to miR-515-5p

miRNA usually inhibits the expression of the target gene by binding to the 3'-UTR of the targeted mRNA to degrade it (23). To explore the target molecules of miR-515-5p, we used the TargetScan algorithm to identify SLC7A11 as a candidate target oncogene (). To further verify the combination of these 2 molecules, we used the dual-luciferase reporter assay in KGN cells. The results showed that the fluorescence intensity of cells transfected with SLC7A11-wt + miR-NC mimic was higher than that of cells transfected with SLC7A11-wt + miR-515-5p mimics. The fluorescence intensity of cells transfected with SLC7A11-mut + miR-515-5p mimics was indifferent from that of cells transfected with SLC7A11-mut + miR-515-5p mimics (). qRT-PCR showed that, in KGN and SVOG cells, the expression level of SLC7A11 in miR-NC mimic transfected cells was significantly higher than that in miR-515-5p mimic transfected cells (). To further verify the binding of SLC7A11 to miR-515-5p, we performed Ago2-RIP detection in KGN cells. The results showed that circRHBG, SLC7A11, and miR-515-5p were mainly enriched in the Ago2 complex (). In addition, knocking down the expression of circRHBG significantly reduced the enrichment of circRHBG in the Ago2 complex, while increased the enrichment of SLC7A11 in the Ago2 complex (), indicating that circRHBG competes with SLC7A11 to bind to miR-515-5p.

Figure 4

circRHBG competes with SLC7A11 to bind to miR-515-5p. (A) Prediction of the binding site of SLC7A11 on miR-515-5p. (B) Luciferase activity of KGN cells transfected with SLC7A11-wt + miR-NC mimics, SLC7A11-wt + miR-515-5p mimics, SLC7A11-mut + miR-NC mimics, and SLC7A11-mut + miR-515-5p mimics. (C) Expression of SLC7A11 in KGN and SVOG cells transfected with miR-NC and miR-515-5p mimics, respectively. (D) RNA immunoprecipitation (RIP) assay detected the relative enrichment of circRHBG, SLC7A11, and miR-515-5p in the anti-Argonaute 2 (Ago2) complex (left); relative enrichment was detected by an Ago2-RIP assay (right). **P<0.01.

Knockdown of circRHBG promotes ferroptosis in PCOS

Many studies have confirmed that the overexpression of SLC7A11 promotes tumor growth by inhibiting iron death (24). qRT-PCR showed that, in KGN and SVOG cells, miR-515-5p mimics significantly increased the expression of miR-515-5p, and PBS-formulated LNA-modified oligonucleotide (LNA-anti-miR) effectively antagonized miR-515-5p in these 2 cells (). CCK-8 assay results showed that, in KGN and SVOG cells, cells with a high expression of miR-515-5p (transfected with miR-515-5p mimics) showed higher cell viability than those with a low expression of miR- 515-5p cells (transfected with LNA-anti-miRs) (). Western blot results also confirmed that interference with circRHBG expression in KGN and SVOG can significantly upregulate the expression of SLC7A11, and simultaneous transfection of si-circRHBG and miR-515-5p LNA can reverse the changes in SLC7A11 expression (). In addition, we found that circRHBG knockdown leads to a decrease in the GSH/GSSG ratio (), and LNA-mediated miR-515-5p silencing can reverse the change in the GSH/GSSG ratio (). The decrease in the GSH/GSSG ratio can lead to the inactivation of GPX4, a key molecule that regulates the ferroptosis pathway, and reduce lipid peroxides (25). Therefore, it can be concluded that circRHBG knockdown promotes ferroptosis in PCOS cells through the circRHBG/miR-515-5p/SLC7A11 axis.

Figure 5

circRHBG knockdown promotes ferroptosis in PCOS. (A) Expression of mir-515-5p in KGN and SVOG cells transfected with miR-NC mimics, miR-515-5p mimics, miR-NC LNA, and miR-515-5p LNA, respectively. (B) Cell Counting Kit-8 assay was used to evaluate the proliferation ability of KGN and SVOG cells. (C) Western blot was performed to determine the expression levels of SLC7A11 and GAPDH. (D,E) GSH/GSSG Quantification Kit was used to determine the glutathione (GSH)/oxidized GSH (GSSG) ratio. **P<0.01. PCOS, polycystic ovary syndrome.

Discussion

Many studies have found that a variety of miRNA and long-chain non-coding RNA (lncRNA) can regulate the progress of PCOS. However, few studies have investigated the role of circRNAs in PCOS (23). So far, four studies have investigated the difference in circRNA between PCOS patients and controls, in which hsa_circ_0083952, hsa_circ_0082709, hsa_circ_0002425, and hsa_circ_0015168, etc. have been confirmed to be dysregulated in PCOS patients (22,26-28). In the present study, we determined the expression profile of circRNA in PCOS and found that the expression of circRHBG was upregulated in PCOS granulosa cells. A series of cell experiments showed that knockdown of circRHBG expression enhanced the proliferation ability of KGN and SVOG cells. These results indicate that circRHBG plays a vital role in the progression of PCOS and could be a biomarker and therapeutic target for PCOS diagnosis.

The mechanism of ceRNA has become a research hotspot in recent years. It refers to the phenomenon that lncRNA, circRNA, and even pseudogenes regulate the expression of miRNA target genes through competition for miRNA binding sites (24). circRNA has been shown to regulate gene expression in a variety of diseases through the ceRNA mechanism (11,25). In the present study, we tested circRNA expression profiles in PCOS granulosa cells and then conducted a bioinformatics analysis. The findings indicated that the research molecule, circRHBG, contained miR-515-5p miRNA response element. We then performed functional tests to confirm the direct relationship between circ0097009 and miR-515-5p. Previous studies have found that miR-515-5p plays a regulatory role in a variety of diseases. Zhang et al. confirmed that miR-515-5p exerts a tumor-suppressor effect by targeting TRIP13 in prostate cancer (29). In addition, miR-515p has been shown to be involved in the epithelial-mesenchymal transition of epithelial cells in the process of smoking-induced airway remodeling (30). Pardo et al. also discussed the role of miR-515-5p/MARK4 in cancer cell migration (31). However, the role of miR-515-5p in the progression of PCOS has not been well studied. In the present study, we found that both circRHBG and SLC7A11 can directly bind to miR-515-5p, which indicates that circRHBG may act as a miR-515-5p sponge to regulate the expression of SLC7A11 through the ceRNA mechanism. Our research has progressively confirmed through multiple experiments that circRHBG acts as the ceRNA of SLC7A11 in PCOS and as a sponge for miR-1261. As mentioned earlier, the bioinformatics analysis showed that both 3'-UTR and circRHBG of SLC7A11 contained miR-515-5p binding sites. Second, the dual-luciferase reporter experiment and MS2bp-MS2bs-based RIP detection confirmed this prediction. Third, knockdown of the circRHBG reduced the expression of SLC7A11. Finally, miR-515-5pLNA reversed the effect of circRHBG knockdown on SLC7A11. These findings indicate that circRHBG acts as a ceRNA to regulate the expression of SLC7A11 through sponge miR-515-5p.

Ferroptosis is a novel type of cell death discovered in recent years. The process of cell death is usually characterized by a large amount of iron accumulation and lipid peroxidation (32). Iron metabolism is involved in endocrine diseases, including PCOS (33,34), while its underlying mechanism remains unclear. SLC7A11, the cystine/glutamate antiporter, functions to transport cysteine for glutathione synthesis and antioxidation, which is an upstream regulator of the ferroptosis pathway and can inhibit ferroptosis when its expression is upregulated (35). Ferroptosis is a recently discovered form of regulated cell death (36). It is characterized by the accumulation of lipid peroxidation products and reactive oxygen species (ROS) produced by iron metabolism, which can be pharmacologically inhibited by iron chelating agents and lipid peroxidation inhibitors (37). Our results showed that circRHBG knockdown reduced the expression level of SLC7A11, indicating that circRHBG is involved in the progression of PCOS and may play a role in regulating SLC7A11 to affect the ferroptosis pathway. In addition, we also found that circRHBG knockdown resulted in a decrease in the GSH/GSSG ratio, which leads to the inactivation of GPX4, a key molecule that regulates the ferroptosis pathway, and reduces lipid peroxides GPX4 (38). Studies have shown that silencing GPX4 makes tumor cells sensitive to ROS-induced cell death (39). Therefore, based on these findings, we conclude that the high expression of circRHBG inhibits iron death in PCOS through the Xc- and GPX4 pathways.

The findings of the present study confirmed the carcinogenic effect of circRHBG in PCOS, which plays a regulatory role in the circRHBG/miR515-5p/SLC7A11 axis. The establishment of the ceRNA network, circRHBG/miR515-5p/SLC7A11, in the PCOS model may provide new diagnostic molecular markers and therapeutic targets for PCOS.

The article’s supplementary files as