Figures
Abstract
Background
The circadian system has a major role in maintaining homeostasis and proper body functions including reproductive capacity. The aim of this study was to examine whether there is an association between genetic variability in the primary clock genes CLOCK and ARNTL and male infertility in humans.
Methodology/Principal Findings
We performed a case-control study, where we searched for an association between polymorphisms of CLOCK and ARNTL genes and male infertility in 961 Slovenian and Serbian Caucasian men. The study group consisted of 517 patients with idiopathic infertility and a control group of 444 fertile men. A statistically significant difference was found in genotype distribution between the two groups in the CLOCK gene: rs11932595 (p = 6·10−5, q = 4·10−4, OR equaled 1.9 with 95% CI 1.4–2.7), rs6811520 (p = 2·10−3, q = 8·10−3, OR = 1.7 with 95% CI 1.2–2.2) and rs6850524 (p = 0.01, q = 0.02, OR = 1.4 with 95% CI 1.1–1.9). Further analyses of haplotypes were consistent with genotyping results.
Citation: Hodžić A, Ristanović M, Zorn B, Tulić C, Maver A, Novaković I, et al. (2013) Genetic Variation in Circadian Rhythm Genes CLOCK and ARNTL as Risk Factor for Male Infertility. PLoS ONE 8(3): e59220. https://doi.org/10.1371/journal.pone.0059220
Editor: Francesc Palau, Instituto de Ciencia de Materiales de Madrid - Instituto de Biomedicina de Valencia, Spain
Received: October 25, 2012; Accepted: February 12, 2013; Published: March 18, 2013
Copyright: © 2013 Peterlin et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by National Research Agency of the Republic of Slovenia, grant no.: P3–0326 (http://www.arrs.gov.si/sl/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Infertility affects about 9% couples [1] and male factors contribute to approximately half of them [2]. About one third of infertility cases are unexplained (idiopathic male infertility), mostly due to our poor understanding of basic molecular mechanisms underlying male fertility [3], [4]. It is estimated that genetic factors are implicated in the pathogenesis of 50% males with idiopathic male infertility [5].
There is a growing body of evidence that circadian rhythms have a major role in maintaining homeostasis and proper body function including reproductive capacity [6]. This is highlighted most obviously in mutant mouse models whereby mutations in the Bmal1 gene significantly reduce fertility in male mice [7]. The circadian clock system as well as primary clock genes, are highly conserved between species. It is known that the circadian system influences testosterone production in humans, showing morning peaks and low levels in the evening [8], [9], [10]. Moreover, serum levels of sex steroids have been associated with genetic variants in circadian rhythm genes [11].
Molecular components of the circadian clock network represent positive and negative transcriptional-translational feedback loops of many genes. CLOCK and ARNTL genes represent the central node in the network, generating a positive loop, and heterodimerizeing and initiating the transcription of other clock genes. Resultant proteins forming the negative feedback loop inhibit the CLOCK and ARNTL transcriptional activity [12].
Therefore, we hypothesized that genetic variability of the CLOCK and ARNTL genes may be associated with male infertility in humans. To test the hypothesis we performed a retrospective case-control genetic association study of polymorphic sites in these two genes on a population of patients with male infertility in comparison with fertile male control population.
Materials and Methods
Ethics Statement
The study was approved by committees in both countries participating in the study: Slovenian National medical ethics committee (reference number: 73/05/12) and by the Ethics committee at Medical faculty at University of Belgrade, Serbia (reference number: 29/II-3). All patients gave informed written consent to participate in the study.
Subjects
Male partners of infertile couples, attending the infertility outpatient clinic of the Department of Obstetrics and Gynecology, University Medical Centre Ljubljana, and the Institute of Human Genetics, and Institute of Urology and Nephrology, Faculty of Medicine, University of Belgrade, participated in the study. The diagnosis of infertility was based on clinical assessment (testicular volume was measured using Prader's orchidometer), semen analysis, follicle stimulating hormone (FSH) testing and histologic evaluation of testicular biopsy specimens. Exclusion criteria were: patients with a history of testicular carcinoma, obstructive azoospermia including congenital bilateral absence of vas deferens, cytogenetic abnormalities and Y chromosome microdeletions.
Clinical, sperm and hormonal characteristics of the 517 infertile patients are summarized in Table 1.
The control group of proven fertility consisted of 444 healthy males, who were fathers of at least one child and who reported no history of infertility. The seminal parameters and biopsy were not assessed for subjects in the control group. All patients and control subjects were Caucasian of the Slavic (Slovene or Serbian) origin and were recruited in a consecutive manner.
Semen Parameters
In the 517 infertile men, semen analysis was performed according to the World Health organization (WHO) criteria [13]. Sperm was assessed in terms of volume, concentration, rapid progressive motility and normal morphology.
Hormonal Parameters
Serum FSH was measured by Microparticle Enzyme Immunoassay (AxSYM System, Abbott Laboratories, Chicago, IL, USA); the reference interval for FSH was 1–8 mIU/ml.
Testicular Biopsy and Histological Evaluation
Testicular biopsy and histological evaluation was performed in 77 patients with diagnosis of non obstructive azoospermia. In all patients testicular biopsy was performed for diagnostic and therapeutic (sperm recovery followed by cryostorage) purposes.
Following unilateral hemiscrototomy, a small incision was made and at least two samples of testicular tissue were taken from each testis. A systematic histological evaluation was performed under light microscopy; more than 100 seminiferous tubules were scored for each patient and results expressed as a relative number of tubules showing Sertoli cells, spermatogonia, spermatocytes, round and elongated spermatids, and spermatozoa. The diagnoses were as follows: Sertoli-cell only syndrome (SCOS) alone 46.8%, Maturation arrest (MA) alone 11.3%, Germinal hypoplasia (GH) alone 7.8%, and following combinations: GH+MA 14.1%, SCOS+MA 11.3%, GH+SCOS+MA 6.3% and GH+SCOS 2.1%.
SNP Selection and Genotyping
Genetic variants in both key circadian rhythm-regulating genes, CLOCK and ARNTL, were genotyped in this study. Eight tagging single nucleotide polymorphisms (SNPs) were chosen from both genes. Of these, four SNPs, rs11932595, rs6811520, rs6850524, and rs13124436 were selected in the CLOCK gene, and four SNPs, rs3789327, rs1481892, rs4757144, and rs12363415 in the ARNTL gene. We based the selection of SNPs on the known genetic linkage in both genes, according to HapMap Phase 3 (http://www.hapmap.org). The SNPs were regarded as proxies for neighbouring SNPs, when their pairwise r2 values exceeded 0.80. The set of most representative tagSNPs for ARNTL and CLOCK gene regions was obtained using the Tagger algorithm available through the Haploview software (Haploview, version 4.2) [14]. To maximize statistical power, only the SNPs with minor allele frequency values exceeding 0.05 were included in the investigated set.
Genomic DNA was isolated from the peripheral blood samples using standard procedures. SNPs genotyping was carried out by real time PCR method performed on 7000 Sequence Detection System (Applied Biosystems, Foster City, CA, USA) using KASPar SNP genotyping chemistry. The PCR reaction mix of 8 µL final volume consisted of 4 µL of DNA sample, 4 µL of Reaction Mix 2X, 0.11 µL Assay Mix, 0.064 µL 50 mM MgCl2, and 3.826 µL H2O. The protocol for PCR amplification was: initial denaturation step at 94°C for 15 minutes, then 20 cycles of denaturation at 94°C for 10 sec, followed by 57°C or 61°C for 5 sec, 72°C for 10 sec, 94°C for 10 sec, 57°C or 61°C for 20 sec, and final extension at 72°C for 40 sec.
The allelic discrimination analysis was performed using SDS Software Version 1.2 (Applied Biosystems, Foster City, CA, USA).
Genotype assignment was performed and interpreted independently by two investigators.
Statistical Analysis
Significances of associations between allelic and genotype frequencies and disease status were analyzed using the Chi-square test (χ2). Odds ratios (OR) and their respective 95% confidence intervals (CI) were also calculated to compare allelic and genotype distribution in patients and control subjects. Analyses were performed using R statistical language (R 2.15.0 for Windows). χ2 goodness-of-fit tests for deviation of genotype distribution from those predicted by Hardy-Weinberg equilibrium were also calculated, providing an additional quality control step of genotyping process. The investigated associations were regarded as significant when they reached p≤0.05.
As multiple SNPs were investigated, appropriate corrections of significance values were also applied using the Benjamini-Hochberg correction method (false-discovery rate – FDR values). We studied the extent of genetic linkage between SNPs using Haploview software version 4.2, which enables calculation of D’ and r2, representing the extent of pairwise linkage between two SNPs.
Haplotype frequencies and haplotype-disease associations were estimated using haplo.stats package for R (http://mayoresearch.mayo.edu/mayo/research/schaid_lab/software.cfm). Functions in haplo.stats package permit the calculation of indirectly measured haplotypes, under the founding assumption that subjects in the study are not related and the linkage phase unknown [15]. Here, the haplo.score function was used to directly ascertain differences in haplotype distributions across the groups of infertile men and control individuals. A global test of association, and per-haplotype association test were performed for both investigated genes. To reduce the effects of multiple testing we excluded all haplotypes with a frequency below 5% from downstream tests.
For power calculations, the pbsize2 function belonging to the gap package for R was utilized, which allows accurate power estimations under a variety of disease inheritance models (available at web address: http://cran.r-project.org/web/packages/gap/index.html).
Results
Genotype frequencies of investigated polymorphisms were in accordance with those predicted by the Hardy-Weinberg equilibrium in the study and in the control group (p<0.05), with the exception of rs12363415, which was excluded from further analyses. Genotype and allelic distribution of the CLOCK and ARNTL polymorphisms of the 517 infertile man and 444 fertile controls are shown in Table 2.
We found a statistically significant difference in the allelic distribution of rs11932595 (p = 6·10−5, q = 4·10−4), rs6811520 (p = 2·10−3, q = 8·10−3) and rs6850524 (p = 0.01, q = 0.02). However, we did not, find any significant association between rs13124436 polymorphism from CLOCK gene and male infertility. Under recessive genotype model, OR estimates ranged between 1.4 and 1.9 for CLOCK gene polymorphisms: rs11932595 (p = 6·10−5, OR = 1.9 with 95% CI 1.4–2.7), rs6811520 (p = 2·10−3, OR = 1.7 with 95% CI 1.2–2.2) and rs6850524 (p = 0.01, OR = 1.4 with 95% CI 1.1–1.9).
We also analyzed the inferred haplotypes in both investigated genes. The frequencies of predicted haplotypes in the study and the control group are presented in Table 3.
A statistically significant difference in haplotype distributions was also confirmed at the CLOCK gene locus when comparing the frequencies of haplotypes TCAG (p = 8·10−5, simulated p value after 10.000 permutations was 5·10−5), CCGG (p = 8·10−3, simulated p value was 8·10−3) and CGGG (p = 0.04, simulated p value equaled to 0.04) between the infertile patients and fertile controls.
The SNPs interrogated in the ARNTL gene, rs3789327, rs1481892 and rs4757144 did not show significant associations of genotype or allelic distribution between the two groups.
Accordingly, we did not find any significant difference comparing the frequencies of 6 most frequent haplotypes for the 4 analyzed SNPs in the ARNTL gene in the study and control groups.
Power analyses were performed to estimate the lower sensitivity threshold of our study to detect the variants characterized by low-to-modest effect sizes. For this purpose, the pbsize2 function of gap package for R was utilized. Calculations showed that our power to detect a significant result in the presence of the actual genotype-phenotype effect with genotype relative risk equal to at least 1.7, was 82.3%, when taking into account the sample size, the significance threshold of 0.05, the prevalence of male infertility in the general population equal to 4.5%, and disease allele frequency of at least 10%, and considering multiplicative model of genetic association.
Discussion
In the case-control study we found evidence of an association between male infertility and gene variants of the CLOCK gene in a sample of 961 men.
The circadian clock is internal timing system which allows an organism to provide environmental changes and adapt to them. Therefore, circadian rhythms manage a large variety of physiological and metabolic functions; any disruption of these rhythms may affect human health. Only a few studies have investigated the role of the circadian system in male fertility. It has been shown that night shift workers have an increased risk of infertility [16] and that infertile night shift workers have increased serotonin levels and decreased sperm quality compared to fertile night shift workers [17]. Serotonin is namely another component of the circadian system, potentially playing a role in human fertility [18] and is necessary for the development of normal spermatogenesis in rats [19]. More data is available regarding the role of the circadian system in male infertility in animals. A mutation of circadian genes influences reproductive fitness in Drosophila melanogaster. Beaver et al. (2002) have shown that clock-mutant males produce significantly fewer progeny, and release smaller quantities of sperm [20], whereas the studies of ARNTL knock-out mice have shown deficiencies in steroidogenesis [7], namely, male ARNTL KO mice had altered levels of reproductive hormones, indicating a defect in testicular Leydig cells. On the other hand, there is evidence that reproductive hormones directly influence the circadian system [21]; dysregulation of the either axis could therefore contribute to the reduced fertility. The circadian system influences testosterone production in humans, showing morning peaks and evening decrease [8], [9], [10]. Moreover, serum levels of sex steroid hormones have been associated with genetic variants in the circadian rhythm genes [11].
In women, genetic variability in the circadian rhythm genes ARNTL and NPAS2 has been suggested to contribute to fertility and seasonality, whereas the genetic variability in the ARNTL gene has been related to a higher number of pregnancies and also to a higher number of miscarriages; polymorphisms in the Npas2 gene have been associated with a decreased number of miscarriages [22].
In addition to its direct effects, the indirect effect of the circadian system on many physiological processes might possibly influence the male reproductive function. For example, the genetic variability in the CLOCK gene has been associated with increased weight and obesity [23], and loss of the ARNTL gene functions was shown to result in development of metabolic syndrome in knock-out rats [24]. Additionally, altered expression in the circadian rhythm genes has been shown to occur in obese males [25]. The negative effect of obesity on male reproductive function has been substantiated in numerous studies [26]–[28].
In conclusion, we provide evidence that genetic variability in the CLOCK gene might be associated with male infertility, and consequently imply the role of the circadian timing system in human reproduction. Further confirmation and mechanistic investigation of circadian system in male infertility is warranted.
Author Contributions
Conceived and designed the experiments: BP. Performed the experiments: AH. Analyzed the data: AH MR AM IN BP. Contributed reagents/materials/analysis tools: MR CT BZ. Wrote the paper: AH BP.
References
- 1. Boivin J, Bunting L, Collins JA, Nygren KG (2007) International estimates of infertility prevalence and treatment-seeking: potential need and demand for infertility medical care. Hum Reprod 22: 1506–1512.
- 2. Agarwal A, Said TM (2003) Role of sperm chromatin abnormalities and DNA damage in male infertility. Hum Reprod Update 9: 331–345.
- 3. Jungwirth A, Giwercman A, Tournaye H, Diemer T, Kopa Z, et al. (2012) European association of urology guidelines on male infertility: the 2012 update. Eur Urol 62: 324–32.
- 4. Matzuk MM, Lamb DJ (2002) Genetic dissection of mammalian fertility pathways. Nat Cell Biol 4 Suppl: s41–49
- 5. Hwang K, Yatsenko AN, Jorgez CJ, Mukherjee S, Nalam RL, et al. (2010) Mendelian genetics of male infertility. Ann N Y Acad Sci 1214: E1–E17.
- 6. Kennaway DJ, Boden MJ, Varcoe TJ (2012) Circadian rhythms and fertility. Mol Cell Endocrinol 349: 56–61.
- 7. Alvarez JD, Hansen A, Ord T, Bebas P, Chappell PE, et al. (2008) The circadian clock protein BMAL1 is necessary for fertility and proper testosterone production in mice. J Biol Rhythms 23: 26–36.
- 8. Cooke RR, McIntosh JE, McIntosh RP (1993) Circadian variation in serum free and non-SHBG-bound testosterone in normal men: measurements, and simulation using a mass action model. Clin Endocrinol (Oxf) 39: 163–171.
- 9. Teo W, McGuigan MR, Newton MJ (2011) The effects of circadian rhythmicity of salivary cortisol and testosterone on maximal isometric force, maximal dynamic force, and power output. J Strength Cond Res 25: 1538–1545.
- 10. Bribiescas RG, Hill KR (2010) Circadian variation in salivary testosterone across age classes in Ache Amerindian males of Paraguay. Am J Hum Biol 22: 216–220.
- 11. Chu LW, Zhu Y, Yu K, Zheng T, Chokkalingam AP, et al. (2008) Correlation between circadian gene variants and serum levels of sex steroids and insulin-like growth factor-I. Cancer Epidemiol Biomarkers Prev 17: 3268–3273.
- 12. Reppert SM, Weaver DR (2002) Coordination of circadian timing in mammals. Nature 418: 935–941.
- 13.
World Health Organization (1999) Laboratory Manual for the Examination of Human Semen and Sperm-Cervical Mucus Interaction, Fourth ed. UK: Cambridge University Press.
- 14. de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, et al. (2005) Efficiency and power in genetic association studies. Nat Genet 37: 1217–1223.
- 15. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA (2002) Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet 70: 425–434.
- 16. El-Helaly M, Awadalla N, Mansour M, El-Biomy Y (2010) Workplace exposures and male infertility - a case-control study. Int J Occup Med Environ Health 23: 331–338.
- 17. Ortiz A, Espino J, Bejarano I, Lozano GM, Monllor F, et al. (2010) The correlation between urinary 5-hydroxyindoleacetic acid and sperm quality in infertile men and rotating shift workers. Reprod Biol Endocrinol 8: 138.
- 18. Segal S, Sadovsky E, Palti Z, Pfeifer Y, Polishuk WZ (1975) Serotonin and 5-hydroxyindoleacetic acid in fertile and subfertile men. Fertil Steril 26: 314–316.
- 19. Aragon MA, Ayala ME, Marin M, Aviles A, Damian-Matsumura P, et al. (2005) Serotoninergic system blockage in the prepubertal rat inhibits spermatogenesis development. Reproduction 129: 717–727.
- 20. Beaver LM, Gvakharia BO, Vollintine TS, Hege DM, Stanewsky R, et al. (2002) Loss of circadian clock function decreases reproductive fitness in males of Drosophila melanogaster. Proc Natl Acad Sci U S A 99: 2134–2139.
- 21. Hagenauer MH, Lee TM (2011) Time for testosterone: the suprachiasmatic nucleus gets sexy. Endocrinology 152: 1727–1730.
- 22. Kovanen L, Saarikoski ST, Aromaa A, Lonnqvist J, Partonen T (2010) ARNTL (BMAL1) and NPAS2 gene variants contribute to fertility and seasonality. PLoS ONE 5: e10007.
- 23. Sookoian S, Gemma C, Gianotti TF, Burgueno A, Castano G, et al. (2008) Genetic variants of Clock transcription factor are associated with individual susceptibility to obesity. Am J Clin Nutr 87: 1606–1615.
- 24. Shimba S, Ogawa T, Hitosugi S, Ichihashi Y, Nakadaira Y, et al. (2011) Deficient of a Clock gene, Brain and Muscle Arnt-Like Protein-1 (BMAL1), induces Dyslipidemia and Ectopic Fat formation. PLos ONE 6: e25231.
- 25. Tahira K, Ueno T, Fukuda N, Aoyama T, Tsunemi A, et al. (2011) Obesity alters the expression profile of clock genes in peripheral blood mononuclear cells. Arch Med Sci 7: 933–940.
- 26. Cabler S, Agarwal A, Flint M, du Plessis SS (2010) Obesity: modern man's fertility nemesis. Asian J Androl 12: 480–489.
- 27. Du Plessis SS, Cabler S, McAlister DA, Sabanegh E, Agarwal A (2010) The effect of obesity on sperm disorders and male infertility. Nat Rev Urol 7: 153–161.
- 28. Hammoud AO, Gibson M, Peterson CM, Hamilton BD, Carrell DT (2006) Obesity and male reproductive potential. J Androl 27: 619–626.