Abstract
Polycystic ovary syndrome (PCOS), characterized by polycystic ovarian morphology, ovarian follicular maturation arrest, and hormonal disorders, affects numerous women in the reproductive age worldwide. A recent study has found that mitochondria are likely to play an essential role in oocyte quality. However, it is still unclear whether oocyte development failure is associated with mitochondria in patients with PCOS. We analyzed the single-cell RNA sequencing data from the previous study, including data from 14 oocytes from 7 healthy fertile women and 20 oocytes from 9 patients with PCOS at the germinal vesicle (GV) stage, metaphase I (MI) stage, and metaphase II (MII) stage. We revealed the transcriptomic dynamics by weighted gene co-expression network analysis (WGCNA) and investigated the differences between stages using PCA and Deseq2 analyses to identify the differential expression genes (DEGs). Gene ontology (GO) was performed using clusterProfiler R package and Metascape. Our results indicated that specific gene modules were related to different stages of oocyte development using WGCNA. Functional enrichment analysis and gene co-expression network analysis found significant enrichment of the mitochondrial regulation genes at the GV stage. PCA (principal component analysis) and differential gene expression analysis suggested that GV was significantly different from the MI and MII stages between the two groups. Further analysis demonstrated that the upregulated differentially expressed genes at the GV stage of patients with PCOS mainly related to mitochondrial function, such as COX6B1, COX8A, COX4l1, and NDUFB9. Meanwhile, these genes tended to be activated at the MII stage in healthy cells, suggesting that some mitochondrial functions may be prematurely activated at the GV stage of PCOS oocytes, whereas this process occurs at the MII stage in healthy oocytes. Collectively, our study showed that aberrant mitochondrial function at the GV stage may contribute to a decline in oocyte quality of PCOS patients.
Introduction
Polycystic ovary syndrome (PCOS) is a common disorder affecting 5–20% of the women of reproductive age worldwide. PCOS is mainly characterized by hyperandrogenism, ovulatory dysfunction, and polycystic ovarian morphology (PCOM) (). The interaction between insulin resistance (), hyperinsulinemia (), and gonadotropin disturbance () can induce PCOS, but the pathogenesis is still unclear. Currently, the treatment strategies for PCOS mainly include medical treatment, ovarian surgery, and in vitro fertilization (IVF) therapy (; ; ; ; ). Women diagnosed with PCOS, who undergo IVF treatment, show a high probability for pregnancy and reach the mean level of the pregnancy rate (). However, several serious problems exist. Besides the high risk of ovarian hyperstimulation syndrome and preeclampsia, oocyte amounts and quality obtained from patients with PCOS undergoing IVF are worse than those from healthy fertile women (; ).
Mitochondria are the center of energy metabolism, and can both supply enough energy for cells and regulate cellular survival and development (; ; ). Accumulating evidence has shown that mitochondria play an essential role in oocyte maturation, influencing oocyte and subsequent embryonic development (; ). Mitochondrial dysfunction affects oocyte maturation, which can consequently decrease female fertility (). It can be inferred that a decline in oocyte quality in patients with PCOS may be related to abnormal mitochondrial function. However, the intracellular alterations and mitochondria-associated mechanism in oocytes of PCOS are still unclear owing to the limitations in patient material and technologies.
In this context, we analyzed the single-cell RNA-seq (scRNA-seq) data of healthy and PCOS oocytes from different stages of its maturation, including germinal vesicle (GV), metaphase I (MI), and metaphase II (MII). Combined with the bioinformatics algorithms, we determined the mitochondrial dynamics in PCOS oocytes and explored whether aberrant mitochondrial function is associated with the failure of oocyte development in patients with PCOS.
Materials and Methods
Data Retrieval
The scRNA-seq data for analyses were obtained from parts of a previous study published by our group, inclusion criteria and ethical statements have been previously described (). Briefly, ovarian stimulation and oocyte retrieval were performed according to the published protocols (). Each picked up oocyte was transferred to a tube containing lysis buffer and stored at −80°C for further analysis. We conducted single-cell RNA-sequencing followed by the published Smart-seq2 protocol () and sequence based on illumina HiSeq 2500.
For calculating the expression of genes, the adaptor contaminants and low-quality based on the scRNA-seq reads were trimmed by Trimmomatic (V0.33) (), then the clean reads were aligned to hg38 reference genome using STAR (version 2.7.1a) (). The number of mapped reads in the BAM file of each sample was counted, then the gene expression was defined using read count (Figure 1A). Additionally, normalized counts Fragments Per Kilobase Million (FPKM) was calculated based on the previous paper (). According to the protocol, the read count was used for DESeq2 to calculated the differentially expressed genes (DEGs), FPKM was used for WGCNA analysis (Figure 1B).
FIGURE 1
Weighted Gene Co-Expression Network Analysis (WGCNA)
To investigate the gene expression network during the maturation process of oocytes, we performed WGCNA (V1.68) () and determined the relationship of gene co-expression in different stages between healthy individuals and individuals with PCOS. The criteria of filtering were clustering height above 15000 or below 10. The power used in WGCNA was set according to the SoftThreshold, and we plotted module detection via dynamic tree cutting. For the demonstration of the relationship between different modules, results were visualized using module and eigengene relation heatmap, and the gene co-expression network was extracted and further processed using MCODE in Cytoscape software (V3.7.2) (). For a distinct and precise representation of the data, genes contained in the significantly relevant modules were extracted and analyzed using Metascape ().
Gene Clustering and Differential Expression Analysis
To explore the clusters in different stages, we used the principal component analysis (PCA) and visualized the dotplots by selecting the top two principal components. Differential gene expression analysis was performed using the DESeq2 package (V1.24.0) (). Briefly, read counts were input to build dds objects, and DEGs were extracted from the dds results. DEGs at different stages were illustrated using a volcano plot according to the EnhancedVolcano package (V1.2.0) (). Genes with a p-value < 0.05 and absolute fold change ≥2 were regarded as differentially expressed genes.
Gene Ontology (GO) and KEGG analysis
To detect the genome-wide expression profiles between healthy and PCOS oocyte, we performed gene set enrichment analysis (GSEA) for hallmark gene sets and GO gene sets separately (; ) using datasets of whole samples in these groups. Gene ontology enrichment analysis was conducted using the clusterProfiler package (V3.12.0) (). Before analysis, we transferred gene names from “symbol” to “entrezid” according to org.Hs.eg.db package (V3.8.2) (). We harvested genes using biological progress (BP) subontology. Terms with the p-value < 0.05 were considered and enriched as significant. Dotplot and gene-concept networks were illustrated using the enrichplot package (V1.4.0) (). Relative expression of selected genes involved in GO terms was calculated using samples filtered by clustering height above 15000 or below 10 and Boxplots were illustrated using the ggplot2 package (V3.2.1) (). Genes involved in particular GO sets were highlighted in the KEGG pathway using clusterProfiler package (V3.12.0) ().
Statistical Analysis
Statistical analyses in Figure 5 were performed by GraphPad Prism (V8.2.1). The FPKM data of selected mitochondria-related genes were used for testing by unpaired, two tailed Mann-Whitney U tests. p-value < 0.05 were regarded as statistically significant and marked with asterisk.
Results
WGCNA Reveals Different Genetic Dynamics of Maturation Between Healthy and PCOS Oocytes
The scRNA-seq data of 6 GV oocytes, 5 MI oocytes, 3 MII oocytes from 7 healthy fertile women and 7 GV oocytes, 4 MI oocytes, 9 MII oocytes from 9 PCOS patients were analyzed. The clinical variables of these 16 donors confirmed the phenotype of PCOS patients, including significant decreased of antral follicle, LH, LH/FSH, testosterone, sex hormone-binding globulin (SHBG) and ovarian volume (p-value < 0.05) (Table 1). For exploring the genetic program macroscopically, we profiled 63129 genes from these 35 samples using WGCNA and 16 gene co-expression modules were identified (Figure 2A). Particularly, we found that the blue module was the most significant and correlated with the GV stage, while the yellow module was the most significant and correlated with the MI stage, and magenta module was the most significant and correlated with the MII stage (Figure 2B). GO functional enrichment of hub genes in these three modules were analyzed using Metascape (). In the blue module, genes were mainly involved in the ncRNA metabolic process, mitochondrion organization and others, while chemotaxis, pid p53 downstream pathway and others were abundant in the yellow module, and reversible hydration of carbon dioxide, cell-cell adhesion and others were highly represented in the magenta module (Figures 2C–E, Supplementary Figure S1, and Supplementary Table S1).
TABLE 1
| Healthy donors (n = 7) | PCOS patients (n = 9) | p-value | |
| Age (years) | 29.7 ± 0.6 | 27.4 ± 1.1 | 0.1266 |
| BMI (Kg/m2) | 18.6 ± 0.4 | 19.7 ± 0.6 | 0.1427 |
| SBP (mmHg) | 119.2 ± 3.9 | 116.3 ± 4.4 | 0.6515 |
| DBP (mmHg) | 77.1 ± 2.5 | 76.4 ± 3.3 | 0.8739 |
| AFC (antral follicle) | 10.1 ± 1.4 | >24 | <0.001(0.000)** |
| FSH (mlU/ml) | 6.1 ± 0.5 | 5.7 ± 0.5 | 0.5516 |
| LH (mIU/ml) | 5.0 ± 0.6 | 8.7 ± 0.6 | <0.001(0.000)** |
| LH/FSH | 3.9 ± 0.3 | 6.2 ± 0.4 | <0.001(0.000)** |
| Insulin (IU) | 10.5 ± 1.9 | 12.7 ± 2.1 | 0.464 |
| Testosterone (nmol/L) | 1.1 ± 0.1 | 1.6 ± 0.2 | <0.05(0.0311)* |
| Androstenedione | 2.2 ± 0.3 | 2.3 ± 0.2 | 0.7234 |
| Sex Hormone-Binding Globulin (SHBG, nmol/L) | 38.0 ± 2.5 | 53.7 ± 5.4 | <0.05(0.0313)* |
| Ovarian volume total | 3809.6 × 408.0 | 5789.8*555.2 | <0.05(0.0165)* |
The clinical variables of healthy donors and PCSO patients.
FIGURE 2
To find the biological pathways which play a key role during the dynamic progression of oocyte development, the gene set enrichment analysis (GESA) was performed and oxidative phosphorylation was significantly enriched in PCOS group using both hallmark gene sets and GO gene sets, suggesting the significant difference in mitochondrial function between healthy and PCOS oocytes (Figure 2F). Furthermore, only hub genes in blue module which related to GV stage enriched with mitochondrial function terms (Figure 2C). Therefore, we further analyzed the genes in co-expression network of blue modules using Cytoscape MCODE algorithm. The hub genes with high MCODE score were related to many important functions during the meiosis and oocyte maturation, including regulation of microtubule and centrosome (RMDN3, CETN2) (; ), DNA repair and replication (POLE) (), oocyte maturation and cleavage (ANAPC7, CSTF1) (; ), nucleic acid binding (ZNF548, DDX19A) (; ). Notably, a subset of genes related to mitochondrial functions was identified, including MRPL22, AIFM1, ADA, ATP5L, NDUFA1, NDUFA2 which are broadly involved in nucleoside triphosphate metabolic processes, mitochondrial oxidation respiratory and energy generation (Figure 2G). Overall, these results suggested that mitochondrial function may have an essential impact in GV stage progression.
Significant Difference in Transcriptome Variations Between Healthy and PCOS Oocytes at the GV Stage
To further explore the discrepancy between healthy and PCOS oocytes in the GV, MI, and MII stages, we clustered genes using PCA. We identified that the clusters were markedly different in GV oocytes between the healthy and PCOS groups, while similar clusters in most of the oocytes at the MI and MII stages between the two groups were not observed (Figure 3A). We collected the highly-differential genes at diverse stages from healthy and PCOS oocytes and visualized the data by volcano plots (Figure 3B and Supplementary Table S2). By setting the standard of p-value < 0.05 and fold change ≥2, a total of 2426, 1298, and 1181 DEGs were identified in the GV, MI, and MII stages, respectively. Prominent DEGs were identified in each stage, for instance, the downregulation of CYP26A1, MTRNR2L1, ELOA and upregulation of FAM53A, PPP1R35, BLM at GV, MI, MII stages in PCOS oocytes, respectively (Figure 3B). It has been reported that CYP26A1, MTRNR2L1, ELOA, and FAM53A were all related with ovarian carcinomas, and CYP26A1 may acts as a meiosis-inhibiting factor and related with the gonadal retinoic acid-degradation (). In addition, DEGs from different stages showed more independence and less overlap between each other (Figure 3C and Supplementary Table S3). These results suggest that different development stages are relatively isolated on the transcriptome level, and variations in the GV stage are the most significant.
FIGURE 3
An Aberrant Status of Mitochondrial Energy Metabolism in PCOS Iocytes
To further investigate the mechanism of differences between healthy and PCOS oocytes at the GV stage, we performed GO analysis using DEGs identified by DEseq2 and found the most of significant enrichment functions of DEGs at the GV stage were mitochondria-related processes, such as ATP synthesis coupled electron transport, mitochondrial respiratory chain complex assembly and similar, however, these terms were not present at the MI and MII stages (Figure 4A). Using GO enrichment analysis, we found that upregulated genes in GV oocytes of patients with PCOS were highly related to mitochondrial energy metabolism (Figure 4B), while upregulated genes in the healthy samples were relevant to oocyte development (Figure 4C), suggesting the mitochondrial function abnormality is likely to be an important reason which lead to the difference between PCOS and healthy oocytes at GV stage.
FIGURE 4
To reveal the dynamics of mitochondria-related genes at different stages, we compared genes which related to mitochondrial function and were variably expressed at the GV stage between healthy and PCOS oocytes using boxplot (Figures 4A, 5). We found that genes such as COX6B1, COX8A, COX4l1, and NDUFB9 that were highly expressed at the GV stage in PCOS oocytes tended to be activated at the MII stage in healthy oocytes (Figure 5). To investigate the specific function of these genes in mitochondria, we further performed KEGG pathway analysis and found that these genes are involved in most parts of the oxidative phosphorylation process (Figure 6). Based on these findings, we suggest that some mitochondrial functions may be prematurely activated at the GV stage of PCOS oocytes, whereas this process occurs at the MII stage in healthy oocytes.
FIGURE 5
FIGURE 6
Discussion
Previous single-cell transcriptome study of PCOS oocyte reported by showed genes related to oocyte meiosis and maturation were differentially downregulated and DNA repair-related genes were upregulated at early stages and these DEGs were the potential reasons caused PCOS oocytes disorder. In the present study, we compared the transcriptomic profiles between healthy and PCOS oocytes and found one of the main differences during oocyte development was mitochondrial function, such as upregulated expression of oxidative phosphorylation related genes using GESA analysis. Through detected specific gene modules and dynamic characteristics of co-expression genes using WGCNA, the mitochondria-related genes were identified as hub genes in the gene module related to GV stage. These results were also confirmed using PCA and differential gene expression analysis, the transcriptomic profile was significantly different at GV stage suggesting that the DEGs at the GV stage are vital during the oocyte development process between women with PCOS and healthy women, and these DEGs are mainly related to mitochondrial functions. Additionally, our results also suggested that mitochondrial function, especially the oxidative phosphorylation process, may be prematurely activated at the GV stage in PCOS oocytes, which differ from the healthy oocytes which activated at MII stage. Overall, we give further insight into the mitochondrial dynamics during oocyte maturation and we find that mitochondrial aberrant activation at GV stage may be one of the important factors contributing to reduced oocyte quality in patients with PCOS.
Polycystic ovary syndrome is one of the most common and severe diseases affecting the female reproductive health, a decline in oocyte quality is more susceptive to occur in PCOS patients, and so far, the pathogenesis of PCOS is still unclear (). The assisted-reproductive technology (ART) presents a way of rescue for PCOS patients, present treatment for PCOS patients includes ovarian stimulation for receiving MII oocyte with subsequent fertilization in vitro (; ), nevertheless, the mature follicular obtained through ovarian stimulation will less than that in women without PCOS for unknown reasons (). The process from GV to MII stages is a key pathway for oocyte development. Although gene abundance have no statistical difference among diverse stages (), dynamic biological process such as differential expression of mitochondria-related genes can be observed. Moreover, the number of mitochondria increases to more than 100,000 in MII oocytes, which is a 1000-fold higher than that in the GV oocytes (; ). Thus, the state of mitochondrial function along with cell death and proliferation can also influence oocyte maturation. As the center of energy metabolism and cell apoptosis in most cell types (; ), mitochondria disorder leads to the failure of human oocyte development (; ).
In this study, we systematically analyzed the changes in the oocytes of patients with PCOS and healthy fertile women at different developmental stages using single-cell transcriptome sequencing. We found that some genes related to gonadal and oocyte development, such as CYP26A1 which is upregulated in oocytes during prophase-I meiotic arrest and downregulate when meiosis resumed (), was downregulated in PCOS oocytes at GV stage. Moreover, We revealed the transcriptomic dynamics during oocyte maturation and found that the mitochondrial function may be prematurely activated since the highly expressed genes at the GV stage of PCOS oocytes show a close relationship with mitochondrial function such as mitochondrial energy metabolism, oxidative phosphorylation process and similar, which differ from the healthy oocytes that these genes were activated at MII stage. It was found that the function of mitochondria is quiescent during follicular development, and energy is mainly provided by the surrounding granular cells (). The premature activation of mitochondria likely leads to the production of harmful metabolites, such as excessive oxygen free radicals, which could induce irreversible damage to the maturating oocyte further resulting in the oocyte retardation (; ).
However, it should be noted that there are some limitations in this study. Since the women with PCOS often have ovulation disorders which mainly cause infertility (), it’s difficult to obtain oocytes from them. Therefore, the number of oocyte samples used in this study is relatively smaller than other oocyte gene expression researches (). In addition, the differential expression of mitochondria-related genes presented in this study was quantified by sequencing and data analysis. Further independent validations such as quantitative RT-PCR and functional researches are needed to comprehensively understand the role of the mitochondrial function in the progress of oocyte maturation.
In conclusion, our findings highlight that the transcriptomic dynamics of PCOS oocytes and aberrant mitochondrial function at the GV stage of oocyte development in patients with PCOS were observed. These results deliver novel meaningful insights into the oocyte quality decline in patients with PCOS.
Statements
Data availability statement
The data is available to access in GEO SRA and the accession number is PRJNA600740.
Author contributions
LQ, BxL, XC, and JX analyzed the data and wrote the manuscript. LQ, BL, QL, and JX performed bioinformatics analysis. WL, BxL, QL, XX, and LW performed the experiments. LQ, ZX, JX, and QZ designed the experiments and revised the manuscript.
Funding
This work was supported by grants from the National Key R&D Program of China (2017YFC1001301 and 2016YFC1000208), the National Natural Science Foundation of China (81771651), the China Postdoctoral Science Foundation (2017M611623), and Huaihua Science and Technology Plan Project (2018N2208).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00396/full#supplementary-material
FIGURE S1GO analysis of hub genes in different modules. (A) GO terms of hub genes in blue module. (B) GO terms of hub genes in yellow modules. (C) GO terms of hub genes in magenta module.
TABLE S1Genes list for three significantly gene modules.
TABLE S2Differentially expressed genes of GV, MI and MII stages between PCOS and healthy group.
TABLE S3Overlap genes of venn diagram.
References
1
AhmadiF.RashidyZ.HaghighiH.AkhoondM.NiknejadiM.HematM.et al (2013). Uterine cavity assessment in infertile women: sensitivity and specificity of three-dimensional hysterosonography versus hysteroscopy.Iran. J. Reprod. Med.11977–982.
2
AzzizR.CarminaE.ChenZ.DunaifA.LavenJ. S.LegroR. S.et al (2016). Polycystic ovary syndrome.Nat. Rev. Dis. Primers2:16057.
3
BabayevE.SeliE. (2015). Oocyte mitochondrial function and reproduction.Curr. Opin. Obstet. Gynecol.27175–181. 10.1097/gco.0000000000000164
4
BligheK.RanaS.LewisM. (2019). EnhancedVolcano: Publication-Ready Volcano Plots with Enhanced Colouring and Labeling. R Package Version 1.4.0.
5
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data.Bioinformatics302114–2120. 10.1093/bioinformatics/btu170
6
ChanD. C. (2006). Mitochondria: dynamic organelles in disease, aging, and development.Cell1251241–1252. 10.1016/j.cell.2006.06.010
7
ChiarattiM. R.GarciaB. M.CarvalhoK. F.MachadoT. S.RibeiroF. K. D.MacabelliC. H. (2018). The role of mitochondria in the female germline: Implications to fertility and inheritance of mitochondrial diseases.Cell Biol. Int.42711–724. 10.1002/cbin.10947
8
ChildT. J.Abdul-JalilA. K.GulekliB.TanS. L. (2001). In vitro maturation and fertilization of oocytes from unstimulated normal ovaries, polycystic ovaries, and women with polycystic ovary syndrome.Fertil. Steril.76936–942. 10.1016/s0015-0282(01)02853-9
9
DarbandiS.DarbandiM.KhorshidH. R. K.SadeghiM. R.AgarwalA.SenguptaP.et al (2017). Ooplasmic transfer in human oocytes: efficacy and concerns in assisted reproduction.Reprod. Biol. Endocrinol.1577.
10
DobinA.DavisC. A.SchlesingerF.DrenkowJ.ZaleskiC.JhaS.et al (2013). STAR: ultrafast universal RNA-seq aligner.Bioinformatics2915–21. 10.1093/bioinformatics/bts635
11
DunaifA. (2006). Insulin resistance in women with polycystic ovary syndrome.Fertil. Steril.86S13–S14.
12
ElkonR.UgaldeA. P.AgamiR. (2013). Alternative cleavage and polyadenylation: extent, regulation and function.Nat. Rev. Genet.14496–506. 10.1038/nrg3482
13
FauserB. C. J. M.ChangJ.AzzizR.LegroR.DewaillyD.FranksS. (2004). Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome (PCOS).Hum. Reprod.1941–47. 10.1093/humrep/deh098
14
FrigolaJ.SabarinathanR.MularoniL.MuinosF.Gonzalez-PerezA.Lopez-BigasN. (2017). Reduced mutation rate in exons due to differential mismatch repair.Nat. Genet.491684–1692. 10.1038/ng.3991
15
GinestetC. (2011). ggplot2: elegant graphics for data analysis.J. R. Stat. Soc. Ser. -Stat. Soc.174245–245.
16
Gomez-SuagaP.PaillussonS.MillerC. C. J. (2017). ER-mitochondria signaling regulates autophagy.Autophagy131250–1251. 10.1080/15548627.2017.1317913
17
HarborneL.FlemingR.LyallH.NormanJ.SattarN. (2003). Descriptive review of the evidence for the use of metformin in polycystic ovary syndrome.Lancet3611894–1901. 10.1016/s0140-6736(03)13493-9
18
HarveyA. J. (2019). Mitochondria in early development: linking the microenvironment, metabolism and the epigenome.Reproduction157R159–R179.
19
HendriksM. L.KetJ. C. F.HompesP. G. A.HomburgR.LambalkC. B. (2007). Why does ovarian surgery in PCOS help? Insight into the endocrine implications of ovarian surgery for ovulation induction in polycystic ovary syndrome.Hum. Reprod. Update13249–264. 10.1093/humupd/dml058
20
HsuM. I.LiouT. H.LiangS. J.SuH. W.WuC. H.HsuC. S. (2009). Inappropriate gonadotropin secretion in polycystic ovary syndrome.Ferti. Steril.911168–1174. 10.1016/j.fertnstert.2008.01.036
21
JansenR. P. S.de BoerK. (1998). The bottleneck: mitochondrial imperatives in oogenesis and ovarian follicular fate.Mol. Cell. Endocrinol.14581–88. 10.1016/s0303-7207(98)00173-7
22
KristensenS. G.PorsS. E.AndersenC. Y. (2017). Improving oocyte quality by transfer of autologous mitochondria from fully grown oocytes.Hum. Reprod.32725–732.
23
LandayM.HuangA.AzzizR. (2009). Degree of hyperinsulinemia, independent of androgen levels, is an important determinant of the severity of hirsutism in PCOS.Fertil. Steril.92643–647. 10.1016/j.fertnstert.2008.06.021
24
LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis.Bmc Bioinformatics9:559. 10.1186/1471-2105-9-559
25
LiJ. N.HuL.LiuY. Y.HuangL.MuY.CaiX. H.et al (2015). DDX19A Senses Viral RNA and Mediates NLRP3-dependent inflammasome activation.J. Immunol.1955732–5749. 10.3389/fendo.2019.00485
26
LiY.LiuH.YuQ.LiuH.HuangT.ZhaoS.et al (2019). Growth hormone promotes in vitro maturation of human oocytes.Front. Endocrinol.10:485.
27
LiuQ. W.LiY. M.FengY.LiuC. J.MaJ. L.LiY. F.et al (2016). Single-cell analysis of differences in transcriptomic profiles of oocytes and cumulus cells at GV, MI, MII stages from PCOS patients.Sci. Rep.6:39638.
28
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.15:550.
29
MahanM. E. (2008). Ovulatory disorders and infertility.J. Reprod Med.53:311.
30
MarcC. (2019). org.Hs.eg.db: Genome Wide Annotation for Human. R package Version 3.8.2.
31
MartinouJ. C.YouleR. J. (2011). Mitochondria in apoptosis: bcl-2 family members and mitochondrial dynamics.Dev. Cell2192–101. 10.1016/j.devcel.2011.06.017
32
May-PanloupP.BoucretL.de la BarcaJ. M. C.Desquiret-DumasV.Ferre-L’HotellierV.MoriniereC.et al (2016). Ovarian ageing: the role of mitochondria in oocytes and follicles.Hum. Reprod. Update22725–743. 10.1093/humupd/dmw028
33
McCartneyC. R.MarshallJ. C. (2016). Polycystic ovary syndrome.N. Engl. J. Med.3751398–1399.
34
MoothaV. K.LindgrenC. M.ErikssonK.-F.SubramanianA.SihagS.LeharJ.et al (2003). PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet.34, 267–273. 10.1038/ng1180
35
NahuisM. J.KoseN.BayramN.van DesselH. J. H. M.BraatD. D. M.HamiltonC. J. C. M.et al (2011). Long-term outcomes in women with polycystic ovary syndrome initially randomized to receive laparoscopic electrocautery of the ovaries or ovulation induction with gonadotrophins.Hum. Reprod.261899–1904. 10.1093/humrep/der141
36
NasriS.HumaraB.AnjomshoaaA.MoradiN.GholipourN.MashjoorS.et al (2019). Early hereditary diffuse gastric cancer (eHDGC) is characterized by subtle genomic instability and active DNA damage response.Pathol. Oncol. Res.25711–721. 10.1007/s12253-018-0547-9
37
PicelliS.BjorklundA. K.FaridaniO. R.SagasserS.WinbergG.SandbergR. (2013). Smart-seq2 for sensitive full-length transcriptome profiling in single cells.Nat. Methods101096–1098. 10.1038/nmeth.2639
38
RequenaA.HerreroJ.LanderasJ.NavarroE.NeyroJ. L.SalvadorC.et al (2008). Use of letrozole in assisted reproduction: a systematic review and meta-analysis.Hum. Reprod. Update14571–582. 10.1093/humupd/dmn033
39
Rodriguez-MariA.CanestroC.BreMillerR. A.CatchenJ. M.YanY. L.PostlethwaitJ. H. (2013). Retinoic acid metabolic genes, meiosis, and gonadal sex differentiation in zebrafish.Plos One8:e73951. 10.1371/journal.pone.0073951
40
SalisburyJ. L.SuinoK. M.BusbyR.SpringettM. (2002). Centrin-2 is required for centriole duplication in mammalian cells.Curr. Biol.121287–1292. 10.1016/s0960-9822(02)01019-9
41
Shalom-PazE.HolzerH.SonW. Y.LevinI.TanS. L.AlmogB. (2012). PCOS patients can benefit from in vitro maturation (IVM) of oocytes.Eur. J. Obstetr. Gynecol. and Reprod. Biol.16553–56. 10.1016/j.ejogrb.2012.07.001
42
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks.Genome Res.132498–2504. 10.1101/gr.1239303
43
SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proceed. Nat. Acad. Sci. U.S.A.102, 15545–15550. 10.1073/pnas.0506580102
44
SpinelliJ. B.HaigisM. C. (2018). The multifaceted contributions of mitochondria to cellular metabolism.Nat. Cell Biol.20745–754. 10.1038/s41556-018-0124-1
45
TangT.GlanvilleJ.OrsiN.BarthJ. H.BalenA. H. (2006). The use of metformin for women with PCOS undergoing IVF treatment.Hum. Reprod.211416–1425. 10.1093/humrep/del025
46
TrapnellC.RobertsA.GoffL.PerteaG.KimD.KelleyD. R.et al (2014). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks (vol 7, pg 562, 2012).Nat. Protoc.92513–2513. 10.1038/nprot1014-2513a
47
Vakifahmetoglu-NorbergH.OuchidaA. T.NorbergE. (2017). The role of mitochondria in metabolism and cell death.Biochem. Biophys. Res. Commun.482426–431. 10.1016/j.bbrc.2016.11.088
48
Van BlerkomJ. (2004). Mitochondria in human oogenesis and preimplantation embryogenesis: engines of metabolism, ionic regulation and developmental competence.Reproduction128269–280. 10.1530/rep.1.00240
49
WangL. Y.WangD. H.ZouX. Y.XuC. M. (2009). Mitochondrial functions on oocytes and preimplantation embryos.J. Zhejiang Univ. Sci. B10483–492. 10.1631/jzus.b0820379
50
WangS.ZhengY. X.LiJ. Y.YuY.ZhangW. Q.SongM. S., et al. (2020). Single-cell transcriptomic atlas of primate ovarian aging.Cell180585.e19–600.e19.
51
WildT.BudzowskaM.HellmuthS.EibesS.KaremoreG.BarisicM., et al. (2018). Deletion of APC7 or APC16 allows proliferation of human cells without the spindle assembly checkpoint.Cell Rep.252317.e3–2328.e5.
52
WitchelS. F.BurghardA. C.TaoR. H.OberfieldS. E. (2019). The diagnosis and treatment of PCOS in adolescents: an update.Curr. Opin. Pediatr.31562–569. 10.1097/mop.0000000000000778
53
YuG. (2019). Enrichplot: Visualization of Functional Enrichment Result. R package Version 1.6.1.
54
YuG. C.WangL. G.HanY. Y.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters.Omics J. Integr. Biol.16284–287. 10.1089/omi.2011.0118
55
ZhouY. Y.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets.Nat. Commun.10:1523.
Summary
Keywords
polycystic ovary syndrome, mitochondrial dynamics, oocyte, single-cell sequence, transcriptomic analysis
Citation
Qi L, Liu B, Chen X, Liu Q, Li W, Lv B, Xu X, Wang L, Zeng Q, Xue J and Xue Z (2020) Single-Cell Transcriptomic Analysis Reveals Mitochondrial Dynamics in Oocytes of Patients With Polycystic Ovary Syndrome. Front. Genet. 11:396. doi: 10.3389/fgene.2020.00396
Received
15 January 2020
Accepted
30 March 2020
Published
30 April 2020
Volume
11 - 2020
Edited by
Ahmed Rebai, Centre of Biotechnology of Sfax, Tunisia
Reviewed by
Said Assou, Université de Montpellier, France; Salma Ben-Salem, Cleveland Clinic, United States
Updates
Copyright
© 2020 Qi, Liu, Chen, Liu, Li, Lv, Xu, Wang, Zeng, Xue and Xue.
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: Zhigang Xue, xuezg@tongji.edu.cnJinfeng Xue, xuejinfeng@me.comQiao Zeng, zengqiao@cellpro.com.cn
†These authors have contributed equally to this work
This article was submitted to Genomic Medicine, a section of the journal Frontiers in Genetics
Disclaimer
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.