Back to Journals » Diabetes, Metabolic Syndrome and Obesity » Volume 13

Identification of Genetic Variants for Female Obesity and Evaluation of the Causal Role of Genetically Defined Obesity in Polycystic Ovarian Syndrome

Authors Ahn Y , Lee H, Cho YS 

Received 11 September 2020

Accepted for publication 21 October 2020

Published 12 November 2020 Volume 2020:13 Pages 4311—4322

DOI https://doi.org/10.2147/DMSO.S281529

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 2

Editor who approved publication: Dr Konstantinos Tziomalos



Yeongseon Ahn,1 Hyejin Lee,2 Yoon Shin Cho1

1Department of Biomedical Science, Hallym University, Chuncheon, Gangwon-do, Republic of Korea; 2Department of Internal Medicine, Ewha Womans University School of Medicine, Seoul, Republic of Korea

Correspondence: Yoon Shin Cho
Department of Biomedical Science, Hallym University, Chuncheon, Gangwon-do 24252, Republic of Korea
Tel +82-33-248-2111
Fax +82-33-256-3420
Email [email protected]

Purpose: Observational studies have demonstrated an increased risk of polycystic ovarian syndrome (PCOS) in obese women. This study aimed to identify genetic variants influencing obesity in females and to evaluate the causal association between genetically defined obesity and PCOS in Korean women.
Methods: Two-stage GWAS was conducted to identify genetic variants influencing obesity traits (such as body mass index [BMI], waist–hip ratio [WHR], and waist circumference [WC]) in Korean women. Two-sample Mendelian randomization (MR) analysis was employed to evaluate the causal effect of variants as genetic instruments for female obesity on PCOS.
Results: Meta-analysis of 9953 females combining discovery (N = 4658) and replication (N = 5295) stages detected four (rs11162584, rs6760543, rs828104, rs56137030), six (rs139702234, rs2341967, rs73059848, rs5020945, rs550532151, rs61971548), and two genetic variants (rs7722169, rs7206790) suggesting a highly significant association (P < 1× 10− 6) with BMI, WHR, and WC, respectively. Of these, an intron variant rs56137030 in FTO achieved genome-wide significant association (P = 3.39× 10− 8) with BMI in females. Using variants for female obesity, their effect on PCOS in 946 cases and 976 controls was evaluated by MR analysis. MR results indicated no significant association between genetically defined obesity and PCOS in Korean women.
Conclusion: This study, for the first time, revealed genetic variants for female obesity in the Korean population and reported no causal association between genetically defined obesity and PCOS in Korean women.

Keywords: female obesity, polycystic ovarian syndrome, causal relation, genome-wide association study, Mendelian randomization

Introduction

Obesity is a medical condition in which excess body fat has accumulated to the extent that it may have a negative effect on health, leading to reduced life expectancy and/or increased health problems. Obesity acts as an important risk factor for various diseases, increasing the risk of high blood pressure, type 2 diabetes, and dyslipidemia in obese people more than in those with adequate weight.1 Body mass index (BMI) and the waist–hip ratio (WHR) or waist circumference (WC) are the commonly used anthropometric measurements to measure general obesity and abdominal (central) obesity, respectively. Environmental and genetic factors influence obesity. Numerous genetic studies have been conducted to gain insight into the genetic basis of obesity.2–7 Genome-wide association studies (GWAS) have identified a large number of loci for BMI,4 WHR,6,8 and WC.6

Obesity in women is associated with alterations in the reproductive cycle, including reduced fertility, as well as an increased risk of polycystic ovarian syndrome (PCOS) and oligo-ovulation or anovulation.9–13 Furthermore, the tendency toward menstrual and ovarian disturbances associated with obesity may predispose women to an increased risk of ovarian, breast, and endometrial cancer.14

PCOS is a common endocrine disorder among women of reproductive age. PCOS is characterized by a variety of phenotypes, such as menstrual dysfunction, hyperandrogenism, and metabolic abnormalities with ethnic differences.15 Accurate diagnosis is difficult because the diagnostic criteria are variable. Three diagnostic criteria have been proposed: the 1990 National Institutes of Health criteria, the 2003 European Society for Human Reproduction and Embryology/American Society for Reproductive Medicine criteria, and the 2006 Androgen Excess Society criteria. The exact pathogenesis of PCOS has not been fully elucidated. PCOS is a complex disorder that is influenced by environmental and genetic factors.16 Early diagnosis and treatment along with weight loss may reduce the risk of long-term complications, such as type 2 diabetes and heart disease.

Mendelian randomization uses genetic variants to determine whether an observational association between a risk factor and an outcome is consistent with a causal effect. Mendelian randomization relies on the natural, random assortment of genetic variants during meiosis yielding a random distribution of genetic variants in a population. Individuals are naturally assigned at birth to inherit a genetic variant that affects or does not affect a risk factor. The Mendelian randomization approach exploits the fact that genotype precedes life events and is therefore not affected by lifestyle, socioeconomic variables, or any factor that follows conception.17 The basic principles of Mendelian randomization can be understood through comparison with randomized clinical trials (RCTs). RCTs are costly and time-consuming and may be impractical to carry out, or there may not be an intervention to test a certain hypothesis, limiting the number of clinical questions that can be answered by an RCT.

In this study, we conducted a meta-analysis of GWASs to identify female-specific genetic variants associated with obesity traits, such as BMI, WHR, and WC, in Korean females. The identified variants for each obesity trait were used as genetic instruments for the subsequent Mendelian randomization analysis to assess the causal relationship between obesity and PCOS.

Methods

Subjects

Subjects for the discovery stage were recruited from the Korean Association Resource study (KARE) cohort, a population-based cohort of 8842 participants. Details of the participant recruitment criteria and the study design have been provided elsewhere.18–20 A total of 4658 females were present in the KARE study cohort (Supplementary Figure 1 and Table 1). Participants included in the studies used for the replication stage included unrelated Korean participants from the Health Examinee shared control study (HEXA) cohort,21 as well as the Cardiovascular disease association study (CAVAS) cohort (formerly RURAL cohort).19 The HEXA cohort consisted of 3695 subjects of which 2048 were females. The CAVAS cohort included two sub-studies of CAVAS1816 and CAVAS3665. The CAVAS1816 study included 1816 subjects, of which 957 were females. A total of 2290 female subjects were selected from the 3665 subjects in the CAVAS3667 study (Supplementary Figure 1 and Table 1). The study design and cohort characteristics of the CAVAS cohort have been described previously.19

Table 1 Clinical Characteristics of Female Subjects in Study Cohorts

Subjects for the PCOS GWAS consisting of 1000 PCOS cases and 1000 controls were recruited from the endocrinology and gynecology clinics at Ewha Woman’s University Hospital for the Mendelian randomization analysis between obesity traits and PCOS. Details of the participant recruitment criteria and the study design have been described previously.22

Genotyping, Imputation, and Quality Control

Single nucleotide polymorphism (SNP) genotyping was carried out for the KARE study subjects using an Affymetrix Genome-wide Human SNP array 5.0. Genotyping using an Affymetrix Genome-Wide Human SNP array 6.0 was carried out for the HEXA and CAVAS1816 study subjects. The CAVAS3667 study subjects were genotyped using the Illumina HumanOmni I-Quad vI array. Samples with a genotyping missing call rate >1% and heterozygosity >30% were excluded from the sample pool. Markers with a missing SNP call rate >5%, with a minor allele frequency (MAF) <0.01, and a Hardy–Weinberg equilibrium (HWE) test P-value <1 × 10−6 were eliminated. Details on genotyping quality control for the genotype data have been described previously.18–21

Subjects from the Ewha Woman’s University Hospital PCOS GWAS were genotyped using the HumanOmni-Quad v I array (Illumina). Samples with a high missing genotype call rate (>2%) or high heterozygosity (>30%) were excluded. Also, subjects whose computed average pairwise identity-by-state value were higher than that estimated from first-degree relatives of Korean sib-pair samples (>0.8) were excluded from analyses. Markers with a high missing gene call rate (>1%), low MAF (<0.05), and a significant deviation from HWE P-value <1 × 10−6 were excluded. The remaining 636,870 SNPs in 1922 samples (976 cases and 946 controls) were used for association in subsequent analyses.22

For the genotype data, including the KARE, HEXA, CAVAS1816, and CAVAS3667 studies, SNP imputation was performed to increase the coverage of common variants employing the minimac3 program with 1000 Genome Phase3 individuals as the imputation reference panel.23 For the genotype data of the Ewha Woman’s University Hospital PCOS study, SNP imputation was performed using the IMPUTE2 program with 1000 Genome phase3 individuals as the imputation reference panel.24 Imputed SNPs of poor imputation quality (Rsq <0.3 for minimac3 imputation and info score <0.5 for IMPUTE2 imputation) were excluded. In addition, imputed SNPs with missing gene call rates >1% (MAF < 0.01, and HWE test P-value <1 × 10−7) were omitted for subsequent analyses.

Phenotype Measurements

Based on the World Health Organization Asia-Pacific Perspective 2000, subjects with BMI ≥ 25 kg/m2 were grouped as obese while those with 18.5 ≤ BMI ≤ 22.9 kg/m2 were considered normal. The WHR diagnostic standard for obesity is 0.85 for women and 0.95 for men (Korean Society for the Study of Obesity). The WC diagnostic standard for obesity is 85 cm (34 inches) for women and 90 cm (36 inches) for men (Korean Society for the Study of Obesity).25

In this study, PCOS patients were diagnosed according to the Rotterdam criteria, which have been described elsewhere.22 Individuals with specific disorders, such as adult-onset congenital adrenal hyperplasia, hyperprolactinemia, and androgen-secreting neoplasia, were excluded from the study. Patients taking medications (eg, steroids, oral contraceptives, metformin, or thiazide diuretics) before starting the study were excluded. Among the regular-cycling volunteers, 1000 women were recruited to serve as the control group. None of the controls had a family history of diabetes or PCOS. Subjects were excluded if they had been on hormonal medication within 3 months of the evaluation or had used other drugs that could affect basal parameter status.

Statistical Analyses

The associations between the genetic variants and each obesity trait (BMI, WHR, and WC) were tested after adjusting for age and participant area. Single variant tests for the imputed genotype data of the four GWAS datasets (KARE, HEXA, CAVAS1816, and CAVAS3667 studies) were performed using the Score test in the rvtests software package.26 A meta-analysis using the METAL program27 was performed combining the four GWAS results based on the inverse-variance weighting method assuming fixed effects with Cochran’s Q test used to assess between-study heterogeneity.28

Genetic variants for the obesity traits identified from the meta-analysis were tested for their association with PCOS from the Ewha Woman’s University Hospital PCOS imputation dataset. Frequentist association tests were conducted for the PCOS case-normal control analysis in an additive model using the SNPTEST program.29

Functional Annotation of Associated Loci

For the functional annotation, obesity trait-associated variants were investigated for their overlap with Expression Quantitative Trait Loci (eQTLs) from the GTEx portal (https://www.gtexportal.org/home/).30 Expressed genes (eGenes) showing nominal association (P-value <10−4) between obesity GWAS variants and their expression in each tissue were identified in the GTEx dataset and assigned as the candidate genes. If a variant is not an eQTL, the gene in which a variant is located or closest to a variant is assigned as the candidate gene.

Mendelian Randomization Analysis

A two-sample Mendelian randomization (MR) analysis was performed to investigate the existence of a causal relationship between obesity and PCOS. Summary association results from the non-overlapping obesity and PCOS sets of individuals were used for the two-sample MR in this study. The odds of PCOS risk were divided by the β coefficients of the levels of obesity traits (BMI, WHR, or WC) to determine ratio estimates for each instrumental variable (IV) (here, genetic variants associated with obesity traits). The effects of the individual genetic instruments were combined using inverse-variance weighted (IVW) analysis, resulting in a weighted mean estimate of the risk of PCOS per 1-standard deviation increase in the levels of obesity traits (BMI, WHR, or WC). The MendelianRandomization package in R statistical software (https://CRAN.R-project.org/package=MendelianRandomization) was used to perform the two-sample MR analysis.31

Results

Identification of Genetic Variants for Obesity Traits in Females

We conducted a two-stage sex-stratified GWAS to discover the genetic loci for obesity traits in females. In stage 1, the discovery stage, SNPs across the whole genome were tested for their association with one of the obesity traits (BMI, WHR, or WC) in 4658 female subjects from the KARE study. SNPs selected from the stage 1 linear analysis (P-value <0.05) after adjusting for area and age were taken forward to stage 2, the replication stage. A total of 5295 subjects were included for the replication analysis of female subjects from the HEXA, CAVAS1816, and CAVAS3667 studies. A meta-analysis was conducted for SNPs that were validated for their association with obesity traits in the replication stage by combining the association results from the discovery and replication stages (Supplementary Figure 1).

Our meta-analysis in female subjects identified four SNPs showing evidence of a suggestive association (P-value <10−6) with BMI. Of those, three SNPs (rs11162584, rs6760543, and rs828104) were female-specific without showing an association with BMI in male subjects. One SNP rs56137030 in FTO showed a genome-wide significant association with BMI in females (P-value = 3.39 X 10−8) but it also showed a nominal association with BMI in males (P-value = 2.00 X 10−2). The meta-analysis for WHR and WC identified six and two suggestively associated SNPs (P-value <10−6), respectively. Of those, four WHR (rs139702234, rs73059848, rs550532151, and rs61971548) and two WC (rs7722169 and rs7206790) associated SNPs were female-specific (Table 2, Figure 1, and Supplementary Figure 2).

Table 2 Obesity Traits (BMI, WHR, and WC) Associated Loci in Females. Discovery Stage Was GWAS for Obesity Traits in Each Sex-Stratified Group of KARE Cohort. Overall Association Results (Pmeta) Were Obtained from Meta-Analyses Combining Discovery (KARE) and Replication (HEXA, CAVAS1816, CAVAS3667) Stages

Figure 1 Regional association plots of rs56137030 (A), rs61971548 (B), and rs7722169 (C), showing the evidence of association in females for BMI, WHR, and WC, respectively. In the top panel of each, each SNP is plotted as a circle along the chromosomal position. Association analysis results represented as −log10P for SNPs (y-axis on left) are shown in a genomic region 400 kb to either side of the lead SNP (shown as a purple diamond). Recombination rates (cM/Mb) within loci are estimated from 1000 Genomes Phase 3 ASN and indicated as blue lines (y-axis on right). The magnitude of pair-wise linkage disequilibrium (LD) between the lead SNP and other SNPs is demonstrated by color, ranging from high (red) to low (blue). In the bottom panels, the locations of known genes are indicated in the region. Genomic positions are based on GRCh37/hg19.

Investigation of Functional Relevance of GWAS Loci to Obesity Traits

If a variant responsible for a GWAS locus also affects gene expression, it is known that the relevant gene could be involved in the biological mechanism of disease pathogenesis.32 In this context, to gain insight into the functional relevance of genetic variants identified from our GWA meta-analysis to female obesity, we investigated their effects on gene expression from the GTEx portal.

SNP rs6760543 showing association for BMI locates between LOC105373352 and TMEM18. The GTEx eQTL dataset indicated that this variant was also associated with ALKAL2 expression in subcutaneous adipose tissue (P-value = 3.50 X 10−5). Another BMI associated variant rs828104 showed association with the expression of GAPVD1 in subcutaneous adipose tissue (P-value = 2.10 X 10−8) and PRPS1P2 in visceral adipose tissue (P-value = 8.50 X 10−5). SNP rs56137030 for BMI also was associated with the gene expression of FTO in skeletal muscle (P-value = 4.10 X 10−7) as well as IRX3 in pancreas (P-value = 4.30 X 10−6) (Table 3).

Table 3 eQTLs of GWAS Variants for Obesity Traits

Among variants for WHR or WC, only eQTL information of rs61971548 was available in the GTEx dataset. This variant showed an association with NID expression in thyroid (P-value = 4.70 X 10−6). For WHR or WC variants that had no significant eQTLs or no available eQTL information in the GTEx dataset, the candidate gene was assigned to be the one where a variant is located or closest to a variant (Table 3).

Evaluation of the Causal Role of Genetically Defined Obesity in PCOS

We conducted two-sample Mendelian randomization (MR) to assess the causal correlation between genetically defined obesity traits and PCOS. The key assumptions of two-sample MR are (1) the IV is causally related to the risk factor; (2) confounding factors of the association between risk factors and outcome should not be related to the IV; and (3) the IV only affects the outcome through its effect on the risk factors. In this study, IVs were genetic variants associated with a risk factor; the risk factor is one of the obesity traits (BMI, WHR, or WC), and the outcome is PCOS. Type 2 diabetes and related traits, such as fasting plasma glucose (FPG) and glycated hemoglobin (HbA1C), were considered confounding factors (Figure 2).

Figure 2 Diagram displaying the components of the Mendelian randomization. Genetic variants as the instrumental variables are associated with biomarker (or exposure), but not with confounding factors as well as with outcome disease. Biomarker is a modifiable risk factor for outcome disease.

Association analyses of individual genetic instruments (SNPs) for the risk factors (obesity traits) with the outcome (PCOS) as well as confounding factors (type 2 diabetes and related traits) demonstrated that 2 of 12 SNPs were associated with confounding factors (rs6760543, PT2D = 3.2 × 10−3, PFPG = 1.7 × 10−2, PHbA1C = 1.4 × 10−2; rs2341967, PT2D = 2.8 × 10−2) (Table 4). To fulfill the key assumptions, we conducted a two-sample MR analysis using nine SNPs after excluding two SNPs showing an association with confounding factors and one SNP (rs550532151) that was not available in the PCOS dataset (Table 4). Our MR results indicated no association between genetically defined obesity and PCOS in Korean women (Table 5 and Figure 3).

Table 4 Association of Individual Genetic Instruments for Obesity Traits with PCOS Risk

Table 5 Mendelian Randomization Results for Obesity Traits on PCOS Risk (Inverse-Variance Weighted)

Figure 3 The results of the Mendelian randomization (MR) analyses between obesity traits (all (A), BMI (B), WHR (C), WC (D), or visceral fat (E)) and PCOS. The x-axis shows the genetic association with exposure (obesity traits). The y-axis shows the genetic association with outcome (PCOS).

Discussion

Abnormalities in women’s health are caused by a variety of factors. An individual’s lifestyle, genetics, hormonal imbalances, and ethnicity all play a role in women’s health. Obesity in women is associated with alterations in the reproductive cycle with a reduction in fertility, as well as an increased risk of PCOS and oligo-ovulation or anovulation. As obesity is strongly influenced by genetic and environmental factors (heritability 40–70%), numerous genetic studies have been conducted in diverse ethnic groups. However, a sex-stratified genetic study for obesity has never been carried out specifically in an East Asian population. In this regard, we conducted a genome-wide association meta-analysis for obesity traits (such as BMI, WHR, and WC) in about 18,000 Korean female subjects.

The meta-analyses identified genetic variants suggesting associations (Pfemale-meta < 10−6) with BMI (four SNPs), WHR (six SNPs), and WC (two SNPs) (Table 2) in female subjects. Of the 12 associated SNPs with one of the obesity traits, two SNPs (rs56137030 and rs7206790) were located in the intron region of the known obesity gene, FTO.33 An association between SNP rs56137030 and BMI was also detected in males (Pmale-meta = 0.02) as well as in all subjects (Pall-meta = 2.58 × 10−8). The other FTO SNP, rs7206790, was female-specific showing an association with WC only in females (Pfemale-meta = 4.14 × 10−7) but not in males (Pmale-meta = 0.65). SNP rs56137030 further showed association with the FTO gene expression (PeQTL = 4.10 × 10−7) in subcutaneous adipose tissue from the GTEx dataset, implying that FTO is the causal gene for BMI (Table 3).

SNPs rs6760543 and rs828104 associated with BMI are located between LOC105373352 and TMEM18, and in the intron region of LOC107987127, respectively. The eQTL information available from the GTEx dataset demonstrated that rs6760543 and rs828104 showed association with the gene expression of ALKAL2 (PeQTL = 3.50 × 10−5) and GAPVD1 (PeQTL = 2.10 × 10−8) in subcutaneous adipose tissue, respectively (Table 3). The coding product of ALKAL2 (ALK And LTK Ligand 2) is a ligand for receptor tyrosine kinases LTK and ALK. It has been known that the stimulation of ALK signaling may be involved in regulation of cell proliferation and transformation.34 The encoded protein of GAPVD1 (GTPase Activating Protein And VPS9 Domains 1) has been known to act as a GTPase-activating protein (GAP) and a guanine nucleotide exchange factor (GEF) and participates in various processes such as endocytosis, insulin receptor internalization or LC2A4/GLUT4 trafficking.35 Based on this notion, it is postulated that ALKAL2 and GAPVD1 are likely functional genes for BMI.

SNP rs139702234 was associated with WHR and is located upstream of PPARG (Peroxisome Proliferator Activated Receptor Gamma). The eQTL information of this SNP was not available in the GTEx dataset. The encoded PPARG protein forms heterodimers with retinoid X receptors, and these heterodimers regulate transcription of various genes. PPARG is known to be involved in adipogenesis and has been implicated in the pathology of numerous diseases, including obesity, diabetes, atherosclerosis, and cancer.36 The association between SNP rs139702234 and WHR was only detected in females in this study. The biological functions of candidate genes of the remaining nine SNPs were not clear with regard to obesity traits (Table 3).

Reproductive disturbances are more common in obese women regardless of the PCOS diagnosis. Obesity is a common finding in women with PCOS and aggravates many of its reproductive and metabolic features. The relationship between PCOS and obesity is complex and not well understood.37 PCOS is the most common endocrine cause of infertility, but the risk of type 2 diabetes, hypertension, dyslipidemia, and cardiovascular disease is much higher than that in control women.38 PCOS is a complex disease for women, ranging from metabolic syndromes to reproductive abnormalities. To gain insight into the causal relationship between obesity and PCOS, we performed an MR analysis using SNPs detected in this study for their association with obesity traits in females.

MR is an analytical method that uses genetic variants as a tool to identify modifiable risk factors affecting individual health. In particular, two-sample MR is a method of estimating the causal association between risk factors and outcomes using two different research samples. Integrating multiple data provides an opportunity to significantly improve statistical power.39 Of the 12 SNPs detected in this study in female subjects, two SNPs (rs6760543 and rs2341967) were removed due to their association with type 2 diabetes, FPG, or HbA1C, which were considered confounding factors for obesity and PCOS (Table 4). As the effect of a genetic variant on the results is through a different pathway than through a risk factor, it is problematic for MR studies in violation of instrumental variables (IVs). SNP rs550532151 was also removed because this SNP was not available in the PCOS dataset. The results of the MR with nine SNPs are shown in Table 5. The MR results revealed no significant association between PCOS and any of the obesity traits (BMI, WHR, or WC). The combined MR results for all nine SNPs were also unrelated.

The obesity traits and PCOS were highly correlated when logistic regression was performed using the epidemiological data from the PCOS study. However, the results of the MR analysis using only genetic factors in this study did not reach significance. This result may be due to reverse causality, a form of confounding that is difficult to account for. It arises if the outcome or preclinical aspects of the disease that lead to the outcome affect the risk factor. People with symptoms of PCOS may participate in more exercise and have a better diet than those without symptoms. This would lead to a negative association between obesity (risk factor) and PCOS (outcome). Another interpretation of the negative association in our MR analyses is vertical pleiotropy in which one factor affects downstream outcomes.40 As one of the downstream consequences of insulin resistance is obesity, vertical pleiotropy may have confounded the causal relationship between obesity traits and PCOS. An MR analysis between genetically defined insulin dysregulation and PCOS would be valuable to pursue in the future. In addition, it may be assumed that the racial differences have been reflected in our negative MR results because Korean PCOS patients are less obese than Caucasian patients.

Ethics Approval and Consent to Participate

This study was approved by the Institutional Review Board (HIRB-2017-030) of Hallym University. Informed consent was obtained from all subjects of the study.

Acknowledgments

All authors acknowledged bioresources from the National Biobank of Korea and the Centers for Disease Control and Prevention, Republic of Korea (KBN-2017-025).

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Funding

This study was supported by the National Research Foundation of Korea (NRF) Grant (NRF-2020R1I1A2075302) and Hallym University Research Fund 2020 (HRF-202002-008).

Disclosure

The authors declare no conflicts of interest.

References

1. Kyrou I, Randeva HS, Tsigos C, et al. Clinical problems caused by obesity. Feingold KR, Anawalt B, Boyce A, editors. In Endotext [Internet]. South Dartmouth (MA): MDText. com, Inc.; 2000

2. Yengo L, Sidorenko J, Kemper KE, et al. Meta-analysis of genome-wide association studies for height and body mass index in approximately 700000 individuals of European ancestry. Hum Mol Genet. 2018;27(20):3641–3649. doi:10.1093/hmg/ddy271

3. Akiyama M, Okada Y, Kanai M, et al. Genome-wide association study identifies 112 new loci for body mass index in the Japanese population. Nat Genet. 2017;49(10):1458–1467. doi:10.1038/ng.3951

4. Locke AE, Kahali B, Berndt SI, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206.

5. Wen W, Cho YS, Zheng W, et al. Meta-analysis identifies common variants associated with body mass index in east Asians. Nat Genet. 2012;44(3):307–311.

6. Wen W, Kato N, Hwang JY, et al. Genome-wide association studies in East Asians identify new loci for waist-hip ratio and waist circumference. Sci Rep. 2016;6(1):17958. doi:10.1038/srep17958

7. Rask-Andersen M, Karlsson T, Ek WE, Johansson A. Genome-wide association study of body fat distribution identifies adiposity loci and sex-specific genetic effects. Nat Commun. 2019;10(1):339. doi:10.1038/s41467-018-08000-4

8. Heid IM, Jackson AU, Randall JC, et al. Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat Genet. 2010;42(11):949–960. doi:10.1038/ng.685

9. Moran LJ, Lombard CB, Lim S, Noakes M, Teede HJ. Polycystic ovary syndrome and weight management. Womens Health (Lond). 2010;6(2):271–283. doi:10.2217/WHE.09.89

10. Rich-Edwards JW, Goldman MB, Willett WC, et al. Adolescent body mass index and infertility caused by ovulatory disorder. Am J Obstet Gynecol. 1994;171(1):171–177. doi:10.1016/0002-9378(94)90465-0

11. Bhattacharya S, Campbell DM, Liston WA, Bhattacharya S. Effect of Body Mass Index on pregnancy outcomes in nulliparous women delivering singleton babies. BMC Public Health. 2007;7(1):168. doi:10.1186/1471-2458-7-168

12. Moran C, Arriaga M, Rodriguez G, Moran S. Obesity differentially affects phenotypes of polycystic ovary syndrome. Int J Endocrinol. 2012;2012:317241. doi:10.1155/2012/317241

13. Zeng X, Xie YJ, Liu YT, Long SL, Mo ZC. Polycystic ovarian syndrome: correlation between hyperandrogenism, insulin resistance and obesity. Clin Chim Acta. 2020;502:214–221. doi:10.1016/j.cca.2019.11.003

14. Templeton A. Obesity and women’s health. Facts Views Vis Obgyn. 2014;6(4):175–176.

15. Engmann L, Jin S, Sun F, et al. Racial and ethnic differences in the polycystic ovary syndrome metabolic phenotype. Am J Obstet Gynecol. 2017;216(5):493e491–493 e413. doi:10.1016/j.ajog.2017.01.003

16. Kosova G, Urbanek M. Genetics of the polycystic ovary syndrome. Mol Cell Endocrinol. 2013;373(1–2):29–38. doi:10.1016/j.mce.2012.10.009

17. McArdle PF, Whitcomb BW, Tanner K, Mitchell BD, Shuldiner AR, Parsa A. Association between bilirubin and cardiovascular disease risk factors: using Mendelian randomization to assess causal inference. BMC Cardiovasc Disord. 2012;12(1):16. doi:10.1186/1471-2261-12-16

18. Shin JH, Lee KM, Shin J, Kang KD, Nho CW, Cho YS. Genetic risk score combining six genetic variants associated with the cellular NRF2 expression levels correlates with type 2 diabetes in the human population. Genes Genomics. 2019;41(5):537–545. doi:10.1007/s13258-019-00791-0

19. Cho YS, Go MJ, Kim YJ, et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat Genet. 2009;41(5):527–534. doi:10.1038/ng.357

20. Kong S, Cho YS. Identification of female-specific genetic variants for metabolic syndrome and its component traits to improve the prediction of metabolic syndrome in females. BMC Med Genet. 2019;20(1):99. doi:10.1186/s12881-019-0830-y

21. Kim YJ, Go MJ, Hu C, et al. Large-scale genome-wide association studies in East Asians identify new genetic loci influencing metabolic traits. Nat Genet. 2011;43(10):990–995. doi:10.1038/ng.939

22. 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

23. Das S, Forer L, Schonherr S, et al. Next-generation genotype imputation service and methods. Nat Genet. 2016;48(10):1284–1287. doi:10.1038/ng.3656

24. Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012;44(8):955–959. doi:10.1038/ng.2354

25. Lee SY, Park HS, Kim DJ, et al. Appropriate waist circumference cutoff points for central obesity in Korean adults. Diabetes Res Clin Pract. 2007;75(1):72–80. doi:10.1016/j.diabres.2006.04.013

26. Zhan X, Hu Y, Li B, Abecasis GR, Liu DJ. RVTESTS: an efficient and comprehensive tool for rare variant association analysis using sequence data. Bioinformatics. 2016;32(9):1423–1426. doi:10.1093/bioinformatics/btw079

27. Willer CJ, Li Y, Abecasis GR. Metal: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–2191. doi:10.1093/bioinformatics/btq340

28. Ioannidis JP, Patsopoulos NA, Evangelou E. Heterogeneity in meta-analyses of genome-wide association investigations. PLoS One. 2007;2(9):e841. doi:10.1371/journal.pone.0000841

29. Marchini J, Howie B, Myers S, McVean G, Donnelly P. A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007;39(7):906–913. doi:10.1038/ng2088

30. GTEx Consortium. Genetic effects on gene expression across human tissues. Nature. 2017;550(7675):204–213.

31. Yavorska OO, Burgess S. Mendelian randomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol. 2017;46(6):1734–1739. doi:10.1093/ije/dyx034

32. Hormozdiari F, van de Bunt M, Segre AV, et al. Colocalization of GWAS and eQTL signals detects target genes. Am J Hum Genet. 2016;99(6):1245–1260. doi:10.1016/j.ajhg.2016.10.003

33. Frayling TM, Timpson NJ, Weedon MN, et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science. 2007;316(5826):889–894. doi:10.1126/science.1141634

34. Reshetnyak AV, Murray PB, Shi X, et al. Augmentor alpha and beta (FAM150) are ligands of the receptor tyrosine kinases ALK and LTK: hierarchy and specificity of ligand-receptor interactions. Proc Natl Acad Sci U S A. 2015;112(52):15862–15867. doi:10.1073/pnas.1520099112

35. Hunker CM, Galvis A, Kruk I, Giambini H, Veisaga ML, Barbieri MA. Rab5-activating protein 6, a novel endosomal protein with a role in endocytosis. Biochem Biophys Res Commun. 2006;340(3):967–975. doi:10.1016/j.bbrc.2005.12.099

36. Wan J, Xiong S, Chao S, et al. PPARgamma gene C161T substitution alters lipid profile in Chinese patients with coronary artery disease and type 2 diabetes mellitus. Cardiovasc Diabetol. 2010;9(1):13. doi:10.1186/1475-2840-9-13

37. Sam S. Obesity and polycystic ovary syndrome. Obes Manag. 2007;3(2):69–73. doi:10.1089/obe.2007.0019

38. Palomba S, Santagni S, Falbo A, La Sala GB. Complications and challenges associated with polycystic ovary syndrome: current perspectives. Int J Womens Health. 2015;7:745–763. doi:10.2147/IJWH.S70314

39. Pierce BL, Burgess S. Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. 2013;178(7):1177–1184. doi:10.1093/aje/kwt084

40. Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601. doi:10.1136/bmj.k601

Creative Commons License © 2020 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.