{"paper_id":"bcd8f1b5-190b-4c36-93e5-b758c334a5ba","body_text":"Polycystic ovary syndrome (PCOS) is a common endocrine disorder and a common cause of female infertility in women of reproductive age. \n 1 \n  One of the main characteristics of the syndrome is granulosa cell (GC) dysfunction, which largely contributes to hyperandrogenism, abnormal follicle development and anti‐Mullerian hormone excess. \n 2 \n  Aberrant gene expression profile was found in GCs of PCOS patients. \n 3 \n  Multiple differential expressed genes were high related to the pathogenesis of PCOS. However, the regulation mechanisms of these genes were largely unknown.\nForkhead Box O3 (FOXO3) plays important roles in diverse cellular processes including apoptosis, metabolism, cell proliferation and cell survival. \n 4 \n  FOXO3 is regulated at several mechanistic levels, such as transcriptional activity, cellular localization, mRNA expression and protein stability. Oxidative stress induces FOXO activation and nuclear translocation by c‐Jun N‐terminal kinase (JNK) or mammalian Ste20‐like kinase 1 (MST1) activation despite phosphorylation by protein kinase B (Akt). \n 5 \n ,  \n 6 \n  After energy deprivation, the increased AMP/ATP ratio leads to AMP‐activated protein kinase (AMPK) activation. \n 7 \n  AMPK activates FOXO3 activity by phosphorylation at six different residues. \n 8 \n  Insulin‐like growth factor‐1 (IGF‐I)/insulin and phosphoinositide 3‐kinase (PI3K)/Akt signalling pathway inactivate FOXOs by phosphorylation resulting in FOXOs nuclear exclusion. \n 9 \n  Zhao et al reported that elevated wnt family member 5A (WNT5a) activates PI3K/Akt signalling in GCs of PCOS patients. \n 10 \n  By contrast, Rice and colleagues showed metabolic insulin resistance in GCs of PCOS patients, suggesting that PI3K/Akt signalling is impaired. \n 11 \n  While many studies have been focused on the regulation of FOXO3 activity by post‐translational modifications, the regulation of  FOXO3  expression is largely unknown.  FOXO3  is one of the differential expressed genes which have higher expression in PCOS. \n 3 \n  Its overexpression in GCs is associated with higher apoptosis in PCOS. \n 12 \n  To date, the factors that up‐regulate and activate FOXO3 in GCs of PCOS patients are still unclear.\nm 6 A is the most prevalent modification of mRNA in higher eukaryotes. \n 13 \n  The modifications are reversible and dynamically regulated by the m 6 A modulators. m 6 A modulators consist of the ‘writers’, the ‘erasers’ and the ‘readers’. Briefly, m 6 A modifications are installed by the ‘writers’ (Methyltransferase like 3 (METTL3) and Methyltransferase like 14 (METTL14)), removed by the ‘erasers’ (Fat mass and obesity‐associated protein (FTO) and alkB homolog 5(ALKBH5)), and recognized by the ‘readers’ (YTH domain‐containing proteins and Eukaryotic initiation factor 3 (eIF3)). \n 14 \n ,  \n 15 \n ,  \n 16 \n ,  \n 17 \n ,  \n 18 \n ,  \n 19 \n ,  \n 20 \n ,  \n 21 \n ,  \n 22 \n  m 6 A modification has diverse biological functions such as nuclear RNA export, RNA splicing, protein translation regulation and RNA decay. In the YTHDF2‐mediated decay pathway, mRNAs with increased m 6 A abundance in 3’‐UTR are down‐regulated due to reduced RNA stability. \n 19\nSeveral human diseases are associated with altered m 6 A modification. m 6 A levels of sperm RNA are elevated in patients with asthenozoospermina. \n 23 \n  Increased m 6 A levels in the granulosa cells were reported in patients with premature ovarian insufficiency (POI), affecting apoptosis and cell proliferation in GCs. \n 24 \n  Aberrant m 6 A modification, through the effects on RNA metabolism, plays critical roles in a variety of cancers. \n 25 \n  However, whether m 6 A modification plays a role in the pathogenesis of PCOS is unknown.\nGiven that m 6 A modification affects the GC function in patients with POI, \n 24 \n  we suggested alteration of the m 6 A profile in the luteinized GCs of PCOS patients, which may account for the dysregulation of certain key genes for PCOS. Here, we showed the differences of m 6 A distribution between the luteinized GCs of normovulatory women and PCOS patients following controlled ovarian hyperstimulation. We identified  FOXO3  mRNA with differential m 6 A peaks, targeted for decay by YTHDF2. Our results indicated that hypomethylated  FOXO3  mRNA caused the dysregulation of  FOXO3  in luteinized GCs from PCOS patients following controlled ovarian hyperstimulation.\n\nForty‐three control patients with tubal factor infertility or male infertility and 36 PCOS patients were recruited in the Reproductive Medicine Center at the First Affiliated Hospital of Wenzhou Medical University between December 2017 and August 2019. Inclusion criteria for control patients were aged between 25 and 35 years, having a regular menstrual cycle, serum testosterone (T) level < 2 nmol/L, 18.5 < body mass index (BMI) < 27, having a normal ovarian reserve and a normal uterus. PCOS was diagnosed according to the Rotterdam revised criteria. \n 26 \n  The inclusion criteria for PCOS patients were aged between 25 and 35 years, having oligo‐ or anovulation, serum testosterone level > 2 nmol/L, 18.5 < BMI < 27, anti‐Müllerian hormone (AMH) > 7 ng/mL and having a normal uterus. The exclusion criteria were smoking, systemic diseases, endometriosis, abnormal serum level of prolactin, dysfunctional thyroid or having previous long‐term medication use.\nThis study was approved by the Ethics Committees of the First Affiliated Hospital of Wenzhou Medical (Approved number: YS2016‐063 and YS2019‐046). Written informed consent was obtained from all the participants.\nPlasma levels of AMH, luteinizing hormone (LH), follicle‐stimulating hormone (FSH), T and estradiol (E2) were collected and measured by chemiluminescence immunoassay (CLIA) between day 3 and day 5 of the menstrual cycle before the controlled ovarian stimulation. Serum fasting glucose levels and fasting insulin levels were measured by an oxidase‐peroxidase method and a CLIA method, respectively. The inter‐assay coefficients of variation were 6.1% for AMH, 7.2% for LH, 5.3% for FSH, 10.0% for T, 8.9% for E2, 4.3% for glucose and 5.4% for insulin. The intra‐assay coefficients of variation were 3.8% for AMH, 5.3% for LH, 4.6% for FSH, 8.1% for T, 6.8% for E2, 2.1% for glucose and 2.9% for insulin. Insulin in follicle fluids was detected by enzyme‐linked immunosorbent assay (ELISA, R&D Systems). HOMA‐IR was calculated by the formula (HOMA‐IR = fasting insulin (mIU/L) × fasting glucose (mmol/L)/ 22.5). \n 27\nUltrasound examination was performed on the 3rd‐5th day of the menstrual cycle with a 7 MHz transvaginal transducer (LOGIC 400, General Electric Medical Systems) to calculate follicle number. The basal antral follicle count (AFC) was assessed as the sum of all follicles of 2‐10 mm in diameter. The polycystic ovary was defined as 12 or more AFC in each ovary.\nAll participants underwent controlled ovarian stimulation with gonadotrophin‐releasing hormone (GnRH) agonist long protocol. Briefly, patients received 0.1 mg/day of GnRH agonist (Decapeptyl, Ferring, Germany) from day 20 of a spontaneous menstrual cycle until the day of human chorionic gonadotrophin (hCG) injection. When pituitary down‐regulation was achieved (usually after 14 days of GnRH agonist injection), ovarian stimulation was initiated with a minimum of 150 IU/day of recombinant human FSH (rhFSH, GonalF, Merck Serono). The dosages of rhFSH were adjusted according to serum E2 levels and follicle growth of the patients. 250 µg of recombinant hCG (Ovidrel, Merck Serono) was administered, when at least three follicles were >17 mm in diameter. Transvaginal aspiration was performed 34‐36 hours later to retrieve oocytes and follicular fluids. The follicular fluids from follicles >14 mm in diameter and without obvious blood contamination were collected. The follicular fluids were centrifuged at 340 g for 8 minutes to pellet the GCs. The cell pellets were resuspended in 1× PBS solution (Gibco), overlaid on 40%/80% gradient solution (PureCeption, SAGE) and centrifuged at 320 g for 20 minutes. GCs in the interface were collected and washed with 1× PBS solution.\nTotal RNA was isolated using the Total RNA Kit II Kit (Omega Bio‐Tek). Polyadenylated RNA was purified from total RNA using the GenElute mRNA Miniprep Kit (Sigma). 500 ng total RNA or polyadenylated RNA was used to determine the RNA m 6 A methylation levels using the EpiQuik m 6 A RNA Methylation Quantification Kit (Epigentek). Briefly, a standard curve was prepared according to the manufacturer's instruction. The RNA samples were coated on the strip wells, followed by incubation with capture antibody. After washes with washing buffer, detection antibody and enhance solution were added separately. Then, the signals were developed by the detection solution. The RNA m 6 A levels were quantified by the absorbance at 450 nm, and the m 6 A contents were calculated based on the standard curve.\nMeRIP was performed using Magna MeRIP m 6 A Kit (Millipore, Billerica, MA, USA) following the manufacturer's instruction. Briefly, 50 µg of total RNA was precipitated using 2.5 volumes of 100% ethanol, one‐tenth volumes of 3 mol/L sodium acetate and 1 mg/mL glycogen. Then, the RNA was resuspended, fragmented and immunoprecipitated with anti‐m 6 A antibody or normal mouse IgG. 10% of fragmented RNA was keep as input. The immunoprecipitated RNA was eluted by competition with m 6 A sodium salt and recovered using miRNeasy Mini Kit (QIAGEN, Hilden, Germany).\nThe recovered RNA was treated with DNase I (0.1 U/µL, Thermo Scientific, Waltham, MA, USA) and reverse‐transcribed by random hexamer primer and RevertAid Reverse Transcriptase (Thermo Scientific). Subsequently, second strand cDNA was synthesized using DNA Polymerase I and RNase H (Thermo Scientific). Sequencing library was constructed using NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB) and purified using AMPure XP system (Beckman Coulter). Library quality was assessed on a Bioanalyzer 2100 (Agilent). The clustering of the index‐coded samples was generated on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3‐cBot‐HS (Illumina). Library preparations were sequenced on a Hiseq 4000 Platform (Illumina).\nReal‐time quantitative RT‐PCR analysis (qRT‐PCR) was carried out with the samples with anti‐m 6 A antibody, the samples with normal IgG, and input samples in triplicates. Ct of the samples with anti‐m 6 A antibody were normalized to input by subtracting the Ct of input from the Ct of IP samples: ∆Ct = Ct IP ‐(Ct input ‐Log2[Input Dlilution Factor]). Then, the per cent of input for each IP sample was calculated: % Input = 2 −ΔCt(normalized IP) . The primers for MeRIP‐qPCR were listed in Table  S1 .\nSequencing data were mapped to the reference genome ( ftp://ftp.ensembl.org/pub/release‐94/fasta/homo_sapiens/dna/ ) using BWA mem (version 0.7.12). m 6 A peaks were identified by peak finding algorithm in MACS2 (version 2.1.0). The threshold of enrichment was set at  q  < 0.05. Differential peak analysis was based on the fold enrichment of peaks between the PCOS group and the control group. When the odds ratio > 2, a differential peak was determined. The statistics of pathway enrichment of genes with differential peaks in KEGG pathway was tested using KOBAS software.\nThe purified luteinized GCs were cultured in high‐glucose Dulbecco's modified Eagle's medium (DMEM, Gibco) containing 10% foetal bovine serum (FBS, Gibco) and 1% penicillin‐streptomycin (Gibco) at 37°C in a humidified incubator under 5% CO 2 . The siRNAs against YTHDF2, METTL3, METTL14, FTO and ALKBH5 were obtained from GenePharma Corporation (China). When the confluence reached 80%, the cells were dissociated and seeded in 24‐well plates (10 5  cells/well). One day later, the cells were transfected with siRNA (GenePharma) targeting negative control, YTHDF2, METTL3, METTL14, FTO or ALKBH5 using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instruction. The sequences of siRNA were listed in Table  S2 .\nFor the GCs collected from follicular fluids, total RNA was extracted using the Total RNA Kit II Kit (Omega Bio‐Tek). For the GCs cultured in 24‐well plates, total RNA in each well was extracted using MicroElute Total RNA Kit (Omega Bio‐Tek). The total RNA was reverse‐transcribed by the PrimeScript RT reagent Kit with gDNA Eraser (Takara) according to the manufacturer's instruction. Real‐time quantitative RT‐PCR analysis (qRT‐PCR) was performed on a Corbett Rotor Gene 6000 Real‐time Cycler using TB Green Premix Ex Taq II (Takara). Each sample was analysed in triplicates to obtain the threshold cycle number.  GAPDH  mRNA was used as an internal control for normalization. Relative expression levels were quantified using comparative 2 ‐∆∆Ct  method and expressed as fold changes relative to the controls. The sequences of the primers were listed in Table  S1 .\nThe total protein was extracted from the cultured GCs using RIPA buffer (Beyotime) containing 1mM PMSF (Beyotime). The concentration of the protein was measured using BCA protein assay kit (Solarbio). The cell extracts were separated on a 12% SDS‐PAGE Gel (ExpressPlus PAGE Gel, GenScript) and transferred to methanol‐activated polyvinylidenefluoride membranes. The membranes were blocked with 5% non‐fat milk powder for 1 hour and incubated with the primary antibody for 2 hour at room temperature. The mouse anti‐FOXO3 (Proteintech, 66428‐1‐Ig) or mouse anti‐Beta‐actin (Proteintech, 60008‐1‐Ig) antibody was used as primary antibody. Then, the membranes were incubated with HRP‐conjugated affinipure goat anti‐mouse IgG (H + L) (Biosharp, BL1001A) for 1 hour at room temperature. The signals were developed with SuperSignal West Pico PLUS Chemiluminescent Substrate (Thermo Scientific).\nThe proximal 3’‐UTR of  FOXO3  was amplified and inserted into 3’‐downstream of the  Renilla  luciferase gene of the psiCHECK TM ‐2 Vector between XhoI and NotI restriction sites. The motif of m 6 A binding was mutated by fusion PCR with a pair of complementary primers containing the A‐to‐T mutation. The primary culture human GCs were transfected with wild‐type m 6 A motif reporter plasmid or mutant m 6 A motif reporter plasmid using Lipofectamine 2000. The  Firefly  and  Renilla  luciferase activities were measured using the Dual‐Luciferase Reporter Assay System (Promega) 48 hours later. The  Renilla  luciferase activities were normalized to  Firefly  luciferase activities.\nData are expressed as the means ± SD. The normality and homogeneity of variance of the data were assessed. Differences between groups were determined by one‐way analysis of variance (ANOVA) followed by the Tukey multiple comparison test using SPSS 16.0 software package. Significance was set at  P  < 0.05.\n\nTo investigate the role of m 6 A modification in the luteinized GCs of PCOS patients, we recruited normovulatory controls and PCOS patients in our reproductive centre (Table  S3 ). The insulin‐glucose parameters indicated that the PCOS patients had insulin resistance in the present study. We examined m 6 A levels in the total RNA of the luteinized GCs collected following ovarian hyperstimulation by RNA m 6 A quantification. m 6 A levels were twofold higher in the luteinized GCs of PCOS patients compared with the controls (Figure  1A ). We further isolated polyadenylated RNA from the total RNA and examined the m 6 A levels. Consistently, the m 6 A levels of polyadenylated RNA were increased in the luteinized GCs of PCOS patients (Figure  1B ).\nQuantification analysis of m 6 A modification in the luteinized GCs of the controls and PCOS patients. The m 6 A content in total RNA (A) and mRNA (B) of GCs were detected by colorimetric assay. Bars represent means ± SD, n = 5. * P  < 0.05 vs the control\nWe conducted MeRIP‐seq to analyse the transcriptome‐wide distribution of m 6 A modification in GCs. 2764 and 3405 m 6 A peaks were identified in the controls and PCOS patients, respectively (Table  S4 ). We selected 8 transcripts with m 6 A peaks to verify the sequencing data by MeRIP‐qPCR. The qRT‐PCR results were in agreement with the sequencing data (Figure  S1 ). We analysed the differential peaks based on the fold enrichment of peaks of the two groups. 1719 and 2195 peaks were distinct in the controls and PCOS patients, respectively, while 996 peaks were overlapped in both groups (Figure  2A  and Table  S5 ).\nDistribution of m 6 A peaks in the luteinized GCs of the controls and PCOS patients. A, Comparison of the m 6 A peaks between the controls and PCOS patients. (B and C) Metagene profiles of enrichment of m 6 A peaks across mRNA transcriptome of the controls (B) and PCOS patients (C). (D and E) Distribution and enrichment of the m 6 A peaks within different gene regions in the GCs of the controls (D) and PCOS patients (E). TSS, transcription start sites; 5’‐UTR, 5’‐untranslated region; CDS, coding sequence; Stop, translation termination sites; 3’‐UTR, 3’‐untranslated region\nSimilar to the other studies, m 6 A peaks were strongly enriched around the stop codon in the controls (Figure  2B,D ). \n 28 \n ,  \n 29 \n ,  \n 30 \n  However, m 6 A peaks in PCOS patients showed a less prominent enrichment around the stop codon and increased locations to the CDS (coding sequence) and the TSS (transcription start sites) regions (Figure  2C,E ).\nWe suggested that m 6 A modification plays a role in dysregulation of PCOS‐associated genes. We examined the expression of  AMH ,  AR ,  FOXO1  and  FOXO3  in our study. Consistent with the previous studies, \n 2 \n ,  \n 3 \n ,  \n 31 \n  the expression of  AMH ,  AR ,  FOXO1  and  FOXO3  were up‐regulated in the luteinized GCs of PCOS patients compared with the controls (Figure  3A ). To screen potential candidates, we compared the expression of the m 6 A methyltransferases and demethylases in the luteinized GCs of the controls and PCOS patients. The expression of  METTL3 ,  METTL14 ,  FTO  and  ALKBH5  were all elevated in PCOS patients compared with the controls (Figure  3B ). We selectively knocked down the expression of  METTL3 ,  METTL14 ,  FTO  or  ALKBH5  in the luteinized GCs of the controls and PCOS patients, respectively. The knockdown efficiency was confirmed by qRT‐PCR analysis (Figure  S2 ). Due to differential m 6 A peaks identified in  FOXO1  and  FOXO3  transcripts, we examined the expression of  FOXO1  and  FOXO3  in these cells. In the controls, depletion of METTL3 or METTL14 increased the expression of  FOXO3 , and silencing of FTO suppressed the expression of  FOXO3 , suggesting m 6 A modification regulated  FOXO3  mRNA degradation (Figure  3C,E ,G). However, depletion of the m 6 A methyltransferases or demethylases did not affect the expression of  FOXO3  in PCOS patients (Figure  3D,F,H ), which indicated that effects of m 6 A modification on  FOXO3  mRNA were cell‐specific. Our results also showed that knockdown of methyltransferases or demethylases did not affect the expression of  FOXO1  in GCs (Figure  3C‐J ). The common target gene  SBK1  was used as a positive control to confirm the results. Collectively, these results demonstrated that selectively knocking down m 6 A methyltransferases or demethylases did not alter the expression of  FOXO3  in the luteinized GCs of PCOS patients, but did so in the controls, suggesting an absence of m 6 A‐regulated transcription of  FOXO3  in the luteinized GCs of PCOS patients.\nInvolvement of m 6 A modification in regulation of  FOXO3  mRNA in the luteinized GCs of the controls but not PCOS patients. A, mRNA levels of PCOS‐associated genes ( AMH ,  AR ,  FOXO1  and  FOXO3 ) in the luteinized GCs from the controls and PCOS patients. B, Expression of the methyltransferases ( METTL3  and  METTL14 ) and demethylases ( FTO  and  ALKBH5 ) in the luteinized GCs of controls and PCOS patients. C‐J, The effects of selectively knockdown of the methyltransferases or demethylases on the expression of  FOXO1 ,  FOXO3  and  SBK1  in the luteinized GCs of the controls and PCOS patients. GCs of the controls and PCOS patients were collected, cultured and subjected to METTL3‐siRNA (C,D), METTL14‐siRNA (E,F), FTO‐siRNA (G,H) or ALKBH5‐siRNA (I,J), respectively. qRT‐PCR analysis was performed to examine the mRNA levels of  FOXO1 ,  FOXO3  and  SBK1 . Data were expressed as fold changes relative to the blank controls. Bars represent means ± SD, n = 5. * P  < 0.05 vs the blank controls\nOur MeRIP‐seq showed a differential m 6 A peak in 3’‐UTR and near the stop codon of the  FOXO3  transcript (Figure  4A ). We performed MeRIP‐qPCR to confirm the differential m 6 A levels on the m 6 A target site (Figure  4B ). To determine how m 6 A modification regulates the  FOXO3  transcript, we analysed the sequence of the m 6 A peak. A putative m 6 A site was identified within the m 6 A peak of 3’‐UTR at position + 2279 (+1 relative to the translation start site) (Figure  4C ). As YTHDF2 is the main binding protein that account for the decay of m 6 A‐modified mRNAs, we next examined the effects of YTHDF2‐knockdown on  FOXO3  expression in human GCs. Depletion of YTHDF2 increased the amount of  FOXO3  transcript, and total protein levels of FOXO3 (Figure  4D,E ). The knockdown efficiency was checked by qRT‐PCR (Figure  S2 ). To assess the functionalities of the m 6 A site in the 3’‐UTR of  FOXO3  transcript, we constructed a reporter plasmid bearing FOXO3‐3’‐UTR with the putative m 6 A site mutated (Figure  4F ). FTO‐siRNA suppress the luciferase activities of the reporter plasmid with the wild‐type m 6 A site (Figure  4G ). However, mutation of the m 6 A site reversed the suppressive effects of FTO‐siRNA on the luciferase activities (Figure  4G ). The presence of the m 6 A readers in human GCs was showed by a previous study. \n 32 \n  To exclude the potential of lacking the m 6 A readers in GCs of PCOS patients, we examined the expression of the m 6 A readers in PCOS. The expression of  YTHDF2  was significantly increased in the luteinized GCs of PCOS patients (Figure  S3 ). These results indicated that m 6 A modification regulated  FOXO3  mRNA decay through the m 6 A site in 3’‐UTR in the luteinized GCs of the controls.\nDifferentially methylated  FOXO3  transcript was regulated via the YTHDF2‐mediated decay pathway in the luteinized GCs of the controls. A, m 6 A peaks distribution across 3’‐UTR of the  FOXO3  transcript. The Y‐axis shows read number. Blue boxes indicate exons, and blue line indicates introns. B, qRT‐PCR analysis of m 6 A peak region and non‐peak region in the 3’‐UTR of the  FOXO3  transcript. Bars represent means ± SD, n = 5. * P  < 0.05 for the differences between the indicated groups. C, Schematic representation of the position of m 6 A motif in the 3’‐UTR of the  FOXO3  transcript. D, The mRNA level of  FOXO3  in the YTHDF2‐knockdown cells. Bars represent means ± SD, n = 5. * P  < 0.05 for the differences between the indicated groups. E, Western blot analysis of FOXO3 total protein level in the YTHDF2‐knockdown cells. F, Wild‐type or m 6 A motif mutant (A‐to‐T mutation) FOXO3‐3’‐UTR fused with  Renilla  luciferase reporter. G, Relative luciferase activity of the FOXO3‐3’‐UTR with wild‐type or mutant m 6 A motif after cotransfection with negative control siRNA or FTO‐siRNA  Renilla  luciferase activity was normalized to firefly luciferase activity. Bars represent means ± SD, n = 3. * P  < 0.05 vs the control siRNA group\n\nIncreasing evidence showed that m 6 A modification play important roles in various physiological processes. \n 17 \n ,  \n 33 \n ,  \n 34 \n ,  \n 35 \n ,  \n 36 \n  Elevated or decreased m 6 A level was associated with various human diseases. \n 23 \n ,  \n 24 \n ,  \n 25 \n  We showed that the level of m 6 A modification was higher in the luteinized GCs of PCOS patients than that in the controls. Aberrant m 6 A levels affect gene expression and biological functions. However, whether human diseases alter m 6 A profiles is largely unknown. Here, we showed that not only elevated m 6 A level but also altered m 6 A profile are associated with PCOS in the luteinized GCs.\nThe site selection of m 6 A modification is largely unknown. miRNAs regulate m 6 A formation at corresponding target sites. \n 37 \n  Altered miRNAs and other associated factors partially explain the differences of m 6 A profiles among cell types. \n 37 \n  Furthermore, m 6 A profile was altered when cells exposed to heat shock, ultraviolet radiation or signalling molecules. \n 38 \n  Acute stress regulates a fraction of genes with m 6 A modification in the cortex. \n 33 \n  The GCs of PCOS patients present different cell characteristics. Cell survival, proliferation rate and responsiveness to FSH were altered in the GCs of PCOS patients. \n 39 \n ,  \n 40 \n  Differentially expressed miRNAs were also identified in the GCs of PCOS patients. \n 41 \n ,  \n 42 \n  Moreover, abnormal serum AMH, androgen, activin A and follistatin were observed in PCOS patients. \n 2 \n ,  \n 43 \n ,  \n 44 \n  All these alterations may account for the altered m 6 A profile in PCOS. We identified 2195 differential m 6 A peaks in the luteinized GCs of PCOS patients compared to those of the controls. KEGG analysis showed that differential peaks were enriched in metabolic pathway (Figure  S4 ). Whether m 6 A modification is involved in the regulation of these transcripts with the differential m 6 A modification needs further investigation.\nInsulin resistance is prevalent in PCOS patients and plays a crucial role in metabolic abnormalities. Although obesity increases risks for insulin resistance, several studies have demonstrated that both obese and non‐obese PCOS patients had impaired insulin signalling. \n 45 \n ,  \n 46 \n  In the present study, the non‐obese PCOS patients had significantly higher insulin 0, insulin 120 and HOMA‐IR compared with the controls. Moreover, the insulin levels in follicular fluid were notably higher in the non‐obese PCOS patients. These indicated that the PCOS patients in our study had insulin resistance. As FOXO proteins mediate the effects of insulin signalling on metabolism, \n 47 \n  it is likely that altered expression of  FOXO3  is involved in the defects of insulin signalling in the luteinized GCs of PCOS patients. Here, we demonstrated an absence of m 6 A‐mediated  FOXO3  mRNA destabilization which altered expression of  FOXO3  in the luteinized GCs of PCOS patients. Our findings have offered a potential mechanism for insulin resistance.\nIt is well known that insulin activates PI3K/Akt signalling to induce the nuclear exclusion of FOXO3. \n 9 \n  However, the effects of PI3K/Akt signalling on GCs of PCOS patients are controversial. Zhao et al showed that WNT5a expression was elevated in GCs of PCOS patients. \n 10 \n  Although they have demonstrated that WNT5a increased inflammation and oxidative stress via PI3K/Akt/NF‐κB signalling in human GCs, they did not confirm that WNT5a activated PI3K/Akt/NF‐κB signalling in GCs of PCOS patients. As insulin resistance is attributed to defects in PI3K/Akt signalling, metabolic insulin resistance in GCs of PCOS patients suggests impaired PI3K/Akt signalling. \n 9 \n ,  \n 11 \n  Thus, it is possible that abnormalities in other signalling cascade rather than PI3K/Akt signalling contributed to chronic inflammation in GCs of PCOS patients. To date, the factors that up‐regulate and activate FOXO3 remain unclear in GCs of PCOS patients. Here, we demonstrated that YTHDF2 mediated  FOXO3  mRNA decay via a putative m 6 A site in the 3’‐UTR in the luteinized GCs of the controls. In contrast to the findings in the controls, although the expression of YTHDF2 was elevated in the luteinized GCs of PCOS patients,  FOXO3  mRNA was hypomethylated and irresponsive to m 6 A modification. Our results showed an absence of m 6 A‐based regulation of FOXO3 expression in the luteinized GCs of PCOS patients. Taken together, we inferred that altered m 6 A modification caused upregulation of FOXO3 in the luteinized GCs of PCOS patients following controlled ovarian hyperstimulation.\nIn summary, the present study provided m 6 A profiles of normovulatory women and PCOS patients in luteinized GCs following controlled ovarian hyperstimulation. We demonstrated that altered m 6 A modification disturbed the regulation of FOXO3 expression in the luteinized GCs of PCOS patients. However, the site selection mechanism of m 6 A in PCOS needs to be explored in future studies.\n\nThe authors confirm that there are no conflicts of interest.\n\nShen Zhang:  Conceptualization (lead); Data curation (lead); Funding acquisition (equal); Investigation (lead); Methodology (lead); Writing‐original draft (lead); Writing‐review & editing (lead).  Wenli Deng:  Conceptualization (supporting); Data curation (supporting); Investigation (supporting); Methodology (supporting); Writing‐original draft (supporting); Writing‐review & editing (supporting).  Qiongyou Liu:  Investigation (supporting); Methodology (supporting); Writing‐review & editing (supporting).  Peiyu Wang:  Data curation (supporting); Investigation (supporting); Methodology (supporting).  Wei Yang:  Investigation (supporting); Methodology (supporting); Writing‐review & editing (supporting).  Wuhua Ni:  Funding acquisition (equal); Project administration (lead).\n\nFig S1\nClick here for additional data file.\nFig S2\nClick here for additional data file.\nFig S3\nClick here for additional data file.\nFig S4\nClick here for additional data file.\nTable S1\nClick here for additional data file.\nTable S2\nClick here for additional data file.\nTable S3\nClick here for additional data file.\nTable S4\nClick here for additional data file.\nTable S5\nClick here for additional data file.","source_license":"CC-BY-4.0","license_restricted":false}