- Short report
- Open Access
Assessing the association of single nucleotide polymorphisms in thyroglobulin gene with age of puberty in bulls
Journal of Animal Science and Technologyvolume 56, Article number: 17 (2014)
Puberty is a stage of sexual development determined by the interaction of many loci and environmental factors. Identification of genes contributing to genetic variation in this character can assist with selection for early pubertal bulls, improving genetic progress in livestock breeding. Thyroid hormones play an important role in sexual development and spermatogenic function. The objective of this study was to evaluate the association between single nucleotide polymorphisms (SNPs) located in thyroglobulin(TG) gene with age of puberty in Angus bulls. Four SNPs were genotyped in 273 animals using SEQUENOM technology and the association between markers and puberty age was analyzed. Results showed a significant association (P < 0.05) between these markers and puberty age estimated at a sperm concentration of 50 million and a progressive motility of 10%. This is the first report of an association of TG polymorphisms with age of puberty in bulls, and results suggest the importance of thyroidal regulation in bovine sexual development and arrival to puberty.
Puberty is a stage of sexual development determined by the interaction of many loci hierarchically arranged in networks, and environmental factors (Ojeda et al., ). Puberty in cattle is an important target for genetic improvement so early prediction using genetic markers is a goal for livestock breeding. The identification of new genes and/or mutations contributing to genetic variation in puberty can assist with the selection for early pubertal bulls, reducing the generation interval and increasing fertility and genetic progress (Johnston et al., , Fortes et al., ).
Thyroid hormones (THs) exert a broad range of effects on metabolism, growth, homeostasis control, and other biological processes (Warner and Mittag, , Mullur et al., ), and show connection with nearly every biological endocrine system. Several studies provide evidence to confirm the role of THs in sexual differentiation and gonadal development (Jannini et al., , Mendis-Handagama and Siril Ariyaratne, , Flood et al., , Duarte-Guterman et al., ). For example, it has been discovered that deiodinases [the enzymes responsible for ioding thyroglobulin (TG) to obtain the active forms triiodothyronine (T3) and thyroxine (T4)] and thyroid receptors (TRs, encoded by trα and trβ genes) are present within gonadal tissues, suggesting that THs must have an action on these organs (Wagner et al., ). The presence of TH machinery in testicular tissues implies that TH axis must regulate aspects of testicular functioning. Indeed, it has been shown that hypo- and hyperthyroid males exhibit testes and sperm dysfunction (Krassas et. al., ). In addition, THs are involved in the regulation of androgen receptors(ARs) expression in testicular tissues through thyroid response elements (TREs) located in the promoter of ar genes(Flood et al., ). Furthermore, THs may also regulate other genes involved in androgen biosynthesis and signaling. For example, THs enhance 5α-reductase expression and activity within the testes, increasing 5α-dihydrotestosterone concentrations (Ram and Waxman, , Duarte-Guterman et al., ). On the other hand, it has been reported that TH-related transcription factors influence the expression of sox9, which induces differentiation of the bipotential cells in the testes into Sertoli cells (Zhou et al., ). Sox9 stimulates the nuclear receptor steroidogenic factor 1 (sf1), which is primarily expressed in Leydig cells and plays an important role in sexual differentiation (Zhao et al., ). Dax1 gene works in parallel with sf1to regulate testicular differentiation and can regulate TH-related gene expression (Sugawara et. al., ).
Beyond the fact that THs have considerable effects on the hypothalamic-pituitary-gonadal axis (HPG) and sexual development, some studies examined the involvement of androgens in TH synthesis and metabolism. For example, gonadotropin-releasing hormone (GnRH) interferes with the hypothalamic-pituitary-thyroid axis (HPT), increasing thyroid stimulating hormone (TSH) secretion. Moreover, ARs were identified in the thyroid gland of different vertebrate species, suggesting that the androgen axis directly regulates TH synthesis and metabolism (Pelletier et al., ). Other works have demonstrated the potential of ARs to regulate TH axis showing that TRs transcriptional levels and distribution within testes are responsive to androgen fluctuations through the presence of androgen response elements (AREs) in the promoter regions of TH-related genes. These evidences demonstrate the existence of a considerable cross-regulation between both axes.
Considering the existing evidence on the biological connection between metabolic status, regulated to a large extent by thyroid hormones, and sexual development, and taking into account the fact that THs play a role in the regulation of sexual function, we decided to evaluate the possible associations between SNPs located in the 3’ flanking region of TG gene with the age of puberty in bull calves.
In order to reach this objective, DNA was extracted from blood samples belonging to 273 Angus bulls using Wizard Genomic kit, following manufacturer instructions (Promega, Madison, WI, USA). Four SNPs (rs378215592, rs110406764, rs109662686, rs109057985) in TG previously reported by Hou et al.,  were genotyped using SEQUENOM platform by GeneSeek Inc. genotyping services (Lincoln, NE, USA). This technology is based on primer-extension reaction that generates allele-specific products with distinct masses, which are then detected through MALDI-TOF mass spectrometry (http://www.sequenom.com/). Detailed information of the studied SNPs is presented in Table 1. Animal samples, phenotypic measurements and estimation of puberty ages used in this work were reported in Liron et al., . The estimated puberty ages were: i. age at 28 cm of scrotal circumference (SC28), and ii. age at sperm concentration 50 million and percentage of progressive motility 10% (C50 - M10). In order to evaluate the linkage disequilibrium (LD) between the three studied SNPs, the haplotypes for each individual were constructed using Phase algorithm (Li and Stephens, ).
The association between haplotype markers of TG gene and the estimated puberty ages was analyzed utilizing MIXED procedure implemented in SAS 9.0 software (SAS Inst. Inc.). The linear mixed model used to analyze the association between puberty age and genotypes was the following:
Where Yijkl = phenotypic observation of the I bull, μ = the overall mean, Si = the fixed effect of ith year, Gj = the fixed effect of jth genotype, Bk = the fixed effect of kth herd, Ol = random effect of lth sire, and eijkI = random error.
Results and discussion
After analyzing the genotyping results, one SNP (rs109662686) was removed given that it exhibited a call rate lower than 57%. The LD analysis indicated that the three remaining SNPs were completely linked (r2 = 1). Only two (TGT and CAG) of the eight possible haplotypes were found, with gene frequencies of 0.81and 0.19 for TGT and CAG haplotypes, respectively. Genotype frequencies values were 0.64 for homozygote TGT, 0.03 for homozygote CAG and 0.33 for heterozygote bulls. The obtained haplotypes for the 273 bulls were tested for association with phenotypic data for the two estimated puberty ages mentioned above. The association analysis showed a significant association (P < 0.05) between the haplotype markers and puberty age estimated at C50 and M10. Homozygote TGT exhibited a mean +/- S.E. age at C50 and M10 of 289.74 ± 8.13 days, while homozygote CAG showed 347.67 ± 22.51 days, resulting in a difference of 57.93 days of age. Heterozygote animals showed a mean +/- S.E. age of 299.65 ± 8.72 days of age. No significant association was found between both haplotypes and age at puberty estimated at SC28 (TGT/TGT = 277.95 ± 23.85 days, TGT/CAG = 281.09 ± 25.82 days and CAG/CAG = 292.00 ± 17.76 days, P > 0.05).
The results obtained here constitute the first report of an association of TG gene polymorphisms with age of puberty in bulls and could be explained by the vast amount of works in which the modulatory influence of THs on male reproduction is demonstrated (Flood et al., , Duarte-Guterman et al., ). The potential of THs in the modulation of male reproductive functions during or preceding puberty was determined to be of such importance at the point that any alteration in their expression and/or concentration has profound effects on male reproduction (Krassas et al., , Weber et al., ). Despite the indirect regulation of sexual maturation by THs through their known roles in development, metabolism, hormonal regulation and other physiological processes, evidence indicates that THs have direct effects on sexual development, reproductive function and associated molecular mechanisms and pathways. Although the specific mechanisms underlying this regulation are not completely established, these direct effects could be exerted mainly through the presence of TH machinery in gonadal tissues. For example, TRs are widely distributed and expressed in different compartments of the testis in mammalian species, which suggests a direct regulatory role for THs in male gonadal development and function (Jannini et al., , Kumar et al., , Wagner et al., ). Studies have also identified deiodinases in the testes of vertebrate species, whose role within testicular functioning in mammalian species has been reviewed by Wagner et al., . Another evidence for the involvement of THs in puberty is the presence of TREs in androgen receptors in the testis and in GnRH and luteinizing hormone receptor (LHR) promoter region (Tsai-Morris et al., ), which demonstrates that THs can directly regulates androgen biosynthesis. Hyper- and hypo-thyroidic conditions alter GnRH concentrations in mammals, consequently affecting LH and FSH production and secretion (Chiao et al., ). Furthermore, several studies have shown a fall in circulating testosterone levels in hypothyroid humans (Kumar et al, ). As we said before, this regulation is bidirectional and there area lot of studies describing the considerable cross-regulation existing between HPT and HPG axes in vertebrates. Interestingly, Fortes et al.,  detected SERPINA7 gene on chromosome X associated with percentage of normal sperm (PNS) in Brahman bulls. This gene codes for thyroxine-binding globulin (TBG), the major TH transport protein in serum. This evidence reinforces the hypothesis that THs play an important role in bovine male fertility.
Despite the broad range of effects in vertebrates and the connections in nearly every biological endocrine system, a lot of studies provide enough evidence to confirm the role of THs in sexual differentiation and gonadal development in mammalian and non-mammalian species. We can affirm that THs influence steroidogenesis and spermatogenesis and as described above there is extensive evidence that links TG gene to testicular development, existing a direct crosstalk between HPG and thyroid hormones axis (Wagner et al., , Nonneman et al., ).
In conclusion, we detected an association between TG polymorphisms and age at puberty at C50 and M10 in male Angus cattle. Our results could contribute to the investigation on regulation of bovine puberty and fertility. THs are essential for normal growth, sexual development and reproductive function, and normal thyroid activity seems to be a requisite for an adequate male reproductive function. However, the knowledge about the interaction between both endocrine axis, HTP and HTG, is still rudimentary and needs further investigation.
Ojeda SR, Dubay C, Lomniczi A, Kaidar G, Matagne V, Sandau US, Dissen GA: Gene networks and the neuroendocrine regulation of puberty. Mol Cell Endocrinol. 2010, 324: 3-11. 10.1016/j.mce.2009.12.003.
Johnston DJ, Barwick SA, Corbet NJ, Fordyce G, Holroyd RG, Williams PJ, Burrow HM: Genetics of heifer puberty in two tropical beef genotypes in northern Australia and associations with heifer and steer-production traits. Anim Prod Sci. 2009, 49: 399-412. 10.1071/EA08276.
Fortes MRS, Lehnert SA, Bolormaa S, Reich C, Fordyce G, Corbet NJ, Whan V, Hawken RJ, Reverter A: Finding genes for economically important traits: Brahman cattle puberty. Anim Prod Sci. 2012, 52: 143-150. 10.1071/AN11165.
Warner A, Mittag J: Thyroid hormone and the central control of homeostasis. J Mol Endocrinol. 2012, 49: 29-35. 10.1530/JME-12-0068.
Mullur R, Liu YY, Brent GA: Thyroid hormone regulation of metabolism. Physiol Rev. 2014, 94 (Suppl2): 355-382.
Jannini EA, Crescenzi A, Rucci N, Screponi E, Carosa E, de Matteis A, Macchia E, d’Amati G, D’Armiento M: Ontogenetic pattern of thyroid hormone receptor expression in the human testis. J Clin Endocrinol Metab. 2000, 85: 3453-3457. 10.1210/jcem.85.9.6803.
Mendis-Handagama SM, Siril Ariyaratne HB: Leydig cells, thyroid hormones and steroidogenesis. Indian J Exp Biol. 2005, 43 (Suppl11): 939-962.
Flood DE, Fernandino JI, Langlois VS: Thyroid hormones in male reproductive development: evidence for direct crosstalk between the androgen and thyroid hormone axes. Gen Comp Endocrinol. 2013, 192: 2-14.
Duarte-Guterman P, Navarro-Martín L, Trudeau VL: Mechanisms of crosstalk between endocrine systems: Regulation of sex steroid hormone synthesis and action by thyroid hormones. Gen Comp Endocrinol. 2014, doi:10.1016/j.ygcen.2014.03.015
Wagner MS, Wajner SM, Maia AL: Is there a role for thyroid hormone on spermatogenesis?. Microsc Res Tech. 2009, 72: 796-808. 10.1002/jemt.20759.
Krassas GE, Papadopoulou F, Tziomalos K, Zeginiadou T, Pontikides N: Hypothyroidism has an adverse effect on human spermatogenesis: A prospective, controlled study. Thyroid. 2008, 18 (Suppl12): 1255-1259.
Ram PA, Waxman DJ: Pretranslational control by thyroid hormone of rat liver steroid 5 alpha-reductase and comparison to the thyroid dependence of two growth hormone-regulated CYP2C mRNAs. J Biol Chem. 1990, 265 (31): 19223-19229.
Zhou RJ, Bonneaud N, Yua CX, Barbara PD, Boizet B, Tibor S, Scherer G, Roeder RG, Poulat F, Berta P: SOX9 interacts with a component of the human thyroid hormone receptor-associated protein complex. Nucleic Acids Res. 2002, 30 (Suppl 14): 3245-3252.
Zhao L, Bakke M, Krimkevich Y, Cushman LJ, Parlow AF, Camper SA, Parker KL: Steroidogenic factor 1 (SF1) is essential for pituitary gonadotrope function. Development. 2001, 128 (Suppl 2): 147-154.
Sugawara T, Kiriakidou M, McAllister JM, Holt JA, Arakane F, Strauss JF: Regulation of expression of the steroidogenic acute regulatory protein (StAR) gene: a central role for steroidogenic factor 1. Steroids. 1997, 62: 5-9. 10.1016/S0039-128X(96)00152-3.
Pelletier G: Localization of androgen and estrogen receptors in rat and primate tissues. Histol Histopathol. 2000, 15: 1261-1270.
Hou GY, Yuan ZR, Zhou HL, Zhang LP, Li JY, Gao X, Wang DJ, Gao HJ, Xu SZ: Association of thyroglobulin gene variants with carcass and meat quality traits in beef cattle. Mol Biol Rep. 2011, 38: 4705-4708. 10.1007/s11033-010-0605-1.
Lirón JP, Prando A, Fernández ME, Ripoli MV, Rogberg-Muñoz A, Goszczynski DE, Posik DM, Peral-García P, Baldo A, Giovambattista G: Association between GNRHR, LHR and IGF1 polymorphisms and timing of puberty in male Angus cattle. BMC Genet. 2012, 5: 13-26.
Li N, Stephens M: Modelling linkage disequilibrium, and identifying recombination hotspots using SNP data. Genetics. 2003, 165: 221-233.
Weber G, Vigone MC, Stroppa L, Chiumello G: Thyroid function and puberty. J Pediatr Endocrinol Metab. 2003, 16 (Suppl 2): 253-257.
Kumar A, Shekhar S, Dhole B: Thyroid and male reproduction. Indian J Endocr Metab. 2014, 18 (Suppl 1): 23-31.
Tsai-Morris CH, Xie XZ, Buczko E, Dufau ML: Promoter and regulatory regions of the rat luteinizing hormone receptor gene. J Biol Chem. 1993, 268 (Suppl 6): 4447-4452.
Chiao YC, Lee HY, Wang SW, Hwang JJ, Chien CH, Huang SW, Lu CC, Chen JJ, Tsai SC, Wang PS: Regulation of thyroid hormones on the production of testosterone in rats. J Cell Biochem. 1999, 73: 554-562. 10.1002/(SICI)1097-4644(19990615)73:4<554::AID-JCB13>3.0.CO;2-L.
Fortes MRS, Reverter A, Hawken RJ, Bolormaa S, Lehnert SA: Candidate Genes Associated with Hormone Levels of Inhibin, Luteinising Hormone, and Insulin-Like Growth Factor 1, Testicular Development, and Sperm Quality in Brahman Bulls. Biol Reprod. 2012, 87 (Suppl 3): 58-
Nonneman D, Rohrer GA, Wise TH, Lunstra DD, Ford JJ: A variant of porcine thyroxine‒binding globulin has reduced affinity for thyroxine and is associated with testis size. Biol Reprod. 2005, 72: 214-220. 10.1095/biolreprod.104.031922.
None of authors of this paper has financial, personal or other relationship with other people or organizations that could inappropriately influence the content of this paper.
MEF and DEG participated in the genetic and statistical studies, drafted the manuscript. AP conceptualized the study and collected the phenotypic data. PPG participated in the design of the study. AB conceptualized the study and collected the phenotypic data. GG and JPL conceptualized and supervised the whole study, analyzed the data, performed the statistical analysis, drafted the manuscript. All authors read and approved the final version of the manuscript.
María Elena Fernández, Daniel Estanislao Goszczynski contributed equally to this work.