| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Obstetrics and Gynecology (R.S.L.), Pennsylvania State University, Hershey, Pennsylvania 17033; Department of Obstetrics and Gynecology and Duke Clinical Research Institute (H.X.B., E.R.M.), Duke University Medical Center, Durham, North Carolina; University of Colorado (W.D.S.), Denver, Colorado; University of Texas Southwestern Medical Center (B.R.C.), Dallas, Texas; Departments of Obstetrics and Gynecology (M.P.D.), Wayne State University, Detroit, Michigan; Baylor College of Medicine (S.A.C.), Houston, Texas; University of Alabama (M.P.S.), Birmingham, Alabama; University of Pennsylvania School of Medicine (C.C., K.G.E., R.S.S.), Philadelphia, Pennsylvania; University of Medicine and Dentistry of New Jersey (P.G.M.), Newark, New Jersey; Stanford University (N.A.C.), Stanford, California; University of Pittsburgh (G.G.G.), Pittsburgh, Pennsylvania; Department of Medicine (J.E.N.), Virginia Commonwealth University School of Medicine, Richmond, Virginia; Department of Obstetrics and Gynecology (L.C.G.), University of California at San Francisco, San Francisco, California; and Reproductive Sciences Branch (P.C.L.), National Institute of Child Health and Human Development, Bethesda, Maryland
Address all correspondence and requests for reprints to: Richard S. Legro, M.D., Department of Obstetrics and Gynecology, Pennsylvania State College of Medicine, M.S. Hershey Medical Center, 500 University Drive, H103, Hershey, Pennsylvania 17033. E-mail: RSL1{at}psu.edu.
| Abstract |
|---|
|
|
|---|
Objective: Our objective was to identify predictive genetic polymorphisms and other determinants of ovulatory response.
Design: This was a substudy of a multicenter randomized clinical trial.
Setting: This study was performed at academic medical centers and their affiliates.
Participants: A total of 312 women with PCOS were included in the study.
Main Outcome Measures: Historical, biometric, biochemical, and genetic parameters were performed.
Results: We found that the C allele of a single nucleotide polymorphism in the STK11 gene (expressed in liver; also known as LKB1) was associated with a significantly decreased chance of ovulation in PCOS women treated with metformin. In an analysis of ovulation per cycle, the adjusted odds ratio (OR) comparing the C/C genotype to the G/G genotype was 0.30 [95% confidence interval (CI) 0.14, 0.66], and the OR for the C/G genotype vs. the G/G genotype was also 0.30 (95% CI 0.14, 0.66). In an analysis of metformin-treated subjects, we found that the percentage of women who ovulated increased with the number of G alleles present: 48% (10 of 21) of C/C women, 67% (32 of 48) of C/G women, and 79% (15 of 19) of G/G women ovulated. We also found that increased frequency of ovulation was associated with lower body mass index (BMI) [adjusted OR of 2.36 (95% CI 1.65, 3.36) and 2.05 (95% CI 1.46, 2.88), respectively, for comparisons of BMI less than 30 vs. BMI equal to or more than 35, BMI 30–34 vs. BMI equal to or more than 35, in the analysis of ovulation per cycle], a lower free androgen index (FAI) [adjusted OR of 1.59 (95% CI 1.17, 2.18) for FAI < 10 vs. FAI
10], and a shorter duration of attempting conception [adjusted OR of 1.63 (95% CI 1.20, 2.21) for < 1.5 vs.
1.5 yr].
Conclusions: We have demonstrated that a polymorphism in STK11, a kinase gene expressed in liver and implicated in metformin action, is associated with ovulatory response to treatment with metformin alone in a prospective randomized trial. The interaction with the effects of changes in modifiable factors (e.g. BMI or FAI) requires further study.
| Introduction |
|---|
|
|
|---|
Restoring ovulation is pivotal to treating many of the patient complaints, including infertility, dysfunctional, or absent uterine bleeding, and a desire to prevent endometrial cancer. Improving ovulatory frequency may also benefit hyperandrogenic stigmata such as hirsutism, acne, and alopecia through normalization of circulating sex steroid levels.
We recently conducted a randomized double-blind trial [Pregnancy in Polycystic Ovary Syndrome (PPCOS)] in which we examined the effects of clomiphene, metformin, and the combination of the two on live birth in women with PCOS but also studied ovulation rates during the study (6). Despite treatment, a substantial proportion of women, 52 of 209 (24.9%) in the clomiphene group, 93 of 208 (44.7%) in the metformin group, and 35 of 209 (16.7%) in the combined group, experienced no documented ovulations over the course of the study.
In a pharmacogenetics substudy, we obtained DNA from participants, and sought to examine the predictive effect of genetic variation in key genes associated with drug response, alone or in combination with other factors that might modify drug action in women. Single nucleotide polymorphisms (SNPs) in selected candidate genes were tested for association with response to treatment with clomiphene and/or metformin based on previous findings of association with drug response or action: estrogen receptor 1 (ESR1, formerly known as ERa); two cytochrome P450 genes (CYP2C9 and CYP2D6); and a serine-threonine kinase gene expressed in the liver, STK11 (formerly known as LKB1).
Several studies have reported an association between response to hormone replacement therapies and the ESR1 PvuII polymorphism (IVS1–401, SNP rs2234693) (7, 8), and variation in two cytochrome P450 genes (CYP2C9 and CYP2D6) has been associated with varied metabolism (9), and response to selective estrogen receptor modulators (SERMs), including tamoxifen (10). Although no specific genetic variation has been associated with response to metformin, a recent study has shown that STK11 is required in the liver to control blood glucose levels (11).
STK11 phosphorylates and activates AMP-activated protein kinase, and is also a tumor suppressor gene implicated in the etiology of the Peutz-Jeghers syndrome (11). The deletion of STK11 (LKB1) in the liver of adult mice resulted in hyperglycemia with increased gluconeogenic and lipogenic gene expression, and its presence was required for metformin efficacy. We also tested a microsatellite marker, D19S884 located in an intron of fibrillin 3 (FBN3), which has been consistently associated with PCOS in family studies (12).
| Subjects and Methods |
|---|
|
|
|---|
Participants
The institutional review board at each center approved the protocol, and all participants gave written informed consent. A total of 312 women infertile women with PCOS participated in this substudy out of 626 who participated in the main randomized study. PCOS was defined as unexplained hyperandrogenic chronic anovulation, and women were randomly assigned to one of three treatment arms by an interactive voice system, and stratified based on study site and prior exposure to either of the study drugs (Fig. 1
) (15). Subjects with other causes of infertility, including male factor and tubal factor, were excluded.
|
We performed a double-blind, double-dummy study. Extended-release metformin (Glucophage XR) plus identical placebo were provided by Bristol-Myers Squibb (New York, NY). Overencapsulated clomiphene citrate tablets (Teva Pharmaceuticals, North Wales, PA) and matching placebo capsules were packaged and tested by a commercial pharmacy supply company (Clinical Trial Services, Inc., Research Triangle Park, NC), specifically for this study. Neither manufacturer had another role in the study.
Baseline laboratory testing was performed with subjects in the fasting state. All specimens were analyzed in a core laboratory, which was blinded to treatment allocation (15). Each subject received a monthly medication package consisting of bottle M (metformin in 500-mg tablets or matching placebo) and blister pack C (clomiphene in 50-mg tablets or matching placebo). Both drugs were begun concurrently. Patients gradually increased their dose of bottle M until reaching the maximum dose of four tablets (two tablets twice a day). Patients took 5 d of tablets from blister pack C beginning on d 3 of menses; this dose was maintained if adequate ovulation was documented. However, in non or poor responders, the dose was increased by one tablet a day on a treatment-cycle basis (either after 5-wk anovulation or after a menses until the maximum dose of three tablets per day was reached).
All patients had progesterone levels measured weekly to biweekly in local laboratories to document ovulation (15). If two consecutive elevated progesterone levels (>5 ng/ml) were noted, a weekly pregnancy test was administered until positive or menses occurred. After the baseline visit, subjects returned each month for a more extensive visit, including a limited physical examination, urine pregnancy test, and repeat fasting blood tests (15). A free androgen index (FAI) was calculated from the following formula: FAI = [total testosterone (nmol/liter)/SHBG (nmol/liter)] x 100. Subjects were treated for up to six cycles or 30 wk. All study medication was discontinued upon a positive pregnancy test. The present study was initiated midway through the PPCOS study, and DNA was obtained from blood independently of the study visit in the main PPCOS study.
Genotyping methods
DNA was prepared from whole blood samples using an AutoGen automated nucleic acid extraction system (AutoGen Inc., Holliston, MA) according to the manufacturers protocols. Genotypes for six SNPs and the microsatellite marker D19S884 described in Table 1
were obtained by methods described in Stewart et al. (12). D19S884 genotypes were determined using an Applied Biosystems 3130 capillary sequencer (Applied Biosystems, Foster City, CA) and Applied Biosystems GeneMapper software. We designate SNPs by gene name_SNP ID (e.g. STK11_rs8111699). The corresponding genotypes for these SNPs were tested for Hardy-Weinberg proportions in our sample, and none deviated significantly.
|
The primary outcome of the trial was the association between SNPs in candidate genes and ovulation rate. Secondary outcomes included the change in response parameters such as FAI, SHBG, insulin, and glucose, which at baseline were associated with an ovulatory response to treatment. All data entry, data management, and analyses were performed at the Data Coordinating Center at the Duke Clinical Research Institute. The study began enrolling patients in November 2002 and completed enrollment in December 2004.
For testing differences between subjects who participated and who did not participate in the genetic substudy,
2 or Fishers exact test was used for categorical variables, and Wilcoxon rank sum tests were used for continuous variables. For descriptive statistics, we examined proportion of ovulated cycles for each individual and compared the rate of at least one ovulated cycles, as well as the mean of the proportion of ovulated cycles between genotypes. The following modeling strategies were used. Generalized estimating equations were used for analysis of ovulation rate per cycle, in which each cycle is an observation unit, to account for correlation of multiple ovulation cycles within a subject. Logistic regression was used for analysis of ovulation rate per patient, in which each patient is an observation unit with an outcome of one if ever ovulated and zero if never ovulated over the course of study. Linear regression was used for modeling the change of response parameters, such as FAI, SHBG, insulin, and glucose. The following steps were used to derive our final models. To examine the interaction of genetic polymorphism with treatment, univariate analyses were performed first with both main and interaction effects of the polymorphism and the treatment group to select genes whose polymorphisms are associated at 0.1 level with the outcome within at least one treatment group. A regression selection procedure was then used to derive a final multivariate model that is adjusted for any potential baseline variables. The genes are retained in the final model if they are associated at the 0.1 level with the outcome within at least one treatment group after adjusting for the significant baseline variables. All analyses were performed with Statistical Analysis System version 8.2 (SAS Institute Inc., Cary, NC).
| Results |
|---|
|
|
|---|
A total of 312 subjects participated in the substudy, whereas 314 subjects did not participate (Fig. 1
). Among nonparticipants, the most common reason for nonparticipation was that that the subjects were unable to be contacted (47.5%; n = 149), followed by other reasons (28.0%; n = 88), not interested (15.3%; n = 48), and not available (11.8%; n = 37). There were no significant differences in treatment groups and baseline variables among subjects who did and did not participate in the substudy (Table 2
).
|
The rate of ovulation per cycle in the metformin group was significantly lower, 30.5% (n = 164 of 537) compared with the groups receiving clomiphene (54.8% or n = 164 of 537, in clomiphene alone and 60.2% or n = 343 of 570 in combined group) (P < 0.001 for both comparisons). Similar results are found for ovulation rate per patient (63.3% or n = 62 of 98 in metformin-only group compared with 85.4% or n = 88 of 103 in clomiphene alone and 86.5% or n = 96 of 111 in combined group; P < 0.001 for both comparisons). There was a significant advantage to combined therapy over clomiphene only. Over the course of the study, the mean (SD) number of ovulations per subject was 2.68 (1.78) in the clomiphene-only group, 1.67 (1.74) in the metformin-only group, and 3.09 (2.04) in the combination group (P < 0.001 for both comparisons).
Predictors of ovulation
We noted within the metformin-only group a significant negative association between ovulation per cycle and the C allele of STK11_rs8111699 compared with the G/G genotype [Table 3
; C/C genotype P = 0.011, and C/G genotype P = 0.002, respectively, for testing that the adjusted odds ratio (OR) = 1] in the analysis of ovulation per cycle. The adjusted OR is 0.30 for both comparisons, whereas the unadjusted OR is 0.41 for C/C vs. C/G and 0.49 for C/G vs. G/G. For descriptive statistics and unadjusted analyses, we also examined individual subjects in the metformin group by STK11 genotype; we note that 48% (10 of 21) of women in the metformin arm with genotype C/C ovulated, 67% (32 of 48) with genotype C/G ovulated, and 79% (15 of 19) with genotype G/G ovulated (Fig. 2
) (P = 0.05 for C/C vs. G/G genotype; and P = 0.33 for C/G vs. G/G genotype). The mean (SD) for the proportion of ovulated cycles is 0.26 (0.35), 0.31 (0.30), and 0.46 (0.34) for C/C, C/G, and G/G groups, respectively (P = 0.05 for C/C vs. G/G genotype; and P = 0.09 for C/G vs. G/G genotype). A similar negative association trend was noted in the clomiphene group compared with the G/G genotype (Table 3
; C/C genotype P = 0.071, and C/G genotype P = 0.12, respectively, for testing that the adjusted OR = 1.0). We also noted a significant positive association with the C/C genotype of ESRI_rs2234693 (the polymorphic Pvu site) compared with the T/T genotype (P = 0.034), but only in the combination therapy group. These trends were maintained when we examined the analysis results for ovulation per patient.
|
|
|
| Discussion |
|---|
|
|
|---|
The association between polymorphisms of STK11 and metformin action has not previously been studied in humans. Our study was conducted prospectively in a large cohort of participants, and examined candidate genes that had been reported to be associated with drug response in other pharmacogenetics studies, with the exception of the FBN3 and STK11, for which there were no prior reports. It focused on an outcome of clinical relevance (i.e. ovulation in women with PCOS) chosen a priori. In addition, we analyzed the association within a model that considered other factors associated with ovulatory failure because it is unlikely that any one common polymorphism acts independently of other environmental and constitutional factors to influence response to the medication. However, the sample size was small, so the power to detect modest effects was not large. The P values reported are nominal, not adjusted for multiple testing.
Our results suggest an association between STK11 genotypes and ovulatory response, with a stepwise increase in ovulation from C/C to C/G to G/G genotype in the metformin group. However, this association was not found in the other two study groups, though a similar trend in the clomiphene group suggests that the mechanism may not involve metformin metabolism. The mechanisms for this treatment-dependent response are beyond the range of this current study but may involve drug-drug interactions in the combined group. Furthermore, we a priori focused on ovulation per cycle rather than ovulation per patient. We note that the same negative association holds up when we examine the outcome as ovulation per patient. One bias toward the identification of nonresponse genotypes in this study is that successful ovulators were more likely to become pregnant and, therefore, to end participation in the treatment phase of the study.
We elected not to use conception as the study outcome because the sample size for this outcome would be significantly smaller than for ovulatory events. Restoring ovulation is a treatment goal for women with PCOS independent of their desire to conceive because it will control bleeding complaints, reduce the risk for endometrial abnormalities such as polyps or hyperplasia, and may lessen peripheral hyperandrogenic stigmata such as hirsutism. We were also interested in the response to these medications.
Other studies have noted an association between the PvuII polymorphism in the ESR1 gene and individual response to hormone replacement therapy in terms of cholesterol levels and other circulating cardiovascular risk parameters (8, 16, 17). We did not examine these parameters in our study but found no association between ovulation and this polymorphism, except a significant positive association with the C/C genotype of ESR1_rs2234693 in the combined group, but no effect was noted in the clomiphene group. Again drug-drug interactions may explain these inconsistent findings. Further study of the relationship between this gene and response to other SERMs, such as tamoxifen, a triphenylethylene derivative similar to clomiphene, and raloxifene, a benzothiophene derivative, that are more widely used for other indications such as breast cancer, is needed.
There is more comprehensive literature examining the association between SERM treatment response and polymorphisms of common SERM metabolizing genes. Tamoxifen is metabolized by CYP2D6 into a more active compound, endoxifen, and inhibition of this enzyme ameliorates tamoxifens benefit in treating women with breast cancer (18). This study also showed a significantly worse relapse-free time and disease free survival among 223 women with breast cancer and the CYP2D6 *4/*4 genotype treated with tamoxifen (10). The metabolism of clomiphene is complex and less well studied than that of tamoxifen, and the CYP450 system is also involved in the metabolism of sex steroids, which change in response to treatment with clomiphene, and this study is unable to explore these relationships.
We modeled the effect of genes and factors such as BMI and FAI to identify predictors of ovulation. The factors that we identified are very similar to those from another large prospective series that sought predictive factors of ovulatory response to clomiphene among 201 women with World Health Organization type II oligomenorrhea, a category mainly consisting of women with PCOS (19). These authors found that the FAI, BMI, cycle history (oligomenorrhea vs. amenorrhea), serum androgen (testosterone and/or androstenedione) levels, and mean ovarian volume assessed by transvaginal sonography were all significantly different (P < 0.01) in responders from those in nonresponders, and FAI was chosen to be the best predictor in univariate analysis. We similarly found that decreasing FAI was a strong predictor of ovulation along with decreasing BMI, a shorter duration of infertility, and a history of prior pregnancy loss (19). These latter two measures may reflect a less severe history of anovulation, analogous to the relative benefit of baseline oligomenorrhea vs. amenorrhea noted by Imani et al. (19). We did not find any clear change in either FAI or glycemic parameters related to an STK11 genotype to explain the association with anovulation in our study. We did not find circulating glucose, insulin, or androgen levels, other than as part of the FAI, or ovarian volume to be a predictor in our model. Differences in our population and diagnostic criteria between studies may also contribute to varying results.
Studies that have examined predictive factors on ovulatory response to metformin have been smaller with mixed results. Moghetti et al. (20), using a logistic regression model, found that higher insulin levels, lower serum androstenedione, and less severe oligomenorrhea in women with PCOS were associated with improved menstrual frequency in response to metformin. Another study found that insulin resistance, based on homeostatic markers, identified improved ovulatory response of metformin therapy (21). However, other studies have not confirmed these predictive factors (22).
We have demonstrated that a polymorphism of a liver-expressed kinase gene (STK11) implicated in metformin action is associated with a poor ovulatory response to metformin alone in a prospective randomized trial. In this study we performed statistical tests for several SNP genotypes at each of several genes, making the approach in part exploratory. We report nominal P values, recognizing that most of these would not remain significant after correction for multiple testing. Nevertheless, the most significant of these, for STK11, are very suggestive, in view of the prior plausibility of a role for this gene. The impact of this SNP needs to be verified in larger sample sizes in other prospective trials and also in patients receiving metformin for the treatment of type 2 diabetes.
| Acknowledgments |
|---|
We thank the Ligand Assay and Analysis Core Laboratory at the University of Virginia Center for Research and Reproduction, under the direction of D. Hasenleider, for their substantial contributions. We also thank Dr. Jerry Strauss at Virginia Commonwealth University for the suggestion to investigate the STK11 gene. Finally, we thank members of the Advisory Board: J. Hogan, F. Howard, M. Schiff, J. Wactawski-Wende, and N. Santoro (Chair); and express special gratitude for their careful and devoted oversight of our trial to the Data and Safety Monitoring Committee: E. Thom, J. Peipert, J. Zhang, P. Cato, C. Henderson, and R. Rebar (Chair).
| Footnotes |
|---|
This work was supported by National Institutes of Health/National Institute of Child Health and Human Development Grants U10 HD27049 (to C.C.), U01 HD38997 (to E.R.M.), U10 HD39005 (to M.P.D.), U10 HD27011 (to S.A.C.), U10 HD33172 (to M.P.S.), U10 HD38988 (to B.R.C.), U10 HD38992 (to R.S.L.), U10 HD38998 (to W.D.S.), U10 HD38999 (to P.G.M.), U54 HD34449 National Cooperative Program in Infertility Research (to R.S.S.), U54-HD29834 University of Virginia Center for Research in Reproduction Ligand Assay and Analysis Core, GCRC Grant MO1RR00056 to the University of Pittsburgh, and MO1RR10732, and construction grant C06 RR016499 to Pennsylvania State University. Glucophage XR and matching placebo were provided by Bristol-Myers Squibb.
Disclosure Summary: R.S.L. has served as a consultant to GlaxoSmithKline and Ferring, and has had lecture fees paid by Serono, meeting support from Abbott, and grant support from Pfizer. J.E.N. has equity ownership/stock options in Bristol Myers Squibb. P.C.L., E.R.M., and H.X.B. had grant support from Tap, H.X.B. also had grant support from Ortho Biotech, and E.R.M. has received research support from Merck and is also a consultant. W.D.S. has grant support from Organon and Wyeth. N.A.C. consults for Organon. P.G.M. has grant support from Ferring and Serono. M.P.D. has grant support from Serono, Tap, Glaxo Smith Klein, and Merck, and has served as a consultant for Tap and Serono. S.A.C., C.C., K.G.E., R.S.S., L.C.G., M.P.S., G.G.G., and B.R.C. have no disclosures.
Present address for P.C.L.: Department of Obstetrics and Gynecology, 100 Trent Drive, 212 Baker House, Duke University Medical Center, Durham, North Carolina 27710.
First Published Online November 13, 2007
1 For a list of members of the Reproductive Medicine Network, see Acknowledgments. ![]()
Abbreviations: BMI, Body mass index; ESR1, estrogen receptor 1; FAI, free androgen index; FBN3, fibrillin 3; OR, odds ratio; PCOS, polycystic ovary syndrome; PPCOS, pregnancy in PCOS; SERM, selective estrogen receptor modulator; SNP, single nucleotide polymorphism.
Received August 3, 2007.
Accepted November 7, 2007.
| References |
|---|
|
|
|---|
variants and the cardiovascular effects of hormone replacement therapy. Pharmacogenomics 4:269–277[CrossRef][Medline]This article has been cited by other articles:
![]() |
B. C. J. M. Fauser and M. J. C. Eijkemans Predicting Pregnancy in Women with Polycystic Ovary Syndrome J. Clin. Endocrinol. Metab., September 1, 2009; 94(9): 3183 - 3184. [Full Text] [PDF] |
||||
![]() |
S. Palomba, A. Falbo, F. Zullo, and F. Orio Jr. Evidence-Based and Potential Benefits of Metformin in the Polycystic Ovary Syndrome: A Comprehensive Review Endocr. Rev., February 1, 2009; 30(1): 1 - 50. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |