{"paper_id":"2de01d82-1774-4983-80d9-ba90a25a797e","body_text":"The effects of estrus cycle on the expression of ovarian biological clock-related genes of Xiang pigs | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article The effects of estrus cycle on the expression of ovarian biological clock-related genes of Xiang pigs Liang-Ting Tang, Shihui Huang, Xi Niu, Sheng Li, Jiafu Wang, Xue-Qin Ran, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-2593341/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract The biological clock has been studied to play a critical role in the reproductive system of various living organisms like swine. To examine the effects of estrus cycle on the expression of ovarian biological clock-related gene in Xiang pig, in this study, we analyzed the expression and alternative splicing of biological clock-related genes in ovary during estrus and diestrus periods. In total, we detected 90 clock-related genes expressed in the ovaries of the Xiang pigs and found 33 clock-related genes differentially expressed between estrous and diestrous stages. We identified 44 differential splicing events from the transcripts of 34 biological clock-related genes. Furthermore, we also found 20 genes including the core clock components, arntl and cry1 were differentially regulated only at AS level and 14 genes, including per1 and clock, were differentially regulated at both expression and AS levels. We also proved that the core clock genes per1, cry1, clock and arntl and the clock-related genes, ppp1cb and ntrk1 were rhythmically expressed in Xiang pig ovaries by RT-qPCR experiments. The results demonstrated that the biological clock in the ovaries of Xiang pigs might play an important role in regulating the ovarian physiological functions by the transcriptional and post-transcriptional regulation. Circadian rhythm Xiang pig Ovary Estrus cycle Alternative splicing Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction Circadian rhythm is the regular change of many functional activities of the biological organism to adapt to the alternation of day and night in the environment. It is endogenous and persists even without external synchronising factors, for example, in constant darkness. Circadian rhythm is controlled by central biological clocks related to hypothalamic suprachiasmatic nucleus (SCN) and peripheral biological clocks distributed in various tissues (Schibler, et al. 2003 ). SCN cells modulate the activity of gonadotropin-releasing hormone (GnRH) neurons by direct innervation or other hypothalamic nucleus, which are necessary for the occurrence of ovulation (Vieyra, et al. 2016 ; Mansano, et al. 2022 ). It is known that the molecular model of mammals circadian clock including human is a multi-feedback loop for transcription and translation including transcriptional regulation and post-transcriptional regulation (Takahashi 2017 ). It consists of positive regulatory components (BMAL1 and CLOCK) and negative regulatory components (PER 1/2 and CRY 1/2) (Reppert and Weaver 2002 ). The positive and the negative regulatory elements also participate in the circadian regulation of clock controlled genes (CCG). The clock-controlled genes act as specific regulators of rhythmic physiological functions in cells and tissues (Sen and Sellix 2016 ). For example, the function of ovary is regulated by the biological clock, including the estrus cycle, luteinizing hormone (LH) levels, ovulation and embryo implantation (Gray, et al. 1978 ; Christian, et al. 2005 ). Ovary is an important reproductive organ of female mammals. In each estrus cycle, the ovaries undergo a series of morphological, hormonal and biochemical changes, which directly affect and/or determine the litter size of female animals (Miao and Luo 2013 ). The realization of ovarian function is a very complex biological process, involving the transcriptional regulation of a large number of genes (Terenina, et al. 2017 ). Alternative splicing (AS) is also an important way of gene expression regulation, which occurs widely in various tissues and cells of almost all organisms (Chacko and Ranganathan 2009 ; Kianianmomeni, et al. 2014 ). Many studies demonstrate that the biological clock existed in the ovary (Fahrenkrug, et al. 2006 ; Karman and Tischkau 2006 ; He, et al. 2007b ; Yoshikawa, et al. 2009 ). The expression of clock genes have been described in the ovary of rat (He, et al. 2007a ), mice (Johnson, et al. 2002 ) and ruminant (Cushman, et al. 2007 ). It indicates that clock genes are involved in the process of follicle formation, ovulation, oocyte maturation and steroid hormone synthesis (Sellix 2015 ). The lack of circadian clock might be related with varieties of reproductive disorders. For instance, mice with per1 and per2 knockout present a significant decrease in reproductive rate due to the irregular and non-cyclical estrus cycle (Zheng, et al. 2001 ; Pilorz and Steinlechner 2008 ). Although mice with bmal1 knockout are capable of ovulation, they exhibit delayed puberty, irregular estrus cycles, and small ovaries and uterus (Alvarez, et al. 2008 ; Boden, et al. 2010 ; Liu, et al. 2014 ). Loss of bmal1 also negatively affects progesterone level (Liu, et al. 2014 ), leading to embryo implantation failure (Boden, et al. 2010 ). Mice mutant with 51 amino acid deletions in the transcription activation domain of clock protein show an irregular estrus cycle, and the normal LH surges disappear on the day of estrus (Miller, et al. 2004 ). In addition, genes of ovarian granulosa cells that are important for follicular development have been shown to be controlled by clock genes, such as LH receptor ( lhcgr ), prostaglandin synthase ( ptgs2 ), steroid-producing enzymes (steroid-producing acute regulatory protein, cyp11a1, aromatase, etc.), and gap junction protein Connexin-43 (Chen, et al. 2012 ; Chu, et al. 2012 ; Chen, et al. 2013 ). Estrus is a sign of reproductive maturity. Moreover, it has been confirmed that the central biological clock located in the SCN was the core of maintaining normal reproductive function and the estrus cycle (Brown-Grant and Raisman 1977 ; Wiegand, et al. 1980 ). However, the function of ovarian biological clock during estrus period is still needed to be clear. The Xiang pig is a small local breed in China with some unique biological characteristics, such as small size, precocious puberty, low birth rate and insignificant estrous behavior (Luo, et al. 2018 ). In order to understand more in-depth views of the influence of estrus on biological clock and further explore the physiological functions of clock genes in porcine ovary, we analyzed the expression of biological clock-related genes and alternative splicing of transcripts in the ovaries of Xiang pigs during estrus and diestrus period based on the whole transcriptome sequencing analysis. Materials And Methods Materials About 40 sows of Xiang pig were selected for observation of estrus performance. They were reared in a standard manner that the pigs drank freely and fed timely under the normal condition. The average light intensity was about 200 lux, the ambient temperature was maintained at 15°C -18°C, and the humidity was in the range of 50% ± 10%. The estrus detection was carried out 2 to 3 times a day. Redness and swelling of the vulva and a positive standing reflex test of the boar were defined as estrus manifestations. Ovary samples of 14 Xiang pig sows were taken, 7 of them were in the diestrus period, and the others were subjected to routine castration surgery on the third day of estrus (the sows showed strong redness and swelling of the vulva). The sampling day was at 8–9 am at March 15 to 30, 2022, with the same light and temperature conditions as previously mentioned, and fasting before ovariectomy surgery. All samples were immediately immersed in liquid nitrogen and stored at -80°C for RNA extraction. The animal procedures were approved by the guidelines of Guizhou University Subcommittee of Experimental Animal Ethics with No. of EAE-GZU-2020-P002. Total RNA extraction Total RNA from 14 ovarian samples were extracted using TRIzol reagent kit (Invitrogen, USA) according to the instructions. The quantity and completeness of total RNA were evaluated by the Bio-analyzer system (Agilent 2100, Agilent Technologies, USA). The total RNA was stored at -80°C for subsequent analysis. The same samples were used for sequencing and PCR analysis. Library construction and sequencing Total of 14 cDNA libraries were constructed using RNA as templates from ovarian tissues. The library construction was carried out according to the standard procedures of BGI (Shenzhen Huada Gene Co., Ltd). After DNase I enzyme treatment, magnetic beads with oligo (dT) were used to isolate mRNA. Under appropriate temperature conditions, the purified mRNA was fragmented by adding fragmentation buffer, and then cDNA was synthesized using random primers and mRNA fragment templates. After 3'end adenylation, the DNA fragments were ligated to the linker adaptors. cDNA fragments of 100 ~ 200 bp in length were selected for PCR amplification to generate a cDNA library. The libraries were sequenced on the Illumina HiSeqTM 2500 sequencing platform, and paired-end sequences with 150 bp in length were generated. Sequencing data analysis The original sequences generated by the Illumina platform were saved in fastq format. The clean reads were obtained by removing low-quality sequences, reads (more than 5% unknown nucleotides), and sequencing adapters. The clean reads should be evaluated through FastaQC software (Q20 < 20%). The reference genome sequence (ssc11.1 version) and annotation files of pig genomes were downloaded from the Ensembl database ( http://www.ensembl.org/index.html ). Used HISAT2 software (v2.1.0) to map the clean reads to the pig reference genome. HISAT2 is an efficient alignment tool for RNA-seq experiments. It uses an indexing framework based on the BWT and Ferragina-Manzini (FM) index algorithms for alignment and data normalization.The Subread featureCounts software (v2.0.0) was used to count the number of reads mapped to the gene. Subread featureCounts applies the read summarization tool to analyze the expression level from RNA-seq, which applies efficient chromosome hash algorithm and feature block technology for statistics and data standardization. Analysis of differential expression The expressed level of each gene were calculated by CPM values (counts per million mapped reads) in all samples, and added 0.001 to CPM value as the gene expression level. Then we used bioconductor software package Limma and DESeq2 to analyze differentially expressed genes (DEGs) between two groups with biological replicates. The FDR ≤ 0.1 and |log2FC|≥0.5 were used as the thresholds for determining the significance of gene expression difference. Identification of differential alternative splicing events The rMATS software (v4.0.2) was used to analyze differentially alternative splicing events between two groups. The rMATS recognizes five alternative splicing events, which are skipped exon (SE), alternative 5’ splice site (A5SS), alternative 3’ splice site (A3SS), mutually exclusive exons (MXE) and retained intron (RI) (Shen, et al. 2014 ). Also, FDR ≤ 0.05 was used as a threshold to detect the differential significance of alternative splicing. Validation The same RNA samples from the ovarian for RNA-seq experiment were used to validate the differential expression by qRT-PCR and to validate the AS events by RT-PCR method. The primers were designed using Primer 6.0, which was listed in table S4. GAPDH and β-actin genes were taken as internal controls. The PCR efficiency of the primers was controlled within 100 ± 10%. All qPCR reactions were performed on the CFX98TM real-time system (Bio-rad co.) according to the manufacturer's instructions and recommended cycling conditions. Three replicate reaction tests were performed using a volume of 10 µL (containing 5 µL of 2× Super Real Pre Mix Plus), 0.3 µL of each of the forward and reverse primers (10 pM / µL) and 1 µL cDNA. The qPCR conditions were as following: initial denaturation at 95°C for 3 minutes, then denaturation at 94°C for 10 seconds in 40 cycles, and annealing and extension at 60°C for 30 seconds. The relative expression level of target gene utilized the method of 2(-ΔΔCt) as reported by Livak et al. (Livak and Schmittgen, 2001). The different level of gene expression between two groups was tested by software SPSS (v21.0) taking the P < 0.05 as threshold of significant difference. The results were presented as mean ± standard deviation. The 8 AS events were randomly selected from 8 genes to validate by RT-PCR method. The RT-PCR annealing temperature was 55°C for 30 seconds, and the extension temperature was 72°C for 60 seconds by 34 cycles. To verify the rhythmic expression of circadian rhythm-related genes in Xiang pig ovary, six genes, per1 , cry1 , clock , arntl , ppp1cb and ntrk1 , were selected to detected quantitatively by RT-qPCR method at six time points (08:00, 12:00, 16:00, 20:00, 24:00, 04:00). Three ovary samples of Xiang pigs were collected at each time point, and RNA was extracted from the ovarian samples for qPCR verification of genes. The amplification conditions were the same as the previous RT-qPCR method. Results Summary of sequencing data The 14 cDNA libraries from ovaries in estrus and diestrus stages were constructed and sequenced respectively, and the sequencing data were filtered and normalized, which were then mapped to the reference genome with HISAT2. It generated 75002160–77100634 read pairs, and 93.57%-96.26% of reads were aligned to the reference genome (Table 1 ). The only matching reads were performed the following statistical analysis. Table 1 RNA-Seq data description Samples Input Read Pairs Both Surviving Forward Only Surviving Reverse Only Surviving Dropped Percentage of mapped reads (%) Diestrus D1 77089626 74382355 (96.49%) 2471254 (3.21%) 164408 (0.21%) 69380 (0.09%) 96.01% D2 77040960 74260107 (96.39%) 2550481 (3.31%) 160145 (0.21%) 70227 (0.09) 95.88% D3 76170080 73520243 (96.52%) 2408511 (3.16%) 167409 (0.22%) 73917 (0.10%) 95.93% D4 77100634 74147772 (96.17%) 2713372 (3.52%) 165857 (0.22%) 73633 (0.10%) 96.26% D5 76860177 74164234 (96.49%) 2449821 (3.19%) 177911 (0.23%) 68211 (0.09) 96.06% D6 76903234 74486606 (96.86%) 2179432 (2.83%) 173980 (0.23%) 63216 (0.08%) 96.11% D7 77054008 74411266 (96.57%) 2408339 (3.13%) 163511 (0.21%) 70892 (0.09%) 95.95% Estrus E1 75005968 74960964 (99.94%) 45004 (0.06%) 0 (0.00%) 0 (0.00%) 95.98% E2 74990854 74960858 (99.96%) 29996 (0.04%) 0 (0.00%) 0 (0.00%) 95.90% E3 75206715 75176632 (99.96%) 30083 (0.04%) 0 (0.00%) 0 (0.00%) 95.93% E4 75187643 75135012 (99.93%) 52631 (0.07%) 31 (0.00%) 0 (0.00%) 96.24% E5 75184086 75154012 (99.96%) 30074 (0.04%) 31 (0.00%) 1 (0.00%) 96.15% E6 76109159 74099877 (97.36%) 1682012 (2.21%) 243549 (0.32%) 83720 (0.11%) 94.13% E7 75002160 73149607 (97.53%) 1545044 (2.06%) 225006 (0.30%) 82502 (0.11%) 93.57% Expression of biological clock-related genes in Xiang pig ovary After mapping and annotating based on the pig reference genome assembly Sscrofa11.1, the RNA-seq data of circadian rhythms and clock genes were picked up from the genome-scale data sets (Table S1). Taking threshold value of 0.1 CPM, 90 rhythms genes were determined, in which 81 were expressed in the estrus ovaries and 83 were detected from the diestrus ovaries. The expression levels of those genes were ranged from 867.2 CPM in the highest abundant transcripts to the cut off value of 0.1 CPM (Table S1). Compared with the diestrus ovary, there were more genes (55.6% in estrus, 42.2% in diestrus) at the medium expression level (10–100 CPM), and fewer genes at the high (> 100 CPM) or low (0.1–10 CPM) expression level in the estrus ovary (Fig. 1 ). The top expressed genes were nono and sfpq with 455.6 and 446.1 CPM in the estrus ovary, and sfpq and top2a with 532.6 and 512.1 CPM in the diestrus ovary, respectively. Analysis of differential expression genes We used R package Limma and DESeq2 to analyze differentially expressed genes (DEGs) between estrus and diestrus ovaries. Also, the P-value ≤ 0.01 and absolute value of log 2 FC ≥ 0.5 were used as the threshold for determining the significant difference of gene expression. The results showed that 33 DEGs in the ovaries between estrus and diestrus cycles, comprising 13 down-regulated genes ( top2a, nrip1, prkdc, setx, ncoa2, pparg, pspc1, per3, sirt1, timeless, ezh2, ppara , and ppp1ca ) and 8 up-regulated genes (atf4, csnk1d, klf10, id2, bhlhe40, id4, ass1 , and ntrk1 ) in at least 2-fold at estrus stage. The fold changes for DEGs were ranged from 0.17 to 26. Furthermore, we found that 5 genes ( ppp1cb, id3, per1, klf9 , and nampt ) were at least 1.5-2 fold upregulated and 7 genes ( ep300, ube3a, dyrk1a, sin3a, clock, creb1 , and magel2 ) were at least 1.5-2 fold down regulated in estrus ovary (Table S2). A lot of the genes were expressed in ovaries at medium or high levels (48.9% and 27.4% ), including the core clock components, arntl, clock, per1/2/3 , and cry1/2 . Detection of differential splicing genes The rMATS was used to analyze differentially alternative splicing genes (DSGs) between estrus and diestrus ovaries. Overall, 44 differential splicing events were identified from the transcripts of 34 circadian rhythms genes (P < 0.05) including 3 A5SS, 12 RI, 14 SE and 15 MXE events (Table S3). There were multiple genes with complex alternative splicing patterns. Of those, nono presented the most abundant differential alternative splicing events, which included four differential splicing types: A5SS, RI, SE, and MXE. Also, ppp1cb harbored three differential splicing types (RI, SE and MXE). Then, ENSSSCG00000000107, per1, clock, ube3a and fbxw11 each contained two differential alternative splicing types. It was found that 20 genes comprised the core clock components, arntl and cry1 were differentially regulated only at AS level. At least 14 genes, such as per1 and clock , were differentially regulated in both expression and AS level in ovaries between estrus and diestrus. Validation of DEGs by RT-qPCR method To validate the results from RNA-seq, we used RT-qPCR method to further examined the differentially expressed genes using the same RNA samples from the ovaries (Fig. 2 ). It confirmed that ppp1cb, csnk1d, per1, klf9, nampt, bhlhe40, id4 , and ass1 were significantly high expressed in estrus samples, which were in accordance with the RNA-seq analysis (Table S2). The expression of atf4 and ntrk1 in estrus and diestrus ovaries is not significantly different. Moreover, 8 AS events, selected randomly from 8 genes, were detected by RT-PCR and consistent with the RNA-Seq results (Figure S1).The primer sequences were shown in Table S4 Validation of rhythmic expression of clock-related genes To verify whether the expression of clock-related genes in Xiang pig ovaries was rhythmic, we further examined the expression patterns of per1, cry1, clock, arntl, ppp1cb and ntrk1 at six time points (08:00, 12:00, 16:00, 20:00, 24:00, 04:00) using RT-qPCR. The samples at each time point were obtained from 3 Xiang pig ovaries collected within one hour. Among them, the expression of per1 and ppp1cb peaked at 24:00 midnight. However, the expression patterns of clock and arntl were opposite, in which the peaks presented at 8:00 am and it fell to the lowest at 24:00. The maximum values was appeared at 20:00 and 24:00 pm for gene cry1 , and it was at 16:00 pm for gene ntrk1 . It illustrated that the core clock genes per1 , cry1 , clock and arntl and the clock-related genes ppp1cb and ntrk1 were rhythmic expressed in Xiang pig ovaries (Fig. 3 ). Discussion The mammalian biological clock system has a multi-level structure, including the main clock located in the suprachiasmatic nucleus and the sub-clocks of peripheral organs and tissues (Schibler, et al. 2003 ). The ovary is a pivotal reproductive organ in female animal, and the expression of the biological clock gene has been found in the ovaries of many organisms (Sellix 2015 ). In this study, we analyzed the expression profile and alternative splicing of the biological clock gene in Xiang pig ovaries in estrous cycle from the RNA-seq data. The study found out a total of 90 rhythms genes that were expressed in ovaries. The expression abundance of these genes was ranging from 0.1 to 532.6 CPM. A lot of the genes were expressed in ovaries at medium or high levels (48.9% and 27.4%), including the core clock components, arntl, clock, per1/2/3 , and cry1/2 . Furthermore, 33 genes were detected to undergo differentially expression and 34 genes were detected to undergo differential alternative splicing between estrous and diestrous ovaries. The DEGs and DSGs related with rhythms might have a connection with the regulation of estrus process in Xiang pig. Transcription and translation of core clock components genes ( clock or npas2, arntl/bmal1 , or arntl2/bmal2, per1/2/3 , and cry1/2 ) play a critical role in rhythm generation process. CLOCK and BMAL1 heterodimerize to activate transcription of circadian target genes including per1/2/3 and cry1/2 . PER and CRY interact and conversely inhibit transcription of bmal1 and clock genes. These genes and their protein products are organized into interlocking positive and negative transcriptional and translational feedback loops, which regulates circadian rhythm generation in the brain suprachiasmatic nucleus (SCN) and peripheral organs (Rosenwasser and Turek 2015 ). In this study, we found that three core clock genes were differentially expressed between estrus and diestrus ovaries. Clock and per3 were down-regulated and per1 was up-regulated in the estrous ovaries. CLOCK plays a key role in maintaining the circadian rhythm and activating downstream elements. Inhibition of CLOCK can inhibit cell growth and increase the rate of apoptosis (Li, et al. 2015 ). In vivo experiments show that female mice is injected with CLOCK-shRNA have fewer oocytes, fewer litters, and a higher rate of apoptosis. The results indicate that clock plays an important role in fertility, and down-regulation of clock leads to cell apoptosis and decrease reproductive capacity (Li, et al. 2015 ). It was proposed that the down-regulation of clock in the ovary at estrus stage might be related to the low litter size of the Xiang pig. The core clock gene period 1 ( per1 ) may be a prolific gene in Drosophila (Cushman, et al. 2007 ). Female mice with the per1 mutation showed a normal number of implantation sites but reduced litter size (Pilorz and Steinlechner 2008 ). Per1 mRNA locates in the secondary oocytes and follicles of ruminants (sheep, cattle) and found that there was no relationship between its transcription level and prolificacy, and this gene did not map to the known QTL region of ovulation rate in cattle (Cushman, et al. 2007 ). Treatment with Progesterone for 1 hour could induce the expression of per1 mRNA in MCF-7 cells (Nakamura, et al. 2010 ). In fact, the up-regulation of per1 was a common feature of many tissues in response to certain types of hormonal stress(For xample: luteinising hormone) (Mansuy, et al. 2009 ; Mendoza, et al. 2011 ; Smarr, et al. 2013 ). This expression pattern indicated that PER1 may be the most sensitive effector in the biological clock system (Zhang, et al. 2019 ). Some studies have found that estrus leads to changes in the expression time and amplitude of the biological clock genes (López-Rodríguez, et al. 2019 ). Rising estrogen cause female animal to go into estrous stage (Ravinder, et al. 2016 ).The estrus in Xiang pig might be speculated to reset the biological clock by up-regulating the per1 gene to response the stimulus of steroid hormone.In this study, per3 was down-regulated in the ovaries of estrus pigs. There is little information on the function of per3 gene in reproduction except for brain development (Noda, et al. 2019 ). In contrast, the expression patterns of clock genes in other pig breeds are much different from that in Xiang pig. In the ovarian follicles of Large White pigs and Mi pigs, the core circadian clock genes were all down-regulated and not differentially expressed between estrus and diestrus periods (Chu, et al. 2017 ). Compared with Large White pigs and Mi pigs, the characteristics of Xiang pig are specific, such as low birth rate and insignificant estrous behavior. It would be interesting to prove whether the expression profiles of core biological clock genes in ovary is a reason for the special reproductive traits such as estrus performance and litter size in different pig breeds. The biological clock system consists of an input pathway, a core oscillator and an output pathway. The post-translational modification of clock proteins is essential to maintain the accuracy and robustness of the evolutionarily conserved circadian clock (Gossan, et al. 2014 ). Post-translational modification and degradation of the clock proteins are key steps in determining the length of the circadian clock cycle (Kennaway, et al. 2015 ). Our research found that the products of most differentially expressed genes were related to the post-translational modification of the core clock protein. For example, PPP1CB and PPP1CA can reduce the phosphorylation of PER2, and affect the nuclear localization of the protein PER2, which may at least partially change the cycle and phase shift characteristics of the biological clock (Schmutz, et al. 2011 ). Also, highlighted DYRK1A is the enzyme responsible for the phosphorylation of CRY2 at Ser557 and plays a key role in regulating the protein level of CRY2 (Kurabayashi, et al. 2010 ). The casein kinases CSNK1D and CSNK1E phosphorylate the PER protein and provide a marker for subsequent degradation (Kennaway, et al. 2015 ). The UBE3A binds and degrades BMAL1 in a ubiquitin ligase-dependent manner (Gossan, et al. 2014 ). Furthermore, MAGEL2 regulates the ubiquitination and stability of CRY1, and changes its nuclear and cytoplasmic distribution (Devos J 2011). Lastly, SIRT1 deacetylates PER2 and BMAL1, thereby participates in biological rhythm regulation (Wang, et al. 2016 ). The previous reports illustrate that the estrus cycle affects the localization and degradation of the core clock protein by changing their post-translational modification. And transcription and post-transcriptional regulation are the basis of clock system component activities, and post-transcriptional mechanisms accounts for more than half of the regulatory network (Koike, et al. 2012 ). The results of this study further strengthened the view that the post-translational mechanism participated in the circadian gene regulatory network. In addition, we detected some differentially expressed genes from Xiang pig ovaries (Table S2) that might affect the expression of core clock genes at the transcriptional level. For example, TOP2A binds to the unique GC-rich open chromatin structure of the bmal1 promoter region, indicating that TOP2 on the bmal1 promoter affects transcription of bmal1 (Ogawa, et al. 2014 ). Also, CBP /P300 and tissue-specific cofactors regulate CLOCK /BMAL1 transcription positively or negatively (Hosoda, et al. 2009 ). Accordingly, ID (DNA binding inhibitor) is an important transcription repressor. Each ID protein contains a helix-loop-helix domain through which it can interact with bHLH protein. However, ID lacks the basic domain that allows binding to the E-box element, which leads to the ID being able to modify the transactivation of clock genes and CCG by interfering with the CLOCK-BMAL1 heterodimer, bHLH orange factor and other bHLH factors (Hou, et al. 2009 ). These data indicate that the estrus cycle promoted the time change of the biological clock in reproductive tissues at the level of transcription and post-translational modifications (Fig. 4 ). In addition, we detected several genes controlled by the circadian rhythm, which were differentially expressed in the ovaries of the two periods including klf9, star, ptgs2 . Other studies have shown that klf9 is a clock output gene, and CLOCK and BMAL1 complexes can bind and activate transcription at the 5' flanking region of klf9 , which is blocked by the co-expression of PER1. Also, KLF9 may play a role in regulating the effects of CLOCK/BMAL1 and in the expression of DBP, other clock and clock output genes, thereby changing the timing and amplitude of the circadian oscillations of gene transcription (Knoedler, et al. 2020 ). The previous studies have found that steroidogenic enzymes (STAR), prostaglandin synthase (PTGS2), which are related to ovarian progesterone synthesis, are clock-controlled genes (CCGs). The E-box enhancer exists in the 5'-flanking region of star, which binds to the CLOCK-BMAL1 heterodimer and activates the transcription of gene star (Nakao, et al. 2007 ). Co-expression of the negative regulators PER and CRY attenuates this activation (Wu, et al. 2011 ). The promoter of ptgs2 has an E-box element and a REV-ERBα/RORα response element (RORE). The secretion of PGF2α can be balanced by the inhibition or stimulation of transcriptional regulation on REV-ERBα and BMAL1/CLOCK, respectively (Isayama, et al. 2014 ). It is thus inferred that the peripheral circadian oscillator, such as ovary, could play an essential role in synchronizing local physiology through regulation of the expression of clock-controlled genes (Mohawk, et al. 2012 ). In this study, star and ptgs2 in the ovaries of the estrus pigs were up-regulated by 6.42 and 5.55 times, respectively. It might be possible that the core clock genes might regulate the production of steroid hormones by controlling the ovarian-specific CCG, thereby affecting the estrus cycle of Xiang pigs. Pre-mRNA splicing is a basic biological process through which introns of nascent RNA are removed and exons are merged to form mature RNA, which is then translated into protein (Lee and Rio 2015 ). Through alternative splicing (AS), different proteins are produced, resulting in a wider variety of cellular functions (Genov, et al. 2019 ). Except for transcription or post-translational controls, alternative splicing also plays a key role in circadian rhythms in a cell-type-specific manner (Shakhmantsir, et al. 2018 ). Recent research reports several cases, in which AS is involved in the regulation of biological clocks in plants, mice, and fruit flies (Sanchez, et al. 2010 ; Petrillo, et al. 2011 ; Hughes, et al. 2012 ; McGlincy, et al. 2012 ; Preußner, et al. 2014 ). In this study, Arntl only underwent differential alternative splicing, while clock and per1 experienced both of differential alternative splicing and differential expression. The differential alternative splicing event occurred in per1 was the retention of introns 14–15. This event would result in the deletion of amino acid residues from positions 545 to 1278 in the encoded protein, and the deletion site was located in the core domain to activate the protein. This intron retention event occurred more frequently in the ovaries of Xiang pigs in the diestrus period, and the active PER1 protein decreased. During the estrus period, the occurrence of the intron retention event decreased, and the active PER1 protein level increased, which further increased the expression of PER1. It indicated that the core clock component PER1 may respond to the changes of sterol hormones in the ovary through transcriptional regulation and post-transcriptional regulation. Clock underwent the differential exon 17 skipping event. When this event occurred, the encoded protein would be deleted amino acid residues from positions 484 to 513. The deleted amino acids were not located in the core domain and the influence on the protein activity was still unclear. These findings suggested that AS might directly change the structure of core clock components to regulate the ovarian biological clock, and the AS regulation of biological clock genes might be independent in the gene expression. However, more evidence needs to be accumulated whether the AS events in these biological clock genes would turn on and turn off the circadian rhythm in pig. Conclusion In this study, we examined 90 clock-related genes expressed in Xiang pigs’ ovaries and found that 33 clock-related genes were differentially expressed between estrus and diestrous stages. We identified 44 differential splicing events from transcripts of 34 circadian clock-related genes. In addition, we also found that 20 genes including core clock components arntl and cry1 were differentially regulated at AS level only, and 14 genes including per1 and clock were differentially regulated at both expression and AS level. We also demonstrated by RT-qPCR experiments that the core clock genes per1 , cry, clock and arntl and the clock-related genes ppp1cb and ntrk1 were rhythmically expressed in ovaries of Xiang pig. The results suggested that the circadian clock in Xiang pig ovaries might play an key role in regulating ovarian physiological functions through transcriptional and post-transcriptional regulation. Declarations Author Contribution LTT, JFW, and XQR conceived and designed the experiments; LTT performed the experiments; LTT and XN analyzed the data; LTT, XN, SL, and SHH contributed to samples and data analysis; and LTT wrote the paper. All authors reviewed the manuscript. Ethics Approval All animal procedures were approved by the Institutional Animal Care and Use Committee of Guizhou University (Approval number GZU-201709). The procedures were conducted in accordance with the National Research Council Guide for the Care and Use of Laboratory Animals. Funding Statement This work was supported by the National Natural Science Foundation of China (31960641, 31672390), the Guizhou Province Science and Technology Innovation Team Building Special QKHPTRC[2019]5615, the Guizhou Province “Hundred” Innovative Talents Project [2016-4012], the Guizhou Agriculture Research program (QKHZC[2017]2585, QKHZC[2017]2587, QKHNY[2013]3073), and the National High Technology Research and Development Program of China (863 Program) [2013AA102503]. The funders monitored the ethics of experiments and provided instruments and financial funds for all of research in the manuscript. Conflicts of Interest All the authors have no conflict of interest with anyone about this manuscript. We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted. No conflict of interest exits in the submission of this manuscript, and manuscript is approved by all authors for publication. The work described has not been published previously, and not been submitted to any other journal, in whole or in part. All the authors listed have approved the manuscript that is enclosed. References Alvarez JD, Hansen A, Ord T, Bebas P, Chappell PE, Giebultowicz JM, Williams C, Moss S, 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 Boden MJ, Varcoe TJ, Voultsios A, Kennaway DJ (2010) Reproductive biology of female Bmal1 null mice. Reproduction 139:1077-1090.https://doi.org/10.1530/REP-09-0523 Brown-Grant K, Raisman G (1977) Abnormalities in reproductive function associated with the destruction of the suprachiasmatic nuclei in female rats. Proc R Soc Lond B Biol Sci 198:279-296. Chacko E, Ranganathan S (2009) Genome-wide analysis of alternative splicing in cow: implications in bovine as a model for human diseases. BMC Genomics 10 Suppl 3:S11. https://doi.org/10.1186/1471-2164-10-S3-S11 Chen H, Chu G, Zhao L, Yamauchi N, Shigeyoshi Y, Hashimoto S, Hattori MA (2012) Rev-erbα regulates circadian rhythms and StAR expression in rat granulosa cells as identified by the agonist GSK4112. Biochem Biophys Res Commun 420:374-379. https://doi.org/10.1016/j.bbrc.2012.02.164 Chen H, Zhao L, Chu G, Kito G, Yamauchi N, Shigeyoshi Y, Hashimoto S, Hattori MA (2013) FSH induces the development of circadian clockwork in rat granulosa cells via a gap junction protein Cx43-dependent pathway. Am J Physiol Endocrinol Metab 304:E566-575. https://doi.org/10.1152/ajpendo.00432.2012 Christian CA, Mobley JL, Moenter SM (2005) Diurnal and estradiol-dependent changes in gonadotropin-releasing hormone neuron firing activity. Proc Natl Acad Sci U S A 102:15682-15687. https://doi.org/10.1073/pnas.0504270102 Chu G, Misawa I, Chen H, Yamauchi N, Shigeyoshi Y, Hashimoto S, Hattori MA (2012) Contribution of FSH and triiodothyronine to the development of circadian clocks during granulosa cell maturation. Am J Physiol Endocrinol Metab 302:E645-653. https://doi.org/10.1152/ajpendo.00470.2011 Chu Q, Zhou B, Xu F, Chen R, Shen C, Liang T, Li Y, Schinckel AP (2017) Genome-wide differential mRNA expression profiles in follicles of two breeds and at two stages of estrus cycle of gilts. Sci Rep 7:5052. https://doi.org/10.1038/s41598-017-04336-x Cushman RA, Allan MF, Jones SA, Rupp GP, Echternkamp SE (2007) Localization of Period 1 mRNA in the ruminant oocyte and investigations of its role in ovarian function. Anim Reprod Sci 99:93-105. https://doi.org/10.1016/j.anireprosci.2006.04.057 Devos J WSV, Wevrick R (2011) Magel2, a Prader-Willi syndrome candidate gene, modulates the activities of circadian rhythm proteins in cultured cells. J Circadian Rhythms 9:12. https://doi.org/10.1186/1740-3391-9-12 Fahrenkrug J, Georg B, Hannibal J, Hindersson P, Gräs S (2006) Diurnal rhythmicity of the clock genes Per1 and Per2 in the rat ovary. Endocrinology 147:3769-3776. https://doi.org/10.1210/en.2006-0305 Genov N, Basti A, Abreu M, Astaburuaga R, Relógio A (2019) A bioinformatic analysis identifies circadian expression of splicing factors and time-dependent alternative splicing events in the HD-MY-Z cell line. Sci Rep 9:11062. https://doi.org/10.1038/s41598-019-47343-w Gossan NC, Zhang F, Guo B, Jin D, Yoshitane H, Yao A, Glossop N, Zhang YQ, Fukada Y, Meng QJ (2014) The E3 ubiquitin ligase UBE3A is an integral component of the molecular circadian clock through regulating the BMAL1 transcription factor. Nucleic Acids Res 42:5765-5775. https://doi.org/10.1093/nar/gku225 Gray GD, Söderstein P, Tallentire D, Davidson JM (1978) Effects of lesions in various structures of the suprachiasmatic-preoptic region on LH regulation and sexual behavior in female rats. Neuroendocrinology 25:174-191. https://doi.org/10.1159/000122739 He PJ, Hirata M, Yamauchi N, Hashimoto S, Hattori MA (2007a) The disruption of circadian clockwork in differentiating cells from rat reproductive tissues as identified by in vitro real-time monitoring system. J Endocrinol 193:413-420. https://doi.org/10.1677/JOE-07-0044 He PJ, Hirata M, Yamauchi N, Hashimoto S, Hattori MA (2007b) Gonadotropic regulation of circadian clockwork in rat granulosa cells. Mol Cell Biochem 302:111-118. https://doi.org/10.1007/s11010-007-9432-7 Hosoda H, Kato K, Asano H, Ito M, Kato H, Iwamoto T, Suzuki A, Masushige S, Kida S (2009) CBP/p300 is a cell type-specific modulator of CLOCK/BMAL1-mediated transcription. Mol Brain 2:34. https://doi.org/10.1186/1756-6606-2-34 Hou TY, Ward SM, Murad JM, Watson NP, Israel MA, Duffield GE (2009) ID2 (inhibitor of DNA binding 2) is a rhythmically expressed transcriptional repressor required for circadian clock output in mouse liver. J Biol Chem 284:31735-31745. https://doi.org/10.1074/jbc.M109.013961 Hughes ME, Grant GR, Paquin C, Qian J, Nitabach MN (2012) Deep sequencing the circadian and diurnal transcriptome of Drosophila brain. Genome Res 22:1266-1281. https://doi.org/10.1101/gr.128876.111 Isayama K, Chen H, Yamauchi N, Hattori MA (2014) REV-ERBalpha inhibits the PTGS2 expression in bovine uterus endometrium stromal and epithelial cells exposed to ovarian steroids. J Reprod Dev 60:362-370. https://doi.org/10.1262/jrd.2014-040 Johnson MH, Lim A, Fernando D, Day ML (2002) Circadian clockwork genes are expressed in the reproductive tract and conceptus of the early pregnant mouse. Reprod Biomed Online 4:140-145. https://doi.org/10.1016/s1472-6483(10)61931-1 Karman BN, Tischkau SA (2006) Circadian clock gene expression in the ovary: Effects of luteinizing hormone. Biol Reprod 75:624-632. https://doi.org/10.1095/biolreprod.106.050732 Kennaway DJ, Varcoe TJ, Voultsios A, Salkeld MD, Rattanatray L, Boden MJ (2015) Acute inhibition of casein kinase 1delta/epsilon rapidly delays peripheral clock gene rhythms. Mol Cell Biochem 398:195-206. https://doi.org/10.1007/s11010-014-2219-8 Kianianmomeni A, Ong CS, Rätsch G, Hallmann A (2014) Genome-wide analysis of alternative splicing in Volvox carteri. BMC Genomics 15:1117. https://doi.org/10.1186/1471-2164-15-1117 Knoedler JR, Avila-Mendoza J, Subramani A, Denver RJ (2020) The Paralogous Kruppel-like Factors 9 and 13 Regulate the Mammalian Cellular Circadian Clock Output Gene Dbp. J Biol Rhythms 35:257-274. https://doi.org/10.1177/0748730420913205 Koike N, Yoo SH, Huang HC, Kumar V, Lee C, Kim TK, Takahashi JS (2012) Transcriptional architecture and chromatin landscape of the core circadian clock in mammals. Science 338:349-354. https://doi.org/10.1126/science.1226339 Kurabayashi N, Hirota T, Sakai M, Sanada K, Fukada Y (2010) DYRK1A and glycogen synthase kinase 3beta, a dual-kinase mechanism directing proteasomal degradation of CRY2 for circadian timekeeping. Mol Cell Biol 30:1757-1768. https://doi.org/10.1128/MCB.01047-09 Lee Y, Rio DC (2015) Mechanisms and Regulation of Alternative Pre-mRNA Splicing. Annu Rev Biochem 84:291-323. https://doi.org/10.1146/annurev-biochem-060614-034316 Li R, Cheng S, Wang Z (2015) Circadian clock gene plays a key role on ovarian cycle and spontaneous abortion. Cell Physiol Biochem 37:911-920. https://doi.org/10.1159/000430218 Liu Y, Johnson BP, Shen AL, Wallisser JA, Krentz KJ, Moran SM, Sullivan R, Glover E, Parlow AF, Drinkwater NR, et al (2014) Loss of BMAL1 in ovarian steroidogenic cells results in implantation failure in female mice. Proc Natl Acad Sci U S A 111:14295-14300. https://doi.org/10.1073/pnas.1209249111 López-Rodríguez D, Franssen D, Sevrin E, Gérard A, Balsat C, Blacher S, Noël A, Parent AS (2019) Persistent vs Transient Alteration of Folliculogenesis and Estrous Cycle After Neonatal vs Adult Exposure to Bisphenol A. Endocrinology 160:2558-2572. https://doi.org/10.1210/en.2019-00505 Luo ZY, Dai XL, Ran XQ, Cen YX, Niu X, Li S, Huang SH, Wang JF (2018) Identification and profile of microRNAs in Xiang pig testes in four different ages detected by Solexa sequencing. Theriogenology 117:61-71. https://doi.org/10.1016/j.theriogenology.2017.06.023 Mansano NDS, Paradela RS, Bohlen TM, Zanardi IM, Chaves FM, Silveira MA, Tavares MR, Donato J, Jr., Frazao R (2022) Vasoactive intestinal peptide exerts an excitatory effect on hypothalamic kisspeptin neurons during estrogen negative feedback. Mol Cell Endocrinol 542:111532. https://doi.org/10.1016/j.mce.2021.111532 Mansuy V, Risold PY, Glauser M, Fraichard A, Pralong FP (2009) Expression of the GABAA receptor associated protein Gec1 is circadian and dependent upon the cellular clock machinery in GnRH secreting GnV-3 cells. Mol Cell Endocrinol 307:68-76. https://doi.org/10.1016/j.mce.2009.02.029 McGlincy NJ, Valomon A, Chesham JE, Maywood ES, Hastings MH, Ule J (2012) Regulation of alternative splicing by the circadian clock and food related cues. Genome Biol 13:R54. https://doi.org/10.1186/gb-2012-13-6-r54 Mendoza J, Lopez-Lopez C, Revel FG, Jeanneau K, Delerue F, Prinssen E, Challet E, Moreau JL, Grundschober C (2011) Dimorphic effects of leptin on the circadian and hypocretinergic systems of mice. J Neuroendocrinol 23:28-38. https://doi.org/10.1111/j.1365-2826.2010.02072.x Miao X, Luo Q (2013) Genome-wide transcriptome analysis between small-tail Han sheep and the Surabaya fur sheep using high-throughput RNA sequencing. Reproduction 145:587-596. https://doi.org/10.1530/REP-12-0507 Miller BH, Olson SL, Turek FW, Levine JE, Horton TH, Takahashi JS (2004) Circadian clock mutation disrupts estrous cyclicity and maintenance of pregnancy. Curr Biol 14:1367-1373. https://doi.org/10.1016/j.cub.2004.07.055 Mohawk JA, Green CB, Takahashi JS (2012) Central and peripheral circadian clocks in mammals. Annu Rev Neurosci 35:445-462. https://doi.org/10.1146/annurev-neuro-060909-153128 Nakamura TJ, Sellix MT, Kudo T, Nakao N, Yoshimura T, Ebihara S, Colwell CS, Block GD (2010) Influence of the estrous cycle on clock gene expression in reproductive tissues: effects of fluctuating ovarian steroid hormone levels. Steroids 75:203-212. https://doi.org/10.1016/j.steroids.2010.01.007 Nakao N, Yasuo S, Nishimura A, Yamamura T, Watanabe T, Anraku T, Okano T, Fukada Y, Sharp PJ, Ebihara S, et al (2007) Circadian clock gene regulation of steroidogenic acute regulatory protein gene expression in preovulatory ovarian follicles. Endocrinology 148:3031-3038. https://doi.org/10.1210/en.2007-0044 Noda M, Iwamoto I, Tabata H, Yamagata T, Ito H, Nagata KI (2019) Role of Per3, a circadian clock gene, in embryonic development of mouse cerebral cortex. Sci Rep 9:5874. https://doi.org/10.1038/s41598-019-42390-9 Ogawa Y, Kawano Y, Yamazaki Y, Onishi Y (2014) Shikonin shortens the circadian period: possible involvement of Top2 inhibition. Biochem Biophys Res Commun 443:339-343. https://doi.org/10.1016/j.bbrc.2013.11.116 Petrillo E, Sanchez SE, Kornblihtt AR, Yanovsky MJ (2011) Alternative splicing adds a new loop to the circadian clock. Commun Integr Biol 4:284-286. https://doi.org/10.4161/cib.4.3.14777 Pilorz V, Steinlechner S (2008) Low reproductive success in Per1 and Per2 mutant mouse females due to accelerated ageing? Reproduction 135:559-568. https://doi.org/10.1530/REP-07-0434 Preußner M, Wilhelmi I, Schultz AS, Finkernagel F, Michel M, Möröy T, Heyd F (2014) Rhythmic U2af26 alternative splicing controls PERIOD1 stability and the circadian clock in mice. Mol Cell 54:651-662. https://doi.org/10.1016/j.molcel.2014.04.015 Ravinder R, Kaipa O, Baddela VS, Singhal Sinha E, Singh P, Nayan V, Velagala CS, Baithalu RK, Onteru SK, Singh D (2016) Saliva ferning, an unorthodox estrus detection method in water buffaloes (Bubalus bubalis). Theriogenology 86:1147-1155. https://doi.org/10.1016/j.theriogenology.2016.04.004 Reppert SM, Weaver DR (2002) Coordination of circadian timing in mammals. Nature 418:935-941. https://doi.org/10.1038/nature00965 Rosenwasser AM, Turek FW (2015) Neurobiology of Circadian Rhythm Regulation. Sleep Med Clin 10:403-412. https://doi.org/10.1016/j.jsmc.2015.08.003 Sanchez SE, Petrillo E, Beckwith EJ, Zhang X, Rugnone ML, Hernando CE, Cuevas JC, Godoy Herz MA, Depetris-Chauvin A, Simpson CG, et al (2010) A methyl transferase links the circadian clock to the regulation of alternative splicing. Nature 468:112-116. https://doi.org/10.1038/nature09470 Schibler U, Ripperger J, Brown SA (2003) Peripheral circadian oscillators in mammals: time and food. J Biol Rhythms 18:250-260. https://doi.org/10.1177/0748730403018003007 Schmutz I, Wendt S, Schnell A, Kramer A, Mansuy IM, Albrecht U (2011) Protein phosphatase 1 (PP1) is a post-translational regulator of the mammalian circadian clock. PLoS One 6:e21325. https://doi.org/10.1371/journal.pone.0021325 Sellix MT (2015) Circadian clock function in the mammalian ovary. J Biol Rhythms 30:7-19. https://doi.org/10.1177/0748730414554222 Sen A, Sellix MT (2016) The Circadian Timing System and Environmental Circadian Disruption: From Follicles to Fertility. Endocrinology 157:3366-3373. https://doi.org/10.1210/en.2016-1450 Shakhmantsir I, Nayak S, Grant GR, Sehgal A (2018) Spliceosome factors target timeless (tim) mRNA to control clock protein accumulation and circadian behavior in Drosophila. Elife 7:e39821. https://doi.org/10.7554/eLife.39821 Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, Zhou Q, Xing Y (2014) rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A 111:E5593-5601. https://doi.org/10.1073/pnas.1419161111 Smarr BL, Gile JJ, de la Iglesia HO (2013) Oestrogen-independent circadian clock gene expression in the anteroventral periventricular nucleus in female rats: possible role as an integrator for circadian and ovarian signals timing the luteinising hormone surge. J Neuroendocrinol 25:1273-1279. https://doi.org/10.1111/jne.12104 Takahashi JS (2017) Transcriptional architecture of the mammalian circadian clock. Nat Rev Genet 18:164-179. https://doi.org/10.1038/nrg.2016.150 Terenina E, Fabre S, Bonnet A, Monniaux D, Robert-Granié C, SanCristobal M, Sarry J, Vignoles F, Gondret F, Monget P, et al (2017) Differentially expressed genes and gene networks involved in pig ovarian follicular atresia. Physiol Genomics 49:67-80. https://doi.org/10.1152/physiolgenomics.00069.2016 Vieyra E, Ramírez DA, Lagunas N, Cárdenas M, Chavira R, Damián-Matsumura P, Trujillo A, Domínguez R, Morales-Ledesma L (2016) Unilaterally blocking the muscarinic receptors in the suprachiasmatic nucleus in proestrus rats prevents pre-ovulatory LH secretion and ovulation. Reproductive biology and endocrinology : RB&E 14:34-34. https://doi.org/10.1186/s12958-016-0168-7 Wang RH, Zhao T, Cui K, Hu G, Chen Q, Chen W, Wang XW, Soto-Gutierrez A, Zhao K, Deng CX (2016) Negative reciprocal regulation between Sirt1 and Per2 modulates the circadian clock and aging. Sci Rep 6:28633. https://doi.org/10.1038/srep28633 Wiegand SJ, Terasawa E, Bridson WE, Goy RW (1980) Effects of discrete lesions of preoptic and suprachiasmatic structures in the female rat. Alterations in the feedback regulation of gonadotropin secretion. Neuroendocrinology 31:147-157. https://doi.org/10.1159/000123066 Wu YG, Bennett J, Talla D, Stocco C (2011) Testosterone, not 5alpha-dihydrotestosterone, stimulates LRH-1 leading to FSH-independent expression of Cyp19 and P450scc in granulosa cells. Mol Endocrinol 25:656-668. https://doi.org/10.1210/me.2010-0367 Yoshikawa T, Sellix M, Pezuk P, Menaker M (2009) Timing of the ovarian circadian clock is regulated by gonadotropins. Endocrinology 150:4338-4347. https://doi.org/10.1210/en.2008-1280 Zhang Y, Meng N, Bao H, Jiang Y, Yang N, Wu K, Wu J, Wang H, Kong S, Zhang Y (2019) Circadian gene PER1 senses progesterone signal during human endometrial decidualization. J Endocrinol. https://doi.org/10.1530/JOE-19-0284 Zheng B, Albrecht U, Kaasik K, Sage M, Lu W, Vaishnav S, Li Q, Sun ZS, Eichele G, Bradley A, et al (2001) Nonredundant roles of the mPer1 and mPer2 genes in the mammalian circadian clock. Cell 105:683-694. https://doi.org/10.1016/s0092-8674(01)00380-4 Additional Declarations No competing interests reported. Supplementary Files supplementaryfile.rar Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {\"props\":{\"pageProps\":{\"initialData\":{\"identity\":\"rs-2593341\",\"acceptedTermsAndConditions\":true,\"allowDirectSubmit\":true,\"archivedVersions\":[],\"articleType\":\"Research Article\",\"associatedPublications\":[],\"authors\":[{\"id\":176975266,\"identity\":\"330a36ff-d6ee-4b33-a7b1-4c3df47112ad\",\"order_by\":0,\"name\":\"Liang-Ting Tang\",\"email\":\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAzElEQVRIiWNgGAWjYBAC+8MMiQ8kDGzk+NkbiNVznOGxgUVBmrFkzwFitZxnfCZR8eFw4oYZCUTqYGxmTpO4YcBsbCD5eOMNhhqbaIJamJnZki1nGLDJmUunFVswHEvLbSCkhY2ZJ/G2hAGPseXsHDMJxobDhLXwMPN/kP5jIJG44eYZIrVIMDMkSUgYGCRuuMFDpBYDZoZkAwmDBGAgA/2SQIxfDPgPAKPyz39gVB7eeONDjQ1hLSjaJRJIUQ7RQqqOUTAKRsEoGBkAADPXO/tLhWSjAAAAAElFTkSuQmCC\",\"orcid\":\"\",\"institution\":\"Guizhou University\",\"correspondingAuthor\":true,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Liang-Ting\",\"middleName\":\"\",\"lastName\":\"Tang\",\"suffix\":\"\"},{\"id\":176975267,\"identity\":\"526ca900-ad8c-497e-8750-916254dfbce3\",\"order_by\":1,\"name\":\"Shihui Huang\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Guizhou University\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Shihui\",\"middleName\":\"\",\"lastName\":\"Huang\",\"suffix\":\"\"},{\"id\":176975269,\"identity\":\"f0767470-2975-4263-ae81-eae0d1405fe0\",\"order_by\":2,\"name\":\"Xi Niu\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Guizhou University\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Xi\",\"middleName\":\"\",\"lastName\":\"Niu\",\"suffix\":\"\"},{\"id\":176975270,\"identity\":\"3a39c71a-02ec-4c87-9b60-524759c66fe9\",\"order_by\":3,\"name\":\"Sheng Li\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Guizhou University\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Sheng\",\"middleName\":\"\",\"lastName\":\"Li\",\"suffix\":\"\"},{\"id\":176975271,\"identity\":\"7f6534cf-1385-4913-98a7-2097ece1d5b7\",\"order_by\":4,\"name\":\"Jiafu Wang\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Guizhou University\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Jiafu\",\"middleName\":\"\",\"lastName\":\"Wang\",\"suffix\":\"\"},{\"id\":176975272,\"identity\":\"3c093f63-b3a3-40eb-a9e2-342e45fc0d1b\",\"order_by\":5,\"name\":\"Xue-Qin Ran\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Guizhou University\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Xue-Qin\",\"middleName\":\"\",\"lastName\":\"Ran\",\"suffix\":\"\"},{\"id\":176975273,\"identity\":\"c5d6eb88-f16a-41f7-8b80-db8ac7a4624b\",\"order_by\":6,\"name\":\"Liang-Ting Tang\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Guizhou Open University\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Liang-Ting\",\"middleName\":\"\",\"lastName\":\"Tang\",\"suffix\":\"\"}],\"badges\":[],\"createdAt\":\"2023-02-16 06:29:26\",\"currentVersionCode\":1,\"declarations\":\"\",\"doi\":\"10.21203/rs.3.rs-2593341/v1\",\"doiUrl\":\"https://doi.org/10.21203/rs.3.rs-2593341/v1\",\"draftVersion\":[],\"editorialEvents\":[],\"editorialNote\":\"\",\"failedWorkflow\":false,\"files\":[{\"id\":33203991,\"identity\":\"26ba6a58-0a93-4b77-a099-c6fa4a092463\",\"added_by\":\"auto\",\"created_at\":\"2023-02-21 02:43:17\",\"extension\":\"png\",\"order_by\":1,\"title\":\"Figure 1\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":10510,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eExpression frequency of circadian rhythm genes in the ovaries of Xiang pigs in estrus and diestrus period \\u003c/strong\\u003e(x-axis: CPM Value; y-axis: Number of genes).\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"1.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-2593341/v1/3610a4860c8ced0a0bfa4509.png\"},{\"id\":33203992,\"identity\":\"f2de874e-7ee2-4a65-a1fe-215a8a588461\",\"added_by\":\"auto\",\"created_at\":\"2023-02-21 02:43:17\",\"extension\":\"png\",\"order_by\":2,\"title\":\"Figure 2\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":178766,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eThe verification of differential expression by RT-qPCR. (\\u003c/strong\\u003eThe gene \\u003cem\\u003eppp1cb, csnk1d, per1, klf9, nampt, bhlhe40, id4, \\u003c/em\\u003eand\\u003cem\\u003e ass1\\u003c/em\\u003ewere significantly high expressed in estrus samples. These data were in accordance with the RNA-seq results and indicated that RNA-Seq results were reliable.\\u003cstrong\\u003e)\\u003c/strong\\u003e\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"2.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-2593341/v1/b58d424ccb62f821bec336c3.png\"},{\"id\":33205204,\"identity\":\"15aa7a7f-0972-40d5-b942-a50767fd8467\",\"added_by\":\"auto\",\"created_at\":\"2023-02-21 02:51:17\",\"extension\":\"png\",\"order_by\":3,\"title\":\"Figure 3\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":403136,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eVerification of rhythmic expression of clock genes RT-qPCR\\u003c/strong\\u003e.RNA was extracted from ovarian samples for qPCR validation of clock-related genes, and data were shown as mean ± SE (n = 3). GraphPad Prism was used to plot line chart, and independent samples t-test and One-way ANOVA were used for statistical analysis of the corresponding test data. P \\u0026lt;0.05 was considered statistically significant.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"3.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-2593341/v1/1677d7e9e99fcd674d521576.png\"},{\"id\":33205203,\"identity\":\"9f2f5746-840c-454b-a609-921ddfefca68\",\"added_by\":\"auto\",\"created_at\":\"2023-02-21 02:51:17\",\"extension\":\"png\",\"order_by\":4,\"title\":\"Figure 4\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":120179,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eThe Effects of estrus cycle on ovarian circadian rhythm genes of Xiang pigs \\u003c/strong\\u003e(This figure showcases the differentially expressed gene products in the estrus and diestrus ovaries are related to the post-translational modification of the core clock protein, which may lead to changes in the expression of core clock genes and clock-controlled genes\\u003cstrong\\u003e).\\u003c/strong\\u003e\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"4.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-2593341/v1/56fea3fe40c423bfd0de6deb.png\"},{\"id\":33355329,\"identity\":\"68775d0a-84cc-4e91-9c2d-439026952bc5\",\"added_by\":\"auto\",\"created_at\":\"2023-02-23 15:59:45\",\"extension\":\"pdf\",\"order_by\":0,\"title\":\"\",\"display\":\"\",\"copyAsset\":false,\"role\":\"manuscript-pdf\",\"size\":945164,\"visible\":true,\"origin\":\"\",\"legend\":\"\",\"description\":\"\",\"filename\":\"manuscript.pdf\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-2593341/v1/012b6593-f888-41ff-bd75-5f123a3c1513.pdf\"},{\"id\":33203990,\"identity\":\"e176ce12-3f57-4ccc-9bcb-1436655ae1cd\",\"added_by\":\"auto\",\"created_at\":\"2023-02-21 02:43:17\",\"extension\":\"rar\",\"order_by\":1,\"title\":\"\",\"display\":\"\",\"copyAsset\":false,\"role\":\"supplement\",\"size\":90308,\"visible\":true,\"origin\":\"\",\"legend\":\"\",\"description\":\"\",\"filename\":\"supplementaryfile.rar\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-2593341/v1/9f03c0a77dabb41e2c3af3ef.rar\"}],\"financialInterests\":\"No competing interests reported.\",\"formattedTitle\":\"The effects of estrus cycle on the expression of ovarian biological clock-related genes of Xiang pigs\",\"fulltext\":[{\"header\":\"Introduction\",\"content\":\"\\u003cp\\u003eCircadian rhythm is the regular change of many functional activities of the biological organism to adapt to the alternation of day and night in the environment. It is endogenous and persists even without external synchronising factors, for example, in constant darkness. Circadian rhythm is controlled by central biological clocks related to hypothalamic suprachiasmatic nucleus (SCN) and peripheral biological clocks distributed in various tissues (Schibler, et al. \\u003cspan citationid=\\\"CR52\\\" class=\\\"CitationRef\\\"\\u003e2003\\u003c/span\\u003e). SCN cells modulate the activity of gonadotropin-releasing hormone (GnRH) neurons by direct innervation or other hypothalamic nucleus, which are necessary for the occurrence of ovulation (Vieyra, et al. \\u003cspan citationid=\\\"CR61\\\" class=\\\"CitationRef\\\"\\u003e2016\\u003c/span\\u003e; Mansano, et al. \\u003cspan citationid=\\\"CR34\\\" class=\\\"CitationRef\\\"\\u003e2022\\u003c/span\\u003e). It is known that the molecular model of mammals circadian clock including human is a multi-feedback loop for transcription and translation including transcriptional regulation and post-transcriptional regulation (Takahashi \\u003cspan citationid=\\\"CR59\\\" class=\\\"CitationRef\\\"\\u003e2017\\u003c/span\\u003e). It consists of positive regulatory components (BMAL1 and CLOCK) and negative regulatory components (PER 1/2 and CRY 1/2) (Reppert and Weaver \\u003cspan citationid=\\\"CR49\\\" class=\\\"CitationRef\\\"\\u003e2002\\u003c/span\\u003e). The positive and the negative regulatory elements also participate in the circadian regulation of clock controlled genes (CCG). The clock-controlled genes act as specific regulators of rhythmic physiological functions in cells and tissues (Sen and Sellix \\u003cspan citationid=\\\"CR55\\\" class=\\\"CitationRef\\\"\\u003e2016\\u003c/span\\u003e). For example, the function of ovary is regulated by the biological clock, including the estrus cycle, luteinizing hormone (LH) levels, ovulation and embryo implantation (Gray, et al. \\u003cspan citationid=\\\"CR15\\\" class=\\\"CitationRef\\\"\\u003e1978\\u003c/span\\u003e; Christian, et al. \\u003cspan citationid=\\\"CR7\\\" class=\\\"CitationRef\\\"\\u003e2005\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003eOvary is an important reproductive organ of female mammals. In each estrus cycle, the ovaries undergo a series of morphological, hormonal and biochemical changes, which directly affect and/or determine the litter size of female animals (Miao and Luo \\u003cspan citationid=\\\"CR38\\\" class=\\\"CitationRef\\\"\\u003e2013\\u003c/span\\u003e). The realization of ovarian function is a very complex biological process, involving the transcriptional regulation of a large number of genes (Terenina, et al. \\u003cspan citationid=\\\"CR60\\\" class=\\\"CitationRef\\\"\\u003e2017\\u003c/span\\u003e). Alternative splicing (AS) is also an important way of gene expression regulation, which occurs widely in various tissues and cells of almost all organisms (Chacko and Ranganathan \\u003cspan citationid=\\\"CR4\\\" class=\\\"CitationRef\\\"\\u003e2009\\u003c/span\\u003e; Kianianmomeni, et al. \\u003cspan citationid=\\\"CR25\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). Many studies demonstrate that the biological clock existed in the ovary (Fahrenkrug, et al. \\u003cspan citationid=\\\"CR12\\\" class=\\\"CitationRef\\\"\\u003e2006\\u003c/span\\u003e; Karman and Tischkau \\u003cspan citationid=\\\"CR23\\\" class=\\\"CitationRef\\\"\\u003e2006\\u003c/span\\u003e; He, et al. \\u003cspan citationid=\\\"CR17\\\" class=\\\"CitationRef\\\"\\u003e2007b\\u003c/span\\u003e; Yoshikawa, et al. \\u003cspan citationid=\\\"CR65\\\" class=\\\"CitationRef\\\"\\u003e2009\\u003c/span\\u003e). The expression of clock genes have been described in the ovary of rat (He, et al. \\u003cspan citationid=\\\"CR16\\\" class=\\\"CitationRef\\\"\\u003e2007a\\u003c/span\\u003e), mice (Johnson, et al. \\u003cspan citationid=\\\"CR22\\\" class=\\\"CitationRef\\\"\\u003e2002\\u003c/span\\u003e) and ruminant (Cushman, et al. \\u003cspan citationid=\\\"CR10\\\" class=\\\"CitationRef\\\"\\u003e2007\\u003c/span\\u003e). It indicates that clock genes are involved in the process of follicle formation, ovulation, oocyte maturation and steroid hormone synthesis (Sellix \\u003cspan citationid=\\\"CR54\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). The lack of circadian clock might be related with varieties of reproductive disorders. For instance, mice with \\u003cem\\u003eper1\\u003c/em\\u003e and \\u003cem\\u003eper2\\u003c/em\\u003e knockout present a significant decrease in reproductive rate due to the irregular and non-cyclical estrus cycle (Zheng, et al. \\u003cspan citationid=\\\"CR67\\\" class=\\\"CitationRef\\\"\\u003e2001\\u003c/span\\u003e; Pilorz and Steinlechner \\u003cspan citationid=\\\"CR46\\\" class=\\\"CitationRef\\\"\\u003e2008\\u003c/span\\u003e). Although mice with \\u003cem\\u003ebmal1\\u003c/em\\u003e knockout are capable of ovulation, they exhibit delayed puberty, irregular estrus cycles, and small ovaries and uterus (Alvarez, et al. \\u003cspan citationid=\\\"CR1\\\" class=\\\"CitationRef\\\"\\u003e2008\\u003c/span\\u003e; Boden, et al. \\u003cspan citationid=\\\"CR2\\\" class=\\\"CitationRef\\\"\\u003e2010\\u003c/span\\u003e; Liu, et al. \\u003cspan citationid=\\\"CR31\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). Loss of \\u003cem\\u003ebmal1\\u003c/em\\u003e also negatively affects progesterone level (Liu, et al. \\u003cspan citationid=\\\"CR31\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e), leading to embryo implantation failure (Boden, et al. \\u003cspan citationid=\\\"CR2\\\" class=\\\"CitationRef\\\"\\u003e2010\\u003c/span\\u003e). Mice mutant with 51 amino acid deletions in the transcription activation domain of clock protein show an irregular estrus cycle, and the normal LH surges disappear on the day of estrus (Miller, et al. \\u003cspan citationid=\\\"CR39\\\" class=\\\"CitationRef\\\"\\u003e2004\\u003c/span\\u003e). In addition, genes of ovarian granulosa cells that are important for follicular development have been shown to be controlled by clock genes, such as LH receptor (\\u003cem\\u003elhcgr\\u003c/em\\u003e), prostaglandin synthase (\\u003cem\\u003eptgs2\\u003c/em\\u003e), steroid-producing enzymes (steroid-producing acute regulatory protein, cyp11a1, aromatase, etc.), and gap junction protein Connexin-43 (Chen, et al. \\u003cspan citationid=\\\"CR5\\\" class=\\\"CitationRef\\\"\\u003e2012\\u003c/span\\u003e; Chu, et al. \\u003cspan citationid=\\\"CR8\\\" class=\\\"CitationRef\\\"\\u003e2012\\u003c/span\\u003e; Chen, et al. \\u003cspan citationid=\\\"CR6\\\" class=\\\"CitationRef\\\"\\u003e2013\\u003c/span\\u003e). Estrus is a sign of reproductive maturity. Moreover, it has been confirmed that the central biological clock located in the SCN was the core of maintaining normal reproductive function and the estrus cycle (Brown-Grant and Raisman \\u003cspan citationid=\\\"CR3\\\" class=\\\"CitationRef\\\"\\u003e1977\\u003c/span\\u003e; Wiegand, et al. \\u003cspan citationid=\\\"CR63\\\" class=\\\"CitationRef\\\"\\u003e1980\\u003c/span\\u003e). However, the function of ovarian biological clock during estrus period is still needed to be clear.\\u003c/p\\u003e \\u003cp\\u003eThe Xiang pig is a small local breed in China with some unique biological characteristics, such as small size, precocious puberty, low birth rate and insignificant estrous behavior (Luo, et al. \\u003cspan citationid=\\\"CR33\\\" class=\\\"CitationRef\\\"\\u003e2018\\u003c/span\\u003e). In order to understand more in-depth views of the influence of estrus on biological clock and further explore the physiological functions of clock genes in porcine ovary, we analyzed the expression of biological clock-related genes and alternative splicing of transcripts in the ovaries of Xiang pigs during estrus and diestrus period based on the whole transcriptome sequencing analysis.\\u003c/p\\u003e\"},{\"header\":\"Materials And Methods\",\"content\":\"\\u003cdiv id=\\\"Sec3\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eMaterials\\u003c/h2\\u003e \\u003cp\\u003eAbout 40 sows of Xiang pig were selected for observation of estrus performance. They were reared in a standard manner that the pigs drank freely and fed timely under the normal condition. The average light intensity was about 200 lux, the ambient temperature was maintained at 15\\u0026deg;C -18\\u0026deg;C, and the humidity was in the range of 50% \\u0026plusmn; 10%. The estrus detection was carried out 2 to 3 times a day. Redness and swelling of the vulva and a positive standing reflex test of the boar were defined as estrus manifestations. Ovary samples of 14 Xiang pig sows were taken, 7 of them were in the diestrus period, and the others were subjected to routine castration surgery on the third day of estrus (the sows showed strong redness and swelling of the vulva). The sampling day was at 8\\u0026ndash;9 am at March 15 to 30, 2022, with the same light and temperature conditions as previously mentioned, and fasting before ovariectomy surgery. All samples were immediately immersed in liquid nitrogen and stored at -80\\u0026deg;C for RNA extraction. The animal procedures were approved by the guidelines of Guizhou University Subcommittee of Experimental Animal Ethics with No. of EAE-GZU-2020-P002.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec4\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eTotal RNA extraction\\u003c/h2\\u003e \\u003cp\\u003eTotal RNA from 14 ovarian samples were extracted using TRIzol reagent kit (Invitrogen, USA) according to the instructions. The quantity and completeness of total RNA were evaluated by the Bio-analyzer system (Agilent 2100, Agilent Technologies, USA). The total RNA was stored at -80\\u0026deg;C for subsequent analysis. The same samples were used for sequencing and PCR analysis.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec5\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eLibrary construction and sequencing\\u003c/h2\\u003e \\u003cp\\u003eTotal of 14 cDNA libraries were constructed using RNA as templates from ovarian tissues. The library construction was carried out according to the standard procedures of BGI (Shenzhen Huada Gene Co., Ltd). After DNase I enzyme treatment, magnetic beads with oligo (dT) were used to isolate mRNA. Under appropriate temperature conditions, the purified mRNA was fragmented by adding fragmentation buffer, and then cDNA was synthesized using random primers and mRNA fragment templates. After 3'end adenylation, the DNA fragments were ligated to the linker adaptors. cDNA fragments of 100\\u0026thinsp;~\\u0026thinsp;200 bp in length were selected for PCR amplification to generate a cDNA library. The libraries were sequenced on the Illumina HiSeqTM 2500 sequencing platform, and paired-end sequences with 150 bp in length were generated.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec6\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eSequencing data analysis\\u003c/h2\\u003e \\u003cp\\u003eThe original sequences generated by the Illumina platform were saved in fastq format. The clean reads were obtained by removing low-quality sequences, reads (more than 5% unknown nucleotides), and sequencing adapters. The clean reads should be evaluated through FastaQC software (Q20\\u0026thinsp;\\u0026lt;\\u0026thinsp;20%). The reference genome sequence (ssc11.1 version) and annotation files of pig genomes were downloaded from the Ensembl database (\\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttp://www.ensembl.org/index.html\\u003c/span\\u003e\\u003cspan address=\\\"http://www.ensembl.org/index.html\\\" targettype=\\\"URL\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e). Used HISAT2 software (v2.1.0) to map the clean reads to the pig reference genome. HISAT2 is an efficient alignment tool for RNA-seq experiments. It uses an indexing framework based on the BWT and Ferragina-Manzini (FM) index algorithms for alignment and data normalization.The Subread featureCounts software (v2.0.0) was used to count the number of reads mapped to the gene. Subread featureCounts applies the read summarization tool to analyze the expression level from RNA-seq, which applies efficient chromosome hash algorithm and feature block technology for statistics and data standardization.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec7\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eAnalysis of differential expression\\u003c/h2\\u003e \\u003cp\\u003eThe expressed level of each gene were calculated by CPM values (counts per million mapped reads) in all samples, and added 0.001 to CPM value as the gene expression level. Then we used bioconductor software package Limma and DESeq2 to analyze differentially expressed genes (DEGs) between two groups with biological replicates. The FDR\\u0026thinsp;\\u0026le;\\u0026thinsp;0.1 and |log2FC|\\u0026ge;0.5 were used as the thresholds for determining the significance of gene expression difference.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec8\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eIdentification of differential alternative splicing events\\u003c/h2\\u003e \\u003cp\\u003eThe rMATS software (v4.0.2) was used to analyze differentially alternative splicing events between two groups. The rMATS recognizes five alternative splicing events, which are skipped exon (SE), alternative 5\\u0026rsquo; splice site (A5SS), alternative 3\\u0026rsquo; splice site (A3SS), mutually exclusive exons (MXE) and retained intron (RI) (Shen, et al. \\u003cspan citationid=\\\"CR57\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). Also, FDR\\u0026thinsp;\\u0026le;\\u0026thinsp;0.05 was used as a threshold to detect the differential significance of alternative splicing.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec9\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eValidation\\u003c/h2\\u003e \\u003cp\\u003eThe same RNA samples from the ovarian for RNA-seq experiment were used to validate the differential expression by qRT-PCR and to validate the AS events by RT-PCR method. The primers were designed using Primer 6.0, which was listed in table S4. GAPDH and β-actin genes were taken as internal controls. The PCR efficiency of the primers was controlled within 100\\u0026thinsp;\\u0026plusmn;\\u0026thinsp;10%. All qPCR reactions were performed on the CFX98TM real-time system (Bio-rad co.) according to the manufacturer's instructions and recommended cycling conditions. Three replicate reaction tests were performed using a volume of 10 \\u0026micro;L (containing 5 \\u0026micro;L of 2\\u0026times; Super Real Pre Mix Plus), 0.3 \\u0026micro;L of each of the forward and reverse primers (10 pM / \\u0026micro;L) and 1 \\u0026micro;L cDNA. The qPCR conditions were as following: initial denaturation at 95\\u0026deg;C for 3 minutes, then denaturation at 94\\u0026deg;C for 10 seconds in 40 cycles, and annealing and extension at 60\\u0026deg;C for 30 seconds. The relative expression level of target gene utilized the method of 2(-ΔΔCt) as reported by Livak et al. (Livak and Schmittgen, 2001). The different level of gene expression between two groups was tested by software SPSS (v21.0) taking the P\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.05 as threshold of significant difference. The results were presented as mean\\u0026thinsp;\\u0026plusmn;\\u0026thinsp;standard deviation. The 8 AS events were randomly selected from 8 genes to validate by RT-PCR method. The RT-PCR annealing temperature was 55\\u0026deg;C for 30 seconds, and the extension temperature was 72\\u0026deg;C for 60 seconds by 34 cycles.\\u003c/p\\u003e \\u003cp\\u003eTo verify the rhythmic expression of circadian rhythm-related genes in Xiang pig ovary, six genes, \\u003cem\\u003eper1\\u003c/em\\u003e, \\u003cem\\u003ecry1\\u003c/em\\u003e, \\u003cem\\u003eclock\\u003c/em\\u003e, \\u003cem\\u003earntl\\u003c/em\\u003e, \\u003cem\\u003eppp1cb\\u003c/em\\u003e and \\u003cem\\u003entrk1\\u003c/em\\u003e, were selected to detected quantitatively by RT-qPCR method at six time points (08:00, 12:00, 16:00, 20:00, 24:00, 04:00). Three ovary samples of Xiang pigs were collected at each time point, and RNA was extracted from the ovarian samples for qPCR verification of genes. The amplification conditions were the same as the previous RT-qPCR method.\\u003c/p\\u003e \\u003c/div\\u003e\"},{\"header\":\"Results\",\"content\":\"\\u003cdiv id=\\\"Sec11\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eSummary of sequencing data\\u003c/h2\\u003e \\u003cp\\u003eThe 14 cDNA libraries from ovaries in estrus and diestrus stages were constructed and sequenced respectively, and the sequencing data were filtered and normalized, which were then mapped to the reference genome with HISAT2. It generated 75002160\\u0026ndash;77100634 read pairs, and 93.57%-96.26% of reads were aligned to the reference genome (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). The only matching reads were performed the following statistical analysis.\\u003c/p\\u003e \\u003cp\\u003e \\u003cdiv class=\\\"gridtable\\\"\\u003e\\u003ctable float=\\\"Yes\\\" id=\\\"Tab1\\\" border=\\\"1\\\"\\u003e \\u003ccaption language=\\\"En\\\"\\u003e \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 1\\u003c/div\\u003e \\u003cdiv class=\\\"CaptionContent\\\"\\u003e \\u003cp\\u003eRNA-Seq data description\\u003c/p\\u003e \\u003c/div\\u003e \\u003c/caption\\u003e \\u003ccolgroup cols=\\\"8\\\"\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c1\\\" colnum=\\\"1\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c2\\\" colnum=\\\"2\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"char\\\" char=\\\".\\\" class=\\\"colspec\\\" colname=\\\"c3\\\" colnum=\\\"3\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"char\\\" char=\\\".\\\" class=\\\"colspec\\\" colname=\\\"c4\\\" colnum=\\\"4\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"char\\\" char=\\\".\\\" class=\\\"colspec\\\" colname=\\\"c5\\\" colnum=\\\"5\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"char\\\" char=\\\".\\\" class=\\\"colspec\\\" colname=\\\"c6\\\" colnum=\\\"6\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"char\\\" char=\\\".\\\" class=\\\"colspec\\\" colname=\\\"c7\\\" colnum=\\\"7\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"char\\\" char=\\\".\\\" class=\\\"colspec\\\" colname=\\\"c8\\\" colnum=\\\"8\\\"\\u003e\\u003c/div\\u003e \\u003cthead\\u003e \\u003ctr\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c1\\\"\\u003e\\u0026nbsp;\\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eSamples\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003eInput\\u0026nbsp;Read \\u0026nbsp;Pairs\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003eBoth\\u0026nbsp; Surviving\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003eForward\\u0026nbsp;Only\\u0026nbsp; Surviving\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003eReverse\\u0026nbsp;Only\\u0026nbsp; Surviving\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003eDropped\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003ePercentage\\u0026nbsp;of\\u0026nbsp; mapped\\u0026nbsp;reads (%)\\u003c/p\\u003e \\u003c/th\\u003e \\u003c/tr\\u003e \\u003c/thead\\u003e \\u003ctbody\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\" morerows=\\\"6\\\" rowspan=\\\"7\\\"\\u003e \\u003cp\\u003eDiestrus\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e77089626\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74382355 (96.49%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2471254 (3.21%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e164408 \\u0026nbsp;(0.21%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e69380 (0.09%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e96.01%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e77040960\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74260107 (96.39%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2550481 (3.31%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e160145\\u0026nbsp; (0.21%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e70227 \\u0026nbsp;(0.09)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e95.88%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e76170080\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e73520243 (96.52%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2408511 (3.16%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e167409 (0.22%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e73917 (0.10%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e95.93%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e77100634\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74147772 (96.17%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2713372 (3.52%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e165857 \\u0026nbsp;(0.22%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e73633 (0.10%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e96.26%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e76860177\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74164234 (96.49%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2449821\\u0026nbsp; (3.19%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e177911 (0.23%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e68211 (0.09)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e96.06%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD6\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e76903234\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74486606 (96.86%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2179432\\u0026nbsp; (2.83%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e173980 (0.23%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e63216 (0.08%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e96.11%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eD7\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e77054008\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74411266 (96.57%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e2408339\\u0026nbsp; (3.13%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e163511 \\u0026nbsp;(0.21%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e70892\\u0026nbsp; (0.09%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e95.95%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\" morerows=\\\"6\\\" rowspan=\\\"7\\\"\\u003e \\u003cp\\u003eEstrus\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e75005968\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74960964 (99.94%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e45004 (0.06%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e95.98%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e74990854\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74960858 (99.96%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e29996 (0.04%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e95.90%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e75206715\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e75176632 (99.96%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e30083 (0.04%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e95.93%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e75187643\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e75135012 (99.93%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e52631 (0.07%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e31 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e0 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e96.24%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e75184086\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e75154012 (99.96%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e30074 (0.04%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e31 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e1 (0.00%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e96.15%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE6\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e76109159\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e74099877 (97.36%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e1682012 (2.21%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e243549 (0.32%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e83720 (0.11%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e94.13%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eE7\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e75002160\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e73149607 (97.53%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e1545044 (2.06%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e225006 (0.30%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e82502 (0.11%)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"char\\\" char=\\\".\\\" colname=\\\"c8\\\"\\u003e \\u003cp\\u003e93.57%\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003c/tbody\\u003e \\u003c/colgroup\\u003e \\u003c/table\\u003e\\u003c/div\\u003e \\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec12\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eExpression of biological clock-related genes in Xiang pig ovary\\u003c/h2\\u003e \\u003cp\\u003eAfter mapping and annotating based on the pig reference genome assembly Sscrofa11.1, the RNA-seq data of circadian rhythms and clock genes were picked up from the genome-scale data sets (Table S1). Taking threshold value of 0.1 CPM, 90 rhythms genes were determined, in which 81 were expressed in the estrus ovaries and 83 were detected from the diestrus ovaries. The expression levels of those genes were ranged from 867.2 CPM in the highest abundant transcripts to the cut off value of 0.1 CPM (Table S1). Compared with the diestrus ovary, there were more genes (55.6% in estrus, 42.2% in diestrus) at the medium expression level (10\\u0026ndash;100 CPM), and fewer genes at the high (\\u0026gt;\\u0026thinsp;100 CPM) or low (0.1\\u0026ndash;10 CPM) expression level in the estrus ovary (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). The top expressed genes were \\u003cem\\u003enono\\u003c/em\\u003e and \\u003cem\\u003esfpq\\u003c/em\\u003e with 455.6 and 446.1 CPM in the estrus ovary, and \\u003cem\\u003esfpq\\u003c/em\\u003e and \\u003cem\\u003etop2a\\u003c/em\\u003e with 532.6 and 512.1 CPM in the diestrus ovary, respectively.\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec13\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eAnalysis of differential expression genes\\u003c/h2\\u003e \\u003cp\\u003eWe used R package Limma and DESeq2 to analyze differentially expressed genes (DEGs) between estrus and diestrus ovaries. Also, the P-value\\u0026thinsp;\\u0026le;\\u0026thinsp;0.01 and absolute value of log\\u003csub\\u003e2\\u003c/sub\\u003eFC\\u0026thinsp;\\u0026ge;\\u0026thinsp;0.5 were used as the threshold for determining the significant difference of gene expression. The results showed that 33 DEGs in the ovaries between estrus and diestrus cycles, comprising 13 down-regulated genes (\\u003cem\\u003etop2a, nrip1, prkdc, setx, ncoa2, pparg, pspc1, per3, sirt1, timeless, ezh2, ppara\\u003c/em\\u003e, and \\u003cem\\u003eppp1ca\\u003c/em\\u003e) and 8 up-regulated genes \\u003cem\\u003e(atf4, csnk1d, klf10, id2, bhlhe40, id4, ass1\\u003c/em\\u003e, and \\u003cem\\u003entrk1\\u003c/em\\u003e) in at least 2-fold at estrus stage. The fold changes for DEGs were ranged from 0.17 to 26. Furthermore, we found that 5 genes (\\u003cem\\u003eppp1cb, id3, per1, klf9\\u003c/em\\u003e, and \\u003cem\\u003enampt\\u003c/em\\u003e) were at least 1.5-2 fold upregulated and 7 genes (\\u003cem\\u003eep300, ube3a, dyrk1a, sin3a, clock, creb1\\u003c/em\\u003e, and \\u003cem\\u003emagel2\\u003c/em\\u003e) were at least 1.5-2 fold down regulated in estrus ovary (Table S2). A lot of the genes were expressed in ovaries at medium or high levels (48.9% and 27.4% ), including the core clock components, \\u003cem\\u003earntl, clock, per1/2/3\\u003c/em\\u003e, and \\u003cem\\u003ecry1/2\\u003c/em\\u003e.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec14\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eDetection of differential splicing genes\\u003c/h2\\u003e \\u003cp\\u003eThe rMATS was used to analyze differentially alternative splicing genes (DSGs) between estrus and diestrus ovaries. Overall, 44 differential splicing events were identified from the transcripts of 34 circadian rhythms genes (P\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.05) including 3 A5SS, 12 RI, 14 SE and 15 MXE events (Table S3). There were multiple genes with complex alternative splicing patterns. Of those, \\u003cem\\u003enono\\u003c/em\\u003e presented the most abundant differential alternative splicing events, which included four differential splicing types: A5SS, RI, SE, and MXE. Also, ppp1cb harbored three differential splicing types (RI, SE and MXE). Then, \\u003cem\\u003eENSSSCG00000000107, per1, clock, ube3a\\u003c/em\\u003e and \\u003cem\\u003efbxw11\\u003c/em\\u003e each contained two differential alternative splicing types. It was found that 20 genes comprised the core clock components, \\u003cem\\u003earntl\\u003c/em\\u003e and \\u003cem\\u003ecry1\\u003c/em\\u003e were differentially regulated only at AS level. At least 14 genes, such as \\u003cem\\u003eper1\\u003c/em\\u003e and \\u003cem\\u003eclock\\u003c/em\\u003e, were differentially regulated in both expression and AS level in ovaries between estrus and diestrus.\\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec15\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eValidation of DEGs by RT-qPCR method\\u003c/h2\\u003e \\u003cp\\u003eTo validate the results from RNA-seq, we used RT-qPCR method to further examined the differentially expressed genes using the same RNA samples from the ovaries (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e). It confirmed that \\u003cem\\u003eppp1cb, csnk1d, per1, klf9, nampt, bhlhe40, id4\\u003c/em\\u003e, and \\u003cem\\u003eass1\\u003c/em\\u003e were significantly high expressed in estrus samples, which were in accordance with the RNA-seq analysis (Table S2). The expression of \\u003cem\\u003eatf4\\u003c/em\\u003e and \\u003cem\\u003entrk1\\u003c/em\\u003e in estrus and diestrus ovaries is not significantly different. Moreover, 8 AS events, selected randomly from 8 genes, were detected by RT-PCR and consistent with the RNA-Seq results (Figure S1).The primer sequences were shown in Table S4\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003c/div\\u003e \\u003cdiv id=\\\"Sec16\\\" class=\\\"Section2\\\"\\u003e \\u003ch2\\u003eValidation of rhythmic expression of clock-related genes\\u003c/h2\\u003e \\u003cp\\u003eTo verify whether the expression of clock-related genes in Xiang pig ovaries was rhythmic, we further examined the expression patterns of \\u003cem\\u003eper1, cry1, clock, arntl, ppp1cb\\u003c/em\\u003e and \\u003cem\\u003entrk1\\u003c/em\\u003e at six time points (08:00, 12:00, 16:00, 20:00, 24:00, 04:00) using RT-qPCR. The samples at each time point were obtained from 3 Xiang pig ovaries collected within one hour. Among them, the expression of \\u003cem\\u003eper1\\u003c/em\\u003e and \\u003cem\\u003eppp1cb\\u003c/em\\u003e peaked at 24:00 midnight. However, the expression patterns of \\u003cem\\u003eclock\\u003c/em\\u003e and \\u003cem\\u003earntl\\u003c/em\\u003e were opposite, in which the peaks presented at 8:00 am and it fell to the lowest at 24:00. The maximum values was appeared at 20:00 and 24:00 pm for gene \\u003cem\\u003ecry1\\u003c/em\\u003e, and it was at 16:00 pm for gene \\u003cem\\u003entrk1\\u003c/em\\u003e. It illustrated that the core clock genes \\u003cem\\u003eper1\\u003c/em\\u003e, \\u003cem\\u003ecry1\\u003c/em\\u003e, \\u003cem\\u003eclock\\u003c/em\\u003e and \\u003cem\\u003earntl\\u003c/em\\u003e and the clock-related genes \\u003cem\\u003eppp1cb\\u003c/em\\u003e and \\u003cem\\u003entrk1\\u003c/em\\u003e were rhythmic expressed in Xiang pig ovaries (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig3\\\" class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003c/div\\u003e\"},{\"header\":\"Discussion\",\"content\":\"\\u003cp\\u003eThe mammalian biological clock system has a multi-level structure, including the main clock located in the suprachiasmatic nucleus and the sub-clocks of peripheral organs and tissues (Schibler, et al. \\u003cspan citationid=\\\"CR52\\\" class=\\\"CitationRef\\\"\\u003e2003\\u003c/span\\u003e). The ovary is a pivotal reproductive organ in female animal, and the expression of the biological clock gene has been found in the ovaries of many organisms (Sellix \\u003cspan citationid=\\\"CR54\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). In this study, we analyzed the expression profile and alternative splicing of the biological clock gene in Xiang pig ovaries in estrous cycle from the RNA-seq data. The study found out a total of 90 rhythms genes that were expressed in ovaries. The expression abundance of these genes was ranging from 0.1 to 532.6 CPM. A lot of the genes were expressed in ovaries at medium or high levels (48.9% and 27.4%), including the core clock components, \\u003cem\\u003earntl, clock, per1/2/3\\u003c/em\\u003e, and \\u003cem\\u003ecry1/2\\u003c/em\\u003e. Furthermore, 33 genes were detected to undergo differentially expression and 34 genes were detected to undergo differential alternative splicing between estrous and diestrous ovaries. The DEGs and DSGs related with rhythms might have a connection with the regulation of estrus process in Xiang pig.\\u003c/p\\u003e \\u003cp\\u003eTranscription and translation of core clock components genes (\\u003cem\\u003eclock\\u003c/em\\u003e or \\u003cem\\u003enpas2, arntl/bmal1\\u003c/em\\u003e, or \\u003cem\\u003earntl2/bmal2, per1/2/3\\u003c/em\\u003e, and \\u003cem\\u003ecry1/2\\u003c/em\\u003e) play a critical role in rhythm generation process. CLOCK and BMAL1 heterodimerize to activate transcription of circadian target genes including \\u003cem\\u003eper1/2/3\\u003c/em\\u003e and \\u003cem\\u003ecry1/2\\u003c/em\\u003e. PER and CRY interact and conversely inhibit transcription of bmal1 and clock genes. These genes and their protein products are organized into interlocking positive and negative transcriptional and translational feedback loops, which regulates circadian rhythm generation in the brain suprachiasmatic nucleus (SCN) and peripheral organs (Rosenwasser and Turek \\u003cspan citationid=\\\"CR50\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). In this study, we found that three core clock genes were differentially expressed between estrus and diestrus ovaries. \\u003cem\\u003eClock\\u003c/em\\u003e and \\u003cem\\u003eper3\\u003c/em\\u003e were down-regulated and \\u003cem\\u003eper1\\u003c/em\\u003e was up-regulated in the estrous ovaries. CLOCK plays a key role in maintaining the circadian rhythm and activating downstream elements. Inhibition of CLOCK can inhibit cell growth and increase the rate of apoptosis (Li, et al. \\u003cspan citationid=\\\"CR30\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). In vivo experiments show that female mice is injected with CLOCK-shRNA have fewer oocytes, fewer litters, and a higher rate of apoptosis. The results indicate that clock plays an important role in fertility, and down-regulation of clock leads to cell apoptosis and decrease reproductive capacity (Li, et al. \\u003cspan citationid=\\\"CR30\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). It was proposed that the down-regulation of clock in the ovary at estrus stage might be related to the low litter size of the Xiang pig. The core clock gene period 1 (\\u003cem\\u003eper1\\u003c/em\\u003e) may be a prolific gene in Drosophila (Cushman, et al. \\u003cspan citationid=\\\"CR10\\\" class=\\\"CitationRef\\\"\\u003e2007\\u003c/span\\u003e). Female mice with the \\u003cem\\u003eper1\\u003c/em\\u003e mutation showed a normal number of implantation sites but reduced litter size (Pilorz and Steinlechner \\u003cspan citationid=\\\"CR46\\\" class=\\\"CitationRef\\\"\\u003e2008\\u003c/span\\u003e). \\u003cem\\u003ePer1\\u003c/em\\u003e mRNA locates in the secondary oocytes and follicles of ruminants (sheep, cattle) and found that there was no relationship between its transcription level and prolificacy, and this gene did not map to the known QTL region of ovulation rate in cattle (Cushman, et al. \\u003cspan citationid=\\\"CR10\\\" class=\\\"CitationRef\\\"\\u003e2007\\u003c/span\\u003e). Treatment with Progesterone for 1 hour could induce the expression of \\u003cem\\u003eper1\\u003c/em\\u003e mRNA in MCF-7 cells (Nakamura, et al. \\u003cspan citationid=\\\"CR41\\\" class=\\\"CitationRef\\\"\\u003e2010\\u003c/span\\u003e). In fact, the up-regulation of \\u003cem\\u003eper1\\u003c/em\\u003e was a common feature of many tissues in response to certain types of hormonal stress(For xample: luteinising hormone) (Mansuy, et al. \\u003cspan citationid=\\\"CR35\\\" class=\\\"CitationRef\\\"\\u003e2009\\u003c/span\\u003e; Mendoza, et al. \\u003cspan citationid=\\\"CR37\\\" class=\\\"CitationRef\\\"\\u003e2011\\u003c/span\\u003e; Smarr, et al. \\u003cspan citationid=\\\"CR58\\\" class=\\\"CitationRef\\\"\\u003e2013\\u003c/span\\u003e). This expression pattern indicated that PER1 may be the most sensitive effector in the biological clock system (Zhang, et al. \\u003cspan citationid=\\\"CR66\\\" class=\\\"CitationRef\\\"\\u003e2019\\u003c/span\\u003e). Some studies have found that estrus leads to changes in the expression time and amplitude of the biological clock genes (L\\u0026oacute;pez-Rodr\\u0026iacute;guez, et al. \\u003cspan citationid=\\\"CR32\\\" class=\\\"CitationRef\\\"\\u003e2019\\u003c/span\\u003e). Rising estrogen cause female animal to go into estrous stage (Ravinder, et al. \\u003cspan citationid=\\\"CR48\\\" class=\\\"CitationRef\\\"\\u003e2016\\u003c/span\\u003e).The estrus in Xiang pig might be speculated to reset the biological clock by up-regulating the \\u003cem\\u003eper1\\u003c/em\\u003e gene to response the stimulus of steroid hormone.In this study, \\u003cem\\u003eper3\\u003c/em\\u003e was down-regulated in the ovaries of estrus pigs. There is little information on the function of \\u003cem\\u003eper3\\u003c/em\\u003e gene in reproduction except for brain development (Noda, et al. \\u003cspan citationid=\\\"CR43\\\" class=\\\"CitationRef\\\"\\u003e2019\\u003c/span\\u003e). In contrast, the expression patterns of clock genes in other pig breeds are much different from that in Xiang pig. In the ovarian follicles of Large White pigs and Mi pigs, the core circadian clock genes were all down-regulated and not differentially expressed between estrus and diestrus periods (Chu, et al. \\u003cspan citationid=\\\"CR9\\\" class=\\\"CitationRef\\\"\\u003e2017\\u003c/span\\u003e). Compared with Large White pigs and Mi pigs, the characteristics of Xiang pig are specific, such as low birth rate and insignificant estrous behavior. It would be interesting to prove whether the expression profiles of core biological clock genes in ovary is a reason for the special reproductive traits such as estrus performance and litter size in different pig breeds.\\u003c/p\\u003e \\u003cp\\u003eThe biological clock system consists of an input pathway, a core oscillator and an output pathway. The post-translational modification of clock proteins is essential to maintain the accuracy and robustness of the evolutionarily conserved circadian clock (Gossan, et al. \\u003cspan citationid=\\\"CR14\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). Post-translational modification and degradation of the clock proteins are key steps in determining the length of the circadian clock cycle (Kennaway, et al. \\u003cspan citationid=\\\"CR24\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). Our research found that the products of most differentially expressed genes were related to the post-translational modification of the core clock protein. For example, PPP1CB and PPP1CA can reduce the phosphorylation of PER2, and affect the nuclear localization of the protein PER2, which may at least partially change the cycle and phase shift characteristics of the biological clock (Schmutz, et al. \\u003cspan citationid=\\\"CR53\\\" class=\\\"CitationRef\\\"\\u003e2011\\u003c/span\\u003e). Also, highlighted DYRK1A is the enzyme responsible for the phosphorylation of CRY2 at Ser557 and plays a key role in regulating the protein level of CRY2 (Kurabayashi, et al. \\u003cspan citationid=\\\"CR28\\\" class=\\\"CitationRef\\\"\\u003e2010\\u003c/span\\u003e). The casein kinases CSNK1D and CSNK1E phosphorylate the PER protein and provide a marker for subsequent degradation (Kennaway, et al. \\u003cspan citationid=\\\"CR24\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). The UBE3A binds and degrades BMAL1 in a ubiquitin ligase-dependent manner (Gossan, et al. \\u003cspan citationid=\\\"CR14\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). Furthermore, MAGEL2 regulates the ubiquitination and stability of CRY1, and changes its nuclear and cytoplasmic distribution (Devos J 2011). Lastly, SIRT1 deacetylates PER2 and BMAL1, thereby participates in biological rhythm regulation (Wang, et al. \\u003cspan citationid=\\\"CR62\\\" class=\\\"CitationRef\\\"\\u003e2016\\u003c/span\\u003e). The previous reports illustrate that the estrus cycle affects the localization and degradation of the core clock protein by changing their post-translational modification. And transcription and post-transcriptional regulation are the basis of clock system component activities, and post-transcriptional mechanisms accounts for more than half of the regulatory network (Koike, et al. \\u003cspan citationid=\\\"CR27\\\" class=\\\"CitationRef\\\"\\u003e2012\\u003c/span\\u003e). The results of this study further strengthened the view that the post-translational mechanism participated in the circadian gene regulatory network. In addition, we detected some differentially expressed genes from Xiang pig ovaries (Table S2) that might affect the expression of core clock genes at the transcriptional level. For example, TOP2A binds to the unique GC-rich open chromatin structure of the \\u003cem\\u003ebmal1\\u003c/em\\u003e promoter region, indicating that TOP2 on the \\u003cem\\u003ebmal1\\u003c/em\\u003e promoter affects transcription of \\u003cem\\u003ebmal1\\u003c/em\\u003e (Ogawa, et al. \\u003cspan citationid=\\\"CR44\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). Also, CBP /P300 and tissue-specific cofactors regulate CLOCK /BMAL1 transcription positively or negatively (Hosoda, et al. \\u003cspan citationid=\\\"CR18\\\" class=\\\"CitationRef\\\"\\u003e2009\\u003c/span\\u003e). Accordingly, ID (DNA binding inhibitor) is an important transcription repressor. Each ID protein contains a helix-loop-helix domain through which it can interact with bHLH protein. However, ID lacks the basic domain that allows binding to the E-box element, which leads to the ID being able to modify the transactivation of clock genes and CCG by interfering with the CLOCK-BMAL1 heterodimer, bHLH orange factor and other bHLH factors (Hou, et al. \\u003cspan citationid=\\\"CR19\\\" class=\\\"CitationRef\\\"\\u003e2009\\u003c/span\\u003e). These data indicate that the estrus cycle promoted the time change of the biological clock in reproductive tissues at the level of transcription and post-translational modifications (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig4\\\" class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003cp\\u003eIn addition, we detected several genes controlled by the circadian rhythm, which were differentially expressed in the ovaries of the two periods including \\u003cem\\u003eklf9, star, ptgs2\\u003c/em\\u003e. Other studies have shown that \\u003cem\\u003eklf9\\u003c/em\\u003e is a clock output gene, and CLOCK and BMAL1 complexes can bind and activate transcription at the 5' flanking region of \\u003cem\\u003eklf9\\u003c/em\\u003e, which is blocked by the co-expression of PER1. Also, KLF9 may play a role in regulating the effects of CLOCK/BMAL1 and in the expression of DBP, other clock and clock output genes, thereby changing the timing and amplitude of the circadian oscillations of gene transcription (Knoedler, et al. \\u003cspan citationid=\\\"CR26\\\" class=\\\"CitationRef\\\"\\u003e2020\\u003c/span\\u003e). The previous studies have found that steroidogenic enzymes (STAR), prostaglandin synthase (PTGS2), which are related to ovarian progesterone synthesis, are clock-controlled genes (CCGs). The E-box enhancer exists in the 5'-flanking region of star, which binds to the CLOCK-BMAL1 heterodimer and activates the transcription of gene \\u003cem\\u003estar\\u003c/em\\u003e (Nakao, et al. \\u003cspan citationid=\\\"CR42\\\" class=\\\"CitationRef\\\"\\u003e2007\\u003c/span\\u003e). Co-expression of the negative regulators PER and CRY attenuates this activation (Wu, et al. \\u003cspan citationid=\\\"CR64\\\" class=\\\"CitationRef\\\"\\u003e2011\\u003c/span\\u003e). The promoter of \\u003cem\\u003eptgs2\\u003c/em\\u003e has an E-box element and a REV-ERBα/RORα response element (RORE). The secretion of PGF2α can be balanced by the inhibition or stimulation of transcriptional regulation on REV-ERBα and BMAL1/CLOCK, respectively (Isayama, et al. \\u003cspan citationid=\\\"CR21\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). It is thus inferred that the peripheral circadian oscillator, such as ovary, could play an essential role in synchronizing local physiology through regulation of the expression of clock-controlled genes (Mohawk, et al. \\u003cspan citationid=\\\"CR40\\\" class=\\\"CitationRef\\\"\\u003e2012\\u003c/span\\u003e). In this study, \\u003cem\\u003estar\\u003c/em\\u003e and \\u003cem\\u003eptgs2\\u003c/em\\u003e in the ovaries of the estrus pigs were up-regulated by 6.42 and 5.55 times, respectively. It might be possible that the core clock genes might regulate the production of steroid hormones by controlling the ovarian-specific CCG, thereby affecting the estrus cycle of Xiang pigs.\\u003c/p\\u003e \\u003cp\\u003ePre-mRNA splicing is a basic biological process through which introns of nascent RNA are removed and exons are merged to form mature RNA, which is then translated into protein (Lee and Rio \\u003cspan citationid=\\\"CR29\\\" class=\\\"CitationRef\\\"\\u003e2015\\u003c/span\\u003e). Through alternative splicing (AS), different proteins are produced, resulting in a wider variety of cellular functions (Genov, et al. \\u003cspan citationid=\\\"CR13\\\" class=\\\"CitationRef\\\"\\u003e2019\\u003c/span\\u003e). Except for transcription or post-translational controls, alternative splicing also plays a key role in circadian rhythms in a cell-type-specific manner (Shakhmantsir, et al. \\u003cspan citationid=\\\"CR56\\\" class=\\\"CitationRef\\\"\\u003e2018\\u003c/span\\u003e). Recent research reports several cases, in which AS is involved in the regulation of biological clocks in plants, mice, and fruit flies (Sanchez, et al. \\u003cspan citationid=\\\"CR51\\\" class=\\\"CitationRef\\\"\\u003e2010\\u003c/span\\u003e; Petrillo, et al. \\u003cspan citationid=\\\"CR45\\\" class=\\\"CitationRef\\\"\\u003e2011\\u003c/span\\u003e; Hughes, et al. \\u003cspan citationid=\\\"CR20\\\" class=\\\"CitationRef\\\"\\u003e2012\\u003c/span\\u003e; McGlincy, et al. \\u003cspan citationid=\\\"CR36\\\" class=\\\"CitationRef\\\"\\u003e2012\\u003c/span\\u003e; Preu\\u0026szlig;ner, et al. \\u003cspan citationid=\\\"CR47\\\" class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e). In this study, \\u003cem\\u003eArntl\\u003c/em\\u003e only underwent differential alternative splicing, while \\u003cem\\u003eclock\\u003c/em\\u003e and \\u003cem\\u003eper1\\u003c/em\\u003e experienced both of differential alternative splicing and differential expression. The differential alternative splicing event occurred in \\u003cem\\u003eper1\\u003c/em\\u003e was the retention of introns 14\\u0026ndash;15. This event would result in the deletion of amino acid residues from positions 545 to 1278 in the encoded protein, and the deletion site was located in the core domain to activate the protein. This intron retention event occurred more frequently in the ovaries of Xiang pigs in the diestrus period, and the active PER1 protein decreased. During the estrus period, the occurrence of the intron retention event decreased, and the active PER1 protein level increased, which further increased the expression of PER1. It indicated that the core clock component PER1 may respond to the changes of sterol hormones in the ovary through transcriptional regulation and post-transcriptional regulation. \\u003cem\\u003eClock\\u003c/em\\u003e underwent the differential exon 17 skipping event. When this event occurred, the encoded protein would be deleted amino acid residues from positions 484 to 513. The deleted amino acids were not located in the core domain and the influence on the protein activity was still unclear. These findings suggested that AS might directly change the structure of core clock components to regulate the ovarian biological clock, and the AS regulation of biological clock genes might be independent in the gene expression. However, more evidence needs to be accumulated whether the AS events in these biological clock genes would turn on and turn off the circadian rhythm in pig.\\u003c/p\\u003e\"},{\"header\":\"Conclusion\",\"content\":\"\\u003cp\\u003eIn this study, we examined 90 clock-related genes expressed in Xiang pigs\\u0026rsquo; ovaries and found that 33 clock-related genes were differentially expressed between estrus and diestrous stages. We identified 44 differential splicing events from transcripts of 34 circadian clock-related genes. In addition, we also found that 20 genes including core clock components \\u003cem\\u003earntl\\u003c/em\\u003e and \\u003cem\\u003ecry1\\u003c/em\\u003e were differentially regulated at AS level only, and 14 genes including \\u003cem\\u003eper1\\u003c/em\\u003e and \\u003cem\\u003eclock\\u003c/em\\u003e were differentially regulated at both expression and AS level. We also demonstrated by RT-qPCR experiments that the core clock genes \\u003cem\\u003eper1\\u003c/em\\u003e, \\u003cem\\u003ecry, clock\\u003c/em\\u003e and \\u003cem\\u003earntl\\u003c/em\\u003e and the clock-related genes \\u003cem\\u003eppp1cb\\u003c/em\\u003e and \\u003cem\\u003entrk1\\u003c/em\\u003e were rhythmically expressed in ovaries of Xiang pig. The results suggested that the circadian clock in Xiang pig ovaries might play an key role in regulating ovarian physiological functions through transcriptional and post-transcriptional regulation.\\u003c/p\\u003e\"},{\"header\":\"Declarations\",\"content\":\"\\u003cp\\u003e\\u003cstrong\\u003eAuthor Contribution\\u003c/strong\\u003e\\u003c/p\\u003e\\n\\u003cp\\u003eLTT, JFW, and XQR conceived and designed the experiments; LTT performed the experiments; LTT and XN analyzed the data; LTT, XN, SL, and SHH contributed to samples and data analysis; and LTT wrote the paper. All authors reviewed the manuscript.\\u003c/p\\u003e\\n\\u003cp\\u003e\\u003cstrong\\u003eEthics Approval\\u003c/strong\\u003e\\u003c/p\\u003e\\n\\u003cp\\u003eAll animal procedures were approved by the Institutional Animal Care and Use Committee of Guizhou University (Approval number GZU-201709). The procedures were conducted in accordance with the National Research Council Guide for the Care and Use of Laboratory Animals.\\u003c/p\\u003e\\n\\u003cp\\u003e\\u003cstrong\\u003eFunding Statement\\u003c/strong\\u003e\\u003c/p\\u003e\\n\\u003cp\\u003eThis work was supported by the\\u0026nbsp;National Natural Science Foundation of China\\u0026nbsp;(31960641, 31672390), the Guizhou Province Science and Technology Innovation Team Building Special QKHPTRC[2019]5615, the Guizhou Province \\u0026ldquo;Hundred\\u0026rdquo; Innovative Talents Project [2016-4012], the Guizhou Agriculture Research program (QKHZC[2017]2585, QKHZC[2017]2587, QKHNY[2013]3073), and the National High Technology Research and Development Program of China (863 Program) [2013AA102503]. The funders monitored the ethics of experiments and provided instruments and financial funds for all of research in the manuscript.\\u003c/p\\u003e\\n\\u003cp\\u003e\\u003cstrong\\u003eConflicts of Interest\\u003c/strong\\u003e\\u003c/p\\u003e\\n\\u003cp\\u003eAll the authors have no conflict of interest with anyone about this manuscript. We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted. No conflict of interest exits in the submission of this manuscript, and manuscript is approved by all authors for publication. The work described has not been published previously, and not been submitted to any other journal, in whole or in part. All the authors listed have approved the manuscript that is enclosed.\\u003c/p\\u003e\\n\\u003cp\\u003e\\u0026nbsp;\\u003c/p\\u003e\"},{\"header\":\"References\",\"content\":\"\\u003col\\u003e\\n\\u003cli\\u003eAlvarez JD, Hansen A, Ord T, Bebas P, Chappell PE, Giebultowicz JM, Williams C, Moss S, 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\\u003c/li\\u003e\\n\\u003cli\\u003eBoden MJ, Varcoe TJ, Voultsios A, Kennaway DJ (2010) Reproductive biology of female Bmal1 null mice. Reproduction 139:1077-1090.https://doi.org/10.1530/REP-09-0523\\u003c/li\\u003e\\n\\u003cli\\u003eBrown-Grant K, Raisman G (1977) Abnormalities in reproductive function associated with the destruction of the suprachiasmatic nuclei in female rats. Proc R Soc Lond B Biol Sci 198:279-296.\\u003c/li\\u003e\\n\\u003cli\\u003eChacko E, Ranganathan S (2009) Genome-wide analysis of alternative splicing in cow: implications in bovine as a model for human diseases. BMC Genomics 10 Suppl 3:S11. https://doi.org/10.1186/1471-2164-10-S3-S11\\u003c/li\\u003e\\n\\u003cli\\u003eChen H, Chu G, Zhao L, Yamauchi N, Shigeyoshi Y, Hashimoto S, Hattori MA (2012) Rev-erb\\u0026alpha; regulates circadian rhythms and StAR expression in rat granulosa cells as identified by the agonist GSK4112. Biochem Biophys Res Commun 420:374-379. https://doi.org/10.1016/j.bbrc.2012.02.164\\u003c/li\\u003e\\n\\u003cli\\u003eChen H, Zhao L, Chu G, Kito G, Yamauchi N, Shigeyoshi Y, Hashimoto S, Hattori MA (2013) FSH induces the development of circadian clockwork in rat granulosa cells via a gap junction protein Cx43-dependent pathway. Am J Physiol Endocrinol Metab 304:E566-575. https://doi.org/10.1152/ajpendo.00432.2012\\u003c/li\\u003e\\n\\u003cli\\u003eChristian CA, Mobley JL, Moenter SM (2005) Diurnal and estradiol-dependent changes in gonadotropin-releasing hormone neuron firing activity. Proc Natl Acad Sci U S A 102:15682-15687. https://doi.org/10.1073/pnas.0504270102\\u003c/li\\u003e\\n\\u003cli\\u003eChu G, Misawa I, Chen H, Yamauchi N, Shigeyoshi Y, Hashimoto S, Hattori MA (2012) Contribution of FSH and triiodothyronine to the development of circadian clocks during granulosa cell maturation. Am J Physiol Endocrinol Metab 302:E645-653. https://doi.org/10.1152/ajpendo.00470.2011\\u003c/li\\u003e\\n\\u003cli\\u003eChu Q, Zhou B, Xu F, Chen R, Shen C, Liang T, Li Y, Schinckel AP (2017) Genome-wide differential mRNA expression profiles in follicles of two breeds and at two stages of estrus cycle of gilts. Sci Rep 7:5052. https://doi.org/10.1038/s41598-017-04336-x\\u003c/li\\u003e\\n\\u003cli\\u003eCushman RA, Allan MF, Jones SA, Rupp GP, Echternkamp SE (2007) Localization of Period 1 mRNA in the ruminant oocyte and investigations of its role in ovarian function. Anim Reprod Sci 99:93-105. https://doi.org/10.1016/j.anireprosci.2006.04.057\\u003c/li\\u003e\\n\\u003cli\\u003eDevos J WSV, Wevrick R (2011) Magel2, a Prader-Willi syndrome candidate gene, modulates the activities of circadian rhythm proteins in cultured cells. J Circadian Rhythms 9:12. https://doi.org/10.1186/1740-3391-9-12\\u003c/li\\u003e\\n\\u003cli\\u003eFahrenkrug J, Georg B, Hannibal J, Hindersson P, Gr\\u0026auml;s S (2006) Diurnal rhythmicity of the clock genes Per1 and Per2 in the rat ovary. Endocrinology 147:3769-3776. https://doi.org/10.1210/en.2006-0305\\u003c/li\\u003e\\n\\u003cli\\u003eGenov N, Basti A, Abreu M, Astaburuaga R, Rel\\u0026oacute;gio A (2019) A bioinformatic analysis identifies circadian expression of splicing factors and time-dependent alternative splicing events in the HD-MY-Z cell line. Sci Rep 9:11062. https://doi.org/10.1038/s41598-019-47343-w\\u003c/li\\u003e\\n\\u003cli\\u003eGossan NC, Zhang F, Guo B, Jin D, Yoshitane H, Yao A, Glossop N, Zhang YQ, Fukada Y, Meng QJ (2014) The E3 ubiquitin ligase UBE3A is an integral component of the molecular circadian clock through regulating the BMAL1 transcription factor. Nucleic Acids Res 42:5765-5775. https://doi.org/10.1093/nar/gku225\\u003c/li\\u003e\\n\\u003cli\\u003eGray GD, S\\u0026ouml;derstein P, Tallentire D, Davidson JM (1978) Effects of lesions in various structures of the suprachiasmatic-preoptic region on LH regulation and sexual behavior in female rats. Neuroendocrinology 25:174-191. https://doi.org/10.1159/000122739\\u003c/li\\u003e\\n\\u003cli\\u003eHe PJ, Hirata M, Yamauchi N, Hashimoto S, Hattori MA (2007a) The disruption of circadian clockwork in differentiating cells from rat reproductive tissues as identified by in vitro real-time monitoring system. J Endocrinol 193:413-420. https://doi.org/10.1677/JOE-07-0044\\u003c/li\\u003e\\n\\u003cli\\u003eHe PJ, Hirata M, Yamauchi N, Hashimoto S, Hattori MA (2007b) Gonadotropic regulation of circadian clockwork in rat granulosa cells. Mol Cell Biochem 302:111-118. https://doi.org/10.1007/s11010-007-9432-7\\u003c/li\\u003e\\n\\u003cli\\u003eHosoda H, Kato K, Asano H, Ito M, Kato H, Iwamoto T, Suzuki A, Masushige S, Kida S (2009) CBP/p300 is a cell type-specific modulator of CLOCK/BMAL1-mediated transcription. Mol Brain 2:34. https://doi.org/10.1186/1756-6606-2-34\\u003c/li\\u003e\\n\\u003cli\\u003eHou TY, Ward SM, Murad JM, Watson NP, Israel MA, Duffield GE (2009) ID2 (inhibitor of DNA binding 2) is a rhythmically expressed transcriptional repressor required for circadian clock output in mouse liver. J Biol Chem 284:31735-31745. https://doi.org/10.1074/jbc.M109.013961\\u003c/li\\u003e\\n\\u003cli\\u003eHughes ME, Grant GR, Paquin C, Qian J, Nitabach MN (2012) Deep sequencing the circadian and diurnal transcriptome of Drosophila brain. Genome Res 22:1266-1281. https://doi.org/10.1101/gr.128876.111\\u003c/li\\u003e\\n\\u003cli\\u003eIsayama K, Chen H, Yamauchi N, Hattori MA (2014) REV-ERBalpha inhibits the PTGS2 expression in bovine uterus endometrium stromal and epithelial cells exposed to ovarian steroids. J Reprod Dev 60:362-370. https://doi.org/10.1262/jrd.2014-040\\u003c/li\\u003e\\n\\u003cli\\u003eJohnson MH, Lim A, Fernando D, Day ML (2002) Circadian clockwork genes are expressed in the reproductive tract and conceptus of the early pregnant mouse. Reprod Biomed Online 4:140-145. https://doi.org/10.1016/s1472-6483(10)61931-1\\u003c/li\\u003e\\n\\u003cli\\u003eKarman BN, Tischkau SA (2006) Circadian clock gene expression in the ovary: Effects of luteinizing hormone. Biol Reprod 75:624-632. https://doi.org/10.1095/biolreprod.106.050732\\u003c/li\\u003e\\n\\u003cli\\u003eKennaway DJ, Varcoe TJ, Voultsios A, Salkeld MD, Rattanatray L, Boden MJ (2015) Acute inhibition of casein kinase 1delta/epsilon rapidly delays peripheral clock gene rhythms. Mol Cell Biochem 398:195-206. https://doi.org/10.1007/s11010-014-2219-8\\u003c/li\\u003e\\n\\u003cli\\u003eKianianmomeni A, Ong CS, R\\u0026auml;tsch G, Hallmann A (2014) Genome-wide analysis of alternative splicing in Volvox carteri. BMC Genomics 15:1117. https://doi.org/10.1186/1471-2164-15-1117\\u003c/li\\u003e\\n\\u003cli\\u003eKnoedler JR, Avila-Mendoza J, Subramani A, Denver RJ (2020) The Paralogous Kruppel-like Factors 9 and 13 Regulate the Mammalian Cellular Circadian Clock Output Gene Dbp. J Biol Rhythms 35:257-274. https://doi.org/10.1177/0748730420913205\\u003c/li\\u003e\\n\\u003cli\\u003eKoike N, Yoo SH, Huang HC, Kumar V, Lee C, Kim TK, Takahashi JS (2012) Transcriptional architecture and chromatin landscape of the core circadian clock in mammals. Science 338:349-354. https://doi.org/10.1126/science.1226339\\u003c/li\\u003e\\n\\u003cli\\u003eKurabayashi N, Hirota T, Sakai M, Sanada K, Fukada Y (2010) DYRK1A and glycogen synthase kinase 3beta, a dual-kinase mechanism directing proteasomal degradation of CRY2 for circadian timekeeping. Mol Cell Biol 30:1757-1768. https://doi.org/10.1128/MCB.01047-09\\u003c/li\\u003e\\n\\u003cli\\u003eLee Y, Rio DC (2015) Mechanisms and Regulation of Alternative Pre-mRNA Splicing. Annu Rev Biochem 84:291-323. https://doi.org/10.1146/annurev-biochem-060614-034316\\u003c/li\\u003e\\n\\u003cli\\u003eLi R, Cheng S, Wang Z (2015) Circadian clock gene plays a key role on ovarian cycle and spontaneous abortion. Cell Physiol Biochem 37:911-920. https://doi.org/10.1159/000430218\\u003c/li\\u003e\\n\\u003cli\\u003eLiu Y, Johnson BP, Shen AL, Wallisser JA, Krentz KJ, Moran SM, Sullivan R, Glover E, Parlow AF, Drinkwater NR, et al (2014) Loss of BMAL1 in ovarian steroidogenic cells results in implantation failure in female mice. Proc Natl Acad Sci U S A 111:14295-14300. https://doi.org/10.1073/pnas.1209249111\\u003c/li\\u003e\\n\\u003cli\\u003eL\\u0026oacute;pez-Rodr\\u0026iacute;guez D, Franssen D, Sevrin E, G\\u0026eacute;rard A, Balsat C, Blacher S, No\\u0026euml;l A, Parent AS (2019) Persistent vs Transient Alteration of Folliculogenesis and Estrous Cycle After Neonatal vs Adult Exposure to Bisphenol A. Endocrinology 160:2558-2572. https://doi.org/10.1210/en.2019-00505\\u003c/li\\u003e\\n\\u003cli\\u003eLuo ZY, Dai XL, Ran XQ, Cen YX, Niu X, Li S, Huang SH, Wang JF (2018) Identification and profile of microRNAs in Xiang pig testes in four different ages detected by Solexa sequencing. Theriogenology 117:61-71. https://doi.org/10.1016/j.theriogenology.2017.06.023\\u003c/li\\u003e\\n\\u003cli\\u003eMansano NDS, Paradela RS, Bohlen TM, Zanardi IM, Chaves FM, Silveira MA, Tavares MR, Donato J, Jr., Frazao R (2022) Vasoactive intestinal peptide exerts an excitatory effect on hypothalamic kisspeptin neurons during estrogen negative feedback. Mol Cell Endocrinol 542:111532. https://doi.org/10.1016/j.mce.2021.111532\\u003c/li\\u003e\\n\\u003cli\\u003eMansuy V, Risold PY, Glauser M, Fraichard A, Pralong FP (2009) Expression of the GABAA receptor associated protein Gec1 is circadian and dependent upon the cellular clock machinery in GnRH secreting GnV-3 cells. Mol Cell Endocrinol 307:68-76. https://doi.org/10.1016/j.mce.2009.02.029\\u003c/li\\u003e\\n\\u003cli\\u003eMcGlincy NJ, Valomon A, Chesham JE, Maywood ES, Hastings MH, Ule J (2012) Regulation of alternative splicing by the circadian clock and food related cues. Genome Biol 13:R54. https://doi.org/10.1186/gb-2012-13-6-r54\\u003c/li\\u003e\\n\\u003cli\\u003eMendoza J, Lopez-Lopez C, Revel FG, Jeanneau K, Delerue F, Prinssen E, Challet E, Moreau JL, Grundschober C (2011) Dimorphic effects of leptin on the circadian and hypocretinergic systems of mice. J Neuroendocrinol 23:28-38. https://doi.org/10.1111/j.1365-2826.2010.02072.x\\u003c/li\\u003e\\n\\u003cli\\u003eMiao X, Luo Q (2013) Genome-wide transcriptome analysis between small-tail Han sheep and the Surabaya fur sheep using high-throughput RNA sequencing. Reproduction 145:587-596. https://doi.org/10.1530/REP-12-0507\\u003c/li\\u003e\\n\\u003cli\\u003eMiller BH, Olson SL, Turek FW, Levine JE, Horton TH, Takahashi JS (2004) Circadian clock mutation disrupts estrous cyclicity and maintenance of pregnancy. Curr Biol 14:1367-1373. https://doi.org/10.1016/j.cub.2004.07.055\\u003c/li\\u003e\\n\\u003cli\\u003eMohawk JA, Green CB, Takahashi JS (2012) Central and peripheral circadian clocks in mammals. Annu Rev Neurosci 35:445-462. https://doi.org/10.1146/annurev-neuro-060909-153128\\u003c/li\\u003e\\n\\u003cli\\u003eNakamura TJ, Sellix MT, Kudo T, Nakao N, Yoshimura T, Ebihara S, Colwell CS, Block GD (2010) Influence of the estrous cycle on clock gene expression in reproductive tissues: effects of fluctuating ovarian steroid hormone levels. Steroids 75:203-212. https://doi.org/10.1016/j.steroids.2010.01.007\\u003c/li\\u003e\\n\\u003cli\\u003eNakao N, Yasuo S, Nishimura A, Yamamura T, Watanabe T, Anraku T, Okano T, Fukada Y, Sharp PJ, Ebihara S, et al (2007) Circadian clock gene regulation of steroidogenic acute regulatory protein gene expression in preovulatory ovarian follicles. Endocrinology 148:3031-3038. https://doi.org/10.1210/en.2007-0044\\u003c/li\\u003e\\n\\u003cli\\u003eNoda M, Iwamoto I, Tabata H, Yamagata T, Ito H, Nagata KI (2019) Role of Per3, a circadian clock gene, in embryonic development of mouse cerebral cortex. Sci Rep 9:5874. https://doi.org/10.1038/s41598-019-42390-9\\u003c/li\\u003e\\n\\u003cli\\u003eOgawa Y, Kawano Y, Yamazaki Y, Onishi Y (2014) Shikonin shortens the circadian period: possible involvement of Top2 inhibition. Biochem Biophys Res Commun 443:339-343. https://doi.org/10.1016/j.bbrc.2013.11.116\\u003c/li\\u003e\\n\\u003cli\\u003ePetrillo E, Sanchez SE, Kornblihtt AR, Yanovsky MJ (2011) Alternative splicing adds a new loop to the circadian clock. Commun Integr Biol 4:284-286. https://doi.org/10.4161/cib.4.3.14777\\u003c/li\\u003e\\n\\u003cli\\u003ePilorz V, Steinlechner S (2008) Low reproductive success in Per1 and Per2 mutant mouse females due to accelerated ageing? Reproduction 135:559-568. https://doi.org/10.1530/REP-07-0434\\u003c/li\\u003e\\n\\u003cli\\u003ePreu\\u0026szlig;ner M, Wilhelmi I, Schultz AS, Finkernagel F, Michel M, M\\u0026ouml;r\\u0026ouml;y T, Heyd F (2014) Rhythmic U2af26 alternative splicing controls PERIOD1 stability and the circadian clock in mice. Mol Cell 54:651-662. https://doi.org/10.1016/j.molcel.2014.04.015\\u003c/li\\u003e\\n\\u003cli\\u003eRavinder R, Kaipa O, Baddela VS, Singhal Sinha E, Singh P, Nayan V, Velagala CS, Baithalu RK, Onteru SK, Singh D (2016) Saliva ferning, an unorthodox estrus detection method in water buffaloes (Bubalus bubalis). Theriogenology 86:1147-1155. https://doi.org/10.1016/j.theriogenology.2016.04.004\\u003c/li\\u003e\\n\\u003cli\\u003eReppert SM, Weaver DR (2002) Coordination of circadian timing in mammals. Nature 418:935-941. https://doi.org/10.1038/nature00965\\u003c/li\\u003e\\n\\u003cli\\u003eRosenwasser AM, Turek FW (2015) Neurobiology of Circadian Rhythm Regulation. Sleep Med Clin 10:403-412. https://doi.org/10.1016/j.jsmc.2015.08.003\\u003c/li\\u003e\\n\\u003cli\\u003eSanchez SE, Petrillo E, Beckwith EJ, Zhang X, Rugnone ML, Hernando CE, Cuevas JC, Godoy Herz MA, Depetris-Chauvin A, Simpson CG, et al (2010) A methyl transferase links the circadian clock to the regulation of alternative splicing. Nature 468:112-116. https://doi.org/10.1038/nature09470\\u003c/li\\u003e\\n\\u003cli\\u003eSchibler U, Ripperger J, Brown SA (2003) Peripheral circadian oscillators in mammals: time and food. J Biol Rhythms 18:250-260. https://doi.org/10.1177/0748730403018003007\\u003c/li\\u003e\\n\\u003cli\\u003eSchmutz I, Wendt S, Schnell A, Kramer A, Mansuy IM, Albrecht U (2011) Protein phosphatase 1 (PP1) is a post-translational regulator of the mammalian circadian clock. PLoS One 6:e21325. https://doi.org/10.1371/journal.pone.0021325\\u003c/li\\u003e\\n\\u003cli\\u003eSellix MT (2015) Circadian clock function in the mammalian ovary. J Biol Rhythms 30:7-19. https://doi.org/10.1177/0748730414554222\\u003c/li\\u003e\\n\\u003cli\\u003eSen A, Sellix MT (2016) The Circadian Timing System and Environmental Circadian Disruption: From Follicles to Fertility. Endocrinology 157:3366-3373. https://doi.org/10.1210/en.2016-1450\\u003c/li\\u003e\\n\\u003cli\\u003eShakhmantsir I, Nayak S, Grant GR, Sehgal A (2018) Spliceosome factors target timeless (tim) mRNA to control clock protein accumulation and circadian behavior in Drosophila. Elife 7:e39821. https://doi.org/10.7554/eLife.39821\\u003c/li\\u003e\\n\\u003cli\\u003eShen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, Zhou Q, Xing Y (2014) rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A 111:E5593-5601. https://doi.org/10.1073/pnas.1419161111\\u003c/li\\u003e\\n\\u003cli\\u003eSmarr BL, Gile JJ, de la Iglesia HO (2013) Oestrogen-independent circadian clock gene expression in the anteroventral periventricular nucleus in female rats: possible role as an integrator for circadian and ovarian signals timing the luteinising hormone surge. J Neuroendocrinol 25:1273-1279. https://doi.org/10.1111/jne.12104\\u003c/li\\u003e\\n\\u003cli\\u003eTakahashi JS (2017) Transcriptional architecture of the mammalian circadian clock. Nat Rev Genet 18:164-179. https://doi.org/10.1038/nrg.2016.150\\u003c/li\\u003e\\n\\u003cli\\u003eTerenina E, Fabre S, Bonnet A, Monniaux D, Robert-Grani\\u0026eacute; C, SanCristobal M, Sarry J, Vignoles F, Gondret F, Monget P, et al (2017) Differentially expressed genes and gene networks involved in pig ovarian follicular atresia. Physiol Genomics 49:67-80. https://doi.org/10.1152/physiolgenomics.00069.2016\\u003c/li\\u003e\\n\\u003cli\\u003eVieyra E, Ram\\u0026iacute;rez DA, Lagunas N, C\\u0026aacute;rdenas M, Chavira R, Dami\\u0026aacute;n-Matsumura P, Trujillo A, Dom\\u0026iacute;nguez R, Morales-Ledesma L (2016) Unilaterally blocking the muscarinic receptors in the suprachiasmatic nucleus in proestrus rats prevents pre-ovulatory LH secretion and ovulation. Reproductive biology and endocrinology : RB\\u0026amp;E 14:34-34. https://doi.org/10.1186/s12958-016-0168-7\\u003c/li\\u003e\\n\\u003cli\\u003eWang RH, Zhao T, Cui K, Hu G, Chen Q, Chen W, Wang XW, Soto-Gutierrez A, Zhao K, Deng CX (2016) Negative reciprocal regulation between Sirt1 and Per2 modulates the circadian clock and aging. Sci Rep 6:28633. https://doi.org/10.1038/srep28633\\u003c/li\\u003e\\n\\u003cli\\u003eWiegand SJ, Terasawa E, Bridson WE, Goy RW (1980) Effects of discrete lesions of preoptic and suprachiasmatic structures in the female rat. Alterations in the feedback regulation of gonadotropin secretion. Neuroendocrinology 31:147-157. https://doi.org/10.1159/000123066\\u003c/li\\u003e\\n\\u003cli\\u003eWu YG, Bennett J, Talla D, Stocco C (2011) Testosterone, not 5alpha-dihydrotestosterone, stimulates LRH-1 leading to FSH-independent expression of Cyp19 and P450scc in granulosa cells. Mol Endocrinol 25:656-668. https://doi.org/10.1210/me.2010-0367\\u003c/li\\u003e\\n\\u003cli\\u003eYoshikawa T, Sellix M, Pezuk P, Menaker M (2009) Timing of the ovarian circadian clock is regulated by gonadotropins. Endocrinology 150:4338-4347. https://doi.org/10.1210/en.2008-1280\\u003c/li\\u003e\\n\\u003cli\\u003eZhang Y, Meng N, Bao H, Jiang Y, Yang N, Wu K, Wu J, Wang H, Kong S, Zhang Y (2019) Circadian gene PER1 senses progesterone signal during human endometrial decidualization. J Endocrinol. https://doi.org/10.1530/JOE-19-0284\\u003c/li\\u003e\\n\\u003cli\\u003eZheng B, Albrecht U, Kaasik K, Sage M, Lu W, Vaishnav S, Li Q, Sun ZS, Eichele G, Bradley A, et al (2001) Nonredundant roles of the mPer1 and mPer2 genes in the mammalian circadian clock. Cell 105:683-694. https://doi.org/10.1016/s0092-8674(01)00380-4\\u003c/li\\u003e\\n\\u003c/ol\\u003e\"}],\"fulltextSource\":\"\",\"fullText\":\"\",\"funders\":[],\"hasAdminPriorityOnWorkflow\":false,\"hasManuscriptDocX\":true,\"hasOptedInToPreprint\":true,\"hasPassedJournalQc\":\"\",\"hasAnyPriority\":false,\"hideJournal\":true,\"highlight\":\"\",\"institution\":\"\",\"isAcceptedByJournal\":false,\"isAuthorSuppliedPdf\":false,\"isDeskRejected\":\"\",\"isHiddenFromSearch\":false,\"isInQc\":false,\"isInWorkflow\":false,\"isPdf\":false,\"isPdfUpToDate\":true,\"isWithdrawnOrRetracted\":false,\"journal\":{\"display\":true,\"email\":\"info@researchsquare.com\",\"identity\":\"researchsquare\",\"isNatureJournal\":false,\"hasQc\":true,\"allowDirectSubmit\":true,\"externalIdentity\":\"\",\"sideBox\":\"\",\"snPcode\":\"\",\"submissionUrl\":\"/submission\",\"title\":\"Research Square\",\"twitterHandle\":\"researchsquare\",\"acdcEnabled\":true,\"dfaEnabled\":false,\"editorialSystem\":\"\",\"reportingPortfolio\":\"\",\"inReviewEnabled\":false,\"inReviewRevisionsEnabled\":true},\"keywords\":\"Circadian rhythm, Xiang pig, Ovary, Estrus cycle, Alternative splicing\",\"lastPublishedDoi\":\"10.21203/rs.3.rs-2593341/v1\",\"lastPublishedDoiUrl\":\"https://doi.org/10.21203/rs.3.rs-2593341/v1\",\"license\":{\"name\":\"CC BY 4.0\",\"url\":\"https://creativecommons.org/licenses/by/4.0/\"},\"manuscriptAbstract\":\"\\u003cp\\u003eThe biological clock has been studied to play a critical role in the reproductive system of various living organisms like swine. To examine the effects of estrus cycle on the expression of ovarian biological clock-related gene in Xiang pig, in this study, we analyzed the expression and alternative splicing of biological clock-related genes in ovary during estrus and diestrus periods. In total, we detected 90 clock-related genes expressed in the ovaries of the Xiang pigs and found 33 clock-related genes differentially expressed between estrous and diestrous stages. We identified 44 differential splicing events from the transcripts of 34 biological clock-related genes. Furthermore, we also found 20 genes including the core clock components, arntl and cry1 were differentially regulated only at AS level and 14 genes, including per1 and clock, were differentially regulated at both expression and AS levels. We also proved that the core clock genes per1, cry1, clock and arntl and the clock-related genes, ppp1cb and ntrk1 were rhythmically expressed in Xiang pig ovaries by RT-qPCR experiments. The results demonstrated that the biological clock in the ovaries of Xiang pigs might play an important role in regulating the ovarian physiological functions by the transcriptional and post-transcriptional regulation.\\u003c/p\\u003e\",\"manuscriptTitle\":\"The effects of estrus cycle on the expression of ovarian biological clock-related genes of Xiang pigs\",\"msid\":\"\",\"msnumber\":\"\",\"nonDraftVersions\":[{\"code\":1,\"date\":\"2023-02-21 02:43:12\",\"doi\":\"10.21203/rs.3.rs-2593341/v1\",\"editorialEvents\":[{\"type\":\"communityComments\",\"content\":0}],\"status\":\"published\",\"journal\":{\"display\":true,\"email\":\"info@researchsquare.com\",\"identity\":\"researchsquare\",\"isNatureJournal\":false,\"hasQc\":true,\"allowDirectSubmit\":true,\"externalIdentity\":\"\",\"sideBox\":\"\",\"snPcode\":\"\",\"submissionUrl\":\"/submission\",\"title\":\"Research Square\",\"twitterHandle\":\"researchsquare\",\"acdcEnabled\":true,\"dfaEnabled\":false,\"editorialSystem\":\"\",\"reportingPortfolio\":\"\",\"inReviewEnabled\":false,\"inReviewRevisionsEnabled\":true}}],\"origin\":\"\",\"ownerIdentity\":\"2b50a6aa-3652-4144-b659-0f7b03ac62e8\",\"owner\":[],\"postedDate\":\"February 21st, 2023\",\"published\":true,\"recentEditorialEvents\":[],\"rejectedJournal\":[],\"revision\":\"\",\"amendment\":\"\",\"status\":\"posted\",\"subjectAreas\":[],\"tags\":[],\"updatedAt\":\"2023-02-23T15:59:36+00:00\",\"versionOfRecord\":[],\"versionCreatedAt\":\"2023-02-21 02:43:12\",\"video\":\"\",\"vorDoi\":\"\",\"vorDoiUrl\":\"\",\"workflowStages\":[]},\"version\":\"v1\",\"identity\":\"rs-2593341\",\"journalConfig\":\"researchsquare\"},\"__N_SSP\":true},\"page\":\"/article/[identity]/[[...version]]\",\"query\":{\"redirect\":\"/article/rs-2593341\",\"identity\":\"rs-2593341\",\"version\":[\"v1\"]},\"buildId\":\"-HB7Z8yhvgn0wM9Nzuekk\",\"isFallback\":false,\"isExperimentalCompile\":false,\"dynamicIds\":[84888],\"gssp\":true,\"scriptLoader\":[]}","source_license":"CC-BY-4.0","license_restricted":false}