Comparison of Expression Patterns of Six Canonical Clock Genes in Year-Round Estrous and Seasonal Estrous Rams

Weihao Chen1, 2, Zhilong Tian1, Lin Ma1, Shangquan Gan3, Wei Sun2, 4* and Mingxing Chu1*

1Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193, China

2College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China

3State Key Laboratory for Sheep Genetic Improvement and Healthy Production, Xinjiang Academy of Agricultural and Reclamation Sciences, Shihezi 832000, China

4Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225009, China

ABSTRACT

Circadian rhythm is a biological rhythm that is related closely to the rhythmic expression of a series of clock genes. A number of studies have revealed the role of circadian rhythms in the estrous mode of mammals. In the present study, the expression patterns of six canonical clock genes (Clock, BMAL1, Cry1, Cry2, Per1 and Per2) were analyzed in year-round estrous rams (Small Tail Han sheep, STH) and seasonal estrous rams (Sunite sheep, SNT). The result showed that all six genes were expressed in brain, cerebellum, hypothalamus, pituitary, testis, epididymis, vas deferens and adrenal gland tissues in both breeds. The expression level of Clock and BMAL1 showed similar trends in the brain, cerebellum, hypothalamus, pituitary, testis and epididymis in both breeds. The expression levels of Clock, BMAL1, and Cry1 were significantly higher in the pituitary tissue of STH rams than in that of SNT rams, whereas the expression level of Cry2 showed the opposite pattern. We speculate that Cry1 and Cry2 may have opposite roles in the circadian rhythm of rams. Moreover, the expression patterns of Cry1/2 and Per1/2 in the pituitary suggested that the CRY and PER proteins may function in the circadian rhythm either as a complex or as individual, Therefore, we concluded that circadian rhythmicity may regulate the estrous mode of rams via clock genes within transcription/translation feedback/feedforward loops. This is the first study to systematically analyze the expression patterns of clock genes in rams.

Article Information

Received 30 April 2019

Revised 30 June 2019

Accepted 11 September 2019

Available online 23 April 2021

(early access)

Published 10 January 2022

Authors’ Contribution

This study was designed by WC, ZT and MC, WC and ZT conducted the experiments and analyzed the data. WC drafted the manuscript. WS, SG, LM and MC helped in preparation of revised manuscript.

Key words

Ram, Circadian rhythm, Seasonal estrous, Clock genes, Tissue expression

DOI: https://dx.doi.org/10.17582/journal.pjz/20190430140456

* Corresponding author: [email protected]; [email protected]

0030-9923/2022/0002-0657 $ 9.00/0

Copyright 2022 Zoological Society of Pakistan

INTRODUCTION

Circadian rhythm, which is one of the basic characteristics of life, refers to biological processes that oscillate with a period of 24 h. Almost all creatures respond to cyclical changes in the external environment such as darkness, light and temperature through signal transduction processes, thus producing their own internal rhythms (Coulon et al., 2016; Lewis et al., 2018). It is widely accepted that the molecular mechanism of circadian rhythm is a feedback-loop model based on transcriptional translation, involving two transcriptional activators (BMAL1 and Clock) and two transcriptional repressors (period (PER) and Cryptochrome (CRY) proteins) (Leloup and Goldbeter, 2003; Preitner et al., 2002). The particular transcriptional feedback loops that are believed to control circadian rhythm are dependent on a small number of canonical clock genes, the expression levels of which are closely related to cycles in behavior and physiology (Janich et al., 2011; Marcheva et al., 2010; Nam et al., 2016).

To date, more than 10 genes have been identified that form the basis of cellular rhythmicity in mammals. These include two clock/cycle-related genes, Clock and BMAL1, two CRY genes, Cry1 and Cry2, and two PER genes, Per1 and Per2. (Angelousi et al., 2019; Moraes et al., 2017; Ye et al., 2018). Recent studies have suggested that Clock and BMAL1 play key roles in the formation of circadian rhythm (Chen et al., 2016; Trott and Menet, 2018). Clock was the first mammalian circadian clock gene to be discovered (King et al., 1997). BMAL1 encodes a functional chaperone of Clock; and their protein products form a dimer that binds to a specific motif, CACGTG (also known as “E-box”), in the promoter regions of Cry1, Cry2, Per1 and Per2. (Zheng et al., 2019).

In mammals, seasonal estrous is one of the biological activities that is regulated by circadian rhythm. In the hypothalamic-pituitary-gonadal axis, seasonal changes in light and temperature are translated into nerve impulses that act on the pineal gland via the suprachiasmatic nucleus (Hastings et al., 2000). This leads to the release of melatonin, which in turn regulates the release of gonadotropin-releasing hormone (Chappell et al., 2003; von Gall et al., 2000), ultimately affecting reproduction.

Sheep (Ovis aries) are a typical seasonal estrous species (Tang et al., 2018). However, in the current research into the molecular mechanism of seasonal estrous mode in sheep, researchers tend to focus mainly on ewes but ignore rams, even though, there also are strong seasonal rhythms related to reproduction in males of many animal species, most notably hamsters (Reiter, 1980) and sheep (Lincoln, 2002). Studies have shown that rams can identify seasonal changes in light, and translate them into molecular signals that affect gonadal function by regulating the secretion of melatonin (Kennaway, 2005). Considering the effects of the clock genes on seasonal estrous in ewes, it is important to explore their potential roles in ram reproduction.

Small Tail Han sheep (STH) and Sunite sheep (SNT) are two Chinese local sheep (Ovis aries) breeds with different estrous modes, year-round estrous and seasonal estrous (mainly in winter and spring), respectively. Both are known for their excellent meat production performance (Tang et al., 2018). Huge difference between the two sheep breeds in estrous modes have resulted in increasing interest in the tissue expression profiles of clock genes in these sheep. In the present study, we compared the tissue expression profiles and mRNA expression levels of six canonical clock genes in eight prolificacy-related tissues between STH and seasonal estrous SNT rams. Our study paves the way for in-depth study of the seasonal estrous mode of rams.

MATERIAL AND METHODS

Selection of experimental sheep and sample collection

Three STH and three SNT rams were supplied by the Yuncheng Breeding Sheep Farm (Yuncheng County, China) and the Sheep and Goat Breeding Farm of Tianjin Institute of Animal Sciences (Tianjin, China), respectively. The six rams were healthy, approximately 2.5 years old and were kept in a sheltered outdoor paddock. Alfalfa hay and concentrate were provided and clear water was available ad libitum. Eight tissues (brain, cerebellum, hypothalamus, pituitary, testis, epididymis, vas deferens and adrenal gland) were collected from each animal. All tissues were snap frozen in liquid nitrogen and then stored at−80 oC to be used for RNA extraction.

All experimental procedures used in the present study were approved by the Science Research Department (in charge of animal welfare issues) of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (IAS-CAAS) (Beijing, China). Ethical approval for the study and all of its protocols was provided by the animal ethics committee of IAS-CAAS (No. IASCAAS-AE-03, 12 December, 2016), which ensured that all efforts were taken to minimize pain and discomfort to the animals while conducting these experiments.

Total RNA extraction and cDNA synthesis

Total RNA was extracted from the eight collected tissues using a total RNA extraction kit for animal tissue (Tiangen, Beijing, China). Trizol (Invitrogen Inc., Carlsbad, CA, USA) was used to dissolve the tissues. Each tissue was homogenized and 50–100 mg samples were used for RNA extraction. The quantity and quality of total RNA were monitored using 1.5% agarose gel electrophoresis (U = 150 V;10 min) and ultraviolet spectrophotometry (UV-1201, Shimadzu, Kyoto, Japan), respectively. The A260/280 ratios (1.8–2.0) of the RNA samples were all 1.9 to 2.0, which showed that the extracted total RNA was of acceptable purity with no contamination or degradation. Therefore, the RNA preparations were deemed fit for use in the follow-up experiments, and so were stored at −80 oC until use.

First strand of cDNA was prepared using a PrimeScriptTM RT Reagent Kit according to the manufacturer’s instructions (TaKaRa Bio Inc., Dalian, China). The PCR thermocycler program was as follows: 37 oC for 15 min, followed by 85 oC for 5 s. The reaction mixture contained 1.0 µL Prime Script RT Enzyme, 1.0 µL random 6-mers, 4.0 µL 5 × Prime Script Buffer (for Real Time), 1.0 µL total RNA and 13 µL RNase-free ddH2O (total volume, 20 μL). Prior to storage at −80 oC, the standard working concentration of cDNA was 200 ng/µL. The quality of cDNA was evaluated by housekeeping gene (RPL-19) amplification, and cDNA were stored at −20 oC until use.

Primer design

Using the Primer Premier software (version 5.0, PREMIER Biosoft Co., Palo Alto, CA, USA), a total of seven primers were designed to amplify different fragments of the ovine Clock, BMAL1, Cry1, Cry2, Per1, Per2 and RPL-19 genes, based on their assembled

 

Table I. Primers used in this study to amplify the six canonical clock genes.

Gene Names	Primer Sequence (5′→3′)	Length (bp)	Tm (°C)	Accession No.
Clock	F: 5′-CAACGCACACATAGGCCTTC-3′ R: 5′-CTATTATGGGTGGTGCCCTGT-3′	181	60	NM_001130932.1
BMAL1	F: 5′-ATTGCAACCGGAAACGCAAG-3′ R: 5′-TGGTGGCACCTCGTAATGTT-3′	288	62	NM_001129734.1
Cry1	F: 5′-ACAGGTGGCGATTTTTGCTT-3′ R: 5′-TCCAGCTTCAGTTGCCAGTT-3′	215	61	NM_001129735.1
Cry2	F: 5′-AGGCTGTTCAAGGAATGGGG-3′ R: 5′-CGTAGGTCTCATCGTGGCTC-3	316	61	NM_001129736.1
Per1	F: 5′- GCCAGACAACCCTTCTACCAGT-3′ R: 5′- GGCTTGCACCTGCTTGACACA-3′	187	61	XM_027974931.1
Per2	F: 5′-TTACGACCACACATTCGCCA-3′ R: 5′-CCCCAGACTGCACGATCTTC-3′	171	61	XM_027967088.1
RPL-19	F: 5′-ATCGCCAATGCCAACTC-3′ R: 5′-CCTTTCGCTTACCTATACC-3′	154	60	XM_012186026.1

 

Abbreviations: bp, base pairs; Tm, melting temperature.

sequences in GenBank. All primers were synthesized by Beijing Tianyi Biotechnology Co., Ltd. (Beijing, China). The housekeeping gene (RPL-19, Genbank: XM_012186026.1) was used as an internal control to normalize the threshold cycle (Ct) values. Primers details are given in Table I.

Quantitative polymerase chain reaction

The expression levels of Clock, BMAL1, Cry1, Cry2, Per1 and Per2 in eight tissues (brain, cerebellum, hypothalamus, pituitary, testis, epididymis, vas deferens and adrenal gland) from SNT and STH rams were measured by quantitative polymerase chain reaction (qPCR).

The qPCR protocol used 20 μL of reaction mixture that contained 10 µL SYBR Premix EX Taq II (TaKaRa Bio Inc., Dalian, China), 0.8 µL each of forward and reverse primer (10 pmol/ul), 6.4 µL RNase-free ddH2O and 2 µL cDNA. Amplifications were performed in triplicate wells using the following PCR thermocycler program: denaturation at 95 °C for 5 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. The dissociation curve was analyzed after amplification. The melting temperature (Tm) peak observed at 85 °C±0.8 on the dissociation curve was used to determine the specificity of the PCR amplification.

Statistical analysis

The 2−ΔCt method (Livak, et al., 2001) was used to process the real-time PCR results. Statistical analyses were carried out using SPSS 19.0 software (IBM, Armonk, NY, USA). Levels of gene expression were analyzed for significant differences by one-way analysis of variance (ANOVA) followed by Fisher’s least significant difference test as a multiple comparison test. All experimental data are presented as mean ± standard error of the mean (SEM). A probability of p ≤ 0.05 was considered statistically significant, and a probability of p ≤ 0.01 was considered to be highly statistically significant.

RESULTS

Expression levels of clock

As shown in Figure 1, Clock was expressed in all the tissues. Brain was found to have the highest Clock expression level, followed by cerebellum and hypothalamic. Expression levels in the hypothalamus, pituitary and epididymis were significantly higher in STH samples compared with SNT samples (p < 0.05).

Expression levels of BMAL1

As shown in Figure 2, BMAL1 was expressed in all the tissues, The highest expression level was found in cerebellum, followed by brain and testis. Expression of BMAL1 in the brain, cerebellum, hypothalamus, pituitary, testis and epididymis were significantly higher in STH samples compared with SNT samples (p < 0.01).

Expression levels of Cry1

As shown in Figure 3, Cry1 was expressed in all the tissues. The highest expression level was found in testis, followed by cerebellum and hypothalamus. Expression levels of Cry1 in the pituitary (p < 0.05), testis (p < 0.01) and adrenal gland (p < 0.05) were significantly higher in SNT samples compared with STH samples.

 

 

Expression levels of Cry2

As shown in Figure 4, Cry2 was expressed in all the tissues. The highest expression level was detected in cerebellum, followed by brain and hypothalamus. Expression levels of Cry2 in the brain (p < 0.05), cerebellum (p < 0.01), hypothalamus (p < 0.01), pituitary (p < 0.05) and adrenal gland (p < 0.05) were significantly higher in STH samples compared with SNT samples.

 

 

Expression levels of Per1

As shown in Figure 5, Per1 was expressed in all the tissues. The highest expression level was found in the brain, followed by adrenal gland and hypothalamus. Expression levels of Per1 in the pituitary and epididymis were significantly higher in SNT samples compare with STH samples (p < 0.01). The expression level of Per1 in the testis was significantly higher in STH samples compare with SNT samples (p < 0.01).

Expression levels of Per2

As shown in Figure 6, Per2 was expressed in all the tissues. The highest expression level was in the brain, followed by cerebellum and hypothalamus. Expression levels of Per2 in the brain, hypothalamus, pituitary, testis and epididymis were significantly higher in SNT samples compared with STH samples (p < 0.05, p < 0.01), whereas the expression levels of Per2 in the vans deferens and adrenal gland were significantly higher in STH samples compared with SNT samples (p < 0.05 and p < 0.01, respectively).

 

 

DISCUSSION

Clock and BMAL1

Circadian rhythm is known to be crucial to the estrous mode of male animals, and members of the canonical clock gene family-Clock and BMAL1-are involved in male reproduction (Peruquetti et al., 2012). Alvarez et al. (2003) reported that down-regulation of Clock in mouse testis leads to a decrease in the fertilization capacity of the corresponding sperm, which highlights the critical function of Clock in spermatogenesis. BMAL1 knock-out male mice are rarely able to successfully fertilize normal female mice. Although the mechanism underlying this effect is not known, it may be related to the influence of BMAL1 on the secretion of testosterone and other steroid hormones (Alvarez et al., 2008). Furthermore, Clock and BMAL1 are involved in the regulation of energy metabolism, lipid metabolism (through a vital role in the differentiation and maturation of fat cells), glucose metabolism (including glucose balance and insulin resistance), and energy balance. (Lamia et al., 2008; Rudic et al., 2004; Shimba et al., 2005).

In sheep, the oscillating pattern of Clock expression is similar to that observed in other mammals (Kennaway et al., 1981). In this study, we found that Clock and BMAL1 were expressed in all the reproductive tissues, which reflects the expression patterns of these genes in humans and mice. This result supports the suggestion that Clock and BMAL1 influence the reproductive activity of rams directly by regulating the expression of related clock genes in the hypothalamic-pituitary-gonadal axis, as well as indirectly influencing ram breeding activities by regulating energy metabolism in other tissues (Pugazhendhi et al., 2019). In addition, the expression patterns of BMAL1 and Clock in the hypothalamus and pituitary (which have key roles in circadian rhythm) are similar to those of ewes (Lincoln et al., 2002), which implies similar roles of these genes in the two sexes.

Many studies have reported the joint effects of Clock and BMAL1 on reproduction (Kondratov et al., 2006; Tamayo et al., 2015; Zhang et al., 2016). Briefly, both genes contain the transcription factor PAS-HLH domain (Kondratov, et al., 2003), which contributes to the formation of Clock: BMAL1 dimers from the gene products. The dimer form activates transcription and translation of PER and CRY (Jung, et al., 2003; Kurbatova et al., 2012; Li et al., 2008), and regulates the frequency of the gonadotropin-releasing hormone pulse generator (Chappell et al., 2003), which in turn affects rhythmic biological activities such as seasonal estrous. In the present study, we found that the expression levels of Clock and BMAL1 in the brain, cerebellum, hypothalamus, pituitary and epididymis were higher in STH than in SNT rams (although this difference was not significant for Clock expression in the brain and cerebellum; p>0.05). Our findings are in agreement with reported expression patterns of Clock and BMAL1 in male black-line hamsters (Liu, 2016), indicating that these genes have similar regulatory functions in rams and other male seasonal estrous animals.

It seems plausible that the long period of estrous in rams may be related to the high expression levels of Clock and BMAL1. However, although expression of BMAL1 in the testis was significantly higher in STH rams compared with SNT rams, there was no significant difference in the expression levels of Clock between STH and SNT testis tissues. One possible explanation is that, in testis, Clock and BMAL1 do not have primary roles that affect ram reproduction. Our findings are in agreement with Morse et al. who found that testis did not exhibit rhythmicity of clock gene expression (Morse et al., 2003). Further studies are needed to investigate the relationship between Clock/BMAL1 and ram reproduction in more depth.

Period and cryptochromes

Besides Clock and BMAL1, the molecular mechanism of seasonal estrous involves the transcriptional activation of Per and Cry by the BMAL1: Clock dimer (Hirayama et al., 2003; Jang et al., 2015; Spoelstra et al., 2014). Cryptochrome is the only negative feedback loop involved in the regulation of circadian rhythm, and this has a strong effect on biological rhythms. In mammals, there are two variants, Cry1 and Cry2, which play opposite roles in the feedback loop (Duong et al., 2011). Cry1 knock-out mice exhibit a short-period circadian rhythm at behavioral and tissue/cell levels, whereas Cry2 knock-out mice exhibit a completely opposite phenotype (van der Horst et al., 1999). The reason behind these opposing phenotypes is still unclear. In mammals, the period gene family has three members, Per1, Per2, and Per3 (Tei et al., 1997). Mutations or deletions in Per1 or Per2 have been shown to cause changes in circadian rhythms, but no evidence of the effect of Per3 on circadian rhythms has been reported so far (Bae et al., 2001). Therefore, Per3 is considered to be unnecessary for circadian rhythm, and so Cry1, Cry2, Per1 and Per2 were selected for investigation of the seasonal estrous of rams in this study.

In the first study of period genes in mice (Tei et al., 1997), it was reported that Per/Cry genes are expressed in mouse heart, brain and testis. Further studies revealed that equivalent genes are expressed in many tissues in human (Hawkins et al., 2008). Our results demonstrate that Per1/2 and Cry1/2 are expressed in all eight of the selected tissues of STH and SNT rams, and are highly expressed in the brain, cerebellum, and hypothalamic-pituitary-gonadal axis. This implies a role of Per/Cry genes in ram reproduction.

Studies have shown that over-expression of Cry in mammalian cells can directly inhibit the transcriptional activity of the BMAL1: Clock dimer, whereas over-expression of Per has little effect on this process (Akashi et al., 2014; Leloup and Goldbeter, 2003). Therefore, we compared the expression levels of Per1, Per2, Cry1, Cry2, BMAL1 and Clock in different tissues, which revealed significantly higher expression of BMAL1 and Clock in the pituitary tissue of STH rams compared with SNT rams. This was also the case for Cry1 and Cry2. Considering the core function of the pituitary gland in circadian rhythm, we speculate that Cry1 and Cry2 may play opposite roles in the circadian rhythm of rams, similar to their reported functions in other mammals. Surprisingly, we found that Per1 and Per2 have similar expression patterns as Clock, BMAL1, and Cry1 in the pituitary tissue, which indicated that PER may be associated with the ram reproduction to some degree. Further research is necessary to confirm this idea.

It is worth noting that the stability of CRY is largely dependent on the presence of PER, they can form a complex to stabilize each other and regulate circadian rhythm (Zhou et al., 2018). We found the expression levels of Cry1, Per1, and Per2 were significantly higher in the pituitary tissue of SNT rams compared with STH rams, whereas, Cry2 showed the opposite pattern. Given the expression of the four genes in the pituitary gland in ewes (Lincoln et al., 2002), we speculate that their regulatory mechanism is similar in rams, namely, PER and CRY bind to form a complex that plays a role in circadian rhythm in the pituitary. However, these four genes did not have similar expression patterns in other tissues, indicating that CRY and PER can work either by forming a dimer or as monomeric proteins.

CONCLUSION

This study describes the expression pattern of six canonical clock genes in year-round estrous (STH) and seasonal estrous (SNT) rams. All six genes were expressed in the eight selected tissues, with high expression levels seen in the brain, cerebellum, hypothalamus, testis and epididymis. Our results suggest that circadian rhythmicity may regulate the estrous mode of rams via clock gene transcription/translation feedback/feedforward loops. However, the specific mechanism remains to be further explored. This is the first study to investigate the tissue-specific expression patterns of the six canonical clock genes in rams, and it provides a foundation for elucidating the molecular mechanisms of the estrous mode in rams.

ACKNOWLEDGEMENTS

This research was funded by the National Natural Science Foundation of China (31772580, 31872333); Earmarked Fund for China Agriculture Research System (CARS-38); Agricultural Science and Technology Innovation Program of China (ASTIP-IAS13); China Agricultural Scientific Research Outstanding Talents and their Innovative Teams Program; China High-level Talents Special Support Plan Scientific and Technological Innovation Leading Talents Program (W02020274); Tianjin Agricultural Science and Technology Achievements Transformation and Popularization Program (201704020); the Projects of Domesticated Animals Platform of the Ministry of Science, Key Research and Development Plan (modern agriculture) in Jiangsu Province (BE2018354); Major new varieties of agricultural projects in Jiangsu Province (PZCZ201739); Jiangsu Agricultural Science and Technology Innovation Fund (CX(18)2003); the Priority Academic Program Development of Jiangsu Higher Education Institutions; Major projects of Natural Science Research of Colleges and Universities in Jiangsu Province (17KJA230001); the project of six peak of talents of Jiangsu Province of China, Postgraduate Research & Practice Innovation Program of Jiangsu Province (SJCX18_0804) and the Graduate International Academic Exchange Fund of Yangzhou University. We thank Amy Phillips, PhD, from Liwen Bianji, Edanz Editing China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Statement of conflicts of Interest

All authors declare no conflicts of interest.

References

Akashi, M., Okamoto, A., Tsuchiya, Y., Todo, T., Nishida, E. and Node, K., 2014. A positive role for PERIOD in mammalian circadian gene expression. Cell Rep., 7: 1056-1064. https://doi.org/10.1016/j.celrep.2014.03.072

Alvarez, J.D., Chen, D., Storer, E. and Sehgal, A., 2003. Non-cyclic and developmental stage-specific expression of circadian clock proteins during murine spermatogenesis. Biol. Reprod., 69:81-91. https://doi.org/10.1095/biolreprod.102.011833

Alvarez, J.D., Hansen, A., Ord, T., Bebas, P., Chappell, P.E., Giebultowicz, J.M., Williams, C., Moss, S. and Sehgal, A., 2008. The circadian clock protein BMAL1 is necessary for fertility and proper testosterone production in mice. J. biol. Rhythms, 23:26-36. https://doi.org/10.1177/0748730407311254

Angelousi, A., Kassi, E., Nasiri-Ansari, N., Randeva, H.S., Kaltsas, G.A. and Chrousos, G.P., 2019. Clock genes and cancer development in particular in endocrine tissues. Endocr. Relat. Cancer, https://doi.org/10.1530/ERC-19-0094

Bae, K., Jin, X., Maywood, E.S., Hastings, M.H., Reppert, S.M. and Weaver, D.R., 2001. Differential functions of mPer1, mPer2, and mPer3 in the SCN circadian clock. Neuron, 30: 525-536. https://doi.org/10.1016/S0896-6273(01)00302-6

Busino, L., Bassermann, F., Maiolica, A., Lee, C., Nolan, P.M., Godinho, S.I., Draetta, G.F. and Pagano, M., 2007. SCFFbxl3 controls the oscillation of the circadian clock by directing the degradation of cryptochrome proteins. Science, 316: 900-904. https://doi.org/10.1126/science.1141194

Chappell, P.E., White, R.S. and Mellon, P.L., 2003. Circadian gene expression regulates pulsatile gonadotropin-releasing hormone (GnRH) secretory patterns in the hypothalamic GnRH-secreting GT1-7 cell line. J. Neurosci., 23: 11202-11213. https://doi.org/10.1523/JNEUROSCI.23-35-11202.2003

Chen, Y., Zhu, D., Yuan, J., Han, Z., Wang, Y., Qian, Z., Hou, X., Wu, T. and Zou, J., 2016. Clock-BMAL1 regulate the cardiac L-type calcium channel subunit CACNA1C through PI3K-Akt signaling pathway. Can. J. Physiol. Pharmacol., 94:1023-1032. https://doi.org/10.1139/cjpp-2015-0398

Coulon, N., Brailly-Tabard, S., Walter, M. and Tordjman, S., 2016. Altered circadian patterns of salivary cortisol in individuals with schizophrenia: A critical literature review. J. Physiol. Paris, 110: 439-447. https://doi.org/10.1016/j.jphysparis.2017.05.002

Duong, H.A., Robles, M.S., Knutti, D. and Weitz, C.J., 2011. A molecular mechanism for circadian clock negative feedback. Science, 332:1436-1439. https://doi.org/10.1126/science.1196766

Hastings, M. and Maywood, E.S., 2000. Circadian clocks in the mammalian brain. Bioessays, 22 :23-31. https://doi.org/10.1002/(SICI)1521-1878(200001)22:1<23::AID-BIES6>3.0.CO;2-Z

Hawkins, G. A., Meyers, D. A., Bleecker, E. R., and Pack, A. I., 2008. Identification of coding polymorphisms in human circadian rhythm genes per1, per2, per3, clock, arntl, cry1, cry2 and timeless in a multi-ethnic screening panel. DNA Seq., 19: 44-49. https://doi.org/10.1080/10425170701322197

Hirayama, J., Nakamura, H., Ishikawa, T., Kobayashi, Y. and Todo, T., 2003. Functional and structural analyses of cryptochrome. Vertebrate CRY regions responsible for interaction with the CLOCK: BMAL1 heterodimer and its nuclear localization. J. biol. Chem., 278: 35620-35628. https://doi.org/10.1074/jbc.M305028200

Jang, A.R., Moravcevic, K., Saez, L., Young, M.W. and Sehgal, A., 2015. Drosophila TIM binds importin alpha1, and acts as an adapter to transport PER to the nucleus. PLoS Genet., 11: e1004974. https://doi.org/10.1371/journal.pgen.1004974

Janich, P., Pascual, G., Merlos-Suarez, A., Batlle, E., Ripperger, J., Albrecht, U., Cheng, H.Y., Obrietan, K., Di Croce, L. and Benitah, S.A., 2011. The circadian molecular clock creates epidermal stem cell heterogeneity. Nature, 480: 209-214. https://doi.org/10.1038/nature10649

Jung, H., Choe, Y., Kim, H., Park, N., Son, G.H., Khang, I. and Kim, K., 2003. Involvement of CLOCK: BMAL1 heterodimer in serum-responsive mPer1 induction. Neuroreport, 14: 15-19. https://doi.org/10.1097/00001756-200301200-00003

Kennaway, D.J., 2005. The role of circadian rhythmicity in reproduction. Human Reprod. Update, 11: 91-101. https://doi.org/10.1093/humupd/dmh054

Kennaway, D.J., Obst, J.M., Dunstan, E.A. and Friesen, H.G., 1981. Ultradian and seasonal rhythms in plasma gonadotropins, prolactin, cortisol and testosterone in pinealectomized rams. Endocrinology, 108: 639-646. https://doi.org/10.1210/endo-108-2-639

King, D.P., Zhao, Y., Sangoram, A.M., Wilsbacher, L.D., Tanaka, M., Antoch, M.P., Steeves, T.D., Vitaterna, M.H., Kornhauser, J.M., Lowrey, P.L., Turek, F.W. and Takahashi, J.S., 1997. Positional cloning of the mouse circadian clock gene. Cell, 89: 641-653. https://doi.org/10.1016/S0092-8674(00)80245-7

Kondratov, R.V., Chernov, M.V., Kondratova, A.A., Gorbacheva, V.Y., Gudkov, A.V. and Antoch, M.P., 2003. BMAL1-dependent circadian oscillation of nuclear CLOCK: Posttranslational events induced by dimerization of transcriptional activators of the mammalian clock system. Genes Dev., 17: 1921-1932. https://doi.org/10.1101/gad.1099503

Kondratov, R.V., Shamanna, R.K., Kondratova, A.A., Gorbacheva, V.Y. and Antoch, M.P., 2006. Dual role of the CLOCK/BMAL1 circadian complex in transcriptional regulation. FASEB J., 20: 530-532. https://doi.org/10.1096/fj.05-5321fje

Kurbatova, I.V., Kolomeichuk, S.N., Topchieva, L.V., Korneva, V.A. and Nemova, N.N., 2012. Expression of the CLOCK, BMAL1, and PER1 circadian genes in human oral mucosa cells as dependent on CLOCK gene polymorphic variants. Dokl. biol. Sci., 446: 323-326. https://doi.org/10.1134/S0012496612050146

Lamia, K.A., Storch, K. and Weitz, C.J., 2008. Physiological significance of a peripheral tissue circadian clock. Proc. natl. Acad. Sci. U.S.A., 105: 15172-15177. https://doi.org/10.1073/pnas.0806717105

Leloup, J.C. and Goldbeter, A., 2003. Toward a detailed computational model for the mammalian circadian clock. Proc. natl. Acad. Sci. U.S.A., 100: 7051-7056. https://doi.org/10.1073/pnas.1132112100

Lewis, P., Korf, H.W., Kuffer, L., Gross, J.V. and Erren, T.C., 2018. Exercise time cues (zeitgebers) for human circadian systems can foster health and improve performance: A systematic review. B. M. J. Open Sport Exerc. Med., 4: e443. https://doi.org/10.1136/bmjsem-2018-000443

Li, R., Yue, J., Zhang, Y., Zhou, L., Hao, W., Yuan, J., Qiang, B., Ding, J.M., Peng, X. and Cao, J.M., 2008. CLOCK/BMAL1 regulates human nocturnin transcription through binding to the E-box of nocturnin promoter. Mol. Cell. Biochem., 317: 169-177. https://doi.org/10.1007/s11010-008-9846-x

Lincoln, G.A., 2002. Neuroendocrine regulation of seasonal gonadotrophin and prolactin rhythms: lessons from the Soay ram model. Reprod. Suppl., 59: 131-147.

Liu, T.Q., 2016. Relationship between expression and reproduction of Bmal1 and Clock in hypothalamus and testis in black hamster. Qufu Normal University. Qufu. (in Chinese).

Livak, K.J. and Schmittgen, T.D., 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta c(t)) method. Methods, 25: 402–408. https://doi.org/10.1006/meth.2001.1262

Marcheva, B., Ramsey, K.M., Buhr, E.D., Kobayashi, Y., Su, H., Ko, C.H., Ivanova, G., Omura, C., Mo, S., Vitaterna, M.H., Lopez, J.P., Philipson, L.H., Bradfield, C.A., Crosby, S.D., JeBailey, L., Wang, X., Takahashi, J.S. and Bass, J., 2010. Disruption of the clock components CLOCK and BMAL1 leads to hypoinsulinaemia and diabetes. Nature, 466: 627-631. https://doi.org/10.1038/nature09253

Moraes, M.N., de Assis, L., Magalhaes-Marques, K.K., Poletini, M.O., de Lima, L. and Castrucci, A., 2017. Melanopsin, a canonical light receptor, mediates thermal activation of clock genes. Sci. Rep., 7: 13977. https://doi.org/10.1038/s41598-017-13939-3

Morse, D., Cermakian, N., Brancorsini, S., Parvinen, M. and Sassone-Corsi, P., 2003. No circadian rhythms in testis: Period1 expression is clock independent and developmentally regulated in the mouse. Mol. Endocrinol., 17: 141–151. https://doi.org/10.1210/me.2002-0184

Nam, D., Yechoor, V.K. and Ma, K., 2016. Molecular clock integration of brown adipose tissue formation and function. Adipocyte, 5: 243-250. https://doi.org/10.1080/21623945.2015.1082015

Peruquetti, R.L., de Mateo, S. and Sassone-Corsi, P., 2012. Circadian proteins CLOCK and BMAL1 in the chromatoid body, a RNA processing granule of male germ cells. PLoS One, 7: e42695. https://doi.org/10.1371/journal.pone.0042695

Preitner, N., Damiola, F., Lopez-Molina, L., Zakany, J., Duboule, D., Albrecht, U. and Schibler, U., 2002. The orphan nuclear receptor REV-ERBalpha controls circadian transcription within the positive limb of the mammalian circadian oscillator. Cell, 110: 251-260. https://doi.org/10.1016/S0092-8674(02)00825-5

Pugazhendhi, A., Kumar, G. and Sivagurunathan, P., 2019. Microbiome involved in anaerobic hydrogen producing granules: A mini review. Biotechnol. Rep. (Amst), 21: e301. https://doi.org/10.1016/j.btre.2018.e00301

Reiter, R.J., 1980. The pineal and its hormones in the control of reproduction in mammals. Endocrinol. Rev., 1: 109-131. https://doi.org/10.1210/edrv-1-2-109

Rudic, R.D., McNamara, P., Curtis, A.M., Boston, R.C., Panda, S., Hogenesch, J.B. and Fitzgerald, G.A., 2004. BMAL1 and CLOCK, two essential components of the circadian clock, are involved in glucose homeostasis. PLoS Biol., 2: e377. https://doi.org/10.1371/journal.pbio.0020377

Schmittgen, T.D. and Livak, K.J., 2008. Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc., 3: 1101-1108. https://doi.org/10.1038/nprot.2008.73

Shimba, S., Ishii, N., Ohta, Y., Ohno, T., Watabe, Y., Hayashi, M., Wada, T., Aoyagi, T. and Tezuka, M., 2005. Brain and muscle Arnt-like protein-1 (BMAL1), a component of the molecular clock, regulates adipogenesis. Proc. natl. Acad. Sci. U. S. A., 102: 12071-12076. https://doi.org/10.1073/pnas.0502383102

Spoelstra, K., Comas, M. and Daan, S., 2014. Compression of daily activity time in mice lacking functional Per or Cry genes. Chronobiol. Int., 31: 645-654. https://doi.org/10.3109/07420528.2014.885529

Tamayo, A.G., Duong, H.A., Robles, M.S., Mann, M. and Weitz, C.J., 2015. Histone monoubiquitination by Clock-Bmal1 complex marks Per1 and Per2 genes for circadian feedback. Nature Struct. mol. Biol., 22: 759-766. https://doi.org/10.1038/nsmb.3076

Tang, J., Hu, W., Di R, Liu, Q., Wang, X., Zhang, X., Zhang, J. and Chu, M., 2018. Expression analysis of the prolific candidate genes, BMPR1B, BMP15, and GDF9 in small Tail Han ewes with three fecundity (FecB Gene) genotypes. Animals, 8: 166. https://doi.org/10.3390/ani8100166

Tei, H., Okamura, H., Shigeyoshi, Y., Fukuhara, C., Ozawa, R., Hirose, M. and Sakaki, Y., 1997. Circadian oscillation of a mammalian homologue of the Drosophila period gene. Nature, 389: 512-516. https://doi.org/10.1038/39086

Trott, A.J. and Menet, J.S., 2018. Regulation of circadian clock transcriptional output by CLOCK: BMAL1. PLoS Genet., 14: e1007156. https://doi.org/10.1371/journal.pgen.1007156

van der Horst, G.T., Muijtjens, M., Kobayashi, K., Takano, R., Kanno, S., Takao, M., de Wit, J., Verkerk, A., Eker, A.P., van Leenen, D., Buijs, R., Bootsma, D., Hoeijmakers, J.H. and Yasui, A., 1999. Mammalian Cry1 and Cry2 are essential for maintenance of circadian rhythms. Nature, 398 :627-630. https://doi.org/10.1038/19323

von Gall, C., Weaver, D.R., Kock, M., Korf, H.W. and Stehle, J.H., 2000. Melatonin limits transcriptional impact of phosphoCREB in the mouse SCN via the Mel1a receptor. Neuroreport, 11:1803-1807. https://doi.org/10.1097/00001756-200006260-00002

Ye, Y., Xiang, Y., Ozguc, F.M., Kim, Y., Liu, C.J., Park, P.K., Hu, Q., Diao, L., Lou, Y., Lin, C., Guo, A.Y., Zhou, B., Wang, L., Chen, Z., Takahashi, J.S., Mills, G.B., Yoo, S.H. and Han, L., 2018. The genomic landscape and pharmacogenomic interactions of clock genes in cancer chronotherapy. Cell Systems, 6: 314-328. https://doi.org/10.1016/j.cels.2018.01.013

Zheng, X., Zhao, X., Zhang, Y., Tan, H., Qiu, B., Ma, T., Zeng, J., Tao, D., Liu, Y., Lu, Y. and Ma, Y., 2019. RAE1 promotes BMAL1 shuttling and regulates degradation and activity of CLOCK: BMAL1 heterodimer. Cell Death Dis., 10: 62. https://doi.org/10.1038/s41419-019-1346-2

Zhou L., Yu Y.M., Sun S.W., Zhang T.Q. and Wang M.H., 2018. Cry 1 regulates the clock gene network and promotes proliferation and migration via the akt/p53/p21 pathway in human osteosarcoma cells. J. Cancer, 9: 2480-2491. https://doi.org/10.7150/jca.25213