==== Front J Pers Med J Pers Med jpm Journal of Personalized Medicine 2075-4426 MDPI 33114509 10.3390/jpm10040193 jpm-10-00193 Article Association of an APBA3 Missense Variant with Risk of Premature Ovarian Failure in the Korean Female Population Park JeongMan 1† Park YoungJoon 1† https://orcid.org/0000-0001-6896-9748Koh Insong 2 https://orcid.org/0000-0003-0541-3528Kim Nam Keun 1 https://orcid.org/0000-0001-7662-7190Baek Kwang-Hyun 1 Yun Bo-Seong 3 Lee Kyung Ju 4 Song Jae Yun 5 Lee Eunil 5 Kwack KyuBum 1* 1 Institute Department of Biomedical Science, College of Life Science, CHA University, Seongnam, Gyeonggi-do 13488, Korea; jungman.park@chauniv.ac.kr (J.P.); yjparkcb@gmail.com (Y.P.); nkkim@cha.ac.kr (N.K.K.); baek@cha.ac.kr (K.-H.B.) 2 Department of Biomedical Informatics, Hanyang University, Seoul 04763, Korea; insong@hanyang.ac.kr 3 Department of Obstetrics and Gynecology, CHA Gangnam Medical Center, CHA University, Seongnam, Gyeonggi-do 13497, Korea; bosungyun@chamc.co.kr 4 Department of Obstetrics and Gynecology, Korea University Medical Center, Seoul 02841, Korea; drlkj52551@korea.ac.kr 5 Department of Preventive Medicine, College of Medicine, Korea University, Seoul 02841, Korea; sjyuni105@gmail.com (J.Y.S.); eunil@korea.ac.kr (E.L.) * Correspondence: kbkwack@cha.ac.kr; Tel.: +82-31-881-7141† Equally contributed to this study as first authors. 26 10 2020 12 2020 10 4 19321 9 2020 22 10 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).Premature ovarian failure (POF) is a complex disease of which the etiology is influenced by numerous genetic variations. Several POF candidate genes have been reported. However, no causal genes with high odds ratio (OR) have yet been discovered. This study included 564 females of Korean ethnicity, comprising 60 patients with POF and 182 controls in the discovery set and 105 patients with POF and 217 controls in the replication set. We conducted genome-wide association analysis to search for novel candidate genes predicted to influence POF development using Axiom Precision Medicine Research Arrays and additive model logistic regression analysis. One statistically significant single nucleotide polymorphism (SNP), rs55941146, which encodes a missense alteration (Val > Gly) in the APBA3 gene, was identified with OR values for association with POF of 13.33 and 4.628 in the discovery and replication sets, respectively. No rs55941146 minor allele homozygotes were present in either cases or controls. The APBA3 protein binds FIH-1 that inhibits hypoxia inducible factor-1α (HIF-1α). HIF-1α contributes to granulosa cell proliferation, which is crucial for ovarian follicle growth, by regulating cell proliferation factors and follicle stimulating hormone-mediated autophagy. Our data demonstrate that APBA3 is a candidate novel causal gene for POF. follicle stimulating hormoneprimary ovarian failuregenome-wide association studyhypoxia-inducible factor 1ovarian folliclecell proliferationautophagy ==== Body 1. Introduction Premature ovarian failure (POF) is an idiopathic, complex disease in which menopause occurs before the age of 40 years [1]. The exact cause of POF remains unknown. However, many factors may contribute to its development, including autoimmune disease, radiation therapy, anti-cancer drugs, chromosomal abnormalities, and mental status [2,3,4]. The diagnostic criteria for POF in premenopausal females include an increase of serum follicle stimulating hormone (FSH) levels (>40 mIU/mL), measured twice in the same month, alongside the presence of amenorrhea for 6 months before the normal age of menopause onset [5]. Without appropriate hormone replacement therapy, women with POF can develop severe health issues, including not only failure of normal ovary function, but also various other problems, including cardiovascular disease, coronary artery disease, and stroke [6,7]. The ovarian follicle is important to oocyte growth and folliculogenesis is a crucial mechanism contributing to female fertility [8,9,10]. There are four stages in folliculogenesis. In the primordial stage, small, dormant primordial follicles are enclosed by a single layer of granulosa cells (GCs) in the ovarian cortex [8,9]. In the primary stage, oocytes and GCs exhibit dramatic growth, and meanwhile, GCs change shape from a flat to cuboid form [8,9]. In the secondary stage, stroma-like theca cells enclose the follicles, and GCs increase until there are nine layers [8,9]. In the tertiary and preovulatory stage, the antrum is fully formed in an FSH dependent manner, and all follicles, except for one, undergo follicular atresia [11,12,13]. GCs fulfill their important roles in the ovarian follicle by interacting with oocytes via gap junctions to provide nutrients and signaling molecules [14]. GCs are also important during follicle development under the influence of FSH and luteinizing hormone, which induce GC proliferation [13,15,16]. Furthermore, FSH promotes the development of preantral follicles and induces an anti-apoptotic process in antral follicles [15,16,17,18]. There have been numerous reports of associations between single nucleotide polymorphisms (SNPs) and POF onset [19,20,21]. Therefore, we tried to find candidates in rare variants using the precision medicine research array (PMRA) chip, which contains low allele frequency markers. By conducting a genome-wide association study (GWAS) to identify novel variants associated with POF in Korean women, we detected an SNP variant that could substantially increase the risk of POF occurrence. 2. Results 2.1. Patients and Kinship Analysis Subjects comprised a total of 564 females of Korean ethnicity. Sixty patients with POF and 182 controls in the discovery set, and 120 patients with POF and 218 controls in the replication set. Patient genomic DNA samples were genotyped using an PMRA. Kinship analysis detected 30,876 associations among 242 individuals in the discovery set, and there were no first and second degree relationships. Kinship analysis was also conducted in the replication set and a total of 322 individuals were left. Our results indicated that none of the study participants were related. 2.2. GWAS Raw data generated by genotyping of the discovery set using Axiom PMRA chips were first filtered by removing markers with missing annotation and deletion/insertion markers (n = 346,668). All individuals passed the individual-level missingness threshold of <0.1. Application of Hardy-Weinberg equilibrium (HWE) threshold of >1.0 × 10−5 to the control group resulted in the exclusion of 7937 markers. A marker-level missingness threshold of >0.001 excluded 210,158 markers because of low genotyping rates, and the minor allele frequency threshold (<0.05) further excluded 255,378 markers. Finally, only autosomal markers were included in our analysis. Hence, 20,046 markers mapped to the X and Y chromosomes or mitochondrial DNA were excluded. Consequently, a total of 625,170 markers were removed, while 277,022 remained. After additive model association logistic regression analysis, 45 SNPs were identified as significantly associated with POF in the discovery set. An identical analysis was then conducted using the replication set, resulting in validation of only one of the 45 SNPs from the discovery set with a significant p-value (2.996 × 10−7) followed by a Bonferroni correction (Table 1). This significant SNP, rs55941146, maps to the APBA3 coding region (A > C) on chromosome 19. In the discovery set, significant data were visualized with Manhattan and quantile-quantile plots (Figure 1, Supplementary Figure S1). Furthermore, the recombination rate of rs55941146 with other SNPs in a range of ±400 kb was analyzed using LocusZoom (Figure 2). Additive model analysis of rs55941146 and POF association, generated odds ratio (OR) values of 13.33 and 4.628 in the discovery and validation sets, respectively (Table 1). The entire process of quality control (QC), applied threshold in each process, and the number of excluded markers are included in the flow chart (Supplementary Figure S2). The minor allele of rs55941146 was only present in the heterozygous form in both patients and controls (Supplementary Table S1). The minor allele frequencies of the associated SNP in other populations, according to 1000 genomes data, was compared with that in the POF group (Supplementary Table S2). The frequency of variation of rs55941146 is 0 in the East Asian group, but 0.4 in Korean POF patients. 2.3. Predicted Influence of rs55941146 on Protein Structure The rs55941146 variant is a missense SNP (A > C) in the APBA3 coding region, which causes a valine to glycine substitution at residue 206 in the 575 amino acid full-length APBA3 protein. Scores (0.01 and 0.999) generated using the Sorting Intolerant From Tolerant (SIFT) and PolyPhen-2 programs, respectively, indicated that the rs55941146 variant is predicted to have a deleterious effect on the APBA3 amino acid sequence, and is likely damaging to the APBA3 three-dimensional structure. 3. Discussion The APBA3 gene encodes amyloid-beta precursor protein binding family A member 3, which is also referred to as mammalian uncoordinated 18-1 (MUNC 18-1) interacting protein 3 (MINT3). APBA3 functions both to modulate processing of the amyloid-beta precursor protein (APP) by binding to its C-terminal domain and regulating factor-inhibiting hypoxia inducible factor-1 (FIH-1) via its N-terminal domain [22,23]. FIH-1 inhibits hypoxia inducible factor-1 (HIF-1), which regulates glucose metabolism under hypoxic conditions [24]. FIH-1 is an asparaginyl hydroxylase enzyme that promotes asparaginyl hydroxylation of the COOH-terminal transactivation domain (CAD) of HIF-1, thereby reducing HIF-1 function [25,26]. Both APBA3 and HIF-1α contain identical domains that compete for binding to FIH-1 [22]. Hence, if APBA3 is bound to FIH-1, the asparaginyl hydroxylase modification of the HIF-1 CAD region mediated by FIH-1 is inhibited [22]. Therefore, inhibition of FIH-1 by APBA3 leads to increased HIF-1 expression. In contrast, if HIF-1α binds to FIH-1, the asparagine in the CAD region of HIF-1α is modified [27], leading to the degradation of HIF-1α by the ubiquitin-proteasome pathway [28]. Under normal conditions, FIH-1 interacts with HIF-1α leading to the HIF-1 degradation, while, during hypoxia, HIF-1α is stabilized and activated due to the interaction of FIH-1 with APBA3 [28]. Gonadotrophins, including FSH, have established roles in stimulating follicle growth and preventing the GC apoptosis associated with follicle atresia [29,30,31]. There is ample evidence supporting the importance of HIF-1α in angiogenesis, cell proliferation, and metabolic conversion from oxidative phosphorylation to glycolysis [32,33,34,35,36]. Inevitably, various stresses, including hypoxia and nutritional stress, occur during follicle growth, which involves intense cell proliferation [36]. During ovarian follicle growth, cell proliferation is promoted by FSH, which stimulates accumulation of HIF-1, and HIF-1α increases in response to treatment with FSH both in vitro and in vivo [29,37]. Hence, HIF-1α is a factor downstream of FSH [38], and FSH also stimulates HIF-1α transcription and translation in ovarian cancer cells [39]. Ovarian follicle atresia, in which immature follicles degenerate, is an important stage of follicle growth [40], triggered by GC apoptosis [41]. The enhancement of autophagy stimulates GC apoptosis [42] and is induced by conditions occurring in primordial follicles, including starvation [43]. The absence of autophagy leads to the accumulation of aging-related catabolic waste products during folliculogenesis [29,44]. Hence, autophagy protects ovarian follicles, and, specifically, oocytes, from abnormal conditions, including starvation, which occur in primordial follicles [43]. In summary, appropriate FSH-mediated autophagy, which is important for removing waste and maintaining metabolic balance, is necessary for ovarian follicle growth and preservation of primordial follicles [29,45]. HIF-1α is important, not only as a factor downstream of FSH-mediated autophagy, but also an inducer of proliferation factors. Knockdown of HIF-1α induces downregulation of proliferation markers, such as cyclin D2 (CCND2) and proliferating cell nuclear antigen (PCNA) [46]. Furthermore, HIF-1 regulates cell proliferation differentially in hypoxia and normoxia [46]. PCNA and CCND2 are proliferation markers in various tissues, including ovarian follicles [47,48], and both markers are used to assess GC proliferation levels both in vitro and in vivo [46,49,50,51]. During ovarian follicle growth, HIF-1α influences the expression levels of numerous factors, including PCNA and CCND2 [48]. The majority of candidate genes identified by GWAS as associated with various diseases have OR values < 1.5, and such genes with relatively low OR values have a weak impact in increasing disease risk [52]. Compared with the majority of reported GWAS findings, our resulting OR was remarkably high. Consequently, we infer that rs55941146 likely has a substantial influence on POF development. In conclusion, we identified a variant with a high OR for association with POF, relative to previously reported candidate genes, which is predicted to have a detrimental impact on the amino acid sequence and tertiary structure of the APBA3 protein. This APBA3 SNP (rs55941146) may influence FSH-mediated autophagy and transcription factor induction by association with HIF-1α in granulosa cells. This would be expected to induce down-regulation of autophagy, and of various transcription factors in pre-antral and antral follicles stages, under hypoxic conditions. Variation in APBA3, which regulates the FIH-1/HIF-1 pathway, may lead to impaired FIH-1 downregulation during hypoxia and consequent inappropriate inactivation of HIF-1. Reduction in the HIF-1α level can lead to suppression of normal levels of autophagy and increased transcriptional activity. Therefore, our data suggest that APBA3 is associated with POF in the Korean female population and represents a new candidate gene for this condition. 4. Materials and Methods 4.1. Patient Recruitment For the discovery set, 242 women were, retrospectively, selected for inclusion in this study from a total of 367 individuals who visited Korea University Anam Medical Center from 2016 to 2019 and Bundang CHA Hospital until 2010. For the replication set, 322 women were, retrospectively, selected for inclusion in the analysis from 338 individuals who gathered by CHA University until 2004. Samples used in this GWAS study was authorized by the Institutional Review Board of Korea University Anam Medical Center (2016AN0216) and the Institutional Review Board of Bundang CHA Hospital (2010123). The samples used in the replication study were authorized by the Institutional Review Board of CHA University (1044308-201310-BR-002-01). 4.2. Whole Genome Genotyping Using the PMRA Chip DNA was extracted from peripheral blood samples from the 580 individuals recruited to the study. SNP array analysis was conducted using Axiom Precision Medicine Research Array (PMRA) 2.0 chips. SNPs (n = 902,380) were initially evaluated in samples from the 242 individuals in the discovery set. SNP data pre-processing was performed using the Affymetrix Power Tool to generate dish quality control (DQC) values to ensure samples were of sufficient quality for analysis. The dish quality control (DQC) threshold value was 0.82, where samples with DQC values < 0.82 were excluded from further genotype analysis. Quality control of the remaining samples was by sex check to identify differences between clinically determined sex and sex predicted based on genotype data. In the following step, samples with call rate values less than the threshold (97%) were removed. Finally, abnormal rates of heterozygosity and plate quality control values were evaluated and relationship tests were performed. 4.3. GWAS Process Korean patients with POF (n = 60) and controls (n = 182) in the discovery set, and patients (n = 105) and controls (n = 217) in the replication set, were included, comprising a total of 564 samples. Data were analyzed using descriptive statistics and additive model logistic regression analysis. For the additive model association analysis, thresholds were as follows: individual-level missingness <0.1, marker-level missingness <0.001, and minor allele frequency >0.05. HWE was applied only in controls, using an empirical p-value threshold of >1.0 × 10−5. Finally, 166,866 SNPs were included in the analysis. 4.4. Statistical Analysis To account for multiple comparisons, p-values were corrected using the Bonferroni method by applying the formula, α = 0.05/N. SNP data analysis was conducted using PLINK v1.07 (Free Software Foundation, Inc. Boston, MA, USA) [53], LocusZoom v0.12 (Department of Biostatistics and Center for Statistical Genetics, University of Michigan, Ann Arbor, MI, USA) [54], and qqman (a package in R, A Language and Environment for Statistical Computing, R Core Team, R Foundation for Statistical Computing, Vienna, Austria [55]. Kinship analysis was conducted using Kinship-based Inference for Genome-wide association studies (KING) [56]. 4.5. Protein Structure Analysis SIFT v6.2.1 (https://sift.bii.a-star.edu.sg/, 180820) was used to analyze the functional consequences of the SNPs identified and protein biological function and PolyPhen-2 v2.2.2. (http://genetics.bwh.harvard.edu/pph2, 180820) was used to predict effects of variants on both the amino acid sequence and protein tertiary structure. Acknowledgments We thank Thomas J. Kwack, for the review and editing of the manuscript. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Materials The following are available online at https://www.mdpi.com/2075-4426/10/4/193/s1, Figure S1. Quantile-quantile plot for GWAS data from 60 patients with premature ovarian failure and 182 controls. Figure S2. Flow chart of data processing in the discovery and replication sets. Table S1. Genotypes of rs55941146 in samples included in this study. Table S2. Comparison of allele frequency in recruited samples and variable population of 1000 genomes. EAS, east Asian; CDX, Chinese Dai in Xishuangbanna, China; CHB, Han Chinese in Beijing, China; CHS, Southern Han Chinese; JPT, Japanese in Tokyo, Japan; KHV, Kinh in Ho Chi Minh City, Vietnam. Click here for additional data file. Author Contributions Conceptualization, K.-H.B., J.Y.S., E.L., and K.K. Data curation, N.K.K., B.-S.Y., K.J.L., and K.K. Formal analysis, J.P., Y.P., and I.K. Funding acquisition, K.-H.B., J.Y.S., E.L., and K.K. Methodology, J.P., Y.P., and I.K. Project administration, K.-H.B., B.-S.Y., K.J.L., J.Y.S., E.L., and K.K. Resources, N.K.K., K.-H.B., B.-S.Y., K.J.L., J.Y.S., and E.L. Software, J.P., Y.P., and I.K. Supervision, I.K., K.-H.B., B.-S.Y., K.J.L., J.Y.S., E.L., and K.K. Validation, Y.P. and N.K.K. Visualization, J.P. Writing—original draft, J.P. and Y.P. Writing—review & editing, K.-H.B., J.Y.S., E.L., and K.K. All authors have read and agreed to the published version of the manuscript. Funding This study was funded by the Korea Ministry of Environment (MOE) as ‘the Environmental Health Action Program (2016001360008)’. Conflicts of Interest There is no conflict of interest. Figure 1 Manhattan plot for genome wide association study (GWAS) data from the Korean female population, showing −log10 (p-values) from GWAS and imputation analysis plotted against the chromosome position. Each color represents a different chromosome. The lower line indicates the suggested association threshold (p = 1.0 × 10−5) while the upper line indicates the genome-wide significance threshold (p = 5.0 × 10−8). Figure 2 Regional association plot for loci significantly associated with premature ovarian failure in the Korean female population: chr19:3354338−4154338 [3354338−4154338] (APBA3). jpm-10-00193-t001_Table 1Table 1 The statistical results of additive model logistic regression analysis in the discovery and replication sets. CHR, chromosome. SNP, single nucleotide polymorphism. BP, physical position. A1, minor allele. TEST, logistic regression test method for the genetic model. OR, odds ratio. STAT, Coefficient t-statistic. P, p-value. Information Discovery Replication CHR SNP BP A1 TEST OR STAT p OR STAT p 1 rs3007782 4,472,736 T ADD 4.287 5.488 4.06 × 10−8 1.061 0.1557 0.8763 1 rs74918369 61,457,396 C ADD 8.071 6.242 4.31 × 10−10 2.133 1.283 0.1994 1 rs115143460 1.9 × 108 A ADD 13.33 7.047 1.83 × 10−12 3.37 × 10−9 0.0009 0.9993 1 rs12026894 2.45 × 108 A ADD 6.483 5.903 3.57 × 10−9 1.33 1.01 0.3125 2 rs12615054 22,223,766 G ADD 7.526 6.069 1.29 × 10−9 0.9204 −0.2382 0.8117 2 rs141292341 37,310,903 T ADD 10.14 6.5 8.05 × 10−11 1.38 0.4903 0.6239 2 rs145263938 87,842,313 T ADD 5.618 5.956 2.58 × 10−9 0.5837 −1.115 0.2648 2 rs12692712 1.65 × 108 C ADD 4.92 5.649 1.61 × 10−8 1.668 1.819 0.06891 2 rs79371157 2.28 × 108 T ADD 9.601 6.56 5.38 × 10−11 0.812 −0.3447 0.7303 3 rs28630998 13,245,725 G ADD 8.74 6.678 2.42 × 10−11 4.312 3.007 0.00264 3 rs2323277 74,110,248 A ADD 11.5 6.684 2.33 × 10−11 1.209 0.4123 0.6801 4 rs141002338 24,750,774 C ADD 10.45 6.536 6.32 × 10−11 0.5075 −1.032 0.3019 4 rs1510746 67,469,201 C ADD 6.439 5.839 5.25 × 10−9 0.6251 −0.9697 0.3322 5 rs144522971 22,316,095 T ADD 11.99 6.853 7.22 × 10−12 2.813 1.338 0.1811 5 rs296486 1.13 × 108 A ADD 8.526 6.41 1.45 × 10−10 2.628 2.506 0.01222 6 rs117716146 16,032,432 T ADD 11.65 6.789 1.13 × 10−11 0.8826 −0.1782 0.8585 6 rs111721931 1.1 × 108 G ADD 13.21 7.047 1.83 × 10−12 1.26 × 10−9 −0.0019 0.9985 6 rs137945470 1.18 × 108 A ADD 14.85 7.153 8.51 × 10−13 2.078 0.5152 0.6064 6 rs11759078 1.53 × 108 G ADD 4.917 5.565 2.62 × 10−8 0.9774 −0.0669 0.9467 6 rs113416075 1.66 × 108 A ADD 7.197 5.911 3.39 × 10−9 0.577 −1.466 0.1428 7 rs11761631 6,932,435 A ADD 3.595 5.276 1.32 × 10−7 0.9426 −0.2598 0.795 7 rs3807170 37,898,032 T ADD 6.065 5.735 9.76 × 10−9 0.721 −0.7404 0.4591 8 rs140245000 55,592,196 A ADD 10.09 6.695 2.16 × 10−11 2.11 1.041 0.298 8 rs137854443 1.45 × 108 C ADD 0.208 −5.769 7.99 × 10−9 NA NA NA 9 rs77096227 8,010,136 C ADD 3.712 5.139 2.76 × 10−7 0.9232 −0.3035 0.7615 9 rs149748677 78,256,026 A ADD 8.312 6.272 3.57 × 10−10 0.5173 −0.5861 0.5578 9 rs146186513 1.31 × 108 G ADD 13.33 7.047 1.83 × 10−12 NA NA NA 9 rs77609276 1.38 × 108 T ADD 8.534 6.201 5.60 × 10−10 0.7848 −0.4485 0.6538 10 rs78014704 1,842,302 T ADD 9.568 6.59 4.39 × 10−11 0.5369 −1.34 0.1803 10 rs34498403 86,799,865 T ADD 8.2 6.455 1.09 × 10−10 1.373 0.3439 0.7309 12 rs10877123 58,881,273 A ADD 3.971 5.498 3.84 × 10−8 1.111 0.4515 0.6516 12 rs74650809 92,513,269 A ADD 3.578 5.176 2.27 × 10−7 1.05 0.2121 0.8321 13 rs56238336 23,029,868 C ADD 4.942 5.737 9.65 × 10−9 1.546 1.448 0.1476 13 rs146581261 1.11 × 108 A ADD 14.3 6.992 2.71 × 10−12 1.682 0.9786 0.3278 14 rs1957293 92,708,417 A ADD 4.209 5.448 5.10 × 10−8 1.068 0.2128 0.8315 15 rs78879506 37,133,840 A ADD 10.6 6.534 6.39 × 10−11 1.125 0.3918 0.6952 15 rs57049930 93,855,283 T ADD 5.072 5.615 1.97 × 10−8 0.9163 −0.2447 0.8067 16 rs77141302 54,570,512 C ADD 4.301 5.826 5.68 × 10−9 1.138 0.4584 0.6466 16 rs149428 58,668,566 A ADD 8.74 6.479 9.21 × 10−11 0.2768 −1.683 0.09247 17 rs146234219 9,038,490 T ADD 13.33 7.047 1.83 × 10−12 3.36 × 10−9 0.0009 0.9993 17 rs17773918 11,399,739 G ADD 8.535 6.425 1.32 × 10−10 3.569 1.718 0.08571 19 rs1715092 3,557,005 A ADD 4.858 5.762 8.31 × 10−9 0.9295 −0.2799 0.7796 19 rs55941146 3,754,338 C ADD 13.33 7.047 1.83 × 10−12 4.75 3.885 0.000102 19 rs918371 48,808,545 T ADD 11.69 6.854 7.17 × 10−12 1.29 × 10−9 −0.0008 0.9993 21 rs16991683 35,800,134 A ADD 4.26 5.47 4.49 × 10−8 0.7194 −1.104 0.2696 ==== Refs References 1. Coulam C.B. Adamson S.C. Annegers J.F. Incidence of premature ovarian failure Obstet. Gynecol. 1986 67 604 10.1097/00006254-198703000-00020 3960433 2. Fortuño C. Labarta E. Genetics of primary ovarian insufficiency: A review J. Assist. Reprod. Genet. 2014 31 1573 1585 10.1007/s10815-014-0342-9 25227694 3. Qin Y. Jiao X. Simpson J.L. Chen Z.-J. Genetics of primary ovarian insufficiency: New developments and opportunities Hum. Reprod. Updat. 2015 21 787 808 10.1093/humupd/dmv036 4. Jiao X. Zhang H. Ke H. Zhang J. Cheng L. Liu Y. Qin Y. Chen Z.-J. Premature Ovarian Insufficiency: Phenotypic Characterization Within Different Etiologies J. Clin. Endocrinol. Metab. 2017 102 2281 2290 10.1210/jc.2016-3960 28368522 5. Goswami D. Conway G.S. Premature Ovarian Failure Horm. Res. Paediatr. 2007 68 196 202 10.1159/000102537 6. Kalantaridou S.N. Naka K.K. Papanikolaou E. Kazakos N. Kravariti M. Calis K.A. Paraskevaidis E.A. Sideris D.A. Tsatsoulis A. Chrousos G.P. Impaired Endothelial Function in Young Women with Premature Ovarian Failure: Normalization with Hormone Therapy J. Clin. Endocrinol. Metab. 2004 89 3907 3913 10.1210/jc.2004-0015 15292326 7. Eshre Guideline Group on POI Webber L. Davies M. Anderson R. Bartlett J. Braat D. Cartwright B. Cifkova R. Keizer-Schrama S.D.M. Hogervorst E. ESHRE Guideline: Management of women with premature ovarian insufficiency Hum. Reprod. 2016 31 926 937 27008889 8. Smitz J. Cortvrindt R. The earliest stages of folliculogenesis in vitro Reproduction 2002 123 185 202 10.1530/rep.0.1230185 11866686 9. Rimon-Dahari N. Yerushalmi-Heinemann L. Alyagor L. Dekel N. Ovarian Folliculogenesis Results Probl. Cell Differ. 2016 58 167 190 27300179 10. Patiño R. Sullivan C.V. Ovarian follicle growth, maturation, and ovulation in teleost fish Fish Physiol. Biochem. 2002 26 57 70 10.1023/A:1023311613987 11. Baerwald A.R. Adams G.P. Pierson R.A. Ovarian antral folliculogenesis during the human menstrual cycle: A review Hum. Reprod. Updat. 2011 18 73 91 10.1093/humupd/dmr039 22068695 12. Johnson A.L. Ovarian follicle selection and granulosa cell differentiation Poult. Sci. 2015 94 781 785 10.3382/ps/peu008 25535403 13. Kaipia A. Hsueh A.J.W. Regulation of ovarian follicle atresia Annu. Rev. Physiol. 1997 59 349 363 10.1146/annurev.physiol.59.1.349 9074768 14. Eppig J.J. Intercommunication between mammalian oocytes and companion somatic cells BioEssays 1991 10.1002/bies.950131105 1772412 15. Alam H. Miyano T. Interaction between growing oocytes and granulosa cells in vitro Reprod. Med. Biol. 2019 19 13 23 10.1002/rmb2.12292 31956281 16. Hsueh A.J. Billig H. Tsafriri A. Ovarian follicle atresia: A hormonally controlled apoptotic process Endocr. Rev. 1994 15 707 724 7705278 17. McGee E. Spears N. Minami S. Hsu S.-Y. Chun S.-Y. Billig H. Hsueh A.J. Preantral ovarian follicles in serum-free culture: Suppression of apoptosis after activation of the cyclic guanosine 3′, 5′-monophosphate pathway and stimulation of growth and differentiation by follicle-stimulating hormone Endocrinology 1997 138 2417 2424 10.1210/endo.138.6.5164 9165031 18. Richards J.S. Hormonal control of gene expression in the ovary Endocr. Rev. 1994 15 725 751 10.1210/edrv-15-6-725 7705279 19. Ma L. Chen Y. Mei S. Liu C. Ma X. Li Y. Jiang Y. Ha L. Xu X. Single nucleotide polymorphisms in premature ovarian failure-associated genes in a Chinese Hui population Mol. Med. Rep. 2015 12 2529 2538 10.3892/mmr.2015.3762 25954833 20. Qin Y. Sun M. You L. Wei D. Sun J. Liang X. Zhang B. Jiang H. Xu J. Chen Z.-J. ESR1, HK3 and BRSK1 gene variants are associated with both age at natural menopause and premature ovarian failure Orphanet J. Rare Dis. 2012 7 5 10.1186/1750-1172-7-5 22248077 21. Knauff E.A.H. Franke L. Van Es M.A. Berg L.H.V.D. Van Der Schouw Y.T. Laven J.S. Lambalk C.B. Hoek A. Goverde A.J. Christin-Maitre S. Genome-wide association study in premature ovarian failure patients suggests ADAMTS19 as a possible candidate gene Hum. Reprod. 2009 24 2372 2378 10.1093/humrep/dep197 19508998 22. Sakamoto T. Seiki M. Mint3 Enhances the Activity of Hypoxia-inducible Factor-1 (HIF-1) in Macrophages by Suppressing the Activity of Factor Inhibiting HIF-1 J. Biol. Chem. 2009 284 30350 30359 10.1074/jbc.M109.019216 19726677 23. Tanahashi H. Tabira T. Genomic organization of the human X11L2 gene (APBA3), a third member of the X11 protein family interacting with Alzheimerʼs β-amyloid precursor protein NeuroReport 1999 10 2575 2578 10.1097/00001756-199908200-00025 10574372 24. Denko N.C. Hypoxia, HIF1 and glucose metabolism in the solid tumour Nat. Rev. Cancer 2008 8 705 713 10.1038/nrc2468 19143055 25. Lando D. Peet D.J. Gorman J.J. Whelan D.A. Whitelaw M.L. Bruick R.K. FIH-1 is an asparaginyl hydroxylase enzyme that regulates the transcriptional activity of hypoxia-inducible factor Genes Dev. 2002 16 1466 1471 10.1101/gad.991402 12080085 26. Lando D. Peet D.J. Whelan D.A. Gorman J.J. Whitelaw M.L. Asparagine Hydroxylation of the HIF Transactivation Domain: A Hypoxic Switch Science 2002 295 858 861 10.1126/science.1068592 11823643 27. Mahon P.C. FIH-1: A novel protein that interacts with HIF-1alpha and VHL to mediate repression of HIF-1 transcriptional activity Genes Dev. 2001 15 2675 2686 10.1101/gad.924501 11641274 28. Maxwell P.H. Wiesener M.S. Chang G.-W. Clifford S.C. Vaux E.C. Cockman M.E. Wykoff C.C. Pugh C.W. Maher E.R. Ratcliffe P.J. The tumour suppressor protein VHL targets hypoxia-inducible factors for oxygen-dependent proteolysis Nat. Cell Biol. 1999 399 271 275 10.1038/20459 29. Zhou J. Yao W. Li C. Wu W. Li Q. Liu H. Administration of follicle-stimulating hormone induces autophagy via upregulation of HIF-1α in mouse granulosa cells Cell Death Dis. 2017 8 e3001 10.1038/cddis.2017.371 28817115 30. Peluso J.J. Steger R.W. Role of FSH in regulating granulosa cell division and follicular atresia in rats Reproduction 1978 54 275 278 10.1530/jrf.0.0540275 722676 31. Chun S.Y. Eisenhauer K.M. Minami S. Billig H. Perlas E. Hsueh A.J. Hormonal regulation of apoptosis in early antral follicles: Follicle-stimulating hormone as a major survival factor Endocrinology 1996 137 1447 1456 10.1210/endo.137.4.8625923 8625923 32. Semenza G.L. HIF-1 and tumor progression: Pathophysiology and therapeutics Trends Mol. Med. 2002 8 S62 S67 10.1016/S1471-4914(02)02317-1 11927290 33. Krishnamachary B. Berg-Dixon S. Kelly B. Agani F. Feldser D. Ferreira G. Iyer N. LaRusch J. Pak B. Taghavi P. Regulation of colon carcinoma cell invasion by hypoxia-inducible factor 1 Cancer Res. 2003 63 1138 1143 12615733 34. Seagroves T.N. Ryan H.E. Lu H. Wouters B.G. Knapp M. Thibault P. Laderoute K. Johnson R.S. Transcription Factor HIF-1 Is a Necessary Mediator of the Pasteur Effect in Mammalian Cells Mol. Cell. Biol. 2001 21 3436 3444 10.1128/MCB.21.10.3436-3444.2001 11313469 35. Tacchini L. Bianchi L. Bernelli-Zazzera A. Cairo G. Transferrin receptor induction by hypoxia. HIF-1-mediated transcriptional activation and cell-specific post-transcriptional regulation J. Biol. Chem. 1999 274 24142 24146 10.1074/jbc.274.34.24142 10446187 36. Heiden M.G.V. Plas D.R. Rathmell J.C. Fox C.J. Harris M.H. Thompson C.B. Growth Factors Can Influence Cell Growth and Survival through Effects on Glucose Metabolism Mol. Cell. Biol. 2001 21 5899 5912 10.1128/MCB.21.17.5899-5912.2001 11486029 37. Alam H. Weck J. Maizels E. Park Y. Lee E.J. Ashcroft M. Hunzicker-Dunn M. Role of the Phosphatidylinositol-3-Kinase and Extracellular Regulated Kinase Pathways in the Induction of Hypoxia-Inducible Factor (HIF)-1 Activity and the HIF-1 Target Vascular Endothelial Growth Factor in Ovarian Granulosa Cells in Response to Follicle-Stimulating Hormone Endocrinology 2009 150 915 928 18845636 38. Alam H. Maizels E.T. Park Y. Ghaey S. Feiger Z.J. Chandel N.S. Hunzicker-Dunn M. Follicle-stimulating hormone activation of hypoxia-inducible factor-1 by the phosphatidylinositol 3-kinase/AKT/Ras homolog enriched in brain (Rheb)/mammalian target of rapamycin (mTOR) pathway is necessary for induction of select protein markers of follicular differentiation J. Biol. Chem. 2004 279 19431 19440 14982927 39. Huang Y. Hua K. Zhou X. Jin H. Chen X. Lü X. Yu Y. Zha X. Feng Y. Activation of the PI3K/AKT pathway mediates FSH-stimulated VEGF expression in ovarian serous cystadenocarcinoma Cell Res. 2008 18 780 791 10.1038/cr.2008.70 18574502 40. Matsuda F. Inoue N. Manabe N. Ohkura S. Follicular Growth and Atresia in Mammalian Ovaries: Regulation by Survival and Death of Granulosa Cells J. Reprod. Dev. 2012 58 44 50 10.1262/jrd.2011-012 22450284 41. Yu Y.S. Sui H.S. Bin Han Z. Li W. Luo M.J. Tan J.-H. Apoptosis in Granulosa cells during follicular atresia: Relationship with steroids and insulin-like growth factors Cell Res. 2004 14 341 346 10.1038/sj.cr.7290234 15353131 42. Choi J. Jo M. Lee E. Choi D. Induction of apoptotic cell death via accumulation of autophagosomes in rat granulosa cells Fertil. Steril. 2011 95 1482 1486 10.1016/j.fertnstert.2010.06.006 20630503 43. Song Z.-H. Yu H.-Y. Wang P. Mao G.-K. Liu W.-X. Li M.-N. Wang H.-N. Shang Y.-L. Liu C. Xu Z.-L. Germ cell-specific Atg7 knockout results in primary ovarian insufficiency in female mice Cell Death Dis. 2015 6 e1589 10.1038/cddis.2014.559 25590799 44. Wu J. Carlock C. Zhou C. Nakae S. Hicks J. Adams H.P. Lou Y. IL-33 is required for disposal of unnecessary cells during ovarian atresia through regulation of autophagy and macrophage migration J. Immunol. 2015 194 2140 2147 10.4049/jimmunol.1402503 25617473 45. Gawriluk T.R. Hale A.N. Flaws J.A. Dillon C.P. Green D.R. Rucker I.E.B. Autophagy is a cell survival program for female germ cells in the murine ovary Reproduction 2011 141 759 765 10.1530/REP-10-0489 21464117 46. Baddela V.S. Sharma A. Michaelis M. Vanselow J. HIF1 driven transcriptional activity regulates steroidogenesis and proliferation of bovine granulosa cells Sci. Rep. 2020 10 3906 10.1038/s41598-020-60935-1 32127571 47. Xu B. Hua J. Zhang Y. Jiang X. Zhang H. Ma T. Zheng W. Sun R. Shen W. Sha J. Proliferating Cell Nuclear Antigen (PCNA) Regulates Primordial Follicle Assembly by Promoting Apoptosis of Oocytes in Fetal and Neonatal Mouse Ovaries PLoS ONE 2011 6 e16046 10.1371/journal.pone.0016046 21253613 48. Hernández-Ochoa I. Barnett-Ringgold K.R. Dehlinger S.L. Gupta R.K. Leslie T.C. Roby K.F. Flaws J.A. The Ability of the Aryl Hydrocarbon Receptor to Regulate Ovarian Follicle Growth and Estradiol Biosynthesis in Mice Depends on Stage of Sexual Maturity1 Biol. Reprod. 2010 83 698 706 10.1095/biolreprod.110.087015 20631400 49. Christenson L.K. Gunewardena S. Hong X. Spitschak M. Baufeld A. Vanselow J. Research Resource: Preovulatory LH Surge Effects on Follicular Theca and Granulosa Transcriptomes Mol. Endocrinol. 2013 27 1153 1171 10.1210/me.2013-1093 23716604 50. Yenuganti V.R. Viergutz T. Vanselow J. Oleic acid induces specific alterations in the morphology, gene expression and steroid hormone production of cultured bovine granulosa cells Gen. Comp. Endocrinol. 2016 232 134 144 10.1016/j.ygcen.2016.04.020 27118706 51. Hubbi M.E. Semenza G.L. Regulation of cell proliferation by hypoxia-inducible factors Am. J. Physiol. Physiol. 2015 309 C775 C782 10.1152/ajpcell.00279.2015 26491052 52. Hodge S.E. Greenberg D.A. How Can We Explain Very Low Odds Ratios in GWAS? I. Polygenic Models Hum. Hered. 2016 81 173 180 10.1159/000454804 28171865 53. Purcell S. Neale B. Todd-Brown K. Thomas L. Ferreira M.A.R. Bender D. Maller J. Sklar P. De Bakker P.I.W. Daly M.J. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses Am. J. Hum. Genet. 2007 81 559 575 10.1086/519795 17701901 54. Pruim R.J. Welch R.P. Sanna S. Teslovich T.M. Chines P.S. Gliedt T.P. Boehnke M. Abecasis G.R. Willer C.J. LocusZoom: Regional visualization of genome-wide association scan results Bioinformatics 2010 26 2336 2337 10.1093/bioinformatics/btq419 20634204 55. Turner S.D. qqman: An R package for visualizing GWAS results using Q-Q and manhattan plots Biorxiv 2018 3 731 10.21105/joss.00731 56. Manichaikul A. Mychaleckyj J.C. Rich S.S. Daly K. Sale M. Chen W.-M. Robust relationship inference in genome-wide association studies Bioinformatics 2010 26 2867 2873 10.1093/bioinformatics/btq559 20926424