Next Article in Journal
Transcriptome Profiling Reveals Role of MicroRNAs and Their Targeted Genes during Adventitious Root Formation in Dark-Pretreated Micro-Shoot Cuttings of Tetraploid Robinia pseudoacacia L.
Previous Article in Journal
Serum Periostin Level and Genetic Polymorphisms Are Associated with Vertebral Fracture in Chinese Postmenopausal Women
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 3
10.3390/genes13030440
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

by
Zhen Wan
†,
Hua Yang
†,
Yu Cai
,
Jianyu Ma
,
Peiyong Cheng
,
Zhibo Wang
,
Feng Wang
and
Yanli Zhang
*
Jiangsu Livestock Embryo Engineering Laboratory, College of Animal Science and Technology, Nanjing Agricultural University, NO. 1 Weigang, Nanjing 210095, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2022, 13(3), 440; https://doi.org/10.3390/genes13030440
Submission received: 4 February 2022 / Revised: 20 February 2022 / Accepted: 25 February 2022 / Published: 27 February 2022
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

1. Introduction

Hu sheep, a famous sheep breed in China, are characterized by high fertility, year-round estrus and fast growth. Therefore, Hu sheep have become an ideal model with which to explore the mechanism of high fecundity in sheep breeding. To date, many studies have shown that sheep fecundity is affected by multiple genes, including mutations in bone morphogenetic protein 15 (BMP15) [1], growth differentiation factor 9 (GDF9) [2], bone morphogenetic protein receptor type 1B (BMPR-1B) [3] and Booroola fecundity (FecB) [4]. Among them, FecB is the first major gene associated with multiple birth traits in sheep, which can increase the ovulation rate and litter size [5,6]. The FecB gene has been found in many prolific Chinese sheep breeds, such as Hu sheep, Duolang sheep and Little Tail Han (STH). It has been confirmed that FecB gene mutation plays a vital role in increasing the ovulation rate and prolificacy in ewes, and the frequency of the FecB allele in Hu sheep is up to 53% [7]. In addition, reproductive traits in less-prolific sheep can be improved by crossbreeding with sheep carrying the FecB gene. Therefore, only Hu sheep with FecB BB genotypes were selected in this study.
Long non-coding RNAs (lncRNAs) are non-coding RNAs with a length greater than 200 nt that are transcribed by RNA polymerase II, which cannot encode proteins [8]. Recent studies have shown that lncRNAs are involved in many important physiological processes in mammals, such as testicular maturation [9], oocyte maturation [10], embryo implantation [11] and gonadotropin secretion [12], which have become a hot spot in genetics research [13]. With the rapid development of RNA-Seq technology in the field of animal genetics and breeding [14,15], many lncRNAs that affect reproduction process have been initially investigated [16,17,18]. For example, Danila Cuomo et al. [19] screened 636 lncRNAs that were mainly enriched with PI3K-Akt and TGF-β signaling pathways, which were closely related to the regulation of follicle development and steroid-hormone biosynthesis. Chunyan Li et al. [20] identified six lncRNAs in the sheep pineal gland between the follicular phase and luteal phase, which could participate in hormone activities to affect sheep reproductive performance via regulation of its target genes. Yang et al. [9] identified some lncRNAs, and their target genes were enriched with male-gonad-development signaling pathways. Furthermore, lncRNAs have been identified as key regulators in the hypothalamic–pituitary–gonadal axis that is associated with reproduction. The previous research on the pituitary mainly focused on pituitary tumors, such as lncRNA H19 [21]. However, the study of lncRNAs in the normal pituitary is still limited. More recently, many lncRNAs were identified in the pituitary gland and found to have potential roles in regulating pituitary function. In mice, MIR205HG can bind to Zbtb20 to regulate the transcriptional activity of PIT1, thereby regulating the secretion of growth hormone and prolactin in the anterior pituitary [22]. LncRNA TCONS_00066406 participates in the development of immature and mature sheep pituitary glands by regulating Hydroxysteroid 17-β Dehydrogenase 12 (HSD17B12) [23]. LncRNA MSTRG.259847.2 was confirmed to be associated with high and low prolificacy of Hu sheep [12]. These studies indicated the presence and role of lncRNA in the pituitary gland. However, studies on the function of sheep lncRNA in the follicular phase and luteal phase are limited.
Therefore, this study focused on analyzing transcriptomics in two stages of the pituitary gland in Hu sheep with FecB BB genotypes: the follicular phase (H-F) and luteal phase (H-L). DE lncRNAs were subjected to gene-function analysis and constructed into a lncRNA–mRNA interaction network. More importantly, we identified the interaction of the candidate lncRNA XR_001039544.4 with its target gene LHB in vitro and explored the potential effect of lncRNA XR_001039544.4 downregulation on pituitary function. This study aimed to explore the candidate regulators of sheep prolificacy at the transcriptional level, which is essential for better understanding the molecular mechanisms by which lncRNAs regulate sheep reproduction at different physiological stages, and for providing insight into other animals.

2. Materials and Methods

2.1. Ethics Statement

The animal study was approved by the Animal Care and Use Committee of Nanjing Agricultural University and was strictly conducted under the Animal Experiments guidelines (Approval ID: SYXK2011-0036).

2.2. Animals and Samples Collection

Six healthy female Hu sheep with FecB BB genotypes (N = 6, litter size = 3) were selected from the Taizhou Sheep Farm (Jiangsu, China). All non-pregnant ewes, approximately 2–3 years old with similar weights, were provided with free access to water and feed under natural lighting.
Before the experiment, all sheep were implanted with vaginal sponge (CIDR) for 11 days, then the vaginal sponge was removed and the removal time was set to 0 h. A total of 6 sheep were divided into follicular-phase (H-F) and luteal-phase (H-L) groups. Then, three sheep were slaughtered at 45 h (follicular phase). The other three sheep were slaughtered on the 9th day after CIDR removal (luteal phase). Pituitary samples were immediately collected and stored at −80 °C for further analysis after slaughtering.

2.3. RNA Extraction, Strand-Specific Library Construction and Sequencing

The total RNA of pituitary samples was extracted using the Trizol reagent (Invitrogen, 15596026, Carlsbad, CA, USA) according to the manufacturer’s protocol. RNA purity, quality and concentration assays, strand-specific library construction and sequencing were performed as previously described [23] and sequenced using Illumina Novaseq6000 by Gene Denovo Biotechnology Co (Guangzhou, China). All sequencing data were outsourced to GENE DENOVO (Guangzhou, China).

2.4. Data Analysis and Transcript Assembly

To obtain high-quality clean reads, filtering of the reads was performed as previously described. Then, the clean reads were mapped to reference genome Oar_ v1.0 by TopHat2 [24] (version 2.1.1, Johns Hopkins University, Baltimore, MD, USA). The reconstruction of transcripts was carried out with Stringtie software (version 1.3.4) and the HISAT2 program (version 2.1.0) [25,26], and the transcripts with ≥200 bp in length and ≥2 exons were filtered. All reconstructed transcripts were aligned with the reference genome to identify novel transcripts.
Two software programs, CNCI [27] (version 2) and CPC [28] (version 0.9-r2) (http://cpc.cbi.pku.edu.cn/ (accessed on 8 June 2021) were used to assess the protein-coding potential of novel transcripts by default parameters. According to the location of lncRNAs in the genome, lncRNAs were divided into five types: bidirectional lncRNAs, intergenic lncRNAs, intronic lncRNAs, antisense lncRNAs and sense-overlapping lncRNAs.

2.5. Analysis of Differential Expression of mRNAs and Long Noncoding RNAs

The DE transcript mRNAs and lncRNAs were analyzed by DESeq2 [29] software between two different groups. The false-discovery rate (FDR) < 0.05 and absolute fold change (FC) ≥ 1.5 were the criteria for identifying differentially expressed genes (DEGs).

2.6. Functional Enrichment Analysis of DEGs

To identify significantly enriched metabolic pathways or signal pathways in DEGs, the calculation of p-values with FDR correction was performed, taking FDR ≤ 0.05 as a threshold. The GO term or pathway with FDR ≤ 0.05 was defined as a significantly enriched pathway in DEGs.

2.7. LncRNA-mRNA Interaction Analysis

We predicted the targeted genes of lncRNAs according to antisense, cis-act and trans-act analysis. The complementary correlation of antisense lncRNA and mRNA was predicted by RNAplex (http://www.tbi.univie.ac.at/RNA/RNAplex.1.html (accessed on 8 June 2021). LncRNAs that cis-regulate with their neighboring genes within 10 kb on the same allele were identified as cis-lncRNAs. The trans-regulation lncRNAs were based on the correlation coefficient between lncRNAs and mRNAs, and values > 0.9 were considered trans-acting. Based on these three interaction relationships, pituitary-function-related genes and potential lncRNAs were filtered and constructed into a gene co-expression network.

2.8. Cell Transfection and qRT-PCR Verification

To validate the potential role of candidate DE lncRNAs on pituitary function, primary pituitary cells were isolated on the basis of previous studies [12]. The siRNA of lncRNA XR_001039544.4 was synthesized by the GenePharma company (Shanghai, China). The sequences of siRNA are listed in Table S1.
Pituitary cells were seeded into 6-well plates for siRNA transfection. When pituitary cells reached 60–70% confluence, the siRNA was mixed with the Lipofectamine 3000 reagent (Invitrogen, No. L3000008, Carlsbad, CA, USA) and added to 6-well plates seeded with pituitary cells, according to the manufacturer’s protocol. The cells were harvested 24 h after transfection for RNA extraction. qPCR was performed as previously described [23]. The primer sequences are presented in Table S1.

2.9. FISH Analysis

Using the lncRNA XR_001039544.4 FISH probe mix (Cy3 labeled) for RNA FISH analysis to assess lncRNA localization in pituitary cells with a FISH kit (GenePharma, F12201/50, Shanghai, China), nuclei were stained with DAPI (Beyotime, C1002, Shanghai, China). The 18S was used as a positive control for the cytoplasm. The specific experimental steps followed the manufacturer’s instructions. Images were captured using an LSM 710 laser scanning confocal microscope.

2.10. ELISA Assay

Cell supernatant was collected, and LH was determined using the ELISA assay (Amresco, DRE-S9371, Shanghai, China) after transfection, according to the manufacturer’s protocol. The sensitivity by this assay is 0.1 mIU/mL. The values of inter- and intra-variation coefficients for the kit were <10%.

2.11. Edu Assay

Cell proliferation was detected using kFluor555 Click-iT EdU Assay Kit (KGA337-500, KeyGEN BioTECH, Jiangsu, China), according to the manufacturer’s instructions.

2.12. Statistical Analysis

All experiments were repeated at least three times. Statistical analysis of data was performed using SPSS software (24.0 Edition, Chicago, IL, USA). The significance of differences between samples was assayed by t-test for two-group comparisons and by one-way ANOVA used for three or multiple groups. The data were expressed as mean + standard error (SEM). p-value < 0.05 (*) and p-value < 0.01 (**) indicate significant and extremely significant differences, respectively.

3. Results

3.1. Overview of Sequencing Data in Sheep Pituitary Gland

A total of 11.9 Gb of clean data were obtained in this study. About 73.4–84.5 million clean reads were obtained per sample of pituitary. The clean data of each sample reached 99% and the percentage of the Q30 base was higher than 90%. All data showed that the sequencing data were highly reliable. In addition, the GC content of the pituitary was 48.33% and 48.71% in the H-F and H-L groups, respectively. The average number of reads for the six samples reached 80,050,123, the ratios of mapped reads and unmapped reads were 95.46% and 4.54%, and more than 80% of the reads were uniquely mapped in both groups. The percentage of reads that mapped to the intron, exon and intergenic region were 61.76%, 23.59%, and 14.65%, respectively (Table S2).

3.2. Identification of LncRNAs and mRNAs in Hu Sheep Pituitary Gland

CPC/CNCI software was used to screen and identify lncRNAs after mapping to the reference genome (Figure 1A). We identified 3259 common lncRNAs, including 1984 lincRNAs (60.9%), 578 antisense lncRNAs (17.7%), 388 sense lncRNAs (11.9%) and 65 intronic lncRNAs (2.0%) (Figure 1B) and 16,651 mRNAs in pituitary gland.
The expression level of lncRNAs was lower than mRNAs regardless of being in the H-F or H-L group (Figure 2A). As shown in Figure 2B, these lncRNAs and mRNAs were randomly assigned to the X-chromosome and 26 autosomes. The distribution trends of lncRNA and mRNA lengths were similar in the pituitary gland. The lncRNA and mRNA transcripts were mainly distributed over 3000 bp (Figure 2C). In addition, the majority of the lncRNAs had two exons while the mRNAs had five exons (Figure 2D). The open-reading-frame (ORF) length of the lncRNA transcripts was mainly concentrated within 100 bp, and the average ORF length of the lncRNA transcripts (about 76 bp on average) was shorter than that of the mRNA transcripts (about 953 bp on average, Figure 2E).

3.3. The Profiling and Verification of DE LncRNA and DEGs of Sheep Pituitary Gland

As shown in Figure 3, we identified 144 DE lncRNA transcripts (62 upregulated and 82 downregulated) and 557 DE mRNA transcripts (150 upregulated and 407 downregulated) between the H-F and H-L groups according to the criteria of FC > 1.5 and FDR < 0.05.
For further validation of the reliability of RNA-Seq data, four DE lncRNAs (MSTRG.5332.1, MSTRG.8834.3, MSTRG.14530.3, XR_003591934.1) and four DE mRNAs transcripts (IGFBP3, TP53INP1, FCF1, FSHB) were randomly selected and verified by using RT-qPCR (Figure 4). The results showed that the PCR results were consistent with the RNA-Seq results, which indicated that the data regarding lncRNA and protein-coding gene-transcript expression levels were reliable.

3.4. GO and KEGG Analysis of DEGs

All DEGs were tested in GO terms and subjected to KEGG enrichment analysis to clarify their underlying biological functions. In our study, 117 upregulated genes and 283 downregulated genes can be mapped to the GO database. The results showed that the most enriched GO term in the three parts were the cellular process, cell, and binding, respectively. As shown in Figure 5A, the top 20 most significantly enriched GO terms related to biological processes, including cellular process, single-organism process metabolic process, biological regulation, regulation of biological process, localization, multicellular-organismal process, response to stimulus, and so forth. In addition, GO terms related to pituitary function were enriched, including reproduction, reproductive process and growth (Table S3).
A total of 89 upregulated genes and 185 downregulated genes were annotated to 265 pathways according to KEGG enrichment analysis. As shown in Figure 5B, the highest enriched KEGG pathways in DEGs were mainly involved in the neuroactive ligand–receptor interaction, the cAMP signaling pathway and cell-adhesion molecules. Notably, some pathways were associated with pituitary function and reproduction, such as cAMP, PI3K-Akt, TGF-β, cGMP-PKG, MAPK, mTOR and serval hormone-related pathways (Figure S1). As shown in Figure 6, 39 downregulated and 22 upregulated genes that are associated with reproduction and pituitary functions were screened to construct the interaction network. The results showed that JUN, PGR, EGR1, DUSP1 and PPARA were the hub genes in the network.

3.5. LncRNA-Gene Interaction Network Construction

To further explore the potential functions of lncRNAs in the sheep pituitary gland, we predicted 80 target genes of DE lncRNAs based on the cis and trans RNA–RNA-interaction principle. A total of 36 DE target genes enriched for sheep-reproduction-related pathways (Table S4), such as growth, ovarian steroidogenesis, reproduction, gonadotropin secretion, and cAMP, PI3K-Akt and MAPK signaling pathways were screened to construct lncRNA–gene co-expression network. One DE lncRNA and one cis-target, as well as 42 DE lncRNAs and 35 trans-targets were involved in this network (Figure 7). XR_001042176.3 was cis-acting with EPHA7 by the sequence complementarity action. Meanwhile, XR_001039544.4 and MSTRG.17247.1 were trans-acting with LHB and PGR, respectively. MSTRG.10561.2 was trans-acting with DLL1, EGR1, GH1, NTRK2 and DUSP1. The downregulated DE gene IL2RG was regulated by one upregulated lncRNA and four downregulated lncRNAs. Moreover, the DE genes DMP1 and PRKACA were simultaneously regulated by four identical lncRNAs.

3.6. Verification and Characterization of XR_001039544.4 in Sheep Pituitary Cells

For validating the interaction between the lncRNAs and their target genes, and for exploring their potential function in pituitary cells, we screened lncRNA XR_001039544.4 for further experimentation, of which the target gene is LHB. First, the lncRNA XR_001039544.4 is an intergenic lncRNA located on chromosome 3. Then, we detected the localization of lncRNA XR_001039544.4 in pituitary cells, and a FISH experiment showed that lncRNA XR_001039544.4 was mainly located in the cytoplasm (Figure 8A). This suggests that lncRNA XR_001039544.4 may play a regulatory role by acting as a ceRNA regulatory mechanism. Further, we verified the expression of XR_001039544.4 and LHB in the pituitary at different periods, and qPCR results showed that lncRNA XR_001039544.4 and its target gene LHB were significantly increased the luteal phase (p < 0.01) (Figure 8C). Meanwhile, we detected the expression of lncRNA XR_001039544.4 in the HPO axis, and qPCR results showed that lncRNA XR_001039544.4 had the highest expression level in the pituitary gland (Figure 8B). Altogether, these data showed that lncRNA XR_001039544.4 could be a candidate lncRNA for exploring the molecular regulatory mechanisms of the luteal phase in the Hu sheep pituitary.

3.7. XR_001039544.4 Knockdown Downregulated the Secretion of LH in Pituitary Cells

To identify the function of lncRNA XR_001039544.4 in LH secretion, we transfected the sheep pituitary cells with lncRNA XR_001039544.4 siRNAs and detected the expression of XR_001039544.4 and LHB. The knockdown efficiency of siRNA1 was the highest among the three siRNAs. The expression level of its target gene, LHB, in the siRNA1-transfected group was significantly lower than in the control groups (Figure 9A,B). Furthermore, the ELISA results showed that the LH level significantly decreased (p < 0.05) in lncRNA XR_001039544.4-knockdown pituitary cells compared to the control group (Figure 9C). Overall, lncRNA XR_001039544.4 could regulate LH secretion by modulating LHB.
In addition, the PCNA mRNA expression level in the lncRNA XR_001039544.4-knockdown group was significantly lower than that in the control group (Figure 9E), and the Edu results showed that the cell-proliferation rate was significantly reduced (p < 0.01) in lncRNA XR_001039544.4-knockdown pituitary cells (Figure 9D,F). Taken together, these data indicated that XR_001039544.4 knockdown might significantly reduce LH secretion in pituitary cells by inhibiting cell proliferation.

4. Discussion

Reproductive traits are important factors affecting sheep production. In recent years, studies had investigated the mechanism of FecB affecting the reproductive performance of sheep [30,31,32,33]. FecB mutations can affect the response of oocytes to FSH and LH [34]. Increasing evidence also indicates that lncRNAs play important roles in sheep reproduction. However, the current research on lncRNA in the sheep pituitary gland is still relatively limited. FSH and LH secreted by the pituitary gland are key factors affecting follicular development and ovulation. Hu sheep are a Chinese indigenous breed with high fertility, and it is of great significance to comprehensively and deeply excavate the genes related to prolificacy. Therefore, we performed different lncRNA-expression profiles of the Hu sheep pituitary at different stages, including the follicular phase and the luteal phase, to explore the potential role of lncRNA in sheep prolificacy.
In this study, we investigated 3259 lncRNAs and 16,651 mRNAs in the sheep pituitary gland. We found that the chromosomal distribution trends of mRNAs and lncRNAs were similar. Notably, the distribution ratio of lncRNAs and mRNAs on chromosomes 1, 2 and 3 were greater than those on other chromosomes, which indicated that these three chromosomes may be closely related to pituitary function. Simultaneously, the length of lncRNAs and mRNAs showed the same trend. Moreover, the exon number, transcript and ORF length of mRNAs and lncRNAs in this study had a similar pattern to those in the immature and mature sheep pituitary gland [23]. In addition, we identified the expression levels of randomly selected lncRNAs, and the results were consistent with the RNA-Seq data. All the results suggested that the identified lncRNAs are reliable in the pituitary gland.
It is well known that the pituitary gland, one of the most important endocrine glands, can produce FSH, LH, GH and prolactin (PRL) hormones and so on. Overall, we screened 144 DE lncRNAs and 577 DE mRNAs during the follicular and luteal phases in the sheep pituitary gland. GO annotation and KEGG enrichment were performed on all transcripts. Among them, we found that hormone related genes FSHB, LHB, PRL, and GH were significantly differentially expressed. By the analysis of the interaction between lncRNAs and mRNAs, we finally found that there was one DE lncRNA that interacted with LHB and seven lncRNAs that interacted with GH. Based on the GO enrichment analysis, 61 DEGs were specifically enriched in reproduction, growth and pituitary function. Smad9, RBP4, ZBTB16, CCNF, ATP8B3, CCT4, ROBO2, PGR and NPPC were enriched in reproduction, which suggests that these genes might be associated with pituitary function in Hu sheep.
Moreover, the KEGG enrichment analysis showed that 23 signaling pathways were related to reproduction, such as the cAMP, GnRH, TNF, TGF-b and MAPK signaling pathway. Among them, the number of genes enriched in the cAMP signaling pathway was the largest. In the pituitary gland, the cAMP pathway plays an essential role in regulating cell proliferation and differentiation [35], as well as hormone secretion and production [36]. The cAMP pathway can also interact with the PI3K/AKT/mTOR, MAPK/ERK and WNT pathways to regulate the normal physiology of the pituitary gland [37]. Our study showed that these DE genes, such as BDNF, PPARA, NTRK2, ADCY9, and WNT5B, were involved in multiple signaling pathways at the same time. Studies have found that BDNF can participate in the regulation of the hypothalamic–pituitary–adrenal axis in mice. Knockout of BDNF in the mouse brain increases the release of adrenocorticotrophic hormone [38]. Meanwhile, BDNF and NTRK2 induces neuronal hyperexcitability in TG neurons and pain hypersensitivity in rats by activating p38 MAPK and inhibiting the PI3K pathway [39]. It also means that these genes might be considered as candidate genes for further study on pituitary function. Interestingly, EGR1, JUN, PGR, and PPARA were at the core of the interaction networks. However, whether there are any further links between these key genes, and how they cooperate with each other remain largely unknown and need further study.
LncRNAs are a heterogeneous type of ncRNAs that have been shown to display crucial regulatory functions in some biological processes [40]. Therefore, we constructed the lncRNA-target-gene-interaction networks. In our study, XR_001042176.3 can cis-regulate EPHA7, which was upregulated at luteal phase. Previous research demonstrated the involvement of EPHA7-EFNA5 signaling in the regulation of the LH and E2 negative-feedback pathways in the hypothalamus [41], highlighting the functional role of EPHA7 in female reproduction [42]. We also found that lncRNA MSTRG.8143.5, which was upregulated at the follicular phase, simultaneously targeted DLL1, EGR1, GH1, DUSP1 and NTRK2 in the sheep pituitary gland. The early-growth-response (EGR1) gene resides within the GnRH transcriptional network and belongs to the primary-response-gene family [43]. The detectable changes in EGR1 transcription occur as early as possible within 1 h of GnRH stimulation [44]. As the EGR1 protein accumulates, it regulates the transcription of secondary-like MAPK phosphatase 2 [45] and tertiary-like LHB response genes [46]. In addition, the JNK and ERK1/2 signaling pathway were involved in the GnRH induction of Egr1PT/mRNA in female-rat pituitary cells. Therefore, EGFR1 and MSTRG.8143.5 may be important for pituitary function at the follicular phase [47].
In pituitary cells, LH is a heterodimeric glycoprotein that shares a common α subunit (α-GSU) with other glycoprotein hormones. However, the β-subunits are receptors and hence are hormone-specific [48]. LH and FSH act on gonadal cells to regulate gonadal growth, differentiation, steroid production and gametogenesis. Studies have found that LHB-knockout female mice are hypogonadal and exhibit decreased levels of serum estradiol and progesterone [49]. In our study, according to the networks, lncRNA XR_001039544.4 was trans-acting on its target gene LHB, which could be considered as a candidate gene to investigate its function in the pituitary. LncRNA XR_001039544.4 had the highest expression in the pituitary gland of the HPO axis. After knockdown of lncRNA XR_001039544.4, the expression of LHB was significantly decreased, indicating that lncRNA XR_001039544.4 might play a regulatory role through LHB. As a key factor in the regulation of gene expression, lncRNAs can produce different modes of action according to the subcellular localization pattern [50,51]. LncRNAs can act as cis or trans in the nucleus to silence or activate specific genes [52]. In contrast, lncRNAs appear in the cytoplasm and can participate in a variety of molecular mechanistic processes, including mRNA stability and translation regulation, and are mediated by protein modifications, miRNA precursors or competing endogenous RNAs [53]. The FISH experiment results showed that lncRNA XR_001039544.4 was mainly located in the cytoplasm, which means that it may act as competing endogenous RNAs (ceRNAs) for miRNAs to regulate LHB expression. By downregulating the expression of the lncRNA XR_001039544.4, the PCNA gene was found to be significantly reduced, indicating that lncRNA XR_001039544.4 could maintain proliferation in pituitary cells. Taken together, the DE lncRNAs identified in this study might cooperate with their target genes and DEGs to regulate pituitary functions.

5. Conclusions

In summary, the transcriptomic pituitary-gland study revealed the differential regulation of lncRNAs and mRNAs related to the follicular phase and luteal phase in Hu sheep with high prolificacy. We screened a set of lncRNAs and genes associated with reproduction and pituitary functions. According to the GO and KEGG databases, the DEGs were annotated with multiple physiological processes related to the pituitary gland such as hormone synthesis and regulation, growth, and the reproductive process. Additionally, we predicted the potential function of DE lncRNAs and constructed the lncRNA-gene-interaction network in this study to provide a valuable resource of candidate lncRNAs in the pituitary gland. Furthermore, these differential mRNA- and lncRNA-expression profiles provide valuable resources for determining the molecular mechanism of sheep prolificacy.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13030440/s1, Figure S1: Functional Enrichment Analysis for pituitary function and reproduction of DE mRNAs. Table S1: The list of primer sequences for lncRNAs and mRNAs. Table S2: Statistical results of reads in each pituitary sample. Table S3 The GO enrichment analysis of DE mRNAs. Table S4: The list of LncRNA–gene interaction (cis and trans RNA–RNA interaction).

Author Contributions

Writing—original manuscript, Z.W.(Zhen Wan) and H.Y.; Data analysis, Z.W.(Zhen Wan); Samples collection, J.M. and Z.W.(Zhibo Wang); Software, Y.C. and P.C.; Experiment design, Y.Z. and F.W.; Writing—review, Z.W.(Zhen Wan), H.Y. and Y.Z. All listed authors have made substantial contributions to the research and the publication of the results. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National key research and development program (No. 2021YFF1000702). National Natural Science Foundation of China (No. 32072737).

Institutional Review Board Statement

All experimental procedures were approved by the Animal Care and Use Committee of Nanjing Agricultural University and were strictly conducted in the Animal Experiments guidelines (Approval ID: SYXK2011-0036; Approval Date: 2011).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data in this study are available, if requested.

Acknowledgments

The authors would like to express the sincere thanks to all the members of Jiangsu Livestock Embryo Engineering Laboratory for sample collection and technical assistance during the preparation of this research.

Conflicts of Interest

The authors have declared that they have no conflict of interest.

References

  1. Estienne, A.; Lahoz, B.; Jarrier, P.; Bodin, L.; Folch, J.; Alabart, J.L.; Fabre, S.; Monniaux, D. BMP15 regulates the inhibin/activin system independently of ovulation rate control in sheep. Reproduction 2017, 153, 395–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wang, F.; Chu, M.; Pan, L.; Wang, X.; He, X.; Zhang, R.; Tao, L.; La, Y.; Ma, L.; Di, R. Polymorphism detection of GDF9 gene and its association with litter size in Luzhong mutton sheep (Ovis aries). Animals 2021, 11, 571. [Google Scholar] [CrossRef] [PubMed]
  3. Abdoli, R.; Zamani, P.; Deljou, A.; Rezvan, H. Association of bmpr-1b and gdf9 genes polymorphisms and secondary protein structure changes with reproduction traits in mehraban ewes. Gene 2013, 524, 296–303. [Google Scholar] [CrossRef] [PubMed]
  4. Wilson, T.; Wu, X.Y.; Juengel, J.L.; Ross, I.K.; Lumsden, J.M.; Lord, E.A.; Dodds, K.G.; Walling, G.A.; McEwan, J.C.; O’Connell, A.R.; et al. Highly prolific booroola sheep have a mutation in the intracellular kinase domain of bone morphogenetic protein IB receptor (ALK-6) that is expressed in both oocytes and granulosa cells. Biol. Reprod. 2001, 64, 1225–1235. [Google Scholar] [CrossRef]
  5. Mulsant, P.; Lecerf, F.; Fabre, S.; Schibler, L.; Monget, P.; Lanneluc, I.; Pisselet, C.; Riquet, J.; Monniaux, D.; Callebaut, I.; et al. Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc. Natl. Acad. Sci. USA 2001, 98, 5104–5109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Qi, M.Y.; Xu, L.Q.; Zhang, J.N.; Li, M.O.; Lu, M.H.; Yao, Y.C. Effect of the Booroola fecundity (FecB) gene on the reproductive performance of ewes under assisted reproduction. Theriogenology 2020, 142, 246–250. [Google Scholar] [CrossRef]
  7. Wang, W.; Liu, S.; Li, F.; Pan, X.; Li, C.; Zhang, X.; Ma, Y.; La, Y.; Xi, R.; Li, T. Polymorphisms of the ovine BMPR-IB, BMP-15 and FSHR and their associations with litter size in two Chinese indigenous sheep breeds. Int. J. Mol. Sci. 2015, 16, 11385–11397. [Google Scholar] [CrossRef] [Green Version]
  8. Wilusz, J.E.; Sunwoo, H.; Spector, D.L. Long noncoding RNAs: Functional surprises from the RNA world. Genes Dev. 2009, 23, 1494–1504. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, H.; Wang, F.; Li, F.; Ren, C.; Pang, J.; Wan, Y.; Wang, Z.; Feng, X.; Zhang, Y. Comprehensive analysis of long noncoding RNA and mRNA expression patterns in sheep testicular maturation. Biol. Reprod. 2018, 99, 650–661. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, X.; Pan, J.; Wu, B.; Teng, X. Construction and analysis of a lncRNA (PWRN2)-mediated ceRNA network reveal its potential roles in oocyte nuclear maturation of patients with PCOS. Reprod. Biol. Endocrinol. 2018, 16, 73. [Google Scholar] [CrossRef]
  11. Xu, H.; Zhou, M.; Cao, Y.; Zhang, D.; Han, M.; Gao, X.; Xu, B.; Zhang, A. Genome-wide analysis of long noncoding RNAs, microRNAs, and mRNAs forming a competing endogenous RNA network in repeated implantation failure. Gene 2019, 720, 144056. [Google Scholar] [CrossRef]
  12. Zheng, J.; Wang, Z.; Yang, H.; Yao, X.; Yang, P.; Ren, C.; Wang, F.; Zhang, Y. Pituitary transcriptomic study reveals the differential regulation of lncRNA s and mRNAs related to prolificacy in different FecB genotyping sheep. Genes 2019, 10, 157. [Google Scholar] [CrossRef] [Green Version]
  13. Spizzo, R.; Almeida, M.I.; Colombatti, A.; Calin, G.A. Long non-coding RNAs and cancer: A new frontier of translational research? Oncogene 2012, 31, 4577–4587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ibeagha-Awemu, E.M.; Li, R.; Dudemaine, P.L.; Do, D.N.; Bissonnette, N. Transcriptome analysis of long non-coding RNA in the bovine mammary gland following dietary supplementation with linseed oil and safflower oil. Int. J. Mol. Sci. 2018, 19, 3610. [Google Scholar] [CrossRef] [Green Version]
  15. Standaert, L.; Adriaens, C.; Radaelli, E.; Van Keymeulen, A.; Blanpain, C.; Hirose, T.; Nakagawa, S.; Marine, J.C. The long noncoding RNA neat1 is required for mammary gland development and lactation. RNA 2014, 20, 1844–1849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Li, Q.; Qiao, J.; Zhang, Z.; Shang, X.; Chu, Z.; Fu, Y.; Chu, M. Identification and analysis of differentially expressed long non-coding RNAs of Chinese Holstein cattle responses to heat stress. Anim. Biotechnol. 2020, 31, 9–16. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, B.; Jiao, B.; Ge, W.; Zhang, X.; Wang, S.; Zhao, H.; Wang, X. Transcriptome sequencing to detect the potential role of long non-coding RNAs in bovine mammary gland during the dry and lactation period. BMC Genom. 2018, 19, 605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zheng, X.; Ning, C.; Zhao, P.; Feng, W.; Jin, Y.; Zhou, L.; Yu, Y.; Liu, J. Integrated analysis of long noncoding RNA and mRNA expression profiles reveals the potential role of long noncoding RNA in different bovine lactation stages. J. Dairy Sci. 2018, 101, 11061–11073. [Google Scholar] [CrossRef] [Green Version]
  19. Amirpour Najafabadi, H.; Khansefid, M.; Mahmoud, G.G.; Zhou, H.; Hickford, J.G.H. Identification of polymorphisms in the oocyte-derived growth differentiation growth factor 9 (GDF9) gene associated with litter size in New Zealand sheep (Ovis aries) breeds. Reprod. Domest. Anim. 2020, 55, 1585–1591. [Google Scholar] [CrossRef]
  20. Li, C.; He, X.; Zhang, Z.; Ren, C.; Chu, M. Pineal gland transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in STH sheep with two FecB genotypes. BMC Genom. Data 2021, 22, 9. [Google Scholar] [CrossRef]
  21. Wu, Z.R.; Yan, L.; Liu, Y.T.; Cao, L.; Guo, Y.H.; Zhang, Y.; Yao, H.; Cai, L.; Shang, H.B.; Rui, W.W.; et al. Inhibition of mTORC1 by lncRNA h19 via disrupting 4E-BP1/raptor interaction in pituitary tumours. Nat. Commun. 2018, 9, 4624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Du, Q.; Hoover, A.R.; Dozmorov, I.; Raj, P.; Khan, S.; Molina, E.; Chang, T.C.; de la Morena, M.T.; Cleaver, O.B.; Mendell, J.T.; et al. MIR205HG is a long noncoding RNA that regulates growth hormone and prolactin production in the anterior pituitary. Dev. Cell 2019, 49, 618–631.e5. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, H.; Ma, J.; Wang, Z.; Yao, X.; Zhao, J.; Zhao, X.; Wang, F.; Zhang, Y. Genome-wide analysis and function prediction of long noncoding RNAs in sheep pituitary gland associated with sexual maturation. Genes 2020, 11, 320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [Green Version]
  26. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. Stringtie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [Green Version]
  27. Sun, L.; Luo, H.; Bu, D.; Zhao, G.; Yu, K.; Zhang, C.; Liu, Y.; Chen, R.; Zhao, Y. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. [Google Scholar] [CrossRef]
  28. Kong, L.; Zhang, Y.; Ye, Z.Q.; Liu, X.Q.; Zhao, S.Q.; Wei, L.; Gao, G. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35, W345–W349. [Google Scholar] [CrossRef]
  29. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  30. Reader, K.L.; Haydon, L.J.; Littlejohn, R.P.; Juengel, J.L.; McNatty, K.P. Booroola BMPR1B mutation alters early follicular development and oocyte ultrastructure in sheep. Reprod. Fertil. Dev. 2012, 24, 353–361. [Google Scholar] [CrossRef]
  31. Goyal, S.; Aggarwal, J.; Dubey, P.K.; Mishra, B.P.; Ghalsasi, P.; Nimbkar, C.; Joshi, B.K.; Kataria, R.S. Expression analysis of genes associated with prolificacy in FecB carrier and noncarrier Indian sheep. Anim. Biotechnol. 2017, 28, 220–227. [Google Scholar] [CrossRef] [PubMed]
  32. McNatty, K.P.; Heath, D.A.; Clark, Z.; Reader, K.; Juengel, J.L.; Pitman, J.L. Ovarian characteristics in sheep with multiple fecundity genes. Reproduction 2017, 153, 233–240. [Google Scholar] [CrossRef]
  33. Tang, J.; Hu, W.; Di, R.; Liu, Q.; Wang, X.; Zhang, X.; Zhang, J.; Chu, M. Expression analysis of the prolific candidate genes, BMPR1B, BMP15, and GDF9 in small tail Han ewes with three fecundity (FecB gene) genotypes. Animals 2018, 8, 166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Cognie, Y.; Benoit, F.; Poulin, N.; Khatir, H.; Driancourt, M.A. Effect of follicle size and of the FecB Booroola gene on oocyte function in sheep. J. Reprod. Fertil. 1998, 112, 379–386. [Google Scholar] [CrossRef]
  35. Hernandez-Ramirez, L.C.; Trivellin, G.; Stratakis, C.A. Cyclic 3′,5′-adenosine monophosphate (cAMP) signaling in the anterior pituitary gland in health and disease. Mol. Cell. Endocrinol. 2018, 463, 72–86. [Google Scholar] [CrossRef] [PubMed]
  36. Kits, K.S.; Mansvelder, H.D. Regulation of exocytosis in neuroendocrine cells: Spatial organization of channels and vesicles, stimulus-secretion coupling, calcium buffers and modulation. Brain Res. Brain Res. Rev. 2000, 33, 78–94. [Google Scholar] [CrossRef]
  37. Formosa, R.; Vassallo, J. cAMP signalling in the normal and tumorigenic pituitary gland. Mol. Cell. Endocrinol. 2014, 392, 37–50. [Google Scholar] [CrossRef] [PubMed]
  38. Naert, G.; Zussy, C.; Tran Van Ba, C.; Chevallier, N.; Tang, Y.P.; Maurice, T.; Givalois, L. Involvement of endogenous brain-derived neurotrophic factor in hypothalamic-pituitary-adrenal axis activity. J. Neuroendocrinol. 2015, 27, 850–860. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, H.; Wei, Y.; Pu, Y.; Jiang, D.; Jiang, X.; Zhang, Y.; Tao, J. Brain-derived neurotrophic factor stimulation of t-type Ca2+ channels in sensory neurons contributes to increased peripheral pain sensitivity. Sci. Signal. 2019, 12, eaaw2300. [Google Scholar] [CrossRef]
  40. Han, D.X.; Sun, X.L.; Fu, Y.; Wang, C.J.; Liu, J.B.; Jiang, H.; Gao, Y.; Chen, C.Z.; Yuan, B.; Zhang, J.B. Identification of long non-coding RNAs in the immature and mature rat anterior pituitary. Sci. Rep. 2017, 7, 17780. [Google Scholar] [CrossRef] [Green Version]
  41. Li, S.; Zhai, J.; Xu, B.; Liu, J.; Chu, W.; Wang, D.; Geng, X.; Chen, Z.J.; Du, Y. Erythropoietin-producing hepatocellular receptor A7 restrains estrogen negative feedback of luteinizing hormone via ephrin a5 in the hypothalamus of female rats. Am. J. Physiol. Endocrinol. Metab. 2020, 319, E81–E90. [Google Scholar] [CrossRef] [PubMed]
  42. Li, S.; Zhai, J.; Liu, J.; Di, F.; Sun, Y.; Li, W.; Chen, Z.J.; Du, Y. Erythropoietin-producing hepatocellular A7 triggering ovulation indicates a potential beneficial role for polycystic ovary syndrome. EBioMedicine 2018, 36, 539–552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gajewska, A.; Herman, A.P.; Wolinska-Witort, E.; Kochman, K.; Zwierzchowski, L. In vivo oestrogenic modulation of Egr1 and Pitx1 gene expression in female rat pituitary gland. J. Mol. Endocrinol. 2014, 53, 355–366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Yuen, T.; Wurmbach, E.; Ebersole, B.J.; Ruf, F.; Pfeffer, R.L.; Sealfon, S.C. Coupling of GnRH concentration and the GnRH receptor-activated gene program. Mol. Endocrinol. 2002, 16, 1145–1153. [Google Scholar] [CrossRef]
  45. Zhang, T.; Wolfe, M.W.; Roberson, M.S. An early growth response protein (Egr) 1 cis-element is required for gonadotropin-releasing hormone-induced mitogen-activated protein kinase phosphatase 2 gene expression. J. Biol. Chem. 2001, 276, 45604–45613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Call, G.B.; Wolfe, M.W. Species differences in GnRH activation of the LHβ promoter: Role of Egr1 and Sp1. Mol. Cell. Endocrinol. 2002, 189, 85–96. [Google Scholar] [CrossRef]
  47. Burger, L.L.; Haisenleder, D.J.; Aylor, K.W.; Marshall, J.C. Regulation of Lhb and Egr1 gene expression by GNRH pulses in rat pituitaries is both c-Jun n-terminal kinase (JNK)- and extracellular signal-regulated kinase (ERK)-dependent. Biol. Reprod. 2009, 81, 1206–1215. [Google Scholar] [CrossRef] [Green Version]
  48. Pierce, J.G.; Parsons, T.F. Glycoprotein hormones: Structure and function. Annu. Rev. Biochem. 1981, 50, 465–495. [Google Scholar] [CrossRef]
  49. Kumar, T.R. Functional analysis of LHβ knockout mice. Mol. Cell. Endocrinol. 2007, 269, 81–84. [Google Scholar] [CrossRef]
  50. Hu, W.; Alvarez-Dominguez, J.R.; Lodish, H.F. Regulation of mammalian cell differentiation by long non-coding RNAs. EMBO Rep. 2012, 13, 971–983. [Google Scholar] [CrossRef] [Green Version]
  51. Guttman, M.; Rinn, J.L. Modular regulatory principles of large non-coding RNAs. Nature 2012, 482, 339–346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Chen, L.L. Linking long noncoding RNA localization and function. Trends Biochem. Sci. 2016, 41, 761–772. [Google Scholar] [CrossRef] [PubMed]
  53. Castro-Oropeza, R.; Melendez-Zajgla, J.; Maldonado, V.; Vazquez-Santillan, K. The emerging role of lncRNAs in the regulation of cancer stem cells. Cell Oncol. 2018, 41, 585–603. [Google Scholar] [CrossRef] [PubMed]
Genes 13 00440 g001 550
Figure 1. LncRNA identification and classification in the pituitary gland. (A) The lncRNAs were identified by two prediction software programs, CNCI and CPC. (B) The number statistics of different types of lncRNAs.
Figure 1. LncRNA identification and classification in the pituitary gland. (A) The lncRNAs were identified by two prediction software programs, CNCI and CPC. (B) The number statistics of different types of lncRNAs.
Genes 13 00440 g001
Genes 13 00440 g002 550
Figure 2. Comparison analysis of lncRNAs and mRNAs in sheep pituitary gland. (A) The boxplot shows the sample expression levels (log10 (FPKM)) of lncRNAs and mRNAs in the H-F and H-L groups. (B) The chromosome distribution of lncRNAs and mRNAs. (C) The length comparison of lncRNAs and mRNAs. (D) Exon content of lncRNAs and mRNAs. (E) Length of the open reading frame (ORF) of lncRNAs and mRNAs.
Figure 2. Comparison analysis of lncRNAs and mRNAs in sheep pituitary gland. (A) The boxplot shows the sample expression levels (log10 (FPKM)) of lncRNAs and mRNAs in the H-F and H-L groups. (B) The chromosome distribution of lncRNAs and mRNAs. (C) The length comparison of lncRNAs and mRNAs. (D) Exon content of lncRNAs and mRNAs. (E) Length of the open reading frame (ORF) of lncRNAs and mRNAs.
Genes 13 00440 g002
Genes 13 00440 g003 550
Figure 3. Comparisons of the number of DE lncRNAs and DE genes (DEGs) in H-F and H-L groups. (A) The hierarchical cluster of DE lncRNAs. (B) Total number of upregulated and downregulated DE lncRNAs in each comparison. (C) Total number of upregulated and downregulated DEGs in each comparison. (D) The hierarchical cluster of DEGs.
Figure 3. Comparisons of the number of DE lncRNAs and DE genes (DEGs) in H-F and H-L groups. (A) The hierarchical cluster of DE lncRNAs. (B) Total number of upregulated and downregulated DE lncRNAs in each comparison. (C) Total number of upregulated and downregulated DEGs in each comparison. (D) The hierarchical cluster of DEGs.
Genes 13 00440 g003
Genes 13 00440 g004 550
Figure 4. The verification of expression level of DE lncRNAs and DE mRNAs in different groups. (A) The relative expression level of four DE mRNAs and lncRNAs in different groups determined by RNA-Seq. (B) The relative expression level of four DE mRNAs and lncRNAs in different groups determined by RT-PCR (* p < 0.05; ** p < 0.01).
Figure 4. The verification of expression level of DE lncRNAs and DE mRNAs in different groups. (A) The relative expression level of four DE mRNAs and lncRNAs in different groups determined by RNA-Seq. (B) The relative expression level of four DE mRNAs and lncRNAs in different groups determined by RT-PCR (* p < 0.05; ** p < 0.01).
Genes 13 00440 g004
Genes 13 00440 g005 550
Figure 5. The top gene ontology (GO) and KEGG enrichment analysis. (A) The top 20 GO enrichment analysis of differential expressed genes (DEGs). (B) The top 20 KEGG enrichment analysis of DEGs.
Figure 5. The top gene ontology (GO) and KEGG enrichment analysis. (A) The top 20 GO enrichment analysis of differential expressed genes (DEGs). (B) The top 20 KEGG enrichment analysis of DEGs.
Genes 13 00440 g005
Genes 13 00440 g006 550
Figure 6. The network of 61 screened DE genes. The network of 61 screened DE genes enriched for Hu sheep pituitary functions and reproduction-related pathways were constructed, the red and green circles represent upregulated and downregulated DE genes, respectively. Node size represents the fold change of a node.
Figure 6. The network of 61 screened DE genes. The network of 61 screened DE genes enriched for Hu sheep pituitary functions and reproduction-related pathways were constructed, the red and green circles represent upregulated and downregulated DE genes, respectively. Node size represents the fold change of a node.
Genes 13 00440 g006
Genes 13 00440 g007 550
Figure 7. 43 DE lncRNAs and their 36 predictably interactive cis- and trans-targeted genes comprised this interactive network. The red and green colors represent up and downregulation, quadrilaterals represent lncRNA, the box represents targeted genes, the straight line and dotted line represent the interaction relationship of trans- and cis-regulation, respectively.
Figure 7. 43 DE lncRNAs and their 36 predictably interactive cis- and trans-targeted genes comprised this interactive network. The red and green colors represent up and downregulation, quadrilaterals represent lncRNA, the box represents targeted genes, the straight line and dotted line represent the interaction relationship of trans- and cis-regulation, respectively.
Genes 13 00440 g007
Genes 13 00440 g008 550
Figure 8. Identification of lncRNA XR_001039544.4. (A) lncRNA XR_001039544.4 localization in Hu sheep pituitary cells as detected by fluorescence in situ hybridization (FISH); scale bar, 20 μm. (B) Expression of lncRNA XR_001039544.4 in the Hu sheep hypothalamus, ovary and pituitary gland (HPO axis). (C,D) Expression of lncRNA XR_001039544.4 and its target gene LHB in different stages of Hu sheep pituitary gland. The data are shown as the mean ± SEM of at least three independent experiments. Statistical significance was analyzed by one-way ANOVA and Student t-test (** p < 0.01).
Figure 8. Identification of lncRNA XR_001039544.4. (A) lncRNA XR_001039544.4 localization in Hu sheep pituitary cells as detected by fluorescence in situ hybridization (FISH); scale bar, 20 μm. (B) Expression of lncRNA XR_001039544.4 in the Hu sheep hypothalamus, ovary and pituitary gland (HPO axis). (C,D) Expression of lncRNA XR_001039544.4 and its target gene LHB in different stages of Hu sheep pituitary gland. The data are shown as the mean ± SEM of at least three independent experiments. Statistical significance was analyzed by one-way ANOVA and Student t-test (** p < 0.01).
Genes 13 00440 g008
Genes 13 00440 g009 550
Figure 9. LncRNA XR_001039544.4 knockdown inhibits the secretion of gonadotropins in pituitary cells. (A,B) mRNA expression of lncRNA XR_001039544.4 and LHB in sheep pituitary cells treated with siRNAs were decreased compared to the negative control (NC), as determined by qPCR. (C) LH secretion was detected by ELISA. (E) mRNA expression of PCNA in lncRNA XR_001039544.4-knockdown pituitary cells. (D,F) Pituitary-cell proliferation was detected using EDU. Scale bars, 100 μm. The data are shown as the mean ± SEM of at least three independent experiments. Statistical significance was analyzed by one-way ANOVA and Student t-test (** p < 0.01).
Figure 9. LncRNA XR_001039544.4 knockdown inhibits the secretion of gonadotropins in pituitary cells. (A,B) mRNA expression of lncRNA XR_001039544.4 and LHB in sheep pituitary cells treated with siRNAs were decreased compared to the negative control (NC), as determined by qPCR. (C) LH secretion was detected by ELISA. (E) mRNA expression of PCNA in lncRNA XR_001039544.4-knockdown pituitary cells. (D,F) Pituitary-cell proliferation was detected using EDU. Scale bars, 100 μm. The data are shown as the mean ± SEM of at least three independent experiments. Statistical significance was analyzed by one-way ANOVA and Student t-test (** p < 0.01).
Genes 13 00440 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wan, Z.; Yang, H.; Cai, Y.; Ma, J.; Cheng, P.; Wang, Z.; Wang, F.; Zhang, Y. Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases. Genes 2022, 13, 440. https://doi.org/10.3390/genes13030440

AMA Style

Wan Z, Yang H, Cai Y, Ma J, Cheng P, Wang Z, Wang F, Zhang Y. Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases. Genes. 2022; 13(3):440. https://doi.org/10.3390/genes13030440

Chicago/Turabian Style

Wan, Zhen, Hua Yang, Yu Cai, Jianyu Ma, Peiyong Cheng, Zhibo Wang, Feng Wang, and Yanli Zhang. 2022. "Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases" Genes 13, no. 3: 440. https://doi.org/10.3390/genes13030440

APA Style

Wan, Z., Yang, H., Cai, Y., Ma, J., Cheng, P., Wang, Z., Wang, F., & Zhang, Y. (2022). Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases. Genes, 13(3), 440. https://doi.org/10.3390/genes13030440

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop