Transcriptomic Changes of Photoperiodic Response in the Hypothalamus Were Identified in Ovariectomized and Estradiol-Treated Sheep

Accurate timing of seasonal changes is an essential ability for an animal’s survival, and the change in the photoperiod is the key factor affecting reproductive seasonality in mammals. Emerging evidence has suggested that multiple hypothalamic genes participate in the photoperiod-induced regulation of reproductive activities in sheep, but the mechanism is still unclear. In this study, we initially examined the plasma level of two major reproductive hormones, namely, follicle-stimulating hormone (FSH) and prolactin (PRL), under different photoperiods in ovariectomized and estradiol-treated (OVX + E₂) sheep using radioimmunoassay (RIA). Of the two hormones, the concentration of PRL significantly increased with the extension of the photoperiod, while FSH showed the opposite trend. Subsequently, an examination of the transcriptomic variation between the short photoperiod (SP) and long photoperiod (LP) was conducted. Differential expression analyses and functional annotation showed that several key genes in the insulin secretion (VAMP2, PRKACB, PRKCG, and PLCB1), GnRH (MAPK13, CGA, CDC42, ATF4, and LHB) pathways, and circadian entrainment (KCNJ5, PER1, GNB2, MTNR1A, and RASD1), as well as numerous lncRNAs, including XR_173257.3, XR_173415.3, XR_001435315.1, XR_001024596.2, and XR_001023464.2, were shown potentially vital for the hypothalamic photoperiodic response. Four of the differentially expressed mRNAs and lncRNAs were validated by qPCR. The constructed mRNA–mRNA interaction networks further revealed that transcripts potentially participated in hypothalamic thyroid hormone synthesis, endocrine resistance, and neuroactive ligand–receptor interactions. The interactome analysis of lncRNAs and their targets implied that XR_173257.3 and its target arylalkylamine N-acetyltransferase (AANAT) and XR_173415.3 and its target TH might participate in the regulation of seasonal reproduction. Together, the changes in reproductive hormones and transcriptome will help to determine the important photoperiod-induced lncRNAs and mRNAs and provide a valuable resource for further research on reproductive seasonality in sheep.

Introduction

Seasonal animals can use changes in the photoperiod to guide seasonal reproduction to ensure the survival of the next generation. Several studies have demonstrated that most mammals have a highly accurate mechanism for photoperiod measurement and show dramatic changes in response to small photoperiod changes (Nishiwaki-Ohkawa and Yoshimura, 2016; Nakayama and Yoshimura, 2018). The annual cycles of mammalian reproductive activities are controlled by the hypothalamic–pituitary–gonadal axis (HPGA), which plays a vital role in reproduction in mammals from fetal development through puberty to sexual maturity (Plant, 2015; Gao, et al., 2018; Hill and Elias, 2018). In mammals, the eye is believed to be the only photoperiod receptor, and the light received by the eyes is transmitted to the pineal gland through the suprachiasmatic nucleus (SCN) and causes melatonin secretion (Lincoln, et al., 2006). As the target of melatonin and an important functional region in the brain, the mediobasal hypothalamus (MBH) can release gonadotropin-releasing hormone (GnRH) to activate the secretion of gonadotropins (luteinizing hormone (LH) and follicle-stimulating hormone (FSH)) and initiate a series of estrus and reproductive activities. In these processes, the change in the photoperiod is the original and key factor affecting hypothalamic function. As previously reported, melatonin can alter hundreds of genes in mammals, including the seasonal expression of thyroid-stimulating hormone β (TSHB) in the pars tuberalis (PT) and MBH, by targeting a special population of PT cells (Wood, et al., 2015; Lomet, et al., 2018). TSH in the PT connects melatonin with hypothalamic triiodothyronine (T3) via the variation of Dio2 concentration in the third ventricle, which may be crucial for photoperiod-induced seasonal timing (Wood and Loudon, 2014; Helfer, et al., 2019). Numerous gene expression patterns in the MBH which are associated with the initiation of photoperiod-induced secretion of LH were also studied in Japanese quail (Nakao, et al., 2008). In addition, insulin signaling was also widely detected in the hypothalamus, and the interruption of insulin signaling can lead to abnormalities in puberty in mice, which has been suggested to be related with puberty (Burcelin, et al., 2003; Qiu, et al., 2013). Moreover, a population of Kiss1 neurons in the arcuate nucleus (ARC) is also indispensable in seasonal reproduction in hamster and sheep (Smith, 2012; Harter, et al., 2018; He, et al., 2018).

In recent years, numerous lncRNAs have been found to play a crucial role in the reproductive regulation in mice (Taylor, et al., 2015), goats (Gao, et al., 2018), rats (Gao, et al., 2018), and sheep (Zhang, et al., 2019; Chen, et al., 2021). The hypothalamus plays an important role of the transmission center in the seasonal reproduction of animals. The role of lncRNAs in modulating the function of the hypothalamus is also well documented. For example, a novel lncRNA, GnRH1 enhancer-derived non-coding RNA (GnRH-E1 RNA), is expressed in the hypothalamus in both mice and rats, and the knockdown of GnRH-E1 RNA resulted in a significant decrease in the expression of GnRH1 in mice (Huang, et al., 2016). A previous study showed that RMST, a novel lncRNA, can direct binding to the transcription factor Sox2 to participate in the ontogeny of GnRH neurons and puberty in a Kallmann syndrome (KS) patient (Stamou, et al., 2019). In addition, our earlier hypothalamic studies found that several candidate lncRNAs, such as MSTRG.26777, MSTRG.105228, and MSTRG.95128, may regulate ovine reproduction by targeting nearby genes (Zhang, et al., 2019). Although the function of lncRNA in the hypothalamus has rarely been studied, we can speculate that it is important to further explore the function of the ovine hypothalamus.

To date, our understanding of how photoperiod induced the changes of the molecular neuroendocrine axis and reproductive seasonality remains limited, the majority of studies in quail, hamsters, and sheep focused on the comparison of differences in long and short photoperiods (LP vs. SP) (16 vs. 8 h light), such as the change of hormonal concentration and key gene or protein expression (Dardente, et al., 2019; Guh, et al., 2019). We lack a comprehensive view of the impact of increasing daylengths on the transcriptome level within the hypothalamus. In the present study, we established an animal model for sheep estrus research based on the ovariectomized (OVX), which has been used for the functional study of the mammalian hypothalamus in rats, mice, Siberian hamster, goats, and sheep (Lomet, et al., 2018). In this model, the serum levels of FSH and PRL can provide reliable information of photoperiodic change. On the other hand, high-throughput sequencing was used to investigate the transcriptomic changes in the hypothalamus. The findings will provide some new information for understanding the genetic basis and molecular mechanisms of the hypothalamus in ovine seasonal reproduction.

Materials and Methods

Animals and Sample Collection

A group of 21 Sunite ewes (35–40 kg, 3 y, clinically normal, and non-pregnant) were selected from Urat Middle Banner, Bayan Nur City, Inner Mongolia Autonomous Region, China and housed in a farm at Tianjin Institute of Animal Sciences, Tianjin (39°N latitude), China. All ewes were fed ad libitum and had free access to water. OVX + E₂ sheep and light control rooms were constructed as previously described (La, et al., 2020; He et al., 2021; Xia, et al., 2021). In brief, estradiol treatment was achieved with an inner diameter of 3.35 mm and an outer diameter of 4.65 mm, packed with 20 mg crystalline 17β-estradiol (Sigma Chemical Co., St. Louis, MO). Implants were inserted into the axillary region and designed to produce the circulated E₂ levels of approximately 3–5 pg/ml for 2 weeks (Smith, et al., 2007). Then, all ewes were maintained in one of three photoperiod-controlled rooms, SP: short photoperiod (8/16 h light–dark), LP: long photoperiod (8/16 h light–dark) and SP–LP: photoperiod from short transfer to long). Along with the photoperiodic treatments, slaughter was performed on SP42 or LP42, SP-LP3 (the third day of SP–LP), SP-LP7, SP-LP15, SP-LP21, and SP-LP42 (n = 3). After slaughtering, the hypothalamus was rapidly removed from the brain, flushed with PBS (pH 7.4), snap-frozen at liquid-nitrogen, and stored at −80°C for subsequent study. Simultaneously, jugular venous blood samples were collected twice weekly (24 times, 12 times for different photoperiods) from three ewes in the SP–LP room. After centrifugation, serum was collected and frozen until assayed for hormones. The hormone preparation and sample collection’s schematic diagrams are shown in Figure 1.

FIGURE 1

FIGURE 1. Hormone preparation and sample collection timeline. All OVX + E₂ ewes were assigned to three experimental groups in three rooms: short photoperiod (SP), long photoperiod (LP), and short photoperiod transfer to long photoperiod (SP-LP). A total of 24 pink triangles represent the day of blood sample collection in SP-LP groups, and seven green circles were the time of hypothalamus sample collection. In addition, the tissue and blood samples were collected at 8 a.m. and 2 p.m., respectively.

Hormone Assays

The iodine (¹²⁵I) human follicle-stimulating hormone radioimmunoassay kit and the iodine (¹²⁵I) prolactin radioimmunoassay kit (Beijing North Biotechnology Research Institute Co., Ltd., Beijing, China) were used for FSH and PRL assays by radioimmunoassay (RIA). All samples from different time points were measured in triplicate. The FSH kit had an assay sensitivity of <1.0 mIU/ml, and the PRL kit sensitivity was ≤40 uIU/ml.

RNA Isolation and Library Preparation and Sequencing

Total RNA was isolated from each sample using TRIzol reagent (Invitrogen, Carlsbad, CA, United States), and 1% agarose gels were used to detect the degradation and contamination prior to subsequent. In addition, RNA’s purity, concentration, and integrity were also detected by a NanoPhotometer^® spectrophotometer (IMPLEN, CA, United States), a Qubit^® RNA Assay Kit in Qubit^® 2.0 Fluorometer (Life Technologies, CA, United States), and the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, United States), respectively. Subsequently, libraries were generated using rRNA-depleted RNA with the NEBNext^® Ultra™ Directional RNA Library Prep Kit for Illumina^® (NEB, Ipswich, MA, United States) according to the manufacturer’s recommendations. Finally, the libraries were sequenced on an Illumina HiSeq 4000 platform, and 150 bp paired-end reads were generated.

Quality Control and Bioinformatics Analysis of Sequenced RNAs

Clean data were obtained by removing reads-containing adapter, reads-containing poly-N sequence, and low-quality reads from raw data. Simultaneously, the Q20, Q30, and GC contents of the clean data were calculated. All the following analyses were based on the clean data. The Ovis aries reference genome and gene model annotation files were downloaded from the genome website directly (Oar_v4.0). The index of the reference genome was built using bowtie2 (v2.2.8) (Sirén, et al., 2014). Reads were aligned to the reference genome using HISAT2 (v2.0.4). The mapped reads of each sample were assembled by StringTie (v1.3.1) (Pertea, et al., 2016).

To improve the statistical validity of the identification, a series of strict screening conditions were employed to identify lncRNAs. First, we removed transcripts that had fewer than two exons, and transcripts of less than 200 bp in length were discarded. The remaining transcripts were BLAST to known sheep lncRNAs using Cuffcompare, and lncRNAs that overlapped with the transcripts in the database were incorporated into the subsequent analysis. Next, transcripts with an FPKM ≥ 0.5 were retained after calculating the expression of each transcript by Cuffdiff. Finally, CNCI (Sun, et al., 2013), CPC (Kong, et al., 2007), and Pfam (El-Gebali, et al., 2019), three databases were used to identify the candidate lncRNAs. The intersections of the results from each program were defined as the novel lncRNA transcripts, and those identified by all three tools were defined as the final candidate lncRNAs for further analysis.

Differential Expression Analysis and Validation of Sequencing Data by qPCR

Cuffdiff (v2.1.1) was used to calculate FPKMs, which represent the gene expression level of both the lncRNAs and coding genes in each sample (Trapnell, et al., 2010). FPKMs were computed by summing the FPKMs of transcripts in each gene group. A criterion of absolute log2 (fold change) > 1 and p ≤ 0.05 between the two groups was used to identify differentially expressed genes in Cuffdiff.

Eight transcripts (four mRNAs and four lncRNAs) were selected for each group to validate the RNA-sequencing results by qRT-PCR. The specific primers are listed in Table 1 cDNA was synthesized from the RNA samples remaining after sequencing with the PrimeScript™ RT reagent Kit (TaKaRa, Dalian, China). qPCR was performed with a SYBR Green assay (TaKaRa, Dalian, China) on a Roche LightCycler 480 (Roche Applied Science, Mannheim, Germany). The qPCR mixture and program has been described in our previous study (Zhang, et al., 2019). All samples were examined in triplicate. Differentially expressed mRNA and lncRNA levels were normalized to β-actin to determine the relative expression using the 2^−ΔΔCt method (Livak and Schmittgen, 2000).

TABLE 1

TABLE 1. Primers of the validation of differentially expressed mRNAs and lncRNAs.

lncRNA Target Gene Prediction and Functional Annotation and Enrichment Analysis

To further elucidate the functions of DE-mRNAs and DE-lncRNAs, we predicted the target genes of the lncRNAs in cis and trans. To understand the functional roles of the target genes of the lncRNAs, we used GOseq R package to implement enrichment analysis. In addition, the differentially expressed protein-coding genes were also analyzed using GO. We employed KOBAS software to detect the enrichment of lncRNA target genes or differentially expressed genes in KEGG pathways. GO terms with -log10 (p-value) < 0.05 and KEGG pathways with p < 0.05 were considered significantly enriched.

Construction of mRNA–mRNA and lncRNA–mRNA Networks

To infer the functions of the photoperiod-based DE-lncRNAs and DE-mRNAs in the ovine hypothalamus, a protein–protein interaction (PPI) network between the protein-coding genes was constructed based on information from STRING. Additionally, we constructed a complementary pair network based on mRNA and mRNA as well as between mRNA and lncRNA by using Cytoscape (V3.6.1).

Statistical Analysis

The statistical evaluation of the experimental results was performed with Student’s t-test using SPSS 20.0 statistical software. All data are expressed as the means with standard error (SE). For hormone profiles, RIA data were analyzed by repeated measures ANOVA, and comparisons of hormones between different photoperiods and qPCR validation were performed by one-way ANOVA. p < 0.05 was considered statistically significant.

Results

Hormone Profile Analysis

RIA was used to assess the serum hormone levels of FSH and PRL in three ewes exposed to the changed photoperiod. The FSH data shown in Figure 2A mostly reflected the impact of the photoperiod, and the serum hormone levels differed significantly between SP and LP (p < 0.05) (Figure 2B). For PRL, the concentration increased with the photoperiodic transition from SP to LP (Figure 2C), and the difference also reached significant level (p < 0.01) (Figure 2D).

FIGURE 2

FIGURE 2. Hormonal profiling of the response to the different photoperiodic treatments. (A) Mean follicle-stimulating hormone (FSH) levels in blood serum of ewes sampled twice weekly throughout the experiment. (B) Mean of FSH for the SP and LP (12 time points each treatment). (C,D) Prolactin (PRL) profiles, the picture interpretation identical to those for FSH. Results were expressed as mean ± SE, * represents p < 0.05, ** represents p < 0.01, *** represents p < 0.001.

Summary of RNA-Sequencing Data

Twenty-one hypothalamus samples from seven groups (SP42, LP42, SP-LP3, SP-LP7, SP-LP15, SP-LP21, and SP-LP42, n = 3) were used to construct 21 cDNA libraries for sequencing. More than 268 million raw reads were generated for each group. The Q30 of the data was not less than 89%, and the GC contents of each library ranged from 49.34 to 55.81%. In addition, 303.38 Gb clean reads were retained and used in the following analysis after discarding reads with adapters or a poly-N content >10% and other low-quality reads. Approximately, 97% of the clean reads from each library were mapped to the sheep reference genome (Oar_v4.0) (Table 2), and 552,061 transcripts were assembled from these 21 libraries, which were assembled using Scripture and Cufflinks. After a rigorous 5-step screening, 11,965 putative non-coding transcripts were retained (Figure 3A), and the intersection of CNCI, CPC, and Pfam yielded the same result (Figure 3B), identifying 9,645 lincRNAs (80.6%), 2,320 antisense lncRNAs (19.4%), and 0 intronic lncRNAs (0.0%).

TABLE 2

TABLE 2. Summary of RNA-sequencing data.

FIGURE 3

FIGURE 3. Screening of the candidate lncRNAs in the ovine hypothalamus transcriptome. (A) 522,061 transcripts were assembled and 11,965 putative non-coding transcripts were retained for the next analysis after five steps filtering. (B) lncRNAs identification was also used in CPC, PFAM, and CNCI.

Comparison of mRNA and mRNA Characteristics

In this study, 13,859 lncRNAs and 44,437 mRNAs were identified. Our results indicated that the length and ORFs of lncRNAs were shorter than those of mRNAs (Figure 4A); the expression of lncRNAs was lower than that of mRNAs (Figure 4D); and lncRNAs tended to contain fewer exons (Figure 4B). These results were consistent with those of our previous studies (Zhang, et al., 2019). In addition, some transcripts of uncertain coding potential (TUCP) were expressed in the hypothalamus (Figure 4C).

FIGURE 4

FIGURE 4. Characterization of candidate lncRNA and mRNA. (A–C) Length, exon number, open reading frame (ORF) length distribution of mRNA, annotated lncRNA, and novel lncRNA. (D) violin plot of expression level about mRNAs, annotated-lncRNAs, and novel-lncRNAs.

Differential Expression Analysis and qPCR Verification

The criteria of |log2 (fold change) |>1 and p ≤ 0.05 were used to identify differentially expressed transcripts. The summary of DE transcripts is shown in Table 3, and the details of the differences in expression can be found in Supplementary Table S1 and Supplementary Table S2. Importantly, the differences between different phases of the photoperiod change need to be verified prior to further analysis. In this study, four DE-mRNAs and four DE-lncRNAs were selected from each compared group to validate the RNA-seq data, and the expression of these RNAs was analyzed by qPCR (Supplementary Figure S1). qPCR results showed that the expression levels of selected DE-mRNAs (PLCB1, PRKACB, SMAD4, BCL2, ATF4, RASD1, CLOCK, CREB1, IGF2, DI O 2, BHLHE41, LHB, ESR1, GNB2, TSHB, FBXL3, and PRL) from the six compared groups were all significantly different (Supplementary Figure S1A). Similarly, the expression of 19 DE-lncRNAs (LNC_004953, LNC_002738, LNC_008255, LNC_008321, XR_173415.3, XR_173257.3, XR_001043724.2, XR_001045234.2, XR_001023464.2, XR_001020627.1, XR_001034881.2, XR_001433797.1, XR_001433798.1, XR_001024596.2, XR_001434242.1, XR_001027187.1, XR_001039609.1, XR_001023520.1, and XR_001434385.1) also reached the level of significant difference (Supplementary Figure S1B), and the expression patterns of these genes and lncRNAs were consistent with the RNA-seq results.

TABLE 3

TABLE 3. Summary of differentially expressed lncRNAs and mRNAs.

Functional Annotation of Differentially Expressed lncRNAs and mRNAs in the Hypothalamus

To understand the functions of DE-mRNAs and DE-lncRNAs in greater depth, DE-mRNAs and DE-lncRNA-colocalized mRNAs were subjected to GO and KEGG pathway analyses. The top five significantly enriched GO terms were in molecular function (MF), cellular component (CC), and biological process (BP),three GO types for DE-mRNA. Compared with SP42, when the photoperiod changes from SP-LP3 to SP-LP42, blinding, RNA binding, and cell cycle items were significantly enriched (Figure 5, data in Supplementary Table S3). For DE-lncRNA, we used their targets to conduct GO enrichment and most of the significantly enriched GO terms were consistent with mRNAs. In addition, most of the significantly enriched GO terms participate in the regulation of biological and cellular processes (Figure 6, data in Supplementary Table S4). KEGG pathway analysis revealed that these differentially expressed mRNA target genes were enriched in several significant items, such as steroid hormone biosynthesis, dopaminergic synapse, and circadian entrainment (Figure 7A, data in Supplementary Table S5), KEGG pathway analyses of differentially expressed lncRNAs were also enriched in pathways mentioned earlier. In addition, insulin secretion, the thyroid hormone signaling pathway, and the Wnt signaling pathway, were also significantly enriched with the photoperiod changes (Figure 7B, data in Supplementary Table S6).

FIGURE 5

FIGURE 5. Top enriched Gene Ontology (GO) terms of differentially expressed mRNA between the short photoperiod and long photoperiod in the hypothalamus of Sunite sheep. The GO terms in (A–F) are from SP42 vs. LP42, SP42 vs. SP-LP3, SP42 vs. SP-LP7, SP42 vs. SP-LP15, SP42 vs. SP-LP21, and SP42 vs. SP-LP42, respectively.

FIGURE 6

FIGURE 6. Top enriched Gene Ontology (GO) terms of differentially expressed lncRNAs between the short photoperiod and long photoperiod in the hypothalamus of Sunite sheep. The GO terms in (A–F) are from SP42 vs. LP42, SP42 vs. SP-LP3, SP42 vs. SP-LP7, SP42 vs. SP-LP15, SP42 vs. SP-LP21, and SP42 vs. SP-LP42, respectively.

FIGURE 7

FIGURE 7. KEGG analysis of differentially expressed mRNAs and lncRNAs targets in hypothalamus about photoperiod responsiveness. (A) Five of the top KEGG enrichment pathway for differentially expressed mRNAs; (B) three to five of the top KEGG enrichment pathway for differentially expressed lncRNA targets. The longitudinal axis represents the enrichment pathways in different compared groups, and the horizontal axis represents the rich factor of these pathways. The spot size and color represent the number of differentially expressed genes and statistical significance of each pathway, respectively.

Furthermore, analyses of DE-mRNAs revealed that the key signaling pathways regulating ovine photoperiodic responsiveness and seasonal reproduction, including the thyroid hormone signaling, insulin secretion, Wnt signaling, circadian entrainment, and GnRH signaling pathways were detected after different lengths of photoperiodic treatment (Table 4). For example, compared with LP42, except for PRKACB and PLCB1, all genes were downregulated in the thyroid hormone signaling pathway, and several genes in circadian entrainment (PER1, RASD1, and PRKACB) were significantly changed with the photoperiod from SP to LP. In addition, with the photoperiod change from SP to LP and gradually from SP-LP3 to SP-LP42, eight photoperiod-related genes, namely, PRKACB, AKT3, ATF4, KCNJ5, CLOCK, PLCB1, SMAD4, and MAPK3, were upregulated in SP42, and 11 other photoperiod-related genes (VAMP2, BTRC, PER1, GNB2, MTNR1A, CDC42, ATF4, LHB, TSHB, ESR1, and RASD1) were downregulated. Information on all lncRNAs and their target genes is listed in Supplementary Table S7. In this study, we searched for target protein-coding genes of the differentially expressed lncRNAs to evaluate the potential regulatory functions of lncRNAs, and a total of 130, 75, 132, 101, 145, and 260 target genes of DE-lncRNA were predicted and classified (Supplementary Table S8). Interestingly, we detected some genes related to animal photoperiodic response and estrous seasonality, such as AANAT, TH, CAMK2G, and KCNJ5, which were in near to XR_173257.3, XR_173415.3, XR_001024596.2, and XR_001023464.2, respectively (Table 5).

TABLE 4

TABLE 4. Summary of differentially expressed mRNAs involved in photoperiodic responsiveness and ovine estrous seasonality.

TABLE 5

TABLE 5. Genomic association between lncRNAs and nearby genes related to photoperiodic responsiveness and ovine estrous seasonality.

mRNA–mRNA and lncRNA–mRNA Interaction Networks in Photoperiodic Responsiveness

Several transcripts were selected to construct mRNA–mRNA and lncRNA–mRNA networks impacted by the change of light from SP to LP (Figure 8). For the mRNA–mRNA interaction network, DE-mRNAs are grouped to insulin secretion (VAMP2, PRKACB, PRKCG, and PLCB1), thyroid hormone synthesis (TSHB, CGA, ATF4, and PRKACA), and estrogen secretion (AKT3, ATF4, MAPK3, and ESR1). DE-mRNAs-regulating hypothalamic hormone secretion and photoperiodic responsiveness were also enriched in GnRH (MAPK13, CAMK2D, CGA, CDC42, ATF4, and LHB), Wnt (SMAD4, RHOA, RAC1, and SFRP5) signaling pathways, and circadian entrainment (PRKACB, KCNJ5, PER1, GNB2, MTNR1A, and RASD1) (Figure 8A). Several DE-lncRNAs related to ovine hypothalamic hormone secretion and photoperiodic responsiveness were also selected to construct an lncRNA–mRNA interaction network (Figure 8B). Of these, DE-lncRNAs were highly correlated with steroid hormone biosynthesis, retinol metabolism, and circadian entrainment. The lncRNAs XR_173257.3 targets AANAT, and is potentially involved in the response to photoperiodic change. XR_173415.3 targets TH and IGF2; interestingly, thyroid hormones, encoded by TH, have been proven to play an important role in mammalian seasonal reproduction. In addition, XR_001435315.1 and XR_001023464.2 probably participate in ovine reproductive activity because of their target genes B4GALNT2 (FecL) and KANJ5.

FIGURE 8

FIGURE 8. Interaction network in the hypothalamus about ovine photoperiod responsiveness and estrous seasonality. (A) mRNA–mRNA interaction network involved in photoperiodic responsiveness, the arrow direction represents the targeting direction; (B) lncRNA–mRNA interaction network related to ovine estrus regulation with the photoperiodic change. Circles and “V” represent mRNAs and lncRNAs, red and green represent upregulated and downregulated transcripts, respectively, and the arrow indicates the direction of the target relationship.

Discussion

Mammals exhibit a sensitive photoperiodic response when daylight length reaches a critical boundary (Nishiwaki-Ohkawa and Yoshimura, 2016). Early reports have presented evidence that the sheep hypothalamus plays a vital role in photoperiod-induced reproductive regulation, and that responsiveness to the photoperiod is a key factor that affects ovine puberty onset (Plant, 2015), the estrous cycle transition, and reproductive hormone regulation (Smith, et al., 2007; Smith, 2012; Li, et al., 2019). In previous study, OVX animals have been studied to analyze the mechanism of hormonal feedback, and the location of the major functional area of hypothalamus was in goats (Muñoz, et al., 2017; Liang et al., 2021) and sheep (Goodman, et al., 2011; Weems, et al., 2017; Zhang, et al., 2019; Dardente, et al., 2019); several key genes related to the photoperiodic response, such as KISS1, RFRP, and DI O 2/DI O 3, were identified in these studies. However, few studies have focused on hypothalamic photoperiod regulation at the whole transcriptome level. To explore the effect of photoperiod changes on reproductive hormone regulation and the regulation of hypothalamic gene expression, we established an OVX + E₂ model in one of the typical seasonal estrus breeds.

Changes in Reproductive Hormones Under Different Photoperiodic Conditions

Several lines of evidence indicate that there is a cyclical rhythm in ovine reproduction, and its activation and quiescence can be clearly seen in gonadotrophin secretion in OVX + E₂ ewes (Lincoln, et al., 2006; Menassol, et al., 2012; Dardente, et al., 2019; Guh, et al., 2019). In this study, FSH and PRL were selected as the indicators of the regulation of gonadotropin secretion by photoperiodic changes in OVX + E₂ ewes. For FSH, the concentration was lowest in the anestrous period and significantly increased in the estrus period; Lincoln (2003) found that the concentration of ovine FSH was very low under LP, which is completely consistent with our findings. This means that the release of FSH was also regulated by photoperiod change or season, and a low concentration of FSH may be related to anestrus in ewes. PRL and FSH display opposite responses to the photoperiod, because they are driven by two different neuroendocrine axes (Dardente, 2012), one or several endocrine factors in the PRL axes have been proved to regulate the photoperiodic rhythms (Lincoln, et al., 2006; Dardente, et al., 2019), our findings also demonstrate that the concentration of PRL increases significantly upon photoperiod extension.

Differentially Expressed mRNAs in Hypothalamic Photoperiodic Response

In the current study, most omics studies of the hypothalamus have focused on differences in gene expression profiles based on the different estrus periods or some important economic traits with significant differences in animals (Gao, et al., 2018; Hou, et al., 2018; Gley, et al., 2019; Zhang, et al., 2019; Liang, et al., 2021). For example, it was found that the seasonal energy balance of Siberian hamsters was affected by the species-specific responsiveness of hypothalamic T3 (Bao, et al., 2019), and that the hypothalamus was also involved in the circadian of metabolism (Cedernaes, et al., 2019) and daylength response (Sáenz De Miera et al., 2017). Recent studies by David emphasized that photoperiodic transitions in endocrine output are regulated by a multistep signaling cascade within the mammalian neuroendocrine system (Hazlerigg, et al., 2018), and the differences of hypothalamic transcriptome levels can provide more information about photoperiodic response in ewes. In this study, compared with SP42, most DE-mRNAs in the SP-LP3 to SP-LP42 groups compared with the SP groups were associated with the signaling of circadian entrainment, insulin secretion, thyroid hormone, and estrogen (Table 4). As previously reported, SCN lesions lead to the alterations in the estrous cycle, sexual behavior, and tonic and phasic secretion of GnRH and gonadotropins in rodents, and it has been proven that Clock and Clock-related genes are the key factors in mammalian reproductive cycles (Silva and Domínguez, 2019). In our studies, the expression of Clock-related genes, such as MTNR1A, PER1, DI O 2, and ESR1, was changed with the artificial photoperiod from SP42 to LP21 and LP42, and was synchronously downregulated with gonadotropins FSH. Importantly, previous studies have demonstrated that these genes were significantly correlated with reproductive seasonality in goats and sheep (Wood, 2018; He, et al., 2019). Astoundingly, hypothalamus mRNAs that were altered in the SP–LP groups (SP-LP7, SP-LP21, and SP-LP42) compared to the SP42 group were all associated with circadian entrainment (Table 4), which was proven to initiate the maturation of human islets and insulin secretion (Alvarez-Dominguez, et al., 2019). Moreover, recent evidence has shown that insulin plays a vital role in the arrangement of reproductive activities through direct actions on insulin receptors; coincidentally, these receptors are present in several areas that are well-known to play a crucial role in reproduction, such as the ARC, ventromedial hypothalamic nucleus (VMH), and preoptic area (POA) (Caraty, et al., 1998; Brüning, et al., 2000; Cernea, et al., 2016). Several genes that were significantly differentially expressed upon photoperiod changes in our studies participate in insulin secretion (VAMP2, PRKACB, and PLCB1) and circadian entrainment (PRKACB, RASD1, and KCNJ5), and the change of previously mentioned key genes in circadian entrainment triggers insulin secretion; subsequently, the receptors in various hypothalamic functional areas are activated to regulate the ovine photoperiod-induced reproductive activities.

LncRNAs and Functional Pathways Related to Photoperiodic Response

At present, there is vast availability of evidence suggesting that the expression of lncRNAs can show high correlations with the expression of neighboring mRNAs. As stated in the introduction, lncRNA GnRH-E1 RNA targets its nearby gene GnRH1 to contribute to the regulation of GnRH neuronal maturation (Huang, et al., 2016). Furthermore, XLOC_1041225 and XLOC_446331 target neighboring IGFBP5, and regulate the onset of puberty in goats (Gao, et al., 2018). The ovine hypothalamic lncRNAs found in the present study had many similar characteristics with those of other mammals, such as shorter and fewer exons, shorter ORFs, and lower expression levels than protein-coding transcripts (Figure 4). As discussed in 4.3, genes involved in circadian entrainment and insulin secretion participate in photoperiod-induced reproduction in sheep. Several lncRNAs target genes were associated with the signaling of circadian entrainment (KCNJ5 and CAMK2G), prolactin (TH and AKT1), steroid hormone biosynthesis (UGT1A1), and tryptophan metabolism (AANAT) (Table 5 and Figure 6). As previously reported, the AANAT gene is an essential enzyme for the rhythmic synthesis of melatonin (Maronde, et al., 2011), and TH can also be connected with AANAT using dopa decarboxylase (DDC). Combined with the earlier studies, our results indicated that the XR_173257.3-AANAT and XR_173415.3-TH networks may be pivotal in rhythmic melatonin synthesis and secretion and subsequent reproductive seasonality in sheep. Moreover, KCNJ5, also known as a G protein-activated inward rectifier potassium channel 4 (GIRK4), has been confirmed to be associated with circadian rhythmicity in humans (Hsiao, et al., 2019). From our results, XR_001023464.2-KCNJ5 synchronous expression may be involved in the ovine circadian rhythm induced by photoperiod changes. Interestingly, we also detected the differentially expressed lncRNA XR_001435315.1, in which the target gene B4GALNT2 was proven to be associated with the ovine litter size in our earlier study (Guo, et al., 2018). All these findings suggest that lncRNAs may be a regulator involved in the photoperiodic response in sheep.

Conclusion

In conclusion, reproductive hormone detection and transcriptome analysis are useful for understanding the effects of photoperiod change on reproductive seasonality in ewes. Our studies, for the first time, identified pivotal changes in reproductive hormone levels and transcriptomic expression with photoperiod changes in OVX + E₂ ewes. FSH and PRL concentrations showed significant changes with the photoperiod. Moreover, we suggest that the transcriptome levels change to coordinate with circadian rhythms. We identified major genes (such as VAMP2, PRKACB, PLCB1, RASD1, and KCNJ5) and lncRNAs (such as XR_173257.3, XR_173415.3, and XR_001023464.2) as candidate transcripts for photoperiod responses, and these factors may interact to regulate hormonal changes and gene expression through the enriched pathways of insulin secretion, circadian entrainment, and thyroid hormone signaling. These findings proved that the application of OVX + E₂ is important for the study of photoperiodic effects in ewes; simultaneously, the differences in reproductive hormones and the transcriptome provide a valuable resource for further study of the molecular mechanisms of the hypothalamic-regulated reproductive seasonality in ewes.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by the Science Research Department of the Institute of Animal Sciences of CAAS, and the approval number is IAS2021-24.

Author Contributions

XH and RD performed the experiments, analyzed data, and wrote the first draft. XG, XC, MZ, XL, and QX contributed to the tissue and serum collection and pretreatment of the sample. XW provided analysis tools and data interpretation. XZ and JZ constructed the OVX model and managed animals. QL and MC contributed to the experimental design and manuscript revision.

Funding

This research was funded by the National Natural Science Foundation of China (32102508, 32172704), the Earmarked Fund for China Agriculture Research System of MOF and MARA (CARS-38), the Central Public-Interest Scientific Institution Basal Research Fund (No.2021-YWF-ZYSQ-14), the Agricultural Science and Technology Innovation Program of China (CAAS-ZDRW202106 and ASTIP-IAS13), the Natural Science Foundation of Tianjin (20JCQNJC00630), the Natural Science Foundation of Jilin Province (20210101376JC), and the China Postdoctoral Science Foundation (2021M703202).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.848144/full#supplementary-material

References

Alvarez-Dominguez, J. R., Donaghey, J., Rasouli, N., Kenty, J. H. R., Helman, A., Charlton, J., et al. (2020). Circadian Entrainment Triggers Maturation of Human In Vitro Islets. Cell Stem Cell 26, 108–122. doi:10.1016/j.stem.2019.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bao, R., Onishi, K. G., Tolla, E., Ebling, F. J. P., Lewis, J. E., Anderson, R. L., et al. (2019). Genome Sequencing and Transcriptome Analyses of the Siberian Hamster Hypothalamus Identify Mechanisms for Seasonal Energy Balance. Proc. Natl. Acad. Sci. USA 116, 13116–13121. doi:10.1073/pnas.1902896116

PubMed Abstract | CrossRef Full Text | Google Scholar

Brüning, J. C., Gautam, D., Burks, D. J., Gillette, J., Schubert, M., Orban, P. C., et al. (2000). Role of Brain Insulin Receptor in Control of Body Weight and Reproduction. Science 289, 2122–2125. doi:10.1126/science.289.5487.2122

PubMed Abstract | CrossRef Full Text | Google Scholar

Burcelin, R., Thorens, B., Glauser, M., Gaillard, R. C., and Pralong, F. P. (2003). Gonadotropin-Releasing Hormone Secretion from Hypothalamic Neurons: Stimulation by Insulin and Potentiation by Leptin. Endocrinology 144, 4484–4491. doi:10.1210/en.2003-0457

PubMed Abstract | CrossRef Full Text | Google Scholar

Caraty, A., Fabre-Nys, C., Delaleu, B., Locatelli, A., Bruneau, G., Karsch, F. J., et al. (1998). Evidence that the Mediobasal Hypothalamus Is the Primary Site of Action of Estradiol in Inducing the Preovulatory Gonadotropin Releasing Hormone Surge in the Ewe1. Endocrinology 139, 1752–1760. doi:10.1210/endo.139.4.5904

PubMed Abstract | CrossRef Full Text | Google Scholar

Cedernaes, J., Waldeck, N., and Bass, J. (2019). Neurogenetic Basis for Circadian Regulation of Metabolism by the Hypothalamus. Genes Dev. 33, 1136–1158. doi:10.1101/gad.328633.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Cernea, M., Phillips, R., Padmanabhan, V., Coolen, L. M., and Lehman, M. N. (2016). Prenatal Testosterone Exposure Decreases Colocalization of Insulin Receptors in Kisspeptin/neurokinin B/dynorphin and Agouti-Related Peptide Neurons of the Adult Ewe. Eur. J. Neurosci. 44, 2557–2568. doi:10.1111/ejn.13373

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Guo, X., He, X., Di, R., Zhang, X., Zhang, J., et al. (2021). Transcriptome Analysis Reveals Differentially Expressed Genes and Long Non-coding RNAs Associated with Fecundity in Sheep Hypothalamus with Different FecB Genotypes. Front. Cel Dev. Biol. 9, 633747. doi:10.3389/fcell.2021.633747

CrossRef Full Text | Google Scholar

Dardente, H., Lomet, D., Chesneau, D., Pellicer‐Rubio, M. T., and Hazlerigg, D. (2019). Discontinuity in the Molecular Neuroendocrine Response to Increasing Daylengths in Ile‐de‐France Ewes: Is Transient Dio2 Induction a Key Feature of Circannual Timing? J. Neuroendocrinol 31, e12775. doi:10.1111/jne.12775

PubMed Abstract | CrossRef Full Text | Google Scholar

Dardente, H. (2012). Melatonin-Dependent Timing of Seasonal Reproduction by the Pars Tuberalis: Pivotal Roles for Long Daylengths and Thyroid Hormones. J. Neuroendocrinol 24, 249–266. doi:10.1111/j.1365-2826.2011.02250.x

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Gebali, S., Mistry, J., Bateman, A., Eddy, S. R., Luciani, A., Potter, S. C., et al. (2019). The Pfam Protein Families Database in 2019. Nucleic Acids Res. 47 (D1), D427–D432. doi:10.1093/nar/gky995

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, X., Ye, J., Yang, C., Luo, L., Liu, Y., Ding, J., et al. (2018). RNA-seq Analysis of lncRNA-Controlled Developmental Gene Expression during Puberty in Goat & Rat. BMC Genet. 19, 19. doi:10.1186/s12863-018-0608-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gley, K., Murani, E., Trakooljul, N., Zebunke, M., Puppe, B., Wimmers, K., et al. (2019). Transcriptome Profiles of Hypothalamus and Adrenal Gland Linked to Haplotype Related to Coping Behavior in Pigs. Sci. Rep. 9, 13038. doi:10.1038/s41598-019-49521-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodman, R. L., Holaskova, I., Nestor, C. C., Connors, J. M., Billings, H. J., Valent, M., et al. (2011). Evidence that the Arcuate Nucleus Is an Important Site of Progesterone Negative Feedback in the Ewe. Endocrinology 152, 3451–3460. doi:10.1210/en.2011-0195

PubMed Abstract | CrossRef Full Text | Google Scholar

Guh, Y.-J., Tamai, T. K., and Yoshimura, T. (2019). The Underlying Mechanisms of Vertebrate Seasonal Reproduction. Proc. Jpn. Acad. Ser. B: Phys. Biol. Sci. 95, 343–357. doi:10.2183/pjab.95.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Wang, X., Liang, B., Di, R., Liu, Q., Hu, W., et al. (2018). Molecular Cloning of the B4GALNT2 Gene and its Single Nucleotide Polymorphisms Association with Litter Size in Small Tail Han Sheep. Animals 8, 160. doi:10.3390/ani8100160

PubMed Abstract | CrossRef Full Text | Google Scholar

Harter, C. J. L., Kavanagh, G. S., and Smith, J. T. (2018). The Role of Kisspeptin Neurons in Reproduction and Metabolism. J. Endocrinol. 238, R173–R183. doi:10.1530/JOE-18-0108

PubMed Abstract | CrossRef Full Text | Google Scholar

Hazlerigg, D., Lomet, D., Lincoln, G., and Dardente, H. (2018). Neuroendocrine Correlates of the Critical Day Length Response in the Soay Sheep. J. Neuroendocrinol 30, e12631. doi:10.1111/jne.12631

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Liu, Q., Li, X., Guo, X., Wang, X., Hu, W., et al. (2018). Molecular Cloning and Epigenetic Change Detection of Kiss1 during Seasonal Reproduction in Chinese Indigenous Sheep. Reprod. Fertil. Dev. 30, 734–743. doi:10.1071/RD17028

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Tao, L., Zhong, Y., Di, R., and Xia, Q. (2021). Photoperiod Induced the Pituitary Differential Regulation of LncRNAs and mRNAs Related to Reproduction in Sheep. PeerJ. 9, 10953. doi:10.7717/peerj.10953

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Zhang, Z., Liu, Q., and Chu, M. (2019). Polymorphisms of the Melatonin Receptor 1A Gene that Affects the Reproductive Seasonality and Litter Size in Small Tail Han Sheep. Reprod. Dom Anim. 54, 1400–1410. doi:10.1111/rda.13538

PubMed Abstract | CrossRef Full Text | Google Scholar

Helfer, G., Barrett, P., and Morgan, P. J. (2019). A Unifying Hypothesis for Control of Body Weight and Reproduction in Seasonally Breeding Mammals. J. Neuroendocrinol 31, e12680. doi:10.1111/jne.12680

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, J. W., and Elias, C. F. (2018). Neuroanatomical Framework of the Metabolic Control of Reproduction. Physiol. Rev. 98, 2349–2380. doi:10.1152/physrev.00033.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, Y., Hu, M., Zhou, H., Li, C., Li, X., Liu, X., et al. (2018). Neuronal Signal Transduction-Involved Genes in Pig Hypothalamus Affect Feed Efficiency as Revealed by Transcriptome Analysis. Biomed. Res. Int. 2018, 1–10. doi:10.1155/2018/5862571

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsiao, Y.-T., Chang, W.-A., Kuo, M.-T., Lo, J., Lin, H.-C., Yen, M.-C., et al. (2019). Systematic Analysis of Transcriptomic Profile of the Effects of Low Dose Atropine Treatment on Scleral Fibroblasts Using Next-Generation Sequencing and Bioinformatics. Int. J. Med. Sci. 16, 1652–1667. doi:10.7150/ijms.38571

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, P. P., Brusman, L. E., Iyer, A. K., Webster, N. J. G., and Mellon, P. L. (2016). A Novel Gonadotropin-Releasing Hormone 1 (Gnrh1) Enhancer-Derived Noncoding RNA Regulates Gnrh1 Gene Expression in GnRH Neuronal Cell Models. PLoS ONE 11, e0158597. doi:10.1371/journal.pone.0158597

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, L., Zhang, Y., Ye, Z.-Q., Liu, X.-Q., Zhao, S.-Q., Wei, L., et al. (2007). CPC: Assess the Protein-Coding Potential of Transcripts Using Sequence Features and Support Vector Machine. Nucleic Acids Res. 35, W345–W349. doi:10.1093/nar/gkm391

PubMed Abstract | CrossRef Full Text | Google Scholar

La, Y., He, X., Zhang, L., Di, R., Wang, X., Gan, S., et al. (2020). Comprehensive Analysis of Differentially Expressed Profiles of mRNA, lncRNA, and circRNA in the Uterus of Seasonal Reproduction Sheep. Genes 11, 301. doi:10.3390/genes11030301

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Lu, S., Ma, Q., Wan, P., Liu, C., Yang, H., et al. (2019). The Comparison of Reproductive Hormone Receptor Expressions of the Sheep Ovary and Hormone Concentrations in Two Chinese Breeds. Reprod. Dom Anim. 54, 892–901. doi:10.1111/rda.13441

CrossRef Full Text | Google Scholar

Liang, C., Han, M., Zhou, Z., Liu, Y., He, X., Jiang, Y., et al. (2021). Hypothalamic Transcriptome Analysis Reveals the Crucial microRNAs and mRNAs Affecting Litter Size in Goats. Front. Vet. Sci. 8, 747100. doi:10.3389/fvets.2021.747100

PubMed Abstract | CrossRef Full Text | Google Scholar

Lincoln, G. A., Clarke, I. J., Hut, R. A., and Hazlerigg, D. G. (2006). Characterizing a Mammalian Circannual Pacemaker. Science 314, 1941–1944. doi:10.1126/science.1132009

PubMed Abstract | CrossRef Full Text | Google Scholar

Lincoln, G. A. (2003). Prolactin Cycles in Sheep under Constant Photoperiod: Evidence that Photorefractoriness Develops within the Pituitary Gland Independently of the Prolactin Output Signal. Biol. Reprod. 69, 1416–1423. doi:10.1095/biolreprod.103.017673

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 25, 402–408. doi:10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Lomet, D., Cognié, J., Chesneau, D., Dubois, E., Hazlerigg, D., and Dardente, H. (2018). The Impact of Thyroid Hormone in Seasonal Breeding Has a Restricted Transcriptional Signature. Cell. Mol. Life Sci. 75, 905–919. doi:10.1007/s00018-017-2667-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Maronde, E., Saade, A., Ackermann, K., Goubran-Botros, H., Pagan, C., Bux, R., et al. (2011). Dynamics in Enzymatic Protein Complexes Offer a Novel Principle for the Regulation of Melatonin Synthesis in the Human Pineal Gland. J. Pineal Res. 51, 145–155. doi:10.1111/j.1600-079X.2011.00880.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Menassol, J.-B., Collet, A., Chesneau, D., Malpaux, B., and Scaramuzzi, R. J. (2012). The Interaction between Photoperiod and Nutrition and its Effects on Seasonal Rhythms of Reproduction in the Ewe1. Biol. Reprod. 86, 52. doi:10.1095/biolreprod.111.092817

PubMed Abstract | CrossRef Full Text | Google Scholar

Muñoz, A. L., Chesneau, D., Hernández, H., Bedos, M., Duarte, G., Vielma, J., et al. (2017). Sexually Active bucks Counterbalance the Seasonal Negative Feedback of Estradiol on LH in Ovariectomized Goats. Domest. Anim. Endocrinol. 60, 42–49. doi:10.1016/j.domaniend.2017.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakao, N., Ono, H., Yamamura, T., Anraku, T., Takagi, T., Higashi, K., et al. (2008). Thyrotrophin in the Pars Tuberalis Triggers Photoperiodic Response. Nature 452, 317–322. doi:10.1038/nature06738

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakayama, T., and Yoshimura, T. (2018). Seasonal Rhythms: the Role of Thyrotropin and Thyroid Hormones. Thyroid 28, 4–10. doi:10.1089/thy.2017.0186

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishiwaki-Ohkawa, T., and Yoshimura, T. (2016). Molecular Basis for Regulating Seasonal Reproduction in Vertebrates. J. Endocrinol. 229, R117–R127. doi:10.1530/JOE-16-0066

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., and Salzberg, S. L. (2016). Transcript-level Expression Analysis of RNA-Seq Experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi:10.1038/nprot.2016.095

PubMed Abstract | CrossRef Full Text | Google Scholar

Plant, T. M. (2015). Neuroendocrine Control of the Onset of Puberty. Front. Neuroendocrinology 38, 73–88. doi:10.1016/j.yfrne.2015.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, X., Dowling, A. R., Marino, J. S., Faulkner, L. D., Bryant, B., Brüning, J. C., et al. (2013). Delayed Puberty but normal Fertility in Mice with Selective Deletion of Insulin Receptors from Kiss1 Cells. Endocrinology 154, 1337–1348. doi:10.1210/en.2012-2056

PubMed Abstract | CrossRef Full Text | Google Scholar

Sáenz de Miera, C., Bothorel, B., Jaeger, C., Simonneaux, V., and Hazlerigg, D. (2017). Maternal Photoperiod Programs Hypothalamic Thyroid Status via the Fetal Pituitary Gland. Proc. Natl. Acad. Sci. USA 114, 8408–8413. doi:10.1073/pnas.1702943114

CrossRef Full Text | Google Scholar

Silva, C.-C., and Domínguez, R. (2019). Clock Control of Mammalian Reproductive Cycles: Looking beyond the Pre-ovulatory Surge of Gonadotropins. Rev. Endocr. Metab. Disord. 21, 149–163. doi:10.1007/s11154-019-09525-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Siren, J., Valimaki, N., and Makinen, V. (2014). Indexing Graphs for Path Queries with Applications in Genome Research. Ieee/acm Trans. Comput. Biol. Bioinf. 11, 375–388. doi:10.1109/TCBB.2013.2297101

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, J. T., Clay, C. M., Caraty, A., and Clarke, I. J. (2007). KiSS-1 Messenger Ribonucleic Acid Expression in the Hypothalamus of the Ewe Is Regulated by Sex Steroids and Season. Endocrinology 148, 1150–1157. doi:10.1210/en.2006-1435

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, J. T. (2012). The Role of Kisspeptin and Gonadotropin Inhibitory Hormone in the Seasonal Regulation of Reproduction in Sheep. Domest. Anim. Endocrinol. 43, 75–84. doi:10.1016/j.domaniend.2011.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Stamou, M., Ng, S.-Y., Brand, H., Wang, H., Plummer, L., Best, L., et al. (2019). A Balanced Translocation in Kallmann Syndrome Implicates a Long Noncoding RNA, RMST, as a GnRH Neuronal Regulator. J. Clin. Endocrinol. Metab. 105, e231–e244. doi:10.1210/clinem/dgz01110.1210/clinem/dgz011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, L., Luo, H., Bu, D., Zhao, G., Yu, K., Zhang, C., et al. (2013). Utilizing Sequence Intrinsic Composition to Classify Protein-Coding and Long Non-coding Transcripts. Nucleic Acids Res. 41, e166. doi:10.1093/nar/gkt646

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, D. H., Chu, E. T.-J., Spektor, R., and Soloway, P. D. (2015). Long Non-coding RNA Regulation of Reproduction and Development. Mol. Reprod. Dev. 82, 932–956. doi:10.1002/mrd.22581

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript Assembly and Quantification by RNA-Seq Reveals Unannotated Transcripts and Isoform Switching during Cell Differentiation. Nat. Biotechnol. 28, 511–515. doi:10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

Weems, P., Smith, J., Clarke, I. J., Coolen, L. M., Goodman, R. L., and Lehman, M. N. (2017). Effects of Season and Estradiol on KNDy Neuron Peptides, Colocalization with D2 Dopamine Receptors, and Dopaminergic Inputs in the Ewe. Endocrinology 158, 831–841. doi:10.1210/en.2016-1830

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, S. H., Christian, H. C., Miedzinska, K., Saer, B. R. C., Johnson, M., Paton, B., et al. (2015). Binary Switching of Calendar Cells in the Pituitary Defines the Phase of the Circannual Cycle in Mammals. Curr. Biol. 25, 2651–2662. doi:10.1016/j.cub.2015.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, S. H. (2018). How Can a Binary Switch within the Pars Tuberalis Control Seasonal Timing of Reproduction? J. Endocrinol. 239, R13–R25. doi:10.1530/JOE-18-0177

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, S., and Loudon, A. (2014). Clocks for All Seasons: Unwinding the Roles and Mechanisms of Circadian and Interval Timers in the Hypothalamus and Pituitary. J. Endocrinol. 222, R39–R59. doi:10.1530/JOE-14-0141

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, Q., Chu, M., He, X., Liu, Q., Zhang, X., Zhang, J., et al. (2021). Identification of Photoperiod-Induced LncRNAs and mRNAs in Pituitary Pars Tuberalis of Sheep. Front. Vet. Sci. 8, 644474. doi:10.3389/fvets.2021.644474

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Tang, J., Di, R., Liu, Q., Wang, X., Gan, S., et al. (2019). Comparative Transcriptomics Reveal Key Sheep (Ovis aries) Hypothalamus LncRNAs that Affect Reproduction. Animals 9, 152. doi:10.3390/ani9040152

PubMed Abstract | CrossRef Full Text | Google Scholar