DOI: https://doi.org/10.21203/rs.3.rs-113865/v1
Background
Our study aimed to identify genes related to precocious puberty expressed in the pituitary of different growth stages goats using suppression subtractive hybridization (SSH) screening. Pituitary tissues at 30 day and 90 day and 180 day growth stages of Jining gray goats and Liaoning cashmere goats were used in this study. To identify differentially expressed genes in the pituitary tissues, mRNA from these tissues was extracted and SSH libraries were constructed for screening (API, EPI and BPI groups).
Results
A total of 60, 49 and 58 differently expression genes were annotation in the database. 222 Gene Ontology (GO) terms and 75 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were matched to those genes. Most differentially expressed genes (DEGs) related to the GO terms “structural constituent of ribosome”, “translation” and “GTP binding” were significantly enriched while numerous DEGs related to the Jak-STAT signaling and oocyte meiosis pathways were also significantly enriched. Candidate genes to be involved in precocious puberty and sexual development were retrieved from the SSH libraries. These related genes were discussed taking into account whether they were expressed in different growth stages of pituitary tissues, and of which them influenced the hypothalamic-pituitary-gonadal (HPG) axis.
Conclusions
Interesting findings about precocious puberty related genes from an evolutionary perspective (such as PRLP0, EIF5A and YWHAH) and for putative future goats breeding applications are reported here. We provide a valuable dataset that will facilitate further research into the reproductive biology of goats.
The cost of reserve ewes occupies a large part of the cost of raising sheep. If the age at puberty could be shortened, the production cost of ewes will be greatly reduced. How puberty is initiated is an enigma that still captivates scientists. The onset of puberty is a complex biological process involving numerous factors under the control of the neuroendocrine pathways that are regulated as part of the hypothalamus-pituitary-gonadal (HPG) axis [1]. The HPG axis is a system comprised of endocrine glands whose function is vital to the regulation of reproduction and associated behaviors [2]. Central neurotransmitters, neurohormones, and environmental cues integrate on the HPG-axis and regulate reproduction and puberty onset [3-5]. The pituitary is an endocrine gland that dynamically regulates peripheral tissues to coordinate fundamental physiological functions such as growth, metabolism, sexual maturity and reproduction. The pituitary regulates these homeostatic processes by interpreting hypothalamic signals and, in response, releases hormones from specialized cells in the anterior pituitary. At the onset of puberty, GnRH stimulates the secretion of the pituitary gonadotropins luteinizing hormone (LH) and follicle-stimulating hormone (FSH). LH and FSH in turn stimulate the ovaries to initiate follicular growth and luteal formation that secretes the sex steroid hormones estrogen (E2) and progesterone (P4) [6].
The timing of puberty in humans and other mammals is strongly influenced by genetic regulation [7]. Currently, mutations in the kisspeptin system, MKRN3, and Neurokinin B (NKB) have been identified in sporadic and familial cases of precocious puberty [8]. The loss of function of MKRN3 results in early puberty, implying an inhibitory role of MKRN3 on GnRH secretion [9]. The discovery of the kisspeptin system as a crucial component for pubertal activation of the hypothalamic-pituitary-gonadal axis occurred in 2003, when loss-of-function mutations of the KISS1R (previously known as GPR54) gene were identified in individuals with isolated hypogonadotropic hypogonadism, establishing KISS1R inactivation as a cause of this disorder [10,11]. The KISS1 positive regulation is believed to be important for the LH surge in females of many mammalian species [12]. NKB belongs to a family of closely related peptides called tachykinins. NKB system is necessary for the activation of the hypothalamic-pituitary-gonadal axis in puberty [13].
Suppression subtractive hybridization (SSH) is an effective method for isolation of specific DNA fragments that can be used to differentiate two closely related species [14]. A key feature of this method is simultaneous normalization and subtraction steps that respectively equalize the abundance of DNA fragments within the target population and exclude sequences common to the two populations being compared (10). In this study, the domestics goats breeds Jining Grey goats (JG) and Liaoning Cashmere goats (LC) at different growth stages (30 day, 90 day and 180 day) were selected as the experiment samples, pituitaries were collected and total RNA was extracted. The differential expression genes library was constructed by SSH, and differential expression sexual precocity related genes were screened. The important candidate genes interaction networks were predicted by STING online database. Several genes were identified to further molecular study to reveal the biology mechanism of goat reproduction.
Construction of cDNA libraries
The results of 1.2% agarose gel electrophoresis of double-stranded cDNA products under different cycles is shown in Fig. 1. Treated group have the brightest bands and the widest range in 23 cycles. Control group have the brightest bands in 33 cycles, but the widest range in 23 cycles.
Construction of SSH libraries
324, 295 and 288 white clones were randomly selected from the API, BPI and EPI libraries respectively. Partial result of bacteria liquid PCR is presented in Fig. 2, most of them are positive clones, and the size of the bands is concentrated at 200 bp to 1000 bp, which conforms to the enzyme digestion effect and meets the requirements of the suppression subtractive hybridization libraries.
Sequencing and sequence analysis
In the three libraries, 265, 217 and 225 positive clones were successfully sequenced in the API, EPI and BPI groups, respectively. Excluding the low-quality sequences, the unmatched sequences, the repeat sequences, and the unannotated protein sequences, 91, 58 and 45 valid ESTs were obtained in API, EPI and BPI groups, respectively (Table S1), respectively. These EST sequences were submitted to NR and UNIPROT database to obtain the annotation information of genes (Fig.3). After clustering analysis of Gene ontology on Level 2, data showed that in the API group cDNA library, within the biological process category, 128 ESTs were classified into 11 categories, comprising: metabolic process (25%), cellular process (25.78%), single-organism process (13.28%), biological regulation (9.38%), developmental process (7.03%), multicellular organismal process (5.47%), response to stimulus (5.47%), signaling (3.13%), cellular component organization or biogenesis (3.13%),localization (1.56%) and biological adhesion(0.78%) (Fig.4 A) . Within the cellular component category, they were classified into 11 categories, comprising: cells (20.77%), cell part (20.77%), organelle (19.13%), macromolecular complex (17.49%), organelle part (12.02%), membrane (4.37%), membrane part (2.19%) etc. (Fig.4 B). Within the molecular function category, they were classified into 6 categories, comprising: binding (38.71%), structural molecule activity (37.1%), catalytic activity (14.52%), enzyme regulator activity (4.84%) channel regulator activity (3.23%) and transporter activity (1.61%) (Fig.4 C). In the EPI group subtracted cDNA library, within the biological process category, 88 ESTs were classified into 10 categories, comprising: cellular process (27.27%), metabolic process (26.14%), single-organism process (14.77%), biological regulation (9.09%), cellular component organization or biogenesis (6.82%), response to stimulus (5.68%), developmental process (3.41%) and localization (3.41%) etc. (Fig.5 A). Regarding the cellular component category, they were classified into 10 categories, including cell (21.85%), cell part (21.85%), macromolecular complex (16.81%), organelle (16.81%), organelle part (11.76%), membrane (6.72%) and membrane part (1.68%) etc. (Fig.4 B). In the molecular function category, they were classified into 4 categories, including binding (39.02%), structural molecule activity (34.15%), catalytic activity (24.39%) and channel regulator activity (2.44%) (Fig.5 C). In the BPI group subtracted cDNA library, within the biological process category, 148 ESTs were classified into 12 categories, comprising: cellular process (24.32%), metabolic process (22.97%), single-organism process (13.51%), biological regulation (10.81%), response to stimulus (8.11%), cellular component organization or biogenesis (4.73%), developmental process (4.73%), multicellular organismal process (4.05%) and localization (2.7%) etc. (Fig.6 A). Regarding the cellular component category, they were classified into 8 categories, including cell part (21.02%), cell (20.38%), organelle (18.47%), macromolecular complex (16.56%), organelle part (11.46%), membrane (7.01%) and membrane part (4.46%) etc. (Fig.6 B). In the molecular function category, they were classified into 6 categories, including structural molecule activity (34.92%), binding (26.98%), catalytic activity (23.81%) and enzyme regulator activity (7.94%) etc. (Fig.6 C).
KEGG Pathway analysis
The differentially expressed annotation genes in the three groups were submitted to KEGG Pathway database (Fig.7). The results revealed that most pathways that differentially expressed known genes involved in were ribosome, PI3K-Akt signaling pathway, RNA transport, measles, legionellosis, insulin signaling pathway and Jak-STAT signaling pathway. The most three pathways that the DEGs in API, EPI, BPI groups participated are ribosome, PI3K-Akt signaling pathway, and RNA transport. Some KEGG pathways in the three groups are related with precocious puberty, for instance, estrogen signaling pathway, cAMP signaling pathway, steroid biosynthesis, oocyte meiosis, neuroactive ligand-receptor interaction. In addition, the DEGs in API groups participate in estrogen signaling pathway. Parts of DEGs in EPI are related with oocyte meiosis and some DEGs in BPI group are related with steroid biosynthesis. Some pathways are related with energy metabolism, including PI3K-Akt signaling pathway, insulin signaling pathway, and mTOR signaling pathway (Table 3).
Integrate analysis of Hypothalamic-pituitary-gonadal (HPG) axis genes
Due to the goat database were not provided in the STRING, the high homology species sheep as the alternatives to analysis the correlation of DEGs. Eighteen candidate genes were related to precocious puberty directly or indirectly (Table 4). Most of the DEGs enriched in the one pathway (shown in Fig.8). In our previous study, the DEGs from hypothalamic and ovary tissues were identified by SSH libraries. Integrate all of the data from the three tissues of hypothalamic-pituitary-gonadal (HPG) axis, the result showed that almost 70% DEGs were enriched in the one pathway (Fig.9).
In this study, SSH technology was used to construct the differential gene libraries of pituitary tissues in JG and LC goats at different growth stages. Through sequencing, sequence alignment and statistics, the SSH library of pituitary contains 707 ESTs, 37.48% of which were expressed in API group and 30.69% were expressed in EPI group and 31.82% were expressed in BPI group. At the same time, the GO annotation of these libraries showed that most of the genes were mainly involved in biological process, cellular component or molecular function such as metabolic process, cellular process, catalytic activity, binding, etc., and some genes were also involved in the developmental process. This indicates that the molecular mechanism of different growth stages affecting the precocious puberty of goats is comprehensive and complex, which might regulate the growth, development, metabolism and development of goats through multiple signaling pathways.
According to the GO annotation of SSH library, 18 precocious puberty-related differential expression genes were screened from three groups, 44.44% of them in API group, 11.11% in EPI group and 44.44% in BPI group. The results indicated that most precocious puberty related genes were shown in the same growth stages of different goat breeding. The influence of growth stages is more important for the goat precocious puberty.
Some of these precocious puberty related differential genes were differentially expressed in the three groups. Among them, eighteen genes were appeared more than twice in three groups including RPS25, GNB2L1,RPS12, GH1,PET100,PLP1,RPL10A,RPL37,RPLP0,RPS19,PKN1,RPL23A,RPS10, RPS14,RPS20.S,RPS7, HSPA8 and SLC25A5. RPLP0 protein (ribosomal protein lateral stalk subunit P0) encodes one out of approximately 80 ribosomal proteins in human, which are involved in protein synthesis and apoptosis processes [15]. In a previous study, Ragni et al. classified the RPLP0 as the most stable reference gene in expression assays performed in mesenchymal stem cell differentiation, osteoblasts precursor cells [16]. In the study interactions of cytotoxic T lymphocytes with tumor cells, RPLP0 is stably expressed in melanoma cells [17]. RPLP0 gene as reference genes are stably expressed in GCs in both controls and polycystic ovarian syndrome patients and can be used for normalization in gene expression profiling by qRT-PCR [18]. Nikishin et al. confirmed the conclusion of RPLP0 gene as the most stable reference genes in the case of vitrified/thawed human ovarian cortical tissue [19].
PKN1 (protein kinase N1) is a stress-responsive kinase and a member of the protein kinase novel (PKN) family also known as protein kinase C-related kinases (PRKs) [20, 21]. The interaction of Rho with PKN1 induces a conformational change in PKN1 leading to activation loop phosphorylation by 3-Phosphoinositide-Dependent Kinase-1(PDK1) on Thr774 in the case of human PKN1 (Thr778 in mouse/rat) which is necessary for the catalytic activation of PKN1/214 and critical for the stability of the protein [22,23]. PKN1 is shown to provide a basal cardioprotective function in the face of I/R injury [24]. In the currently study found that PKN1 as a novel key player in fine-tuning the balance between axonal outgrowth and presynaptic differentiation in the parallel fiber–forming (PF-forming) cerebellar granule cells (Cgcs) [25]. In neurons PKN1 is the most abundant isoform and has been implicated in a variety of functions including cytoskeletal organization and neuronal differentiation [26-29]. In the precocious puberty related studies found that PKN1 gene modulated TGFβ and EGF dependent regulation of the responses of HEC-A-1 endometrial cancer cells proliferation, migration, and invasiveness, and therefore is a component of the network signaling downstream of TGFβ and EGF [30]. Attarha et al. used proteomics, systems biology, and immunohistochemistry to explore protein expression in human endometrial tumours, identified PKN1 may be considered as predictive biomarkers of endometrial cancer [31].
SLC25A5 (solute carrier family 25 member 5) facilitated the transport of molecules involved in the urea and citric acid cycles, oxidative phosphorylation, DNA maintenance and iron metabolism processes [32, 33]. In the dysregulated expression of androgen metabolism genes and genetic analysis in hypospadias, the data provided direct evidence that SLC25A5 may contribute the genetic etiology of hypospadias [34]. In the bovine reproduction study, SLC25A5 gene were found high expressed in the SSH library of bovine blastocyst [35]. RPL37 (ribosomal protein L37), as a component of the 60S subunit of ribosomes and belongs to the L37E family of ribosomal proteins. In the Zebra Finch Song System study, RPL37 might all influence sexual differentiation, perhaps with the hormone and proteins interacting, such that an appropriate balance is required for normal development [36]. PRL37 protein as a histotype-specific prognostic biomarker was found in the Ovarian Carcinomas Using Immunohistochemistry [37]. These genes above mentioned were potential molecular markers for the study of goat precocious puberty.
In our experiment, 18 DEGs were related to precocious puberty directly or indirectly. Among them, RPL10A, RPL23A, RPL24, RPL34, RPL37, RPLP0, RPS3A , RPS6 and RPS7 are all common protein family ribosomal protein, they play an important role in the development of ovarian and breast cancers related to precocious puberty [38,39]. And EIF5A, HNRNPDL, HSPA8, PKN1, PRKAR2B, SLC40A1, SLC25A5, SNRPD2 and YWHAH, play key roles in oocyte meiosis, steroid biosynthesis and conduction of reproduction process [40–42].
Precocious puberty is regulated via hormones of the hypothalamic-pituitary-gonadal (HPG) axis. In our study, we analysis of the hypothalamic-pituitary-gonadal (HPG) axis related genes by the online tool STRING. We found that most of the DEGs from the three tissues enriched in one pathway, the EIF5A and YWHAH gene play an important role in the network (Fig.9).
eIF5A (eukaryotic translation initiation factor 5A) is a highly conserved 17-kDa protein that is expressed ubiquitously in all cells. eIF5A participates in diverse cellular processes, including protein translation [43-45], nucleocytoplasmic transport of RNA [46], cell proliferation [47], inflammation, and apoptosis[48,49]. Stimulation of follicular growth by follicle-stimulating hormone (FSH) is associated with increased expression of luteinizing hormone/choriogonadotropin receptors (LHCGRs) in granulosa cells [50, 51]. EIF5A was identified as one of the proteins that interact with LRBP and lead to the degradation of LHCGR mRNA by facilitating the transport of LHCGR mRNA-LRBP complex to P bodies for degradation [52-54]. In the currently, a series studies were carried out to determine the mechanism by which inhibition of hypusination of EIF5A causes an increase in LHCGR mRNA expression [55]. The EIF5A gene usually expressed in ovary tissue, in our study the expression of EIF5A gene were detected in the pituitary tissue. The result indicated that this gene might as an important gene expressed in the hypothalamic-pituitary-gonadal (HPG) axis to regulate the goats precocious puberty.
YWHAH (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein eta) is an isoform of YWHA protein. Numerous studies implicate YWHA as a critical regulator of the cell cycle in meiotic and mitotic cells as well as other cellular processes [56, 57]. The results indicate that oocyte-specific or global elimination of YWHAH protein does not result in abnormal fertility, oocyte maturation or development. Breeding and development of pups was normal in the absence of YWHAH or YWHAE in females with oocyte-specific knockout of these genes. Global inactivation of Ywhah in female mice does not appear to alter oogenesis, oocyte maturation and early development [58]. The previously noted that YWHAH may play an important role in meiotic spindle formation by employing antisense morpholino knockdown approaches [59].
The effect of differential expressed genes from pituitary tissues at different growth stages may be related to the process of precocious puberty. However, the exact pathway remains unknown, and there are few reports on how these differential genes affect the precocious puberty of goats. In the future, more function analysis will be needed to further explore the effects of these genes on precocious puberty in goats.
In summary, the expression of many genes in pituitary tissues from JG and LC goats at different growth stages significantly affected precocious puberty, including 18 related genes. Many genes related with oocyte meiosis and steroid biosynthesis may have significant roles in the puberty onset, of which, PRLP0, EIF5A and YWHAH play a regulatory role in the precocious puberty of goats, thus will be our candidates for further research.
Ethical statement
All experimental procedures mentioned in the present study were approved by the Science Research Department (in charge of animal welfare issue) of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (IAS-CAAS) (Beijing, China). Ethical approval was provided by the animal ethics committee of IASCAAS (No. IAS2019-63).
Animals and tissue collection
All the Jining Gray (JG) goats (Jining Grey goats conservation base, Jiaxiang County, Shandong Province, People’s Republic of China) and Liaoning Cashmere (LC) goats (Liaoning Cashmere Goat Breeding Center, Liaoyang City, Liaoning Province, People’s Republic of China)displayed in Table 1 in this study were housed in open sheepfolds and under the same nutrition condition. The pituitary tissues were collected from each goat, sacrificed after anesthesia with 3% pentobarbital sodium salt injection (20mg/kg body weight) (Merck, Darmstadt, Germany), and preserved in RNA later RNA Stabilization reagent (Qiagen, Hilden, Germany) and kept at -20°C until RNA isolation.
poly A+ RNA purification and concentration
According to the instruction of Dynabeads mRNA DIRECT™ Kit (Invitrogen, Inc. USA), the poly A+ RNA was isolated from pituitary. After mixed in equal quantity from three goats in each group, the poly A+RNA concentrated using RNA clean & concentration-5 mRNA (Zymo Research, Orange, CA, USA) to a suitable concentration.
Suppression subtractive hybridization (SSH)
The concentrated poly A+ RNA from the pituitaries of JG goats was tester and that of LC goats was driver. As shown in Table 2, the poly A+ RNA from the pituitaries of 30-JG was hybridized with that of 30-LC goats and this hybridization was named as A group. The hybridizations between 90-JG and 180-LC, 90-JG and 90-LC were names as B group and E group, respectively.
The regents for SSH were from PCR-Select cDNA Subtraction Kit (Clontech, palo alto, CA) and the procedure were carried out according to the protocols. Simply, the double strand cDNA (dscDNA) were synthesized from 2 μg poly A+ RNA of tester and driver respectively and then were digested by Rsa Ⅰ restriction endonuclease to obtain shorter, blunt-ends dscDNA fragments, required for adaptor ligation and optimal for subtraction. After analysis of Rsa Ⅰ digestion, the digested blunt-ends of tester cDNA, but not adaptor cDNA, were divided into two parts and ligated with two different cDNA adaptors (adaptor 1 5′-CTAATACGACTCACTATAGGGCTCGAGCGGCCGCCCGGGCAGGT-3′ and adaptor 2R 5′-CTAATACGACTCACTATAGGGCAGCGTGGTCGCGGCCGAGGT-3′), respectively. After ligation, the ligation efficiency analysis was performed by PCR experiment with GAPDH primers (GAPDHF: 5′- AGGCTGGGGCTCACTTGAAG-3′, GAPDHR: ATGGCGTGGACAGTGGTCAT-3′) (goat GAPDH mRNA sequence-GenBank: AJ431207) and PCR primer I (5′-CTAATACGACTCACTATAGGGC-3′) using the Advantage cDNA PCR kit (Clontech, palo alto, CA).
Two cycles of hybridization were followed after the ligation. In the first hybridization, 1.5 μL Rsa I-Digested Driver cDNA and 1.5μl Adaptor-Ligated Tester were hybridized in 1.0 μL 4×Hybridization buffer at 68°C for 8 hours. For the second hybridization process, 1.0 μL fresh denatured driver cDNA, 1.0 μL 4×Hybridization Buffer and 2.0 μL ddH2O were then mixed with the two samples from the first hybridization simultaneously and incubated at 68°C overnight.
The final hybridization solution (also called the subtracted library) was employed as a template to amplify the differentially expressed sequences in the tester population by using a set of PCR primer1 and was followed by nested PCR primers (Nested PCR primer 1 5′-TCGAGCGGCCGCCCGGGCAGGT-3′, Nested PCR primer2R 5′-AGCGTGGTCGCGGCCGAGGT-3′), which does not exponentially amplify the non-adaptor (derived from driver cDNA), cDNA with the one adaptor on either end (derived from tester cDNA hybridized with driver cDNA), or cDNA with the same adaptor on both ends (derived from relatively abundant tester cDNA). For the primary PCR, 1 μL sample was added to 24 μL PCR master mix prepared using the reagents supplied in the kit, and cycling conditions commenced as follows: 75 °C for 5 min to extend the adaptors; 94 °C for 25 sec; and 27 cycles at 94 °C for 10 sec, 66 °C for 30 sec, and 72 °C for 1.5 min. Amplified products were diluted 10-fold in sterile water and 1 μL of diluted primary PCR products were added to 24 μL of secondary PCR master mix containing nested primers, 1 and 2R, to ensure specific amplification of double-stranded templates containing both adaptors. Secondary PCR was performed 12 cycles at 94 °C for 10 sec, 68 °C for 30 sec, and 72 °C for 1.5 min. Primary and secondary PCR products were analyzed on a 2% agarose gel. And now the PCR mixture is enriched for differentially expressed cDNA in tester sample.
Cloning, sequencing and sequence analysis
3 μL of purified PCR products was ligated into pGEM-T Easy Vector (Promega, Madison, WI, USA) and transformed into the competent E. coli DH5α. Positive mono clones were detected by PCR and sequenced by Invitrogen Corporation.
The vector and adaptor sequence of sequences obtained by sequencing were remove using the UniVec database (http://www.ncbi.nlm.nih.gov/VecScreen/) to get clean sequence.
The differentially expressed ESTs were classified into known genes, known ESTs and unknown ESTs according to homologies analyzed by BLAST (www.ncbi.nlm.nih.gov/blast/blast.cgi). The differentially expressed genes (DEGs) function predicted in the Cluster of Orthologous Groups of proteins database (COG) (http://www.ncbi.nlm.nih.gov/COG/) and were classified according to the COG database. Then the DEGs were submitted to the KEGG pathway database (http://www.genome.jp/kegg/pathway.html) to analysis which pathway they participate in. The protein-protein interactions were predicted in the STRING (Search Tool for Retrieval of Interacting Genes/Proteins) net (http://string.embl.de/).
SSH: suppression subtractive hybridization
GO: Gene Ontology
KEGG: Kyoto Encyclopedia of Genes and Genomes
DEGs: differentially expressed genes
HPG: hypothalamic-pituitary-gonadal
LH: luteinizing hormone
FSH: follicle-stimulating hormone
E2: estrogen
P4: progesterone
NKB: Neurokinin B
JG: Jining Grey goats
LC: Liaoning Cashmere goats
COG: Cluster of Orthologous Groups of proteins database
PKN: protein kinase novel
PDK1: 3-Phosphoinositide-Dependent Kinase-1
Cgcs: cerebellar granule cells
Acknowledgement
Not Applicable.
Funding
This research was funded by the National Natural Science Foundation of China (30871773, 31201778), the Earmarked Fund for China Agriculture Research System (CARS-38), the Agricultural Science and Technology Innovation Program of China (ASTIP-IAS13), Shandong Provincial Natural Science Foundation (ZR2017BC054).
Author information
Affiliation
Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, China Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193, China
Yufang Liu&Mingxing Chu
College of Life Sciences and Food Engineering ,Hebei University of Engineering, Handan 056021, China
Yufang Liu
College of Agriculture, Liaocheng University, Liaocheng 252059, China
Guiling Cao& Yujing Xie
Contributions
Conceived and designed the experiments: GLC, MXC; Performed the experiments: GLC, YJX, YFL; Analyzed the data: GLC, YJX, YFL; Contributed reagents: MXC; Acquisition of data: GLC, YJX, YFL; YFL and MXC drafted the manuscript. YFL, GLC and MXC contributed in writing the paper; YFL, GLC and MXC substantively revised the paper; Study initiation: GLC, YJX, YFL and MXC; Read and approved the final version of the paper: all co-authors.
Corresponding author
Correspondence to Mingxing Chu.
Ethics declarations
Ethics approval and consent to participate
All experimental procedures mentioned in the present study were approved by the Science Research Department (in charge of animal welfare issue) of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (IAS-CAAS) (Beijing, China). Ethical approval was provided by the animal ethics committee of IASCAAS (No. IAS2019-63).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and analyzed during the current study available from the corresponding author on reasonable request.
Table 1 The detail information of experiment samples
Goat breed |
Age |
Abbreviation |
Number |
Development stage |
Jining Grey goats |
30-day-old |
30-JG |
3 |
Juvenile |
90-day-old |
90-JG |
3 |
Puberty |
|
Liaoning Cashmere goats |
30-day-old |
30-LC |
3 |
Juvenile |
90-day-old |
90-LC |
3 |
Juvenile |
|
180-day-old |
180-LC |
3 |
Puberty |
Table 2 Suppression subtractive hybridization groups. JG: Jining Grey goats; LC: Liaoning Cashmere goats.
Hybridization Groups |
Tester |
Driver |
A |
30-JG |
30-LC |
B |
90-JG |
180-LC |
E |
90-JG |
90-LC |
Table 3 The genes enriched in different KEGG function in the three SSH groups
KOG function |
AO |
EO |
BO |
PI3K-Akt signaling pathway |
GH1 |
LOC102322233 GH1 |
LOC102322233 bPKN1 RPS6 |
Jak-STAT signaling pathway |
GH1 |
LOC102322233 GH1 |
LOC102322233 |
MAPK signaling pathway |
HSPA8 |
|
|
Wnt signaling pathway |
LOC102250913 LOC102324853 |
|
|
Estrogen signaling pathway |
HSPA8 |
|
|
cGMP-PKG signaling pathway |
LOC102250913 |
SLC25A5 |
Slc25a5 |
mTOR signaling pathway |
LOC100073610 |
|
RPS6 |
TGF-beta signaling pathway |
LOC102250913 |
|
|
Insulin signaling pathway |
|
|
PRKAR2B EXOC7 RPS6 |
Steroid biosynthesis |
|
|
SC5D |
Oocyte meiosis |
|
YWHAH |
|
Table 4 Precocious puberty related important DEGs among three groups
Gene name |
Description |
Reference |
EIF5A (eukaryotic translation initiation factor 5A) |
mRNA-binding protein involved in translation elongation |
[54-55,60] |
HNRNPDL(heterogeneous nuclear ribonucleoprotein D like) |
Involved in regulation of gene expression |
[61] |
HSPA8(heat shock protein family A (Hsp70) member 8) |
A chaperone, and binds to nascent polypeptides to facilitate correct folding |
[42] |
PKN1 (protein kinase N1) |
mediate insulin signals to the actin cytoskeleton |
[30-31, 62-64] |
PRKAR2B (protein kinase cAMP-dependent type II regulatory subunit beta) |
important for a variety of cellular functions |
[65,66] |
RPL10A (ribosomal protein L10a) |
the organelles that catalyze protein synthesis |
[67] |
RPL23A (ribosomal protein L23a) |
one of the target molecules involved in mediating growth inhibition by interferon |
[38] |
RPL24 (ribosomal protein L24) |
The protein belongs to the L24E family of ribosomal proteins. |
[68,69] |
RPL34 (ribosomal protein L34) |
Overexpression of this gene has been observed in some cancer cells. |
[70] |
RPL37 (ribosomal protein L37) |
The protein contains a C2C2-type zinc finger-like motif. |
[37] |
RPLP0 (ribosomal protein lateral stalk subunit P0) |
It is a neutral phosphoprotein |
[18-19,71-74] |
RPS3A (ribosomal protein S3A) |
The protein belongs to the S3AE family of ribosomal proteins. |
[75,76] |
RPS6 (ribosomal protein S6) |
It is the major substrate of protein kinases in the ribosome |
[77] |
RPS7 (ribosomal protein S7) |
Exhibits mRNA binding activity |
[78] |
Slc25a5 (solute carrier family 25 member 5) |
Suppressed expression of this gene has been shown to induce apoptosis and inhibit tumor growth. |
[35] |
SLC40A1 (solute carrier family 40 member 1) |
The product of this gene functions as a gated pore that translocates ADP from the cytoplasm into the mitochondrial matrix. |
[41] |
SNRPD2 (small nuclear ribonucleoprotein D2 polypeptide) |
It is required for pre-mRNA splicing and small nuclear ribonucleoprotein biogenesis. |
[79] |
YWHAH (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein eta) |
This gene contains a 7 bp repeat sequence in its 5' UTR, and changes in the number of this repeat have been associated with early-onset schizophrenia and psychotic bipolar disorder. |
[56-59] |