Molecular Vision 2011; 17:2297-2310 <http://www.molvis.org/molvis/v17/a250>
Received 28 April 2011 | Accepted 20 August 2011 | Published 25 August 2011

Polymorphisms in steroidogenesis genes, sex steroid levels, and high myopia in the Taiwanese population

Zoe Tzu-Yi Chen,1,2,3 I-Jong Wang,3,4 Ya-Tang Liao,5 Yung-Feng Shih,3 Luke Long-Kuang Lin3

1Department of Ophthalmology, Taipei City Hospital Zhongxing Branch, Taipei, Taiwan R.O.C; 2College of Medicine, National Taiwan University, Taipei, Taiwan R.O.C; 3Department of Ophthalmology, National Taiwan University Hospital, Taipei, Taiwan R.O.C; 4Graduate Institute of Clinical Medical Science, China Medical University, Taichung, Taiwan R.O.C; 5Genomics Research Center, Academia Sinica, Taiwan R.O.C

Correspondence to: I-Jong Wang, M.D., Ph.D., No. 7, Chung-Shan S. Rd., Taipei, Taiwan R.O.C.; Phone: 886-2-23123456 ext 2131; FAX: 886-2-23412875; email: ijong@ms8.hinet.net

Abstract

Purpose: To evaluate the relationship among single nucleotide polymorphisms (SNPs) in steroidogenesis enzyme genes, serum levels of sex steroids, and high myopia in Taiwanese male and female populations.

Methods: A campus-based sample of 283 cases (145 males and 138 females) with high myopia and 280 controls (144 males and 136 females) with low myopia or emmetropia was studied. Estradiol, progesterone, and testosterone levels were determined using enzyme-linked immunosorbent assay kits. We genotyped six SNPs within five steroidogenesis enzyme genes (17 alpha-hydroxylase/17,20 lyase [CYP17A1], 3 beta-hydroxysteroid dehydrogenase [HSD3B1], 17 beta-hydroxysteroid dehydrogenase 1 [HSD17B1], steroid-5-alpha-reductase, alpha polypeptide 2 [SRD5A2], and aromatase [CYP19A1]) using polymerase chain reaction–restriction fragment length polymorphism methods. Student’s t-tests, χ2 tests, logistic regression, multifactor dimensionality reduction (MDR) methods, and ANOVA were used to determine significance.

Results: An MDR analysis corroborated the synergistic genotype association and demonstrated that synergistic interaction between rs6203 (HSD3B1), rs10046 (CYP19A1), and sex might confer susceptibility to high myopia (p=0.019). In both male and female subjects, levels of testosterone were significantly higher in cases than in controls; in male subjects, the levels of estradiol were significantly higher and those of progesterone were significantly lower in cases (all p-values <0.001). The rs605059 (HSD17B1), with sex-gene interaction, showed association with estradiol levels in males (p=0.035) and testosterone levels in females (p=0.027).

Conclusions: Testosterone levels correlate with high myopia, and interaction of steroidogenesis enzyme genes and sex may be a modulating factor in sex hormone metabolism and high-myopia risk.

Introduction

Myopia is the most common eye disorder in Taiwan and around the world. A refractive error in excess of −6 diopters (D), defined as high myopia, is also termed “pathological myopia” due to its associated potential complications, including cataracts, glaucoma, macular degeneration, and retinal detachment, which can lead to blindness. Many environmental factors [1,2] have been reported to be associated with myopia. Population, family, and twin studies also have provided evidence for a genetic component to high myopia [3-5]. According to a nationwide survey in Taiwan in 2000, girls have a higher prevalence and more severe myopia progression than boys [6]. Several other epidemiological studies [7,8] have also shown similar results. In addition, the female sex was reported to be independently associated with faster myopia progression [9]. Because of sex differences, a hormonal hypothesis has been postulated by a review article examining myopia in opposite-sex twins [10]. In animal models, sex differences were also demonstrated in chick eye growth and experimental myopia [11]. It has also been shown that myopia progresses quickly in childhood, especially around puberty [6-8,12]. These studies have demonstrated that sex may be a modulating factor affecting myopia and have suggested that the differential expression of sex hormones may be as well. Knowledge about sex differences in the progression of myopia and about the pathophysiology underlying those differences may enhance the accuracy and effectiveness of clinical assessment and treatment of myopia. With a remarkable prevalence of high myopia and a homogenous genetic background, the Taiwanese population would be the best candidate for genetic studies to find myopia susceptibility genes, and approaches focusing on high myopia may minimize the impact of environmental confounders.

Thinning of the sclera is a key feature of human myopia development. Modulated by retinoscleral signals, the scleral remodeling mechanism is common [13-16] to different paradigms of myopia induction. It involves the regulation of numerous gene products (mRNAs and proteins) [17-23] in the extracellular matrix, such as collagens, proteoglycans, matrix metalloproteinases (MMPs), and tissue inhibitors of metalloproteinases (TIMMPs). Sex steroids are present in all tissues because they circulate through the blood; however, their effects can be observed only in cells armed with the corresponding receptors [24]. Studies have demonstrated the presence of androgen, estrogen, and/or progesterone receptors in various ocular tissues, such as the lacrimal gland, meibomian gland, lid, conjunctiva, cornea, iris/ciliary body, lens, retina/uvea, retina/choroid, and retinal pigment epithelial cells [24-28]. Nevertheless, both the levels and activities of these receptor proteins are regulated by hormones in the blood, suggesting that the sex hormone axis may play a role in the pathogenesis of certain ocular diseases.

Progesterone and estradiol are the major naturally occurring human hormones in the groups of progestogens and estrogens, respectively. Testosterone, a steroid hormone from the androgen group, is the principal male sex hormone. Steroids have been reported to change the biochemical characteristics of the iris/ciliary body, aqueous outflow pathway, and sclera in the rabbit eye [29]. In humans, estrogen is a modulating factor of the biomechanical properties of the cornea, which may contribute to the development of keratectasia after corneal refractive surgery [30]. Changes in cornea thickness due to fluctuations in estrogen levels are known to occur during the menstrual cycle [31] and pregnancy [32]. Corneal curvature also changes during pregnancy [33]. Besides, steroidogenesis enzymes and sex steroids have been proved to be related to ocular physiology and diseases [24,34-42]. Hormonal balance has been hypothesized to be related to the weakening of the collagenous tissue of the sclera [43]. However, no studies have been published to date regarding the influence of these sex hormones on the sclera or clinical myopia. Further research on this will be helpful to assess and manage myopia.

Sex steroids such as estrogen and dehydroepiandrosterone (DHEA, belonging to the androgen group) have been reported to regulate MMPs. Estrogen is known to upregulate MMP-2 and/or MMP-9 in mouse and human cells [26,44-46]. Depression of MMP-2 activity, on the other hand, has been shown in the rat polycystic ovary (PCO) induced by DHEA [47]. Because MMP-2 upregulation was suggested to take part in scleral remodeling in myopia development [17,23], and because sex hormones might implicate MMP-2 in the pathogenesis of myopia, our hypothesis is that genetic polymorphisms in steroidogenesis enzyme genes may affect circulating sex steroid levels and thus influence the risk of high myopia.

The progestogens, estrogens, and androgens are the three major classes of sex steroids and are derivatives of cholesterol (Figure 1; progestogens in orange boxes, estrogens in yellow boxes, and androgens in blue boxes). The enzymes involved in the biogenesis of sex hormones include 17 alpha-hydroxylase/17,20 lyase (CYP17A1), 3 beta-hydroxysteroid dehydrogenase (HSD3B1), 17 beta-hydroxysteroid dehydrogenase 1 (HSD17B1), steroid-5-alpha-reductase, alpha polypeptide 2 (SRD5A2), and aromatase (CYP19A1). Several notable nonsynonymous SNPs or common SNPs with functional effects have been reported in these genes, including rs743572 (CYP17A1), rs6203 (HSD3B1), rs1047303 (HSD3B1), rs605059 (HSD17B1), rs523349 (SRD5A2), and rs10046 (CYP19A1) [48-59]. rs743572 (CYP17A1) is located in the 5′ UTR and rs10046 (CYP19A1) is located in the 3′ UTR, the rest of these SNPs are all in exon areas.

Herein, we report the first case-control study to assess the relationship among genetic variations of steroidogenesis enzymes, sex hormone levels (including progesterone, estradiol, and testosterone), and high myopia.

Methods

Ethics statement

The study was approved by the Ethics Committees of the National Taiwan University Hospital, Taipei, Taiwan and followed the tenets of the Declaration of Helsinki. All participants were informed about the nature of the study and signed a written consent form before blood withdrawal.

Subjects and clinical assessment

A total of 283 unrelated high-myopia cases (145 males and 138 females) with refractive errors in excess of −6 D in both eyes and 280 controls (144 males and 136 females) with refractive errors of −2 D to 1 D in either eye were recruited from the National Taiwan University for this campus-based control study. No participant had ammetropia of more than 2 D or known ocular diseases/insults that could predispose to myopia. All subjects were of Han Chinese origin and were similar in social background. Ocular examination included objective autorefraction and keratometry (KR-8100; Topcon, Tokyo, Japan) and slit-lamp evaluation.

Measurement of hormone levels

Blood samples were centrifuged at 805× g for 15 min at 4 °C, and serum was separated and stored at −30 °C until analysis. Hormone levels were determined using an enzyme-linked immunosorbent assay (ELISA) kit (Diagnostic Systems Laboratories, San Antonio, TX) according to the manufacturer's instructions. The Steroid ELISA is based on the principle of competitive binding between steroids in the test specimen and steroid-horseradish peroxidase (HRP) conjugate for a constant amount of rabbit anti-steroid (progesterone, estradiol, and testosterone). In the incubation, goat anti-rabbit IgG-coated wells were incubated with standards (25 μl for progesterone and estradiol, 10 μl for testosterone), controls, subject samples (25 μl for progesterone and estradiol, 10 μl for testosterone), 100 μl steroid-HRP conjugate reagent, and 50 μl rabbit anti-steroid hormone reagent at room temperature (18–25 °C) for 90 min. During the incubation, a fixed amount of HRP-labeled steroid hormone competed with the endogenous steroid hormone in the standard, sample, or quality control serum for a fixed number of binding sites of the specific steroid hormone antibody. Thus, the amount of steroid hormone peroxidase conjugate immunologically bound to the well progressively decreased as the concentration of steroid hormone in the specimen increased. Unbound steroid hormone peroxidase conjugate was then removed, and the wells were washed. Next, a solution of H2O2/TMB reagent (100 μl TMB for progesterone and testosterone, 200 μl H2O2/TMB only for estradiol) was added and incubated at room temperature for 20 min, resulting in the development of blue color. The color development was stopped with the addition of stop solution (100 μl 1N HCl for progesterone and testosterone, 100 μl 3N HCl only for estradiol) and the absorbance was measured spectrophotometrically at 450 nm. Standards were assayed in duplicate, and unknowns were assayed in triplicate. The results were expressed as ng/ml or pg/ml.

SNP genotyping

Genomic DNA was extracted from participants’ peripheral leukocytes using an extraction kit (Qiagen GmbH, Hilden, Germany). To analyze polymorphisms in each gene, we used polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) methods to directly genotype the polymorphism. The primers for amplifying 50 ng of genomic DNA were as follows: 5′-GCT AGG GTA AGC AGC AAG AG-3′ and 5′-CAA GGT GAA GAT CAG GGT AG-3′ for rs743572 (CYP17A1), 5′-GGC TTC CTG CTG GAA ATA GTG A-3′ and 5′-GCA AAC CAC ATG GTC TTT CCT C-3′ for rs6203 (HSD3B1), 5′-CAG GCC AAT TTA CAC CTA TCG-3′ and 5′-TCA AAC TAT GTG AAG GAA TGG A-3′ for rs1047303 (HSD3B1), 5′-CCT TCC ACC GCT TCT ACC AAT A-3′ and 5′-CAG GGA CCA CAC AGA CCC AG-3′ for rs605059 (HSD17B1), 5′-AGG ATT TGG AGC TGC TGA GAG C-3′ and 5′-GAC GGT ACT TCT GGG CCT CTT CTA C-3′ for rs523349 (SRD5A2), and 5′-ACT ACT GAT GAG AAA TGC TCC AGA CT-3′ and 5′-GGT TCC TTA TAG GTA CTT TCA GCC A-3′ for rs10046 (CYP19A1). The PCR reactions were performed in 25-μl aliquots containing 20 pmole of each primer, 1× reaction buffer, 100 µM deoxynucleotide triphosphates, and 1 unit of Taq polymerase; a PCR system 9700 thermocycler (Applied Biosystems, Foster City, CA) was used for cycling. The temperatures used during PCR were as follows: for rs743572, 95 °C for 2 min, followed by 30 cycles of 95 °C for 1 min, 57 °C for 1 min, and 72 °C for 1 min, with a final extension at 72 °C for 10 min; for rs6203 and rs605059, 95 °C for 2 min, followed by 30 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 30 s, with a final extension at 72 °C for 10 min; for rs1047303, 95 °C for 2 min, followed by 40 cycles of 95 °C for 30 s, 50 °C for 30 s, and 72 °C for 30 s, with a final extension at 72 °C for 10 min; and for rs523349 and rs10046, 95 °C for 2 min, followed by 30 cycles of 95 °C for 30 s, 52 °C for 30 s, and 72 °C for 30 s, with a final extension at 72 °C for 10 min. Aliquots of PCR products were digested with the following restriction enzymes and conditions: MspA1I for rs743572 at 37 °C for 3 h, BglI for rs6203 at 37 °C overnight, DrdI for rs1047303 at 37 °C overnight, BtsI for rs605059 at 55 °C overnight, HpyCH4V for rs523349 at 37 °C overnight, and BsrCI for rs10046 at 65 °C overnight, in reaction mixtures containing 10 μl of PCR product, digestion buffer, BSA (BSA), and 5 units of restriction enzyme (New England Biolabs, Ipswich, MA). The digested fragments were visualized on 4% agarose gel with ethidium bromide staining to identify the base pair change. To exclude the possibility of introducing PCR error and miscutting of the restriction enzyme, all steps were repeated twice. For equivocal results, we double-checked with either repeating the experiments or directly sequencing the PCR products using dye terminator chemistry (BigDye Terminator version 3.1 on a model 3730 Genetic Analyzer; Applied Biosystems).

Statistical analyses

SAS version 9.1 (SAS, Inc., Cary, NC) was used for all statistical comparisons. Comparisons of age, refractive status, and hormone levels between the case and control groups were performed with Student’s t-tests. Genotypes were assessed for Hardy–Weinberg equilibrium (HWE), and allele frequencies were compared between cases and controls using χ2 tests. Odds ratios (ORs) for specific genotypes and 95% confidence intervals (CIs) were calculated using logistic regression analyses with the wild type genotypes as references. Synergistic analysis was evaluated by pairwise modification of genotypic associations through logistic regression. Global and genotypic-specific p-values were reported. The multifactor dimensionality reduction (MDR) method [60,61] with version 1.1.0 of the open-source MDR software package (Dartmouth Medical School, Hanover, NH) was used to search for the most potential gene-gene interaction models. MDR tests all possible multilocus genotype combinations and reports a combination that provides the best classification [62]. The final models with relatively highest prediction accuracy for single to four-way interactions are selected for cross-validation. This method, which is a nonparametric and model-free approach, combines cross-validation and permutation-testing procedures to minimize false positives caused by multiple testing and has performed well across many genetic simulation scenarios where purely epistatic relationships existed between status and variables without main effects [63]. True-positive models possess estimated testing accuracy of 50% or more cross-validation consistency. ANOVA (ANOVA) including six polymorphisms and sex was performed with the general linear models (GLM) procedure to determine the sources of variation for hormone levels, and F tests computed from the type III sums of squares reported the effect of an independent variable after adjustment for all other independent variables included in the model. Two-way interaction of each polymorphism and sex was also analyzed. Another ANOVA compared between-group hormone levels for different genotypes of steroidogenesis enzyme genes. Data from males and females were also analyzed separately.

Results

Subjects

Our study consisted of 283 high-myopia cases (145 males and 138 females; average age 18.8±0.8 years [mean±standard deviation]; average spherical equivalent −8.2±1.7 D for right eyes, −8.0±1.8 D for left eyes) and 280 emmetropia or low-myopia controls (144 males and 136 femals; average age 19.4±2.4 years [mean±standard deviation]; average spherical equivalent −0.8±0.8 D for right eyes, −0.7±0.8 D for left eyes; Table 1).

Sex steroid levels

The level of progesterone was significantly lower and levels of estradiol and testosterone were significantly higher in the cases than in the controls for male subjects (all p-values <0.001; Table 1). The level of testosterone was also significantly higher in the cases than in the controls for female subjects (p<0.001; Table 1).

Genotype association

Although no statistically significant deviation from HWE was noted for the controls, the genotypes of all SNPs between the cases and controls were not statistically different (all p-values >0.05; Table 2). Pairwise stratified analyses conducted separately among male and female subjects for all six SNPs (rs743572 [CYP17A1], rs6203 [HSD3B1], rs1047303 [HSD3B1], rs605059 [HSD17B1], rs523349 [SRD5A2], and rs10046 [CYP19A1] were indicated as SNP1–6, respectively) were used to evaluate synergistic modification of association with logistic regression. All the significant results are shown in Table 3. Among male study subjects, SNP6-SNP2 and SNP3-SNP4 combinations showed significant modification of genotypic risks for myopia in both global and specific tests. In the SNP6-AG subgroup, the TC genotype in SNP2 was associated with a high risk of myopia (OR=2.38 [1.17–4.83], p-value=0.017); in the SNP4-GA subgroup, SNP3-AC was associated with a higher risk (OR=4.32 [1.37–13.61], p-value=0.013). Meanwhile, these combinations were related to a protective effect among female subjects (OR=0.40 [0.19–0.81], p-value=0.011 and OR=0.57 [0.21–1.50], p-value=0.253, respectively). Among female subjects, SNP2-TC was associated with a significant risk (OR=2.88 [1.05–7.88], p-value=0.04) in the SNP6-AA subgroup, and SNP2-CC was associated with a protective effect (OR=0.26 [0.08–0.88], p-value=0.03) in the SNP6-AG subgroup; however, there were no significant association among male subjects (OR=0.96 [0.39–2.37], p-value=0.929 and OR=0.42 [0.12–1.46], p-value=0.173, respectively). Table 4 summarizes the best interaction models of seven independent variables (including these six SNPs and sex) in relation to high myopia obtained from the MDR analysis. Consistent with synergistic genotypic association, among all two-factor combinations, the model composed of rs6203 (HSD3B1) and rs10046 (CYP19A1) was the best attribute for predicting high-myopia risk, though without statistical significance (testing accuracy, 50.6%; p=0.422;cross-validation consistency, 7/10; p=0.429). The best interaction model of all combinations was a three-factor model composed of rs6203 (HSD3B1), rs10046 (CYP19A1), and sex (testing accuracy, 57.5%; p=0.019; cross-validation consistency, 10/10; p=0.07).

For all hormones tested, sex as an independent variable was significantly associated with variation, with p<0.001 in all cases (Table 5). The results of ANOVA that included a gene-sex interaction term as an independent variable found that the rs605059 (HSD17B1) genotype-sex interaction term was significant (Table 5). In the male subjects, genetic variants of rs605059 (HSD17B1) showed between-group differences in estradiol levels (p=0.035; Table 6). In the female subjects, genetic variants of rs605059 (HSD17B1) showed between-group differences in testosterone levels (p=0.027; Table 7).

Our results reveal that synergistic interaction of rs6203 (HSD3B1), rs10046 (CYP19A1), and sex may confer susceptibility to high myopia; testosterone levels for either sex correlate significantly with high myopia; and rs605059 (HSD17B1) interacting with sex-related factors may be associated with circulating sex hormone levels in high myopia.

Discussion

To our knowledge, this is the first study that describes the relationship between polymorphisms in steroidogenesis enzyme genes and high-myopia risk and their correlations with serum sex hormone levels. The principle findings of this study are the identification of (1) the combination of rs6203 (HSD3B1), rs10046 (CYP19A1), and sex may confer susceptibility to high myopia; (2) significant differences in all progesterone, estradiol, and testosterone levels between high-myopia cases and controls in males and significant differences in testosterone levels between high-myopia cases and controls in females; (3) an association between rs605059 (HSD17B1) and serum estradiol levels in males and serum testosterone levels in females; and (4) sexual dimorphism of these associations, which suggests a sex-gene interaction.

SNP information was obtained from reviewing papers [48-59] and from GenBank. We validated the presence of six common SNPs in our study. A comparison with four populations from the HapMap database (743572, 6203, 1047303, 605059, 523349, and 10046; accessed October 1, 2009) revealed significant differences in allele frequencies between Taiwanese controls and both Utah, United States, residents with European ancestry (CEU) and the Yoruba population in Ibadan, Nigeria (YRI; p<0.001; Table 8). The between-population differences noted here may help explain the well known differences in serum hormone levels and the prevalence of hormone-related diseases between Asian and Caucasian populations [64,65]. A lack of deviation from HWE among controls (Table 2) indicated that population stratification or genotyping errors are not likely. As the first study investigating an association between polymorphisms in steroidogenesis enzyme genes and high-myopia risk, this study did not observe statistically significant associations in any single variant. However, it is conceivable that particular combinations of multiple variants may confer an enhanced susceptibility to high myopia. Both stratified analysis of genotypic association (parametric) and MDR analysis (non-parametric) confirmed a synergistic interaction among rs6203 (HSD3B1), rs10046 (CYP19A1), and sex. The HSD3B gene family has two genes and five pseudogenes, all of which map to chromosome 1p13 [66-68]. The HSD3B1 gene encodes the type I enzyme, which is exclusively expressed in the placenta and peripheral tissues, such as prostate, breast, and skin. The HSD3B2 gene encodes the type II enzyme, which is predominantly expressed in classical steroidogenic tissues, namely the adrenals, testis, and ovary [68-71]. 3β-hydrosteroid dehydrogenases (HSD3B1) catalyze several reactions in the androgen pathway (Figure 1), leading to production of androstenedione and testosterone. They also catalyze the production of progesterone from pregnenolone. Cytochrome P-450 19A1 (CYP19A1), located on chromosome 15, codes for aromatase, which is involved in the rate-limiting step in estrogen production. It is mainly expressed in the ovarian granulose cells; however, aromatase is also expressed in several other tissues, including fat, placenta, breast, central nervous system, skin, and bone. Rosmond et al. [51] proposed that genetic variation in HSD3B1 might lead to an elevation in plasma aldosterone (one of the mineralocorticoids, which are downstream products of progestogens), with a resultant increase in intravascular volume and hypertension, and they did prove the rs6203 silent substitution in exon 4 of HSD3B1 was associated with elevated blood pressure in a population of Swedish men. Thus, it is possible that genetic variations in these steroidogenesis enzyme genes might also influence sex steroids and the pathogenesis of steroid-related diseases. Our results also support the associations between genetic variations of steroidogenesis enzyme genes and high myopia.

This current study showed that high-myopia cases have significantly higher levels of testosterone in both sexes, suggesting that higher testosterone levels for either sex are associated with an increased risk of high myopia. Additionally, in males, serum levels of estradiol are significantly higher and serum levels of progesterone are significantly lower in the high-myopia group. The possibility that steroid levels are associated with an increased rate of axial elongation requires further elucidation by a longitudinal study. Previously, Balacco et al. [72] reported higher levels of steroids (testosterone) in high-myopia cases (in excess of −10 D). Plasma cortisol levels have been shown to be lower in myopia cases than in normal subjects [73], although the results were not highly significant and the study size was small. In another case-control study, no significant differences were noted in serum levels of cortisol, testosterone, or estradiol between high myopes and non-myopes (p>0.1) [74], but the study size was also small. The number of subjects for our study is relatively large, and the results are significant (all p-values <0.001). Modulated by retinoscleral signals, the scleral remodeling mechanism is intrinsic to myopia [13]. The existence of sex steroid receptors in the retina [24,26,28] makes an effect of sex hormones on retinoscleral signals plausible.

Estrogen has been shown to upregulate MMP-2 and/or MMP-9 in human ocular cells [26,45], and MMP-2 upregulation was suggested as part of the scleral remodeling process in myopia development [17,23]. The finding of higher estradiol levels in male high-myopia cases is consistent with our hypothesis that differential sex hormone levels may play a role in myopia pathogenesis by regulating MMP-2.

Because steroidogenesis enzymes catalyze the biosynthesis of sex hormones, we examined whether polymorphisms of these enzymes influence the serum levels of sex hormones. Our data suggested that for males, rs605059 (HSD17B1) A allele is associated with a significant increase in the level of estradiol and a non-significant increase in the level of testosterone; for females, it is associated with a non-significant decrease in the level of estradiol and a significant decrease in the level of testosterone. This sexual dimorphism could be the result of a sex-gene interaction we demonstrated by ANOVA analysis (Table 5), although these results of female subjects should be interpreted with caution because of limitations to the present association study. The effect of polymorphisms in estrogen-metabolizing genes on circulating estrogen levels would be most pronounced in women with no other sources of estrogens because peripheral conversion of androgens in the adipose tissue or hormone replacement therapy (HRT) may mask the genetic effect. In addition, fluctuations in progesterone and estradiol levels during the menstrual cycle or pregnancy may obscure a modest genetic effect. Hence, a future study with a larger sample size and more considerable lifestyle information will help match controls to cases regarding various confounders—such as age, body mass index (BMI), pregnancy or menopausal status, HRT usage, day of period, time of day, and fasting status at blood draw—and determine these associations in females. Evidence proves that sexual dimorphism of the muscarinic system exists in the human brain and bladder [75-79]. In rat models, sex steroids selectively modulate muscarinic receptors in certain brain tissues (e.g., hypothalamus, adenohypophysis, superior cervical ganglia, or anterior pituitary gland) in a region-specific manner but not in the retina [80-85]. The muscarinic receptors and hormone receptors might be in different steps of a common pathway in myopia. Previous studies that have examined the rs605059 (HSD17B1) variant in relation to estrogens found no association [86-88]. The enzymes 17 beta-hydroxysteroid dehydrogenase 1 (HSD17B1) and 17 beta-hydroxysteroid dehydrogenase 2 (HSD17B2) are responsible for the reduction and oxidation of estrogens to more- and less-active forms, respectively. HSD17B1 is responsible for the reduction of estrone to estradiol, the more bioactive estrogen; HSD17B2 is responsible for the oxidation of estradiol to estrone. In addition to activating estrone to estradiol, human HSD17B1 enhances androgen activity in vivo [89]. The functional significance of rs605059 (HSD17B1) is unclear, but it appears to have little effect on the catalytic or immunological properties of this enzyme [90] and is indicated by sequence homology analysis (though with low confidence) that it is unlikely to affect function [91].

We observed associations between certain combinations of polymorphisms at steroidogenesis enzyme genes (rs6203 [HSD3B1], rs10046 [CYP19A1], and sex) and high-myopia risk; meanwhile, there is evidence of correlations between rs605059 (HSD17B1)–sex interaction and sex hormone levels. Sex hormones may play only a modulating rather than mediating role in the formation of high myopia; thus, the genotype associations with hormone levels and high-myopia traits are not parallel, and the small genetic effect of genotype-disease association cannot be detected unless a stratified synergistic association analysis is applied. The complex underlying genetic mechanism, especially in complex diseases, warrants consideration of environment-gene and gene-gene interactions. Our study did clearly reveal that a sex-gene interaction plays a significant role in sex hormone levels and high-myopia traits as well. A larger sample size and more comprehensive confounder controlling will further help elucidate the associations.

In conclusion, these findings reflect the unique genetic characteristics of the Han Chinese population. We found significant correlations between sex steroid levels and high myopia, and the synergistic interaction of sex and steroidogenesis enzyme genes might confer susceptibility to high myopia and influence sex steroid levels. Further studies are required to confirm this relationship. Sex-gene interaction is crucial for understanding sex hormone metabolism and possible high-myopia risk.

Acknowledgments

This study was supported by the National Science Council, Taipei, Taiwan (NSC95–2314-B002–153-MY2) and the Department of Health, Taipei City Government, Taipei, Taiwan (96002–62–093).

References

  1. Wong TY, Foster PJ, Hee J, Ng TP, Tielsch JM, Chew SJ, Johnson GJ, Seah SK. Prevalence and risk factors for refractive errors in adult Chinese in Singapore. Invest Ophthalmol Vis Sci. 2000; 41:2486-94. [PMID: 10937558]
  2. Mutti DO, Mitchell GL, Moeschberger ML, Jones LA, Zadnik K. Parental myopia, near work, school achievement, and children's refractive error. Invest Ophthalmol Vis Sci. 2002; 43:3633-40. [PMID: 12454029]
  3. Goldschmidt E. The importance of heredity and environment in the etiology of low myopia. Acta Ophthalmol (Copenh). 1981; 59:759-62. [PMID: 7315228]
  4. Edmund J, Goldschmidt E. Concordance for myopia and discordance for optic disk cupping in a pair of monozygotic twins. Klin Monatsbl Augenheilkd. 1988; 193:645-6. [PMID: 3221661]
  5. Pacella R, McLellan J, Grice K, Del Bono EA, Wiggs JL, Gwiazda JE. Role of genetic factors in the etiology of juvenile-onset myopia based on a longitudinal study of refractive error. Optom Vis Sci. 1999; 76:381-6. [PMID: 10416932]
  6. Lin LL, Shih YF, Hsiao CK, Chen CJ, Lee LA, Hung PT. Epidemiologic study of the prevalence and severity of myopia among schoolchildren in Taiwan in 2000. J Formos Med Assoc. 2001; 100:684-91. [PMID: 11760374]
  7. Zhao J, Pan X, Sui R, Munoz SR, Sperduto RD, Ellwein LB. Refractive Error Study in Children: results from Shunyi District, China. Am J Ophthalmol. 2000; 129:427-35. [PMID: 10764849]
  8. Goh PP, Abqariyah Y, Pokharel GP, Ellwein LB. Refractive error and visual impairment in school-age children in Gombak District, Malaysia. Ophthalmology. 2005; 112:678-85. [PMID: 15808262]
  9. Hyman L, Gwiazda J, Hussein M, Norton TT, Wang Y, Marsh-Tootle W, Everett D. Relationship of age, sex, and ethnicity with myopia progression and axial elongation in the correction of myopia evaluation trial. Arch Ophthalmol. 2005; 123:977-87. [PMID: 16009841]
  10. Miller EM. Reported myopia in opposite sex twins: a hormonal hypothesis. Optom Vis Sci. 1995; 72:34-6. [PMID: 7731654]
  11. Zhu X, Lin T, Stone RA, Laties AM. Sex differences in chick eye growth and experimental myopia. Exp Eye Res. 1995; 61:173-9. [PMID: 7556481]
  12. Saw SM, Tong L, Chua WH, Chia KS, Koh D, Tan DT, Katz J. Incidence and progression of myopia in Singaporean school children. Invest Ophthalmol Vis Sci. 2005; 46:51-7. [PMID: 15623754]
  13. McBrien NA, Gentle A. Role of the sclera in the development and pathological complications of myopia. Prog Retin Eye Res. 2003; 22:307-38. [PMID: 12852489]
  14. Phillips JR, Khalaj M, McBrien NA. Induced myopia associated with increased scleral creep in chick and tree shrew eyes. Invest Ophthalmol Vis Sci. 2000; 41:2028-34. [PMID: 10892839]
  15. Rada JA, Shelton S, Norton TT. The sclera and myopia. Exp Eye Res. 2006; 82:185-200. [PMID: 16202407]
  16. Siegwart JT, , Jr Norton TT. Regulation of the mechanical properties of tree shrew sclera by the visual environment. Vision Res. 1999; 39:387-407. [PMID: 10326144]
  17. Guggenheim JA, McBrien NA. Form-deprivation myopia induces activation of scleral matrix metalloproteinase-2 in tree shrew. Invest Ophthalmol Vis Sci. 1996; 37:1380-95. [PMID: 8641841]
  18. Jobling AI, Nguyen M, Gentle A, McBrien NA. Isoform-specific changes in scleral transforming growth factor-beta expression and the regulation of collagen synthesis during myopia progression. J Biol Chem. 2004; 279:18121-6. [PMID: 14752095]
  19. McBrien NA, Lawlor P, Gentle A. Scleral remodeling during the development of and recovery from axial myopia in the tree shrew. Invest Ophthalmol Vis Sci. 2000; 41:3713-9. [PMID: 11053267]
  20. Moring AG, Baker JR, Norton TT. Modulation of glycosaminoglycan levels in tree shrew sclera during lens-induced myopia development and recovery. Invest Ophthalmol Vis Sci. 2007; 48:2947-56. [PMID: 17591859]
  21. Rada JA, Nickla DL, Troilo D. Decreased proteoglycan synthesis associated with form deprivation myopia in mature primate eyes. Invest Ophthalmol Vis Sci. 2000; 41:2050-8. [PMID: 10892842]
  22. Siegwart JT, , Jr Norton TT. The time course of changes in mRNA levels in tree shrew sclera during induced myopia and recovery. Invest Ophthalmol Vis Sci. 2002; 43:2067-75. [PMID: 12091398]
  23. Siegwart JT, , Jr Norton TT. Selective regulation of MMP and TIMP mRNA levels in tree shrew sclera during minus lens compensation and recovery. Invest Ophthalmol Vis Sci. 2005; 46:3484-92. [PMID: 16186323]
  24. Wickham LA, Gao J, Toda I, Rocha EM, Ono M, Sullivan DA. Identification of androgen, estrogen and progesterone receptor mRNAs in the eye. Acta Ophthalmol Scand. 2000; 78:146-53. [PMID: 10794246]
  25. Paimela T, Ryhanen T, Mannermaa E, Ojala J, Kalesnykas G, Salminen A, Kaarniranta K. The effect of 17beta-estradiol on IL-6 secretion and NF-kappaB DNA-binding activity in human retinal pigment epithelial cells. Immunol Lett. 2007; 110:139-44. [PMID: 17532054]
  26. Marin-Castaño ME, Elliot SJ, Potier M, Karl M, Striker LJ, Striker GE, Csaky KG, Cousins SW. Regulation of estrogen receptors and MMP-2 expression by estrogens in human retinal pigment epithelium. Invest Ophthalmol Vis Sci. 2003; 44:50-9. [PMID: 12506055]
  27. Suzuki T, Kinoshita Y, Tachibana M, Matsushima Y, Kobayashi Y, Adachi W, Sotozono C, Kinoshita S. Expression of sex steroid hormone receptors in human cornea. Curr Eye Res. 2001; 22:28-33. [PMID: 11402376]
  28. Ogueta SB, Schwartz SD, Yamashita CK, Farber DB. Estrogen receptor in the human eye: influence of gender and age on gene expression. Invest Ophthalmol Vis Sci. 1999; 40:1906-11. [PMID: 10440242]
  29. Knepper PA, Collins JA, Frederick R. Effects of dexamethasone, progesterone, and testosterone on IOP and GAGs in the rabbit eye. Invest Ophthalmol Vis Sci. 1985; 26:1093-100. [PMID: 3926712]
  30. Spoerl E, Zubaty V, Raiskup-Wolf F, Pillunat LE. Oestrogen-induced changes in biomechanics in the cornea as a possible reason for keratectasia. Br J Ophthalmol. 2007; 91:1547-50. [PMID: 17591666]
  31. Giuffrè G, Di RL, Fiorino F, Bubella DM, Lodato G. Variations in central corneal thickness during the menstrual cycle in women. Cornea. 2007; 26:144-6. [PMID: 17251801]
  32. Weinreb RN, Lu A, Beeson C. Maternal corneal thickness during pregnancy. Am J Ophthalmol. 1988; 105:258-60. [PMID: 3344782]
  33. Park SB, Lindahl KJ, Temnycky GO, Aquavella JV. The effect of pregnancy on corneal curvature. CLAO J. 1992; 18:256-9. [PMID: 1424062]
  34. Hamad A, Kluk M, Fox J, Park M, Turner JE. The effects of aromatase inhibitors and selective estrogen receptor modulators on eye development in the zebrafish (Danio rerio). Curr Eye Res. 2007; 32:819-27. [PMID: 17963101]
  35. Salyer DL, Lund TD, Fleming DE, Lephart ED, Horvath TL. Sexual dimorphism and aromatase in the rat retina. Brain Res Dev Brain Res. 2001; 126:131-6. [PMID: 11172896]
  36. Suzuki T, Sullivan DA. Estrogen stimulation of proinflammatory cytokine and matrix metalloproteinase gene expression in human corneal epithelial cells. Cornea. 2005; 24:1004-9. [PMID: 16227852]
  37. Yu X, Tang Y, Li F, Frank MB, Huang H, Dozmorov I, Zhu Y, Centola M, Cao W. Protection against hydrogen peroxide-induced cell death in cultured human retinal pigment epithelial cells by 17beta-estradiol: a differential gene expression profile. Mech Ageing Dev. 2005; 126:1135-45. [PMID: 16029884]
  38. Smith JA, Vitale S, Reed GF, Grieshaber SA, Goodman LA, Vanderhoof VH, Calis KA, Nelson LM. Dry eye signs and symptoms in women with premature ovarian failure. Arch Ophthalmol. 2004; 122:151-6. [PMID: 14769589]
  39. Bigsby RM, Cardenas H, Caperell-Grant A, Grubbs CJ. Protective effects of estrogen in a rat model of age-related cataracts. Proc Natl Acad Sci USA. 1999; 96:9328-32. [PMID: 10430942]
  40. Sator MO, Akramian J, Joura EA, Nessmann A, Wedrich A, Gruber D, Metka M, Huber JC. Reduction of intraocular pressure in a glaucoma patient undergoing hormone replacement therapy. Maturitas. 1998; 29:93-5. [PMID: 9643522]
  41. Haffner SM, Klein R, Dunn JF, Moss SE, Klein BE. Increased testosterone in type I diabetic subjects with severe retinopathy. Ophthalmology. 1990; 97:1270-4. [PMID: 2147054]
  42. Gupta PD, Kalariya N, Nagpal K, Vasavada A. Pinteraction of sex steroid hormones with the eye. Cell Mol Biol (Noisy-le-grand). 2002; 48 Online Pub:OL379-86. [PMID: 12643455]
  43. Balacco-Gabrieli C. The etiopathogenesis of degenerative myopia. Ann Ophthalmol. 1983; 15:312-4. [PMID: 6229208]
  44. Potier M, Elliot SJ, Tack I, Lenz O, Striker GE, Striker LJ, Karl M. Expression and regulation of estrogen receptors in mesangial cells: influence on matrix metalloproteinase-9. J Am Soc Nephrol. 2001; 12:241-51. [PMID: 11158214]
  45. Suzuki T, Sullivan DA. Estrogen stimulation of proinflammatory cytokine and matrix metalloproteinase gene expression in human corneal epithelial cells. Cornea. 2005; 24:1004-9. [PMID: 16227852]
  46. Puistola U, Westerlund A, Kauppila A, Turpeenniemi-Hujanen T. Regulation of 72-kd type IV collagenase-matrix metalloproteinase-2 by estradiol and gonadotropin-releasing hormone agonist in human granulosa-lutein cells. Fertil Steril. 1995; 64:81-7. [PMID: 7789584]
  47. Henmi H, Endo T, Nagasawa K, Hayashi T, Chida M, Akutagawa N, Iwasaki M, Kitajima Y, Kiya T, Nishikawa A, Manase K, Kudo R. Lysyl oxidase and MMP-2 expression in dehydroepiandrosterone-induced polycystic ovary in rats. Biol Reprod. 2001; 64:157-62. [PMID: 11133670]
  48. Feigelson HS, Shames LS, Pike MC, Coetzee GA, Stanczyk FZ, Henderson BE. Cytochrome P450c17alpha gene (CYP17) polymorphism is associated with serum estrogen and progesterone concentrations. Cancer Res. 1998; 58:585-7. [PMID: 9485002]
  49. Zmuda JM, Cauley JA, Kuller LH, Ferrell RE. A common promotor variant in the cytochrome P450c17alpha (CYP17) gene is associated with bioavailability testosterone levels and bone size in men. J Bone Miner Res. 2001; 16:911-7. [PMID: 11341336]
  50. Chang BL, Zheng SL, Hawkins GA, Isaacs SD, Wiley KE, Turner A, Carpten JD, Bleecker ER, Walsh PC, Trent JM, Meyers DA, Isaacs WB, Xu J. Joint effect of HSD3B1 and HSD3B2 genes is associated with hereditary and sporadic prostate cancer susceptibility. Cancer Res. 2002; 62:1784-9. [PMID: 11912155]
  51. Rosmond R, Chagnon M, Bouchard C, Bjorntorp P. Polymorphism in exon 4 of the human 3 beta-hydroxysteroid dehydrogenase type I gene (HSD3B1) and blood pressure. Biochem Biophys Res Commun. 2002; 293:629-32. [PMID: 12054649]
  52. Kraft P, Pharoah P, Chanock SJ, Albanes D, Kolonel LN, Hayes RB, Altshuler D, Andriole G, Berg C, Boeing H, Burtt NP, Bueno-de-Mesquita B, Calle EE, Cann H, Canzian F, Chen YC, Crawford DE, Dunning AM, Feigelson HS, Freedman ML, Gaziano JM, Giovannucci E, Gonzalez CA, Haiman CA, Hallmans G, Henderson BE, Hirschhorn JN, Hunter DJ, Kaaks R, Key T, Le Marchand L, Ma J, Overvad K, Palli D, Pike MC, Riboli E, Rodriguez C, Setiawan WV, Stampfer MJ, Stram DO, Thomas G, Thun MJ, Travis R, Trichopoulou A, Virtamo J, Wacholder S. Genetic variation in the HSD17B1 gene and risk of prostate cancer. PLoS Genet. 2005; 1:e68 [PMID: 16311626]
  53. Tsuchiya M, Nakao H, Katoh T, Sasaki H, Hiroshima M, Tanaka T, Matsunaga T, Hanaoka T, Tsugane S, Ikenoue T. Association between endometriosis and genetic polymorphisms of the estradiol-synthesizing enzyme genes HSD17B1 and CYP19. Hum Reprod. 2005; 20:974-8. [PMID: 15640252]
  54. Wu AH, Seow A, Arakawa K, Van Den Berg D, Lee HP, Yu MC. HSD17B1 and CYP17 polymorphisms and breast cancer risk among Chinese women in Singapore. Int J Cancer. 2003; 104:450-7. [PMID: 12584742]
  55. Makridakis N, Ross RK, Pike MC, Chang L, Stanczyk FZ, Kolonel LN, Shi CY, Yu MC, Henderson BE, Reichardt JK. A prevalent missense substitution that modulates activity of prostatic steroid 5alpha-reductase. Cancer Res. 1997; 57:1020-2. [PMID: 9067262]
  56. Ntais C, Polycarpou A, Ioannidis JP. SRD5A2 gene polymorphisms and the risk of prostate cancer: a meta-analysis. Cancer Epidemiol Biomarkers Prev. 2003; 12:618-24. [PMID: 12869400]
  57. Roberts RO, Bergstralh EJ, Farmer SA, Jacobson DJ, McGree ME, Hebbring SJ, Cunningham JM, Anderson SA, Thibodeau SN, Lieber MM, Jacobsen SJ. Polymorphisms in the 5alpha reductase type 2 gene and urologic measures of BPH. Prostate. 2005; 62:380-7. [PMID: 15389785]
  58. Bethke L, Webb E, Sellick G, Rudd M, Penegar S, Withey L, Qureshi M, Houlston R. Polymorphisms in the cytochrome P450 genes CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1, CYP19A1 and colorectal cancer risk. BMC Cancer. 2007; 7:123 [PMID: 17615053]
  59. Haiman CA, Hankinson SE, Spiegelman D, Brown M, Hunter DJ. No association between a single nucleotide polymorphism in CYP19 and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2002; 11:215-6. [PMID: 11867511]
  60. Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH. Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001; 69:138-47. [PMID: 11404819]
  61. Hahn LW, Ritchie MD, Moore JH. Multifactor dimensionality reduction software for detecting gene-gene and gene-environment interactions. Bioinformatics. 2003; 19:376-82. [PMID: 12584123]
  62. Moore JH. Computational analysis of gene-gene interactions using multifactor dimensionality reduction. Expert Rev Mol Diagn. 2004; 4:795-803. [PMID: 15525222]
  63. Ritchie MD, Hahn LW, Moore JH. Power of multifactor dimensionality reduction for detecting gene-gene interactions in the presence of genotyping error, missing data, phenocopy, and genetic heterogeneity. Genet Epidemiol. 2003; 24:150-7. [PMID: 12548676]
  64. Shimizu H, Ross RK, Bernstein L, Pike MC, Henderson BE. Serum oestrogen levels in postmenopausal women: comparison of American whites and Japanese in Japan. Br J Cancer. 1990; 62:451-3. [PMID: 2206953]
  65. Bernstein L, Yuan JM, Ross RK, Pike MC, Hanisch R, Lobo R, Stanczyk F, Gao YT, Henderson BE. Serum hormone levels in pre-menopausal Chinese women in Shanghai and white women in Los Angeles: results from two breast cancer case-control studies. Cancer Causes Control. 1990; 1:51-8. [PMID: 2102277]
  66. McBride MW, McVie AJ, Burridge SM, Brintnell B, Craig N, Wallace AM, Wilson RH, Varley J, Sutcliffe RG. Cloning, expression, and physical mapping of the 3beta-hydroxysteroid dehydrogenase gene cluster (HSD3BP1–HSD3BP5) in human. Genomics. 1999; 61:277-84. [PMID: 10552929]
  67. Luu-The V, Lachance Y, Labrie C, Leblanc G, Thomas JL, Strickler RC, Labrie F. Full length cDNA structure and deduced amino acid sequence of human 3 beta-hydroxy-5-ene steroid dehydrogenase. Mol Endocrinol. 1989; 3:1310-2. [PMID: 2779585]
  68. Rhéaume E, Lachance Y, Zhao HF, Breton N, Dumont M, de Launoit Y, Trudel C, Luu-The V, Simard J, Labrie F. Structure and expression of a new complementary DNA encoding the almost exclusive 3 beta-hydroxysteroid dehydrogenase/delta 5-delta 4-isomerase in human adrenals and gonads. Mol Endocrinol. 1991; 5:1147-57. [PMID: 1944309]
  69. Simard J, Durocher F, Mebarki F, Turgeon C, Sanchez R, Labrie Y, Couet J, Trudel C, Rheaume E, Morel Y, Luu-The V, Labrie F. Molecular biology and genetics of the 3 beta-hydroxysteroid dehydrogenase/delta5-delta4 isomerase gene family. J Endocrinol. 1996; 150Suppl:S189-207. [PMID: 8943802]
  70. Gingras S, Simard J. Induction of 3beta-hydroxysteroid dehydrogenase/isomerase type 1 expression by interleukin-4 in human normal prostate epithelial cells, immortalized keratinocytes, colon, and cervix cancer cell lines. Endocrinology. 1999; 140:4573-84. [PMID: 10499513]
  71. Gingras S, Moriggl R, Groner B, Simard J. Induction of 3beta-hydroxysteroid dehydrogenase/delta5-delta4 isomerase type 1 gene transcription in human breast cancer cell lines and in normal mammary epithelial cells by interleukin-4 and interleukin-13. Mol Endocrinol. 1999; 13:66-81. [PMID: 9892013]
  72. Balacco C, Santoro G, Santoro M, Abaticchio G, Vavalle G, Giorgino R. Plasmatic and urinary steroids in high myopia. Note 1-Testosterone, 17-OH, 17-KS. Boll Soc Ital Biol Sper. 1978; 54:975-7. [PMID: 743414]
  73. Ojha RK, Singh R, Maurya O, Agrawal JK. Myopia and plasma cortisol. Indian J Ophthalmol. 1989; 37:91-3. [PMID: 2583791]
  74. Yu CY, Shih YF, Lin LL, Hung PT, Hou PK. Study of β-Estradiol, Cortisol and Testosterone on High Myopes. Tzu-Chi medical Journal. 1994; 6:179-83.
  75. Arrighi N, Bodei S, Peroni A, Mirabella G, Zani D, Simeone C, Cunico SC, Spano P, Sigala S. Detection of muscarinic receptor subtypes in human urinary bladder mucosa: age and gender-dependent modifications. Neurourol Urodyn. 2008; 27:421-8. [PMID: 17924537]
  76. Fragkouli A, Stamatakis A, Zographos E, Pachnis V, Stylianopoulou F. Sexually dimorphic effects of the Lhx7 null mutation on forebrain cholinergic function. Neuroscience. 2006; 137:1153-64. [PMID: 16338089]
  77. Rhodes ME, O'Toole SM, Wright SL, Czambel RK, Rubin RT. Sexual diergism in rat hypothalamic-pituitary-adrenal axis responses to cholinergic stimulation and antagonism. Brain Res Bull. 2001; 54:101-13. [PMID: 11226719]
  78. Chrissobolis S, Sobey CG. Influence of gender on K+-induced cerebral vasodilatation. Stroke. 2004; 35:747-52. [PMID: 14764930]
  79. Yoshida T, Kuwabara Y, Sasaki M, Fukumura T, Ichimiya A, Takita M, Ogomori K, Ichiya Y, Masuda K. Sex-related differences in the muscarinic acetylcholinergic receptor in the healthy human brain–a positron emission tomography study. Ann Nucl Med. 2000; 14:97-101. [PMID: 10830526]
  80. Sharif NA. Chronic prolactin, gonadal and thyroid hormone treatments in vivo alter levels of TRH and muscarinic receptors in male and female rat tissues. Brain Res. 1988; 449:364-8. [PMID: 2840173]
  81. Sokolovsky M, Egozi Y, Avissar S. Molecular regulation of receptors: interaction of beta-estradiol and progesterone with the muscarinic system. Proc Natl Acad Sci USA. 1981; 78:5554-8. [PMID: 6946492]
  82. Olsen KL, Edwards E, Schechter N, Whalen RE. Muscarinic receptors in preoptic area and hypothalamus: effects of cyclicity, sex and estrogen treatment. Brain Res. 1988; 448:223-9. [PMID: 3378147]
  83. Al-Dahan MI, Thomas PJ. Contrasting effects of testicular and ovarian steroids upon muscarinic binding sites in the brain. Pharmacology. 1987; 34:250-8. [PMID: 3615568]
  84. Dohanich GP, Witcher JA, Weaver DR, Clemens LG. Alteration of muscarinic binding in specific brain areas following estrogen treatment. Brain Res. 1982; 241:347-50. [PMID: 7104717]
  85. Egozi Y, Avissar S, Sokolovsky M. Muscarinic mechanisms and sex hormone secretion in rat adenohypophysis and preoptic area. Neuroendocrinology. 1982; 35:93-7. [PMID: 7133321]
  86. Setiawan VW, Hankinson SE, Colditz GA, Hunter DJ, De VI. HSD17B1 gene polymorphisms and risk of endometrial and breast cancer. Cancer Epidemiol Biomarkers Prev. 2004; 13:213-9. [PMID: 14973105]
  87. Travis RC, Churchman M, Edwards SA, Smith G, Verkasalo PK, Wolf CR, Wolf H, Key TJ. No association of polymorphisms in CYP17, CYP19, and HSD17–B1 with plasma estradiol concentrations in 1,090 British women. Cancer Epidemiol Biomarkers Prev. 2004; 13:2282-4. [PMID: 15598795]
  88. Dunning AM, Dowsett M, Healey CS, Tee L, Luben RN, Folkerd E, Novik KL, Kelemen L, Ogata S, Pharoah PD, Easton DF, Day NE, Ponder BA. Polymorphisms associated with circulating sex hormone levels in postmenopausal women. J Natl Cancer Inst. 2004; 96:936-45. [PMID: 15199113]
  89. Saloniemi T, Lamminen T, Huhtinen K, Welsh M, Saunders P, Kujari H, Poutanen M. Activation of androgens by hydroxysteroid (17beta) dehydrogenase 1 in vivo as a cause of prenatal masculinization and ovarian benign serous cystadenomas. Mol Endocrinol. 2007; 21:2627-36. [PMID: 17666583]
  90. Puranen TJ, Poutanen MH, Peltoketo HE, Vihko PT, Vihko RK. Site-directed mutagenesis of the putative active site of human 17 beta-hydroxysteroid dehydrogenase type 1. Biochem J. 1994; 304:289-93. [PMID: 7998947]
  91. Johnson MM, Houck J, Chen C. Screening for deleterious nonsynonymous single-nucleotide polymorphisms in genes involved in steroid hormone metabolism and response. Cancer Epidemiol Biomarkers Prev. 2005; 14:1326-9. [PMID: 15894696]