| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Genomic Medicine, Faculty of Medicine (M.G., L.H., N.G., A.B., S.H., M.I.M.) and Institute of Reproductive and Developmental Biology (M.G., S.F.), Imperial College, Hammersmith Campus, London W12 0NN, United Kingdom; Departments of Epidemiology and Public Health (U.S., M.-R.J.) and Reproductive Endocrinology (K.R., M.J.G., S.F.), Imperial College, St. Marys Campus, London W2 1PG, United Kingdom; Wellcome Trust Centre for Human Genetics (B.L.P., C.J.G., M.I.M.) and Oxford Centre for Diabetes, Endocrinology and Metabolism (B.L.P., A.B., C.J.G., M.I.M.), Churchill Hospital, Oxford OX3 7LJ, United Kingdom; Department of Obstetrics and Gynaecology (G.S.C.), Bloomsbury Site, University College, London WC1E 6BT, United Kingdom; and Departments of Clinical Chemistry (A.R., S.T.), Obstetrics and Gynecology (S.T., H.M., A.-L.H.), and Public Health Science and General Practice (M.-R.J., S.T., A.P.), University of Oulu, Oulu, Finland 90014
Address all correspondence and requests for reprints to: Prof. Mark McCarthy, Robert Turner Professor of Diabetes, Oxford Centre for Diabetes, Endocrinology & Metabolism, Churchill Hospital Site, Old Road, Headington, Oxford OX3 7LJ, United Kingdom. E-mail: Mark.mccarthy{at}drl.ox.ac.uk.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
PCOS displays strong familial aggregation (5), and the etiology is likely to be multifactorial, individual susceptibility being determined by the action of multiple genetic and environmental risk factors. A wide range of potential candidate genes have been studied in the search for PCOS-susceptibility variants, but few of these have been studied in multiple populations, and even fewer have shown convincing replication (6, 7).
Some of the strongest evidence presented to date has related to a pentanucleotide repeat 5' to the CYP11A gene. The CYP11A gene is a credible biological candidate because it encodes P450scc (also known as cytochrome P450scc or cholesterol side-chain cleavage), the enzyme catalyzing the initial and rate-limiting step in adrenal and ovarian androgen biosynthesis. Thecal cells cultured from PCOS subjects show increased androgen production (8) and, in most studies, increased CYP11A mRNA expression (9, 10, 11) compared with samples from control women. Maintenance of this cellular phenotype on long-term culture (12) strengthens the hypothesis that CYP11A overexpression is a stable property of PCOS theca cells that contributes to the excess androgen production that characterizes this condition.
In the search for genetic variants that might be associated with CYP11A overexpression, attention has focused on a pentanucleotide repeat approximately 500 bp 5' to the translation start site (D15S520). Gharani et al. (13) reported evidence for linkage to this region in 20 multiplex pedigrees segregating PCOS (LOD,
1.74), and in an analysis of 148 PCOS subjects and 59 controls, associations between pentanucleotide length (specifically absence of the 216-bp, four-repeat allele) and both symptomatic PCOS and raised testosterone levels were found. Subsequent studies have, to varying extents, corroborated these findings. Urbanek et al. (14) found nominal evidence for excess allele sharing in the region in an analysis of 28 multiplex sibships, as well as a nonsignificant reduction in transmission (47%) of the four-repeat allele to PCOS offspring in 163 parent-offspring trios. Diamanti-Kandarakis et al. (15) studied 170 Greek subjects (80 with PCOS) and found a significant increase in the frequency of four-repeat-lacking genotypes in PCOS subjects (26% vs. 13% in controls) and reproduced the previously observed association with testosterone. Most recently, Daneshmand et al. (10) also reported a reduction in the four-repeat allele frequency in PCOS (14% in 51 PCOS women vs. 44% in 280 controls), which was associated with a striking compensatory difference in the frequency of the nine-repeat allele (56% in PCOS women vs. 7% in controls). The only published study that has failed to produce any evidence for a relationship between pentanucleotide genotype and phenotype, a Spanish study of 92 hirsute women (only 34 of whom had PCOS) and 33 controls, was also the smallest and, consequently, likely to have been underpowered (16).
On the basis of these data, CYP11A is generally considered one of the few well-supported examples of a PCOS susceptibility effect (6, 7). However, the sample sizes included in the studies described have been rather modest. Increasing concern over the poor performance of association studies in the identification and characterization of complex trait susceptibility effects (17, 18, 19) prompted us to reevaluate the relationship between CYP11A promoter variation and PCOS disease status, polycystic ovary (PCO)-related symptoms, and serum testosterone levels.
In this study, we have extended genotyping of the CYP11A pentanucleotide, together with a neighboring tetranucleotide repeat element not previously studied, to almost 2500 individuals, which is around twice the combined size of the previous studies. Our analyses of these adequately sized and well-characterized populations indicate that these variants do not contribute significantly to observed phenotypic variation in either serum testosterone levels or PCOS disease status.
| Subjects and Methods |
|---|
|
|
|---|
Women with PCOS (n = 583) were ascertained from infertility and endocrine clinics, predominantly at St. Marys and the Middlesex Hospitals in London. The criteria for diagnosis of PCOS were as previously described (5, 20) and, in line with the revised (2003) consensus on the diagnostic criteria for this syndrome (21), rely on the combination of clinical symptoms, ultrasonographic examination, and biochemical data. The essential criterion for the diagnosis of PCOS in this study was the presence of polycystic ovaries on ultrasound in a patient presenting with anovulation (oligomenorrhea or amenorrhea) and/or hyperandrogenism. A diagnosis of hyperandrogenism required clinical (presence of hirsutism or acne) and/or biochemical [serum total testosterone > 0.78 ng/ml (2.7 nmol/liter)] features. Other potential endocrine and neoplastic causes of hyperandrogenemia were excluded by appropriate tests (21). From this wider patient cohort of diverse ethnic origin, all of whom were typed for the variants of interest, we restrict the analyses reported here to the 371 women with exclusively British/Irish origin (based on three generations of reported ethnicity and birthplace). For 141 of these cases (designated probands), DNA was available from both parents, such that parent-offspring trios could be generated; the remaining 230 patients without parental DNA are designated PCOS cases. Family relationships in the trios were confirmed by genotyping with a panel of highly polymorphic microsatellite markers (combined heterozygosity, >99%). Quantitative trait analyses were performed in all PCOS subjects (i.e. probands and cases combined) but excluded any woman taking oral contraceptive treatment, pregnant at the time of sampling, or with previous oophorectomy. A small number of women (4%, 14 of 371 women) who had received metformin were included in these analyses; their exclusion did not alter the findings.
As population control subjects for case-control analyses, we used DNA (n = 350) from a random UK population of blood donors available from the Centre for Applied Microbiology and Research (Salisbury, UK), used extensively for such purposes. All control subjects were of British UK origin. Metabolic and hormonal status was not known for these individuals. Use of population controls (rather than subjects known not to have PCOS) is generally associated with only a modest reduction in power, which can readily be compensated by increasing the number of control subjects typed (19).
Finnish subjects
For the cohort study, samples were available from the Northern Finland Birth Cohort of 1966, which is based on 12,058 live-born babies born in the northernmost two provinces of Finland arising from pregnancies with expected dates of delivery during calendar year 1966. Of 5926 cohort members (3088 women) who had undergone clinical review and DNA sampling at age 31, we genotyped 1589 for this study (22). These included 527 women who had reported symptoms of hirsutism and/or oligo/amenorrhea on the 31-year review questionnaire (symptomatic women) and 1062 randomly selected women (cohort controls) from the same cohort without such symptoms (generating an approximate 1:2 ratio of cases to controls, all matched for age given the birth cohort design). No woman who was receiving oral contraceptives at the time of metabolic sampling was included in either group. The cohort cases have been shown as a group to have biochemical features consistent with PCOS (e.g. hyperandrogenemia and insulin resistance) (22), but only a proportion have proven PCOS. Ovarian ultrasound scan (USS) examination of these symptomatic women is ongoing, with 67 of 178 so far examined having PCO status confirmed on USS. These 67 individuals are referred to as USS-confirmed cohort cases. For the quantitative trait analyses, which were the main emphasis of analysis in the Finnish subjects, the group size was reduced to 1435 (475 cases and 960 controls) by exclusion of 138 women who were pregnant at the time of the blood sampling and 16 in whom no testosterone measures were available. The single subject who was taking metformin at the time of recruitment is included.
Clinical features of these sample sets are provided in Table 1
. All clinical investigations were conducted in accordance with the guidelines in the Declaration of Helsinki, and the study was approved by the relevant ethics committees in UK and Finland. All subjects provided fully informed consent.
|
In the UK women with PCOS and all the Finnish subjects, fasting samples were used for the analysis of testosterone levels. For the UK subjects, testosterone was measured with an in-house RIA using ether extraction and dextran-coated charcoal separation. Genomic DNA was extracted from whole blood using the PUREGENE kit (Gentra Systems, Minneapolis, MN) as described previously (23). For the Finnish subjects, serum testosterone was determined by automated chemiluminescence system (Ciba-Corning ACS-180; Diamond Diagnostics, Holliston, MA).
Genotyping
The pentanucleotide D15S520 was genotyped by PCR amplification in a 15-µl volume containing 30 ng genomic DNA, 10 pmol of deoxynucleotide triphosphates, 1.0 mM MgCl2, 1U Taq polymerase (Bioline, London, UK), and 8 µM of each of the forward (FAM-labeled 5'-CCTGGGCAACAAGAATGAAA for the UK subjects and FAM-labeled 5'-GGTGAAACTGTGCCATTGC for the Finns) and reverse (5'-GCCCAAGAGAAGGTCAGAGC for the UK subjects and 5'-GTTTGGGGGAAATGAGGGGC for the Finns) primers. Details of PCR conditions are available from the authors. Sequence analysis of the CYP11A gene revealed a second (tetranucleotide) microsatellite lying in intron 1 approximately 18.4 kb upstream of exon 2. This tetranucleotide lies immediately upstream of a novel exon that is expressed as the 5' end of an alternate CYP11A transcript (Gaasenbeek, M., M. I. McCarthy, S. Halford, and S. Franks, unpublished observations). We elected to type this microsatellite, which lies approximately 1500 bp from D15S520, because sequence context suggested a potential functional role analogous to that of the pentanucleotide with respect to the usual promoter. This microsatellite was amplified in a separate reaction (primers: forward, HEX-labeled 5'-AGTAGCTGGGATTATAGACCC; reverse, 5'-GGAGACTGGTGAGGCTAAGT; 1.5 mM MgCl2). Resulting PCR products were diluted 1:4 in water, and 0.5 µl of each PCR reaction was mixed with formamide and 450 bp ROX size standards (Applied Biosystems, Warrington, UK), followed by separation and sizing on an ABI PRISM 3700 DNA sequencer (Applied Biosystems). Genotypes were derived using Genotyper software (Applied Biosystems). When possible, Mendel-checking methods were used to confirm genotyping accuracy. In addition, approximately 10% of samples were re-genotyped. From these data, we estimate that the overall genotyping error is less than 1%. For convenience, alleles are designated by prefix P or T, reflecting pentanucleotide or tetranucleotide locus, respectively, and a number reflecting the number of repeat units (e.g. P4 includes four repeats at the pentanucleotide). We did not genotype either D15S169 or D15S519 [markers typed in some of the previous studies (13, 14)] either because of distance from the gene (D15S169) or uncertain location (D15S519 does not map to the current genome assembly).
Statistical analysis
Linkage disequilibrium relationships between the two microsatellites were determined in UK and Finnish control subjects using the LDMAX option within the GOLD (Graphical Overview of Linkage Disequilibrium) package (24), which implements the expectation-maximization (EM) algorithm. Between-group comparisons of genotype and allele frequencies used standard contingency table methods implemented in STATXACT (Cytel, Cambridge, MA), which allows exact probabilities (or Monte-Carlo estimates thereof) to be generated. Haplotype frequencies were derived by standard maximum likelihood methods implemented in TRANSMIT (25), with group comparisons assessed by likelihood ratio testing and permutation methods (at least 1000 permutations). In the UK trios, family-based association tests (25, 26) were undertaken using TRANSMIT; alleles and haplotypes with a frequency less than 2% were pooled for their respective analyses. Quantitative trait analyses of serum testosterone levels were performed by ANOVA on log-transformed values, using SPSS (version 10; SPSS Inc., Chicago, IL) and SAS (version 8.2; SAS Institute, Inc., Cary, NC), with and without adjustment for age. For analyses of testosterone levels, all women taking oral contraceptive medication at the time of sampling were excluded, as were women who were pregnant or had undergone oophorectomy.
Power calculations
The UK samples (case control + trios combined) provide approximately 85% power (at P = 0.05) to detect effect sizes (odds ratios for PCOS,
1.41.5) significantly smaller than those implied by some earlier reports (odds ratio, >2) (13, 15). The largest group studied for the relationship between CYP11A variation and testosterone levels (the Finnish cohort controls) provided more than 90% power to detect (at P = 0.05) differences between P4-containing and non-P4-containing genotypes of approximately 0.33 of a SD. Previous studies of the pentanucleotide (13, 15) have suggested effect sizes at least twice this.
| Results |
|---|
|
|
|---|
58% in controls). The 169-bp (P4) allele is equivalent to the 216-bp allele designated in some previous papers, the difference reflecting use of differing PCR primers (13). Seven alleles were observed at the tetranucleotide ranging from 276 bp (T6) to 304 bp (T13), with the 284-bp (T8) modal (
82% in controls). In the Finnish controls, allele distributions were similar, with the modal alleles again being P4 (69%) and T8 (81%).
In the UK population controls, the two microsatellites were in strong linkage disequilibrium (D',
0.86; P < 0.001); the most marked departures from equilibrium were positive associations between the P4 and T8 alleles (56.5% observed vs. 47.6% expected) and between P6 and T12 (11.3% vs. 3.1%) and the complementary negative associations between P6 and T8 (10.8% vs. 21.7%) and P4 and T12 (0.4% vs. 6.8%). In the Finnish controls, the extent and direction of linkage disequilibrium relationships between alleles at the two loci were similar (D',
0.87; P < 0.0001).
At the pentanucleotide, comparisons between UK PCOS cases and population controls revealed significant differences in genotype frequencies (Monte-Carlo estimate, P = 0.029; 99% confidence interval, 0.0270.030), which were largely attributable to an excess of P4 homozygotes in the cases (47.2% vs. 34.8% in controls; Table 2
). The same was true of allele frequencies (exact P = 0.005, reflecting P4 allele frequencies of 67.6% in cases and 58.2% in controls; Table 2
). At the tetranucleotide, there were no significant differences in genotype frequencies (P = 0.139), but allele frequency comparisons indicated an excess of T8 alleles in cases (P = 0.0097; 88.4% vs. 82.4% in controls). However, global comparison of haplotype frequencies (by maximum likelihood testing) revealed no significant differences between cases and controls (exact P = 0.15), although, as expected given the allele frequencies, the P4-T8 haplotype was more frequent in cases (65.5% vs. 56.5% in controls, P = 0.005 uncorrected).
|
2 = 4.15, df = 4, P = 0.39; P4 transmissions,
2 = 0.25, df = 1, P = 0.62), for the T8 or any other allele at the tetranucleotide (global analysis,
2 = 5.26, df = 5, P = 0.39; T8 transmissions,
2 = 0.007, df = 1, P = 0.93), or for estimates of haplotype transmission (global analysis,
2 = 4.41, df = 5, P = 0.49; P4-T8 haplotype,
2 = 0.017, df = 1, P = 0.90).
We also analyzed the relationship between microsatellite genotypes and testosterone levels. In the UK dataset, intermediate trait information was unavailable for the population controls, so these analyses were performed in all PCOS subjects (i.e. a combination of the cases and the trio probands) not receiving oral contraception. There was no association between genotype at either locus and testosterone levels (pentanucleotide, P = 0.98, Table 3
; tetranucleotide, P = 0.87). In particular, given previous findings with regard to the P4 allele (13), we compared testosterone levels in subjects with one or more P4 allele [P4/, equivalent to 216+ in previous papers (13, 15)] with all other subjects and found no appreciable difference [P4/: geometric mean, 0.61 ng/ml (2.10 nmol/liter); SD range, 0.380.96 ng/ml (1.333.32 nmol/liter); non-P4: geometric mean, 0.55 ng/ml (1.92 nmol/liter); SD range, 0.370.82 ng/ml (1.292.85 nmol/liter); P = 0.29). Reanalysis after adjusting for subject age had no significant impact on these findings.
|
| Discussion |
|---|
|
|
|---|
These findings are clearly in conflict with several previous studies (10, 14, 15) that have appeared to corroborate the initial findings (13) that variation at the CYP11A pentanucleotide is both linked to and (with specific reference to the absence of the four-repeat allele) associated with both symptomatic PCOS and raised testosterone levels.
We can discount several possible explanations for these discrepant findings. First, we are confident that our failure to detect association was not due to genotyping error; the low rate of discrepancies observed on cryptic duplicate genotyping and the consistent pattern of linkage disequilibrium relationships between the microsatellites across multiple datasets provide strong reassurance in this regard. Second, we can be confident, given the size of the resources typed, that the failure to detect associations with effect sizes similar to those previously reported was not due to lack of power. Third, the explanation is unlikely to lie in between-cohort differences in linkage disequilibrium relationships (for example, between the pentanucleotide and some as yet unidentified etiological variant); the close similarities in haplotype structure observed between the Finnish and UK data sets suggest that there are no profound differences in linkage disequilibrium structure, at least among populations of European origin (in whom the majority of the previous studies were conducted).
Our data suggest that the strength of, indeed the existence of, any association between CYP11A promoter variation and androgen-related phenotypes has been overestimated in previous studies through a combination of small sample size and data exploration conducted without due allowance for the associated inflation of the nominal type 1 error rate. This may have been compounded by the effects of any publication bias that has favored studies generating positive results. These factors are increasingly being recognized as contributing to the inconsistency that has plagued association studies, particularly in the context of multifactorial traits (17, 18, 19).
A further explanation for the discrepant findings may lie in the possibility of ethnic stratification (26). In this article, we have restricted the reported analyses to individuals of European origin. However, in PCOS cases of non-European origin, we have found marked differences in microsatellite allele frequencies. For example, in both Afro-Caribbean and South Asian subjects, the P6, rather than the P4, allele was the commonest allele at the pentanucleotide repeat, with frequencies of 46.6% and 48.6%, respectively. This suggests that studies of this locus may be particularly susceptible to spurious associations arising from the failure to match cases and control subjects adequately for ethnic composition.
Concerns over the rather poor performance of association studies have led to promulgation of recommendations designed to improve study design (17, 18, 19, 27, 28). Several of these have been incorporated in the current study, including: use of well-sized data sets; attempts to seek replication of any nominally significant findings, particularly those arising from data exploration; use of both case-control and family-based association methods; and analyses based on both intermediate phenotypes and discrete disease categories.
It is important to point out, of course, that the present study does not allow us to conclude that these variants have absolutely no effect on PCOS status or on testosterone levels. However, the size of the study means that any such effect must be extremely small and certainly much smaller than the effect size implied by the previous positive findings. We are currently undertaking functional studies of these microsatellites to determine whether they have any impact on CYP11A expression. Equally, the present study does not allow us to exclude the possibility that CYP11A variants other than those studied contribute to PCOS susceptibility and population variation in testosterone levels. Extensive resequencing and genotyping will be needed to test this hypothesis. However, as discussed above, it seems unlikely that the previously reported associations at the pentanucleotide repeat were a reflection of linkage disequilibrium with etiological variants elsewhere in the gene. Otherwise, given the apparent similarities in haplotype structure, at least as judged from the relationship between the microsatellites in these two quite distinct European populations, we should have expected to detect the same effects in this study.
The aim of this study was to test whether the previously reported associations with the CYP11A pentanucleotide held up to reanalysis in larger datasets. We conclude that they do not and that the major sites of genomic variation underlying susceptibility to PCOS and population variation in testosterone levels remain to be identified.
| Acknowledgments |
|---|
| Footnotes |
|---|
Abbreviations: PCO, Polycystic ovary; PCOS, PCO syndrome; USS, ultrasound scan.
Received September 18, 2003.
Accepted February 1, 2004.
| References |
|---|
|
|
|---|
-hydroxylase/C(1720) lyase promoters. Fertil Steril 77:274280[CrossRef][Medline]
This article has been cited by other articles:
![]() |
M. Simoni, C.B. Tempfer, B. Destenaves, and B.C.J.M. Fauser Functional genetic polymorphisms and female reproductive disorders: Part I: polycystic ovary syndrome and ovarian response Hum. Reprod. Update, September 1, 2008; 14(5): 459 - 484. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. C. Sakoda, C. Blackston, J. A. Doherty, R. M. Ray, M. G. Lin, H. Stalsberg, D. L. Gao, Z. Feng, D. B. Thomas, and C. Chen Polymorphisms in Steroid Hormone Biosynthesis Genes and Risk of Breast Cancer and Fibrocystic Breast Conditions in Chinese Women Cancer Epidemiol. Biomarkers Prev., May 1, 2008; 17(5): 1066 - 1073. [Abstract] [Full Text] [PDF] |
||||
![]() |
The ESHRE Capri Workshop Group Genetic aspects of female reproduction Hum. Reprod. Update, April 2, 2008; (2008) dmn009v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. O. Goodarzi, M. R. Jones, H. J. Antoine, M. Pall, Y.-D. I. Chen, and R. Azziz Nonreplication of the Type 5 17 -Hydroxysteroid Dehydrogenase Gene Association with Polycystic Ovary Syndrome J. Clin. Endocrinol. Metab., January 1, 2008; 93(1): 300 - 303. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. H. Olson, E. V. Bandera, and I. Orlow Variants in Estrogen Biosynthesis Genes, Sex Steroid Hormone Levels, and Endometrial Cancer: A HuGE Review Am. J. Epidemiol., February 1, 2007; 165(3): 235 - 245. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. W. Setiawan, I. Cheng, D. O. Stram, E. Giorgi, M. C. Pike, D. Van Den Berg, L. Pooler, N. P. Burtt, L. Le Marchand, D. Altshuler, et al. A Systematic Assessment of Common Genetic Variation in CYP11A and Risk of Breast Cancer Cancer Res., December 15, 2006; 66(24): 12019 - 12025. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Vink, S. Sadrzadeh, C. B. Lambalk, and D. I. Boomsma Heritability of Polycystic Ovary Syndrome in a Dutch Twin-Family Study J. Clin. Endocrinol. Metab., June 1, 2006; 91(6): 2100 - 2104. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Xita and A. Tsatsoulis Fetal Programming of Polycystic Ovary Syndrome by Androgen Excess: Evidence from Experimental, Clinical, and Genetic Association Studies J. Clin. Endocrinol. Metab., May 1, 2006; 91(5): 1660 - 1666. [Abstract] [Full Text] [PDF] |
||||
![]() |
D.H. Abbott, D.K. Barnett, C.M. Bruns, and D.A. Dumesic Androgen excess fetal programming of female reproduction: a developmental aetiology for polycystic ovary syndrome? Hum. Reprod. Update, July 1, 2005; 11(4): 357 - 374. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. F. Escobar-Morreale, M. Luque-Ramirez, and J. L. San Millan The Molecular-Genetic Basis of Functional Hyperandrogenism and the Polycystic Ovary Syndrome Endocr. Rev., April 1, 2005; 26(2): 251 - 282. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Ehrmann Polycystic Ovary Syndrome N. Engl. J. Med., March 24, 2005; 352(12): 1223 - 1236. [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 |