Transcriptome-wide association study identified candidate genes associated with polycystic ovary syndrome
=========================================================================================================
* Yingzhou Shi
* Changrong Li
* Hongjie Yu
* Xinyu Chen
* Ping Shi
* Chao Xu
## Abstract
Polycystic ovary syndrome (PCOS) is one of the most common metabolism and endocrine disorder affecting women of reproductive age. However, the genetic mechanisms of PCOS remain largely elusive. A transcriptome-wide association study(TWAS) was conducted to identify the genes associated with PCOS using gene expression references of the whole blood and ovary. Then mRNA expression profiling analysis was performed to identify common genes with TWAS. Genes detected by TWAS were subjected to Gene Ontology(GO), pathway enrichment and protein-protein analysis. TWAS identified 422 genes for the whole blood and 177 genes for ovary, such as EML6 (TWAS.P = 6.33E-05), GRAP (TWAS.P = 0.000211), MAFTRR(TWAS.P = 1.21E-05). Functional enrichment analysis detected biological processes including ncRNA metabolic process, microtubule-related pathways. The present study identified a number of potential PCOS genes and biological processes. These results may provide insight into the genetic basis of PCOS.
Keyword
* Transcriptome-wide association study
* polycystic ovary syndrome
## Introduction
Polycystic ovary syndrome (PCOS), characterized by hyperandrogenism, irregular anovulation, and polycystic ovarian changes, is one of the most common reproductive metabolic disorders affecting reproductive-aged women. It is estimated that PCOS affect 6–20% of women1, and approximately 80% of those PCOS are infertile2. Although studies have suggested that obesity, age at menopause, insulin resistance/hyperinsulinemia, sex hormone-binding globulin (SHBG) levels, and depression may be risk factors for PCOS, its potential etiologies remain uncertain3.
Genetic factors are thought to be important in the pathogenesis of PCOS. For example, researchers have found that PCOS is twice as common in monozygotic twins as in dizygotic twins and conclude that genetic factors account for about 79% of the liability of PCOS4. The prevalence of PCOS is also higher among the mothers and sisters of affected women, suggesting a phenomenon of familial aggregation5. Several genome-wide association studies(GWAS) have identified susceptibility loci for PCOS6–10. However, the biological significance of their findings remains largely unclear as close to 90% of GWAS-identified SNPs are in the noncoding region. Expression quantitative trait loci (eQTL) are genomic loci that are associated with gene transcript levels. Gene expression is a critical intermediate phenotype between genes and diseases. Such integration of GWAS and eQTL is effective in identifying candidate genes associated with diseases.
The availability and cost of specimens, as well as intrinsic factors, small effects, make large-scale expression-trait associations challenging. To address these issues, transcriptome-wide association studies (TWAS) were developed, which combine gene expression data with large-scale genome-wide association studies. First, datasets (e.g. GTEx) with matched expression data and genotype were used for building a prediction model. Then GWAS data are combined with these models to assess the association between predicted gene expression and disease risk. A growing number of researchers are using TWASs to identify disease-associated genetic loci11–13.
In the present study, we conduct TWAS and mRNA expression profiling analysis to identify genes that are associated with PCOS. Initially, TWAS was performed to identify genes associated with PCOS risk. Then mRNA expression profiling analysis was carried out to identify common genes between TWAS. For the genes identified by TWAS, we conducted an analysis of gene ontology and pathway enrichment to explore their biological functions and related pathways.
## Method
### GWAS summary data for PCOS
We obtained PCOS GWAS summary data from a GWAS meta-analysis of seven cohort studies that included 10074 patients with PCOS and 103164 controls of European ancestry7. The cohort studies include Rotterdam(cases: 1184, controls:5799), UK (cases: 670,controls: 1379), EGCUT(cases:157,controls 2807), deCODE(cases: 658,controls: 6774), Chicago (cases: 984, controls: 2963), Boston(cases: 485,controls: 407), 23andMe(cases: 5184, controls: 82759). Diagnosis of PCOS was made based on National Institutes of Health criteria14, Rotterdam criteria15, or self-report.
### TWAS analysis of PCOS
TWAS analysis was conducted using FUSION software. Pre-computed gene expression weights combined with PCOS summary statistics were used to calculate the association between single gene and PCOS. We derived the gene expression weights for whole blood and the ovary from the FUSION website ([http://gusevlab.org/projects/fusion](http://gusevlab.org/projects/fusion)). In FUSION, expression weights were calculated using five linear models, including the best linear unbiased prediction model (BLUP), Bayesian sparse linear mixed model (BSLMM), least absolute shrinkage and selection operator (LASSO), Elastic Net, and top single nucleotidepolymorphisms (SNPs). In the present study, a TWAS *P*-value was calculated for each gene within the whole blood and ovary for European populations.
### mRNA Expression profiles of PCOS
A genome-wide mRNA expression profile of PCOS was used to identify differentially expressed genes (DEGs). The PCOS expression data were downloaded from the Gene Expression Omnibus (GEO) Datasets GSE34526 ([https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE34526](https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE34526))16. Gonadotropin-stimulated granulosa cells were obtained from 7 patients with PCOS and 3 healthy individuals. GEO2R was used to screen for differentially expressed genes17. Genes with a folding change (FC) >2 or <0.5 (|log2FC |> 1) and p-value < 0.01 were considered to be differentially expressed genes.
### Functional Enrichment
Metascape ([https://metascape.org/gp/index.html](https://metascape.org/gp/index.html)) was used to perform gene ontology and pathway enrichment analyses of the genes identified by TWAS. Metascape analysis is based on databases, including Gene Ontology, KEGG, Reactome Gene Sets, and WikiPathways. A hypergeometric test is used to identify enhanced ontology terms, and the p-value is corrected for false discovery rates. Molecular Complex Detection (MCODE) algorithm was used to extract densely connected regions18.
## Results
### TWAS results of PCOS
A total of genes were identified by TWAS in blood and ovary tissue.
The TWAS identified 422 genes with a p value < 0.05 for the whole blood, such as EML6 (TWAS.P = 6.33E-05), GRAP (TWAS.P = 0.000211), PTGER2(TWAS.P= 0.000403)(Supplementary Table1 and Table2). For the ovary tissue, TWAS identified 177 genes with a p value < 0.05, such as MAFTRR(TWAS.P = 1.21E-05), NEIL2 (TWAS.P = 3.22E-05), ZNF555(TWAS.P= 0.000364). 52 genes with p value < 0.05 were detected in both whole blood and ovary TWAS analysis (Supplementary Table3). Table 1 and Figure1 show the top 20 genes associated with PCOS as determined by TWAS.
View this table:
[Table 1](http://medrxiv.org/content/early/2023/02/10/2023.02.09.23285661/T1)
Table 1
Top 20 genes identified by TWAS analysis
[Figure1](http://medrxiv.org/content/early/2023/02/10/2023.02.09.23285661/F1)
Figure1
Manhattan plot of TWAS-identified top20 genes associated with PCOS
The x-axis represents the physical location of each gene, while the y-axis represents the Z score of association between that gene and PCOS.
### Genes shared between TWAS and mRNA Expression Profiling
The TWAS analysis and the mRNA expression analysis detected 18 common genes, such as CAMK1D, COLEC12 and PLCB2.(Table 2). A Venn diagram of TWAS versus common genes identified by mRNA expression profiling are shown in Figure2.
View this table:
[Table 2](http://medrxiv.org/content/early/2023/02/10/2023.02.09.23285661/T2)
Table 2
Common gene identified by TWAS and mRNA expression profiling
[Figure2](http://medrxiv.org/content/early/2023/02/10/2023.02.09.23285661/F2)
Figure2
Venn diagram of gene identified by TWAS and mRNA expression profiling
### Gene Enrichment Analysis of genes identified by TWAS
In order to obtain an overall impression of the characteristics of these genes identified by TWAS, a total 547 genes were submitted to Metascape performing GO and pathway enrichment, and protein-protein interaction analysis.
GO and pathway enrichment analysis showed that TWAS-identified genes were enriched in ncRNA metabolic process, rRNA metabolic process, transport along microtubule and so on(Figure 3). MCODE enrichment analysis based on PPI enrichment for TWAS-identified genes showed 11 PPI modules(Figure 4). The top enrichment term of MCODE1 included 13 proteins distrusted in three major clusters: “Metabolism of RNA”, “Major pathway of rRNA processing in the nucleolus and cytosol”, “rRNA processing in the nucleus and cytosol”. MCODE2, which consists of 7 proteins was mostly enriched in “Interleukin-20 family signaling”, “JAK-STAT signaling pathway”, “Regulatory circuits of the STAT3 signaling pathway”.
[Figure 3](http://medrxiv.org/content/early/2023/02/10/2023.02.09.23285661/F3)
Figure 3
GO enrichment analysis of the genes identified by TWAS
[Figure 4](http://medrxiv.org/content/early/2023/02/10/2023.02.09.23285661/F4)
Figure 4
Protein-protein interaction network of genes identified by TWAS
## Discussion
The TWAS method is a powerful tool that utilizes genetic variation and gene expression to identify genes whose cis-regulated expression is associated with complex traits. In present study, we performed TWAS to detect candidate genes closely associated with PCOS. Further integrative analysis of TWAS and mRNA expression profiling identified 18 common genes. To the best of our knowledge, this is the first TWAS integrated with mRNA expression profiling analysis for PCOS. The results of our study provide new insight into the genetic mechanism behind PCOS. Our TWAS identified several genes that have been previously reported in GWAS and post-GWAS, such as RPS2619, NEIL28, and ZNF55619 as well as novel genes that have never been reported.
NEIL2 is is the gene encoding a member of the Fpg / Nei DNA glycosylase family. The encoded enzyme is involved in DNA repair mainly by cleaving oxidatively damaged bases and introducing DNA strand breaks through its base cleavage enzyme activity.
The NEIL2 gene encodes NEI endonuclease VIII-like 2, one of the DNA glycosylases involved in DNA repair20. DNA damage has been demonstrated to play a role in the pathogenesis of PCOS21. Further experiments are needed to determine the relationship between NEIL2 and PCOS.
EML6 is the only member of the EML family that is preferentially expressed by oocytes in the mouse ovary24. There is a colocalization of EML6 with meiotic spindles, and this specialized localization plays a critical role in maintaining spindle integrity and homologous chromosome segregation.
RPS26, located on 12q13.2, is a gene encoding a ribosomal protein, which is a component of the 40S subunit and belongs to the ribosomal protein S26E family22. The knockdown of Rps26 in mouse oocytes has previously been shown to result in delayed oocyte development and even the cessation of chromatin structural transition in oocytes, resulting in premature ovarian failure (POF)23.
GO and pathway enrichment indicated that ncRNA metabolic and microtubule-related pathways were enriched. The non-coding RNAs (ncRNAs) are functional RNAs that are not coded for protein production but serve as important regulators of a wide range of biological events25. Several studies have indicated that ncRNAs play a role in the occurrence and development of PCOS. The expression of sncRNAs was significantly different in serum, granulosa cells (GCs), follicular fluid (FF), and other tissues between PCOS patients and the general population26,27. Accordingly, our results further confirm that ncRNA may play a crucial role in the development of PCOS.
Microtubules, which are polymers containing α- and β-tubulin heterodimers, make up the oocyte spindle28. An abnormal spindle can affect the development, differentiation, and growth of a germ cell directly. Recent research has identified several mutations of spindle-related proteins that lead to infertility. An example of this is the TUBB8 gene, which encodes an isoform of beta-tubulin that is specific to humans and primates. TUBB8 mutations in women lead to a range of spindle defects during oocyte meiosis29. The enriched microtubule-related pathways may suggest that microtubules contribute to infertility in women with PCOS.
## Conclusion
In summary, this TWAS study identified novel and susceptible genes and potential pathways for PCOS. Further studies are needed to explore the potential mechanism.
## Data Availability
https://www.repository.cam.ac.uk/handle/1810/283491
[https://www.repository.cam.ac.uk/handle/1810/283491](https://www.repository.cam.ac.uk/handle/1810/283491)
## Acknowledgements
We would like to thank the participants and the investigators for generating the publicly available dataset. The scientific calculations in this paper have been done on the HPC Cloud Platform of Shandong University.
* Received February 9, 2023.
* Revision received February 9, 2023.
* Accepted February 10, 2023.
* © 2023, Posted by Cold Spring Harbor Laboratory
This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)
## Reference
1. 1.Goodarzi MO, Dumesic DA, Chazenbalk G, Azziz R. Polycystic ovary syndrome: etiology, pathogenesis and diagnosis. Nat Rev Endocrinol. 2011;7(4):219–231. doi:10.1038/nrendo.2010.217
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrendo.2010.217&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=21263450&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
2. 2.Melo AS, Ferriani RA, Navarro PA. Treatment of infertility in women with polycystic ovary syndrome: approach to clinical practice. Clinics (Sao Paulo). 2015;70(11):765–769. doi:10.6061/clinics/2015(11)09
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.6061/clinics/2015(11)09&link_type=DOI)
3. 3.Day FR, Hinds DA, Tung JY, et al. Causal mechanisms and balancing selection inferred from genetic associations with polycystic ovary syndrome. Nat Commun. 2015;6:8464. doi:10.1038/ncomms9464
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms9464&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26416764&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
4. 4.Vink JM, Sadrzadeh S, Lambalk CB, Boomsma DI. Heritability of polycystic ovary syndrome in a Dutch twin-family study. J Clin Endocrinol Metab. 2006;91(6):2100–2104. doi:2016092613242900031
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2005-1494&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=16219714&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238095000018&link_type=ISI)
5. 5.Legro RS, Driscoll D, Strauss JF, Fox J, Dunaif A. Evidence for a genetic basis for hyperandrogenemia in polycystic ovary syndrome. Proc Natl Acad Sci U S A. 1998;95(25):14956–14960. doi:10.1073/pnas.95.25.14956
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiOTUvMjUvMTQ5NTYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMi8xMC8yMDIzLjAyLjA5LjIzMjg1NjYxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)
6. 6.Lee H, Oh JY, Sung YA, et al. Genome-wide association study identified new susceptibility loci for polycystic ovary syndrome. Hum Reprod. 2015;30(3):723–731. doi:10.1093/humrep/deu352
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/humrep/deu352&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=25574032&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
7. 7.Day F, Karaderi T, Jones MR, et al. Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria. PLoS Genet. 2018;14(12):e1007813. doi:10.1371/journal.pgen.1007813
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1007813&link_type=DOI)
8. 8.Hayes MG, Urbanek M, Ehrmann DA, et al. Genome-wide association of polycystic ovary syndrome implicates alterations in gonadotropin secretion in European ancestry populations. Nat Commun. 2015;6:7502. doi:10.1038/ncomms8502
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms8502&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26284813&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
9. 9.Shi Y, Zhao H, Shi Y, et al. Genome-wide association study identifies eight new risk loci for polycystic ovary syndrome. Nat Genet. 2012;44(9):1020–1025. doi:10.1038/ng.2384
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2384&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=22885925&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
10. 10.Chen ZJ, Zhao H, He L, et al. Genome-wide association study identifies susceptibility loci for polycystic ovary syndrome on chromosome 2p16.3, 2p21 and 9q33.3. Nat Genet. 2011;43(1):55–59. doi:10.1038/ng.732
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.732&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=21151128&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285683500016&link_type=ISI)
11. 11.Xu J, Zeng Y, Si H, et al. Integrating transcriptome-wide association study and mRNA expression profile identified candidate genes related to hand osteoarthritis. Arthritis Res Ther. 2021;23(1):81. doi:10.1186/s13075-021-02458-2
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13075-021-02458-2&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
12. 12.Doke T, Huang S, Qiu C, et al. Transcriptome-wide association analysis identifies DACH1 as a kidney disease risk gene that contributes to fibrosis. J Clin Invest. 2021;131(10):141801. doi:10.1172/JCI141801
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI141801&link_type=DOI)
13. 13.Yao S, Zhang X, Zou SC, et al. A transcriptome-wide association study identifies susceptibility genes for Parkinson’s disease. NPJ Parkinsons Dis. 2021;7:79. doi:10.1038/s41531-021-00221-7
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41531-021-00221-7&link_type=DOI)
14. 14.Carmina E. Diagnosis of polycystic ovary syndrome: from NIH criteria to ESHRE-ASRM guidelines. Minerva Ginecol. 2004;56(1):1–6.
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=14973405&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
15. 15.Rotterdam ESHRE/ASRM-Sponsored PCOS consensus workshop group. Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome (PCOS). Hum Reprod. 2004;19(1):41–47. doi:10.1093/humrep/deh098
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/humrep/deh098&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=14688154&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187827100007&link_type=ISI)
16. 16.Kaur S, Archer KJ, Devi MG, Kriplani A, Strauss JF, Singh R. Differential gene expression in granulosa cells from polycystic ovary syndrome patients with and without insulin resistance: identification of susceptibility gene sets through network analysis. J Clin Endocrinol Metab. 2012;97(10):E2016–2021. doi:2016092613074200069
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2011-3441&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=22904171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309664400023&link_type=ISI)
17. 17.Barrett T, Wilhite SE, Ledoux P, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013;41(Database issue):D991–995. doi:10.1093/nar/gks1193
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gks1193&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=23193258&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312893300141&link_type=ISI)
18. 18.Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4(1):2. doi:10.1186/1471-2105-4-2
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-4-2&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=12525261&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
19. 19.Sun Q, Gao Y, Yang J, Lu J, Feng W, Yang W. Mendelian Randomization Analysis Identified Potential Genes Pleiotropically Associated with Polycystic Ovary Syndrome. Reprod Sci. 2022;29(3):1028–1037. doi:10.1007/s43032-021-00776-z
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s43032-021-00776-z&link_type=DOI)
20. 20.Bandaru V, Sunkara S, Wallace SS, Bond JP. A novel human DNA glycosylase that removes oxidative DNA damage and is homologous to Escherichia coli endonuclease VIII. DNA Repair (Amst). 2002;1(7):517–529. doi:10.1016/s1568-7864(02)00036-8
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1568-7864(02)00036-8&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=12509226&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
21. 21.Dincer Y, Akcay T, Erdem T, Ilker Saygili E, Gundogdu S. DNA damage, DNA susceptibility to oxidation and glutathione level in women with polycystic ovary syndrome. Scandinavian Journal of Clinical and Laboratory Investigation. 2005;65(8):721–728. doi:10.1080/00365510500375263
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00365510500375263&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=16509054&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233617800010&link_type=ISI)
22. 22.Filipenko ML, Vinichenko NA, Karpova GG, Mertvetsov NP, Amaldi F. Isolation, structural analysis and mapping of the functional gene of human ribosomal protein S26. Gene. 1998;211(2):287–292. doi:10.1016/s0378-1119(98)00108-5
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0378-1119(98)00108-5&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=9602156&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
23. 23.Liu XM, Yan MQ, Ji SY, et al. Loss of oocyte Rps26 in mice arrests oocyte growth and causes premature ovarian failure. Cell Death Dis. 2018;9(12):1144. doi:10.1038/s41419-018-1196-3
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41419-018-1196-3&link_type=DOI)
24. 24.Yin H, Hou X, Zhang T, Shi L, Su YQ. Participation of EML6 in the regulation of oocyte meiotic progression in mice. J Biomed Res. 2020;34(1):44–53. doi:10.7555/JBR.33.20190014
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7555/JBR.33.20190014&link_type=DOI)
25. 25.Mattick JS, Makunin IV. Non-coding RNA. Hum Mol Genet. 2006;15 Spec No 1:R17-29. doi:10.1093/hmg/ddl046
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddl046&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=16651366&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237320500003&link_type=ISI)
26. 26.Naji M, Aleyasin A, Nekoonam S, Arefian E, Mahdian R, Amidi F. Differential Expression of miR-93 and miR-21 in Granulosa Cells and Follicular Fluid of Polycystic Ovary Syndrome Associating with Different Phenotypes. Sci Rep. 2017;7(1):14671. doi:10.1038/s41598-017-13250-1
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-13250-1&link_type=DOI)
27. 27.Murri M, Insenser M, Fernández-Durán E, San-Millán JL, Escobar-Morreale HF. Effects of polycystic ovary syndrome (PCOS), sex hormones, and obesity on circulating miRNA-21, miRNA-27b, miRNA-103, and miRNA-155 expression. J Clin Endocrinol Metab. 2013;98(11):E1835–1844. doi:2016092613073400868
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2013-2218&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=24037889&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000327180700020&link_type=ISI)
28. 28.Goshima G, Kimura A. New look inside the spindle: microtubule-dependent microtubule generation within the spindle. Curr Opin Cell Biol. 2010;22(1):44–49. doi:10.1016/j.ceb.2009.11.012
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ceb.2009.11.012&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=20022736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275219200007&link_type=ISI)
29. 29.Feng R, Sang Q, Kuang Y, et al. Mutations in TUBB8 and Human Oocyte Meiotic Arrest. N Engl J Med. 2016;374(3):223–232. doi:10.1056/NEJMoa1510791
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1510791&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26789871&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F10%2F2023.02.09.23285661.atom)