Screening of Puberty Related Differentially Expressed Genes from Ovary Tissues of Jining Grey Goat Based on Suppression Subtractive Hybridization

Yufang Liu1,2, Guiling Cao3, Yujing Xie3 and Mingxing Chu1*

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

2College of Life Sciences and Food Engineering, Hebei University of Engineering, Handan 056021, China.

3College of Agriculture, Liaocheng University, Liaocheng 252059, China.

Yufang Liu and Guiling Cao made equal contributors to this article.

ABSTRACT

To explore the molecular mechanism of the puberty effects on the reproduction of goats, goats at different growth stages from two breeds (Jining Grey (JG) goats and Liaoning Cashmere (LC) goats) were divided into three groups with three replicates (juvenile stage (JG 30d vs. LC 30d, AO group), puberty (JG 90d vs. LC 180d, BO group) and the same age control of puberty (JG 90d vs. LC 90d, EO group). Ovary tissues were taken from goats on 30 days, 90 days and 180 days, total RNA was extracted, differentially expressed (DE) gene libraries were constructed by suppression subtractive hybridization (SSH) technology, and puberty-related differential genes were screened. A total of 184 differentially expressed genes were screened, including 57, 68 and 59 in the AO, BO and EO groups, respectively. There are 22 differentially expressed genes were directly associated to the ovary development. According to the gene function analysis in KOG database, the known genes were clustered into 10 subdivisions in A group, 8 in E group and 9 in B group under three categories, mainly in posttranslational modification, protein turnover, chaperones, energy production and conversion, and transcription process. Pathway analysis in KEGG pathway database of the known genes revealed that the five pathways that most differentially expressed genes involving: GnRH signaling pathway, oxytocin signaling pathway, melanogenesis, thyroid hormone synthesis, and insulin secretion. Online tool was used to predict the function of DEGs, the results showed that the CYP11A1 gene might regulate the GDF9 and other prolificacy-associated markers in ovary organs of goats, thereby affecting the development of puberty. The expression of various genes and the transduction of hormone secretion signals, improving the reproduction of goats.

Article Information

Received 20 July 2021

Revised 05 September 2021

Accepted 21 September 2021

Available online 10 June 2022

(early access)

Published 05 June 2023

Authors’ Contribution

Conceptualization, MXC. Data analysis, YJX and GLC. Investigation, GLC, YJX, YFL and MXC. Resources, GLC and YFL. Writing original draft preparation, YFL and GLC. Writing review and editing, YFL, GLC and MXC. Supervision, MXC. All authors have read and agreed to the published version of the manuscript.

Key words

Puberty, Ovary, Goat, Suppression subtractive hybridization, Differentially expressed genes

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

* Corresponding author: [email protected]

0030-9923/2023/0004-1637 $ 9.00/0

Copyright 2023 by the authors. Licensee Zoological Society of Pakistan.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

INTRODUCTION

Puberty is the process by which animals maturate into an adult capable of sexual reproduction, characterized by the maturation of gametogenesis, secretion of gonadal hormones. In human, onset of puberty occurs after reactivation of the hypothalamic GnRH secretory system (De Sanctis et al., 2019). At puberty, the pulsatile GnRH secretion, and the subsequent episodic pituitary gonadotropin secretion, necessary for normal gonadal development and function, is triggered by the activation of the GnRH pulse generator (Dulka et al., 2020). Where, the GnRH pulsatile secretion is dependent on the coordinated action of the scattered GnRH neurons, controlled by stimulatory and inhibitory inputs, such as neurokinin B, kisspeptins and gamma aminobutyric acid (Plant, 2015). In sheep, the pubertal onset is marked by a decrease in estradiol-negative feedback, resulting in an increased frequency of GnRH, and subsequently LH, pulses and the increased estradiol production that in turn induces the subsequent GnRH/LH surge and first ovulation (Foster and Hileman, 2015). The immature female lamb is capable of generating requisite high-frequency GnRH/LH pulses, but this pattern is not produced because sensitivity to estradiol feedback inhibition is high (Moenter et al., 1991). Removal of steroid-negative feedback by ovariectomy (OVX) in the pre-pubertal ewe lamb results in an elevation of GnRH/LH secretion, which can readily be reduced with subcutaneous administration of E2 (OVX+E lambs). This E2-induced suppression reflects an inhibition of GnRH/LH pulse frequency that persists until the time when pubertal onset would normally occur (Foster and Hileman, 2015).

In summary, the mechanism of pubertal onset is complex and thought to be associated with environmental factors, neuroendocrine factors, genetic factors and their interactions (Smith, 2001). Understanding how puberty is timed is critical for developing practical strategies to synchronize lambing with market demand and for potentially decreasing the generation interval in farm animals, thereby increasing the overall production of food and fiber. The Jining Gray (JG) goat is a sexual precocious breed that reaches puberty at 60-90 days and the Liaoning Cashmere (LC) goat is a late-puberty breed that reaches puberty at 150-210 days. Our previous study showed that Jining Grey (JG) goats and Liaoning Cashmere (LC) goats provide the best materials to identify factors that regulate the maturation of the hypothalamic-pituitary-gonadal (HPG) axis (Cao et al., 2015). In this study, we selected the domestics JG and LC goats at different growth stages (30, 90 and 180 days) as the experimental animals. Suppression subtractive hybridization (SSH) is an effective method for isolation of specific DNA fragments that can be used to differentiate two closely related species (Diatchenko et al., 1996; Rebrikov et al., 2004). 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 (Gurskaya et al., 1996). We have therefore used the SSH method to isolate puberty-related genes in goat ovaries. These results increase our understanding of the physiological processes involved in goat reproduction. In particular, the differently expression genes in the puberty may play a role in the modulation of the ovary development.

MATERIALS AND METHODS

Ethics approval

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. IASCAAS-AE-03).

Animals and tissue collection

All the Jining Grey (JG) goats and Liaoning Cashmere (LC) goats displayed in Table I in this study were housed in open sheepfolds and under the same nutrition condition. The ovaries 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.

 

Table I. The detail information of experiment samples.

Goat breed	Age	Abbreviation	Number	Development stage
Jining gray 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

 

RNA extraction and concentration of cDNA libraries

Total RNA of ovaries was isolated using Trizol Reagent (Invitrogen, Inc. USA) according to the manufacturer’s instructions and the quantity and quality were determined by electrophoresis and spectrometry. The equal quantity total RNA from three goats in each group were mixed together and used for poly A+ RNA purification using the Oligotex mRNA Midi Kit (Qiagen, Hilden, Germany). After that, the poly A+ RNA was concentrated using RNA clean and concentration-5 mRNA (Zymo Research, Orange, CA, USA) to a suitable concentration.

Construction of SSH libraries

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) was synthesized from 2 μg poly A+ RNA of tester and driver respectively and then were digested by Rsa I restriction endonuclease to obtain shorter, blunt-ends dscDNA fragments which were required for adaptor ligation and optimal for subtraction. After analysis of Rsa I digestion, the digested blunt-ends of tester cDNA were divided into two parts and ligated with two different cDNA adaptors (adaptor 1 ATAGGGCTCGAGCGGCCGCCCGGGCAGGT-3′ and adaptor 2R 5′- CTAATACGACTCACTATAGGGCAGCGTGGTCGCGGCCGAGGT- 3′), respectively and adaptor will not be ligated to the driver dscDNA. 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 bioinformatics analysis

3 μL of PCR products purified using DNA clean and concentration-5 (Zymo Research, Orange, CA, USA) was ligated into pGEM-T Easy Vector (Promega, Madison, WI, USA) and transformed into the competent E. coli DH5α. Most white clones were picked up and cultured for PCR and sequencing. 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/).

RESULTS

Construction of cDNA libraries

The results of 1.2% agarose gel electrophoresis of part of double-stranded cDNA products under different cycles are shown in Figure 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. As can be seen from the figure, the optimal cycles for both treated group and control group was 23.

 

Table II. The SSH groups arrangement.

Hybridization Groups	Tester	Driver
A	30-JG	30-LC
B	90-JG	180-LC
E	90-JG	90-LC

 

Table III. The numbers of differentially expressed ESTs in ovary from A, E and B group.

	A group	E group	B group
Known gene	62 (57.4%)	74 (69.8%)	127 (84.1%)
Known EST	45 (41.7%)	31 (29.2%)	23(15.2%)
Unknown EST	1(0.9%)	1(1.0%)	1(0.07%)
Total	108	106	151

 

*The number in the bracket represents the percent.

 

Homologies analysis of differentially expressed ESTs

108, 106 and 151 differential expressed ESTs were obtained in A group, E group and B group, respectively after the positive clones were sequenced. Total 365 EST were submitted to NR, KEGG and UNIPROT database for sequence alignment. The homology with goat, sheep, cattle, human, pig, mouse and rat were recorded and the ESTs were classified. The ESTs were classified into known gene, known EST and unknown EST, as shown in Table III, and the percentage of known gene of juvenile group is the least. In group A, the 62 ESTs of known genes represent 57 transcripts of genes, the 74 ESTs of known genes represent 59 transcripts of genes in group E, and the 127 ESTs of known genes represent 68 transcripts of genes in B group.

GO and KOG classification

Upon GO functional enrichment, 124, 68, and 55 specific GO terms in biological process, molecular function and cellular component were identified in the three groups, respectively (Fig. 2). Most DEGs in the three groups are enriched in cell, organelle, and cell part in cellular component and metabolic process and cellular process in biological process. In the cellular component classification, most DEGs are enriched in cell, cell part and organelle. In molecular function classification, most DEGs are enriched in binding, structural molecule activity and catalytic activity and metabolic process, cellular process, single-organism process are the top three subdivisions with DEGs in the biological process. The known genes in the three groups were submitted into KOG database. According to the gene function analysis in the database, the known genes were grouped into 13 divisions in A group, 15 in E group and 14 in B group (Fig. 3).

 

 

Table IV. The genes enriched in different KOG function in the three SSH groups.

KOG function	AO	EO	BO
Energy production and conversion	ND4, ATP6, CYTB	ATP5H, NDUFA9, ND4, COX1, ND4L, ATP6, ATP5A1	ND4, COX2, COX1, ND2, CYTB, ATP6, LDHB, ATP5B
RNA processing and modification	RPL4, DROSHA, NOP58	FUBP1, RPL4, SRSF11	RPL4
Past-transcription modification
Amino acid transport and metabolism	IDH3B, SHMT2
Nucleotide transport and metabolism	IMPDH2
Lipid transport and metabolism	TECR, PLA2G7
Transcription process	EEF1B2, IKBKAP, ASCC1, NCOR1
Cytoskeleton	ACTG2, PNN, ARPC5L		ACTB, MYLK, ARPC2,RPS15A
Cell motility		Thymosin beta-b	Thymosin beta-b
Inorganic ion transport		FTH1	ATP1B3, FTH1
Extracellular structures		SPARC	COL3A1, LAMA4
Signal transduction	GNB2L1, FBLN1	GNB2L1, RGS2, LTBP1, PDLIM3	GNAS, GNB2L1
Posttranslational modification, protein turnover, chaperones	OSGEP, POMP, PSMD6, HSP90AA1	YWHAB, C3, YWHAQ, HSP90B1, GDI1, HSC71, HSP90AA1, PPIB	GSTM2, PPIB, RPN1, HSPA8, C3
Intracellular trafficking, secretion, and vesicular transport	TEMD10	CLTA, TM9SF3, COPA, TIMM8B, KDELR3	ARL4A

 

KEGG pathway analysis

The differentially expressed known genes in the three groups were submitted to KEGG Pathway database (Fig. 4 and Table V). The results revealed that most pathways that differentially expressed known genes involved in were metabolic pathways, Parkinson’s disease, oxidative phosphorylation, Huntington’s disease, Alzheimer’s disease and GnRH signaling pathway. The most three pathways that the DEGs in AO, EO, BO groups participated are ribosome, oxidative phosphorylation, and Parkinson’s disease. Some KEGG pathways in the three groups are related with reproduction, for instance, ovarian steroidogenesis, estrogen signaling pathway, oocyte meiosis, GnRH signaling pathway, melanogenesis. In addition, the DEGs in EO and BO groups participate in progesterone-mediated oocyte maturation. Parts of DEGs in EO are related with oocyte meiosis and some DEGs in BO group are related with oxytocin signaling pathway and melanogenesis. Some pathways are related with energy metabolism, including fatty acid metabolism, glycolysis/gluconeogenesis, PI3K-Akt signaling pathway, insulin signaling pathway, and mTOR signaling pathway.

 

Table V. Parts of KEGG pathways in the three SSH groups.

AO	EO	BO
Ovarian steroidogenesis		Ovarian steroidogenesis
Estrogen signaling pathway	Estrogen signaling pathway	Estrogen signaling pathway
Steroid hormone biosynthesis	Oocyte meiosis	Oxytocin signaling pathway
	Progesterone-mediated oocyte maturation	Progesterone-mediated oocyte maturation
		GnRH signaling pathway
		Melanogenesis
		Circadian entrainment
Fatty acid metabolism	Galactose metabolism	Glycolysis / Gluconeogenesis
Biosynthesis of unsaturated fatty acids	Amino sugar and nucleotide sugar metabolism	Carbohydrate digestion and absorption
Fatty acid elongation	Glycolysis / Gluconeogenesis	N-Glycan biosynthesis
	Starch and sucrose metabolism
	Pentose phosphate pathway
		Insulin secretion
		Insulin signaling pathway
MAPK signaling pathway		mTOR signaling pathway
PI3K-Akt signaling pathway	PI3K-Akt signaling pathway	PI3K-Akt signaling pathway
	Glutathione metabolism	Glutathione metabolism

 

 

The pathways that DEGs in AO group participated are related with fat metabolism, such as fatty acid metabolism and biosynthesis of unsaturated fatty acids, while, that in EO and BO groups are related with glucose metabolism, for instance, glycolysis/gluconeogenesis, carbohydrate digestion and absorption, pentose phosphate pathway. Carbohydrate metabolism provide much energy to body and energy obtained from metabolism is usually stored temporarily within cells in the form of ATP.

Protein-protein interaction analysis

The protein interaction of DEGs in A, E and B group analyzed in STRING database were showed in Figure 5. The most obvious interaction network in the three groups is the structure of ribosome/protein translation, especially in the BO group. more proteins are need in the ovary of JG goats than that in the ovary of LC goats, even all are pubertal goats. The numbers of DEGs related with the in oxidative phosphorylation EO and BO group are more than that in AO group. The oxidative phosphorylation is the metabolic pathway in which the mitochondria in cells use their structure, enzymes, and energy released by the oxidation of nutrients to reform ATP which is the molecule that supplies energy to metabolism. This pathway is a highly efficient way of releasing energy. So, these may indicate that more genes related with energy production and storage is expressed in JG goat than that in LC goat in these stages.

Twenty-two candidate genes were related to ovary development directly or indirectly (Table VI). The STRING database was used to predict the relationship between candidate genes and GDF9 and other prolificacy-associated markers in goat, cattle, mouse and human. The results showed that CYP11A1 gene was interacted with GDF9 and BMP15 genes whatever in the four species (Fig. 6).

 

DISCUSSION

In this study, SSH technology was used to construct the differential gene libraries of different growth stages of ovaries in two breed goats. Through sequencing, the SSH library of AO group contains 108 ESTs, 52.8% of which were known genes. There were 106 ESTs in the SSH library of EO group, 55.66% of which were known genes. There were 151 ESTs in the SSH library of BO group, 45.03%

 

Table VI. Puberty-related important DEGs among three groups.

Gene name	Description	References
CYP11A1 (cytochrome P450 family 11 subfamily A member 1)	Cleavage of cholesterol to pregnenolone, the precursor of most steroid hormones	Terry et al., 2010; Wickenheisser et al., 2012; Meng-Chun et al., 2018; Moravek et al., 2016
GNAS	Involved as modulators or transducers in various transmembrane signaling systems	Nicolas et al., 2013; Hirofumi Ando et al., 2016
BCAP31 (B cell receptor associated protein 31)	Involved in the anterograde transport of membrane proteins from the endoplasmic reticulum to the Golgi and in caspase 8-mediated apoptosis	Preethi et al., 2015
HSP90B1 (heat shock protein 90 beta family member 1)	Expression of this protein is associated with a variety of pathogenic states, including tumor formation.	Li et al., 2016; Christophe et al., 2011
GAPDH (glyceraldehyde-3-phosphate dehydrogenase)	Catalyzes an important energy-yielding step in carbohydrate metabolism	Liu et al., 2020
mt-co1 (mitochondrially encoded cytochrome c oxidase I)	Cooperate to transfer electrons derived from NADH and succinate to molecular oxygen	Zhen et al., 2015; Ghaffari et al., 2015
NDUFB6 (NADH: ubiquinone oxidoreductase subunit B6)	Transfers electrons from NADH to the respiratory chain	Schutt et al., 2019
SPARC (secreted protein acidic and cysteine rich)	Promote tumor cell invasion	Tumbarello et al., 2016; John et al., 2019
LAMA4 (laminin subunit alpha 4)	Related to cell adhesion, differentiation, migration, signaling, neurite outgrowth and metastasis	Liu et al., 2020
Mylk (myosin light chain kinase)	Facilitate myosin interaction with actin filaments to produce contractile activity	Wallace et al., 2019
COL3A1 (collagen type III alpha 1 chain)	Mutations in this gene are associated with Ehlers-Danlos syndrome types IV, and with aortic and arterial aneurysms.	Engqvist et al., 2019
YBX1 (Y-box binding protein 1)	Aberrant expression of the gene is associated with cancer proliferation in numerous tissues.	Campbell et al., 2017; Lv, 2020
PGM1 (phosphoglucomutase 1)	Mutations in this gene cause glycogen storage disease type 14.	Balakrishnan et al., 2019
COX1 (cytochrome c oxidase subunit I)	Cooperate to transfer electrons derived from NADH and succinate to molecular oxygen	Barnard et al., 2018
EIF3F (eukaryotic translation initiation factor 3 subunit F)	Specifically targets and initiates translation of a subset of mRNAs involved in cell proliferation	Rafael et al., 2018
SLITRK2 (SLIT and NTRK like family member 2)	Involved in synaptogenesis and promotes excitatory synapse differentiation	Chen et al., 2020
ZFX (zinc finger protein X-linked)	As a transcriptional regulator for self-renewal of stem cell types	Song et al., 2020
JMJD6 (Jumonji domain containing 6, arginine demethylase and lysine hydroxylase)	Identified as a putative phosphatidylserine receptor involved in phagocytosis of apoptotic cells	Siu et al., 2009
RGS2 (regulator of G protein signaling 2)	Acts as a mediator of myeloid differentiation and may play a role in leukemogenesis	Jiang et al., 2016; Matsuo et al., 2013
Ikbkap (inhibitor of kappa light polypeptide gene enhancer in B-cells, kinase complex-associated protein)	Involved in neurogenesis.	Yang et al., 2019
bzw1 (basic leucine zipper and W2 domains 1)	Enhances histone H4 gene transcription but does not seem to bind DNA directly.	Xu et al., 2020; Liu et al., 2018
rpl4-b (60S ribosomal protein L4-B)	Translates the encoded message by selecting cognate aminoacyl-transfer RNA (tRNA) molecules	Kolkova et al., 2013

 

of which were known genes. At the same time, the KOG annotation of these libraries showed that the most DEGs are enriched in the translation, ribosomal structure and biogenesis, especially in the B group (59.5%). In addition, 11.77%, 14.7% and 14.87% DEGs are related with energy production and conversion. These may indicate that more genes expressed in ovaries used for energy production and conversion in juvenile and pubertal JG goats than that in LC goats. The percentages of DEGs enriched in each KOG function in the three groups are different. Compared with group B (0.83%), more genes related with RNA processing and modification expressed in the ovaries of juvenile JG goat (A group, 11.77%) and pubertal JG goat (E group, 8.11%) group, which may indicate that the biology process in ovaries of juvenile and pubertal JG goats need more transcripts participating in past-transcription modification. 3.92%, 1.96% and 4% of the DEGs in A group are related with Amino acid, nucleotide and lipid transport and metabolism respectively, but no in E and B group. 9.8% of the DEGs in group A, but only 0.82% in B group and no in E group, take part in transcription process. And the percentage of DEGs related with cytoskeleton in A (9.8%) group is higher than that in B (4.13%) and E (4.05%) group. These may indicate that more active metabolism occurred in the juvenile JG goats’ ovaries.

Also, the percentages of DEGs enriched in some KOG stem are higher in B and E group but lower or no in A group. Some DEGs related with cell motility are in B and E group but not in A group, maybe because that these DEGs have some roles in the follicle motility and development in the pubertal ovaries. 2.7%, 2.7% and 3.3%, 2.47% of the DEGs in E and B group are related with inorganic ion transport and metabolism, and but extracellular structures, respectively, but no in A group, and whether the related inorganic ions and the extracellular structures correlating with ovaries function need further research. The percentages of the DEGs related with signal transduction in E (6.76%) and B (6.61%) group are higher than that in A (3.92%) group. The percentage of the DEGs related with Posttranslational modification, protein turnover, chaperones, and intracellular trafficking, secretion, and vesicular transport in E (14.87%, 8.11%) group is higher than that in A (7.84%, 1.98%) and B (4.96%, 0.83%) group, respectively. These genes may have some undiscovered roles in ovaries among the onset of puberty, or some of them are unique to JG goats for their high prolificacy.

The KEGG pathway analysis of annotation genes showed that the leptin-insulin signal pathway and PI3K-Akt signaling pathway are found in the three groups. Besides this one, insulin secretion, insulin signaling pathway and mTOR signaling pathway are found related with DEGs in B group. These DEGs related with energy may have some undiscovered roles in the onset of puberty. Also, some DEGs in B and E groups take part in the glutathione metabolism pathway. Glutathione plays a key role in the onset of puberty and oocyte growth. The distribution of glutathione in the hypothalamus was increased to the maximum at puberty and decreased to the adult level thereafter. Administration of glutathione to immature female rats brings about a preponement of puberty. The pituitary, ovaries and uterine mass was increased significantly and the pituitary FSH levels were increased in ovariectomized-steroid primed rats (Vali Pasha et al., 1989; Vali-Pasha, 2007). The glutathione concentration in oocytes from prepubertal gilts was significantly higher than that in cyclic oocytes (Pawlak et al., 2015). Oocytes with higher glutathione concentration have more developmental advantage, such as the higher rate of male pronucleus formation, larger size and proper pattern of mitochondria distribution, and higher blastocyst yield (De Matos and Furnas, 2000; Wu et al., 2007; Jiao et al., 2013; Abazari-Kia et al., 2014). The DEGs in BO and EO group involved in the glutathione metabolism, but no in juvenile group, may indicated that the glutathione metabolism may have some effect on the high fecundity and/or precocious puberty of Jing Grey goats. These data are consistent with the possibility that the kisspeptin/GnRH network is intact before puberty but merely inhibited by heightened sensitivity to estradiol negative feedback. In the immature female sheep, preovulatory gonadotropin surges do not occur spontaneously. Although the surge system is capable of function from a very early age, it remains dormant until puberty because of inadequate tonic (pulsatile) LH secretion. Pulses of LH, along with FSH, are ultimately responsible for the production of estradiol by the ovarian follicle.

According to the information of frequencies in the SSH cDNA libraries, the functional information and location in the protein interaction, some important candidate genes were selected as the targets for further study. More differential expressed ESTs and more known gene were in group B than that in group A and E. Many genes appeared more than once in group B and E, and are most related with translation process. For example, EEF1A gene appeared 9 and 6 times in the ovaries of the B and E group; guanine nucleotide-binding protein, alpha-stimulating activity polypeptide 1 (GNAS) gene appeared 5 times in the B group but none in A and B group. NADH dehydrogenase subunit 4(ND4) gene appeared 3 and 4 times in B and E group. In B group, the genes that appeared more than once are more related with ribosomal protein.

GNAS gene, appeared 5 times in the BO group, participated in GnRH signaling pathway. The GNAS gene encodes the α-subunit of the stimulatory guanine nucleotide-binding protein (Gsα), which transduces signals from a G protein-coupled receptor. GNAS has been shown to elevate intracellular cAMP levels by stimulating adenylyl cyclase, which provokes cellular proliferation through the protein kinase A-ERK signal pathway (Landis et al., 1989). Several studies reported that some McCune-Albright Syndrome (MAS) patients with mutation in GNAS gene displayed precocious puberty (Mariot et al., 2011). GNAS gene also have important role in oocyte mature. Gs-alpha activity in the oocyte is required to maintain meiotic arrest within the ovarian follicle (Mehlmann et al., 2006; Mehlmann, 2013). KEGG pathway analysis revealed that GNAS gene takes part in the GnRH signaling pathway in our results. So, GNAS may be an important gene for the onset of puberty.

Secreted protein acidic and rich in cysteine (SPARC) was originally identified as a bone-specific protein. SPARC is a secreted, multi-domain protein, containing an amino-terminal acidic domain that binds hydroxyapatite and calcium ions, a follistatin-like domain containing multiple cysteine residues, and a carboxy-terminal extracellular calcium-binding (EC) domain containing two EF-hand motifs (Sasaki et al., 1998). Exogenous SPARC can promote apoptosis in ovarian cancer cells and elevated SPARC expression has been shown to occur in the activated stroma surrounding ovarian tumors (Brown et al., 1999; Paley et al., 2000), leading to the suggestion that it may be modulating the ovarian tumor microenvironment through regulation of matrix metalloproteinases, inflammation, and pro-migratory cytokines (Said et al., 2007). So, SPARC gene may be an important gene for the onset of puberty in goat.

CYP11A1 encoded the cholesterol side chain cleavage enzyme (P450scc), and cleaves the cholesterol side chain, converting cholesterol to pregenolone, the precursor of androgens, estrogens, and progesterone (Hu et al., 2002; Terry et al., 2010). The pentanucleotide repeat in the 5’ regulatory region of CYP11A1 has been shown to be associated with ovarian androgen excess and polycystic ovary syndrome (PCOS) (Pusalkar et al., 2009). Elevated steady state CYP11A1 mRNA abundance in PCOS cells results from increased transactivation of the CYP11A1 promoter and increased CYP11A1 mRNA stability (Wickenheisser et al., 2012). Meng-Chun Shih et al., found that the promoter mutation of CYP11A1 did not decrease life span and caused no defect in mouse reproduction (Meng-Chun et al., 2018). Most of studies showed that CYP11A1 gene associated with ovary disease and development, which gene might quite be an important candidate gene related to puberty.

CONCLUSION

Overall, several biological processes were selected by SSH, and of which energy production and conversion may be an initiating factor for pubertal onset in the goat. Many genes related with energy production and conversion and GnRH signal pathway may have significant roles in the puberty onset, of which, CYP11A1, GNAS and SPARC will be our candidates for further research.

Data availability

The original data of the paper are available from the corresponding author upon request.

Financial support

This paper has been supported 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).

Statement of conflict of interest

The authors have declared no conflict of interest.

REFERENCES

Abazari-Kia, A.H., Mohammadi-Sangcheshmeh, A., Dehghani-Mohammadabadi, M., Jamshidi-Adegani, F., Veshkini, A., Zhandi, M., Cinar, M.U. and Salehi, M., 2014. Intracellular glutathione content, developmental competence and expression of apoptosis-related genes associated with G6PDH activity in goat oocyte. J. Assist. Reprod. Genet., 31: 313-321. https://doi.org/10.1007/s10815-013-0159-y

Balakrishnan, B., Verheijen, J., Lupo, A., Raymond, K., Turgeon, C.T., Yang, Y., Carter, K.L., Whitehead, K.J., Kozicz, T., Morava, E. and Lai, K., 2019. A novel phosphoglucomutase-deficient mouse model reveals aberrant glycosylation and early embryonic lethality. J. Inher. Metabol. Dis., 42: 998-1007. https://doi.org/10.1002/jimd.12110

Barnard, M.E., Hecht, J.L., Rice, M.S., Gupta, M., Harris, H.R., Heather Eliassen, A., Rosner, B.A., Terry, K.L, and Tworoger, S.S., 2018. Anti-inflammatory drug use and ovarian cancer risk by COX1/COX2 expression and infiltration of tumor-associated macrophages. Cancer Epidemiol. Biomark. Prevent., 27: 1509-1517. https://doi.org/10.1158/1055-9965.EPI-18-0346

Brown, T.J., Shaw, P.A., Karp, X., Huynh, M.H., Begley, H. and Ringuette, M.J., 1999. Activation of SPARC expression in reactive stroma associated with human epithelial ovarian cancer. Gynecol. Oncol., 75: 25-33. https://doi.org/10.1006/gyno.1999.5552

Campbell, T.M., Castro, M.A.A., de Oliveira, K.G., Ponder, B.A.J. and Meyer, K.B., 2017. ERα binding by transcription factors NFIB and YBX1 enables FGFR2 signaling to modulate estrogen responsiveness in breast cancer. Cancer Res., 78: 410-421. https://doi.org/10.1158/0008-5472.CAN-17-1153

Cao, G.L., Feng, T., Chu, M.X., Di, R., Zhang, Y.L., Huang, D.W., Liu, Q.Y., Hu, W.P. and Wang, X.Y., 2015. Subtraction suppressive hybridisation analysis of differentially expressed genes associated with puberty in the goat hypothalamus, reproduction. Fertil. Dev., 28: 1781-1787. https://doi.org/10.1071/RD14434

Chen, M., Zeng, J., Chen, S.Y., Li, J.J., Wu, H.J., Dong, X.H., Lei, Y., Zhi, X.L. and Yao, L.Q., 2020. SPTBN1 suppresses the progression of epithelial ovarian cancer via SOCS3-mediated blockade of the JAK/STAT3 signaling pathway. Aging, 12: 10896-10911. https://doi.org/10.18632/aging.103303

Christophe, A., Florent, L.M., Colette, C., Zihai, L. and Elisabeth, S.C., 2011. Oocyte–targeted deletion reveals that Hsp90b1 is needed for the completion of first mitosis in mouse zygotes. PLoS One, 6: e17109. https://doi.org/10.1371/journal.pone.0017109

De Matos, D.G. and Furnus, C.C., 2000. The importance of having high glutathione (GSH) level after bovine in vitro maturation on embryo development: Effect of beta-mercaptoethanol, cysteine and cystine. Theriogenology, 53: 761-771. https://doi.org/10.1016/S0093-691X(99)00278-2

De Sanctis, V., Soliman, A.T., Di Maio, S., Soliman, N. and Elsedfy, H., 2019. Long-term effects and significant adverse drug reactions (ADRs) associated with the use of gonadotropin-releasing hormone analogs (GnRHa) for central precocious puberty: A brief review of literature. Acta bio-med. Atenei Parmensis, 90: 345–359.

Diatchenko, L., Lau, Y.F., Campbell, A.P., Chenchik, A., Moqadam, F., Huang, B., Lukyanov, S., Lukyanov, K., Gurskaya, N., Sverdlov, E.D. and Siebert, P.D., 1996. Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc. natl. Acad. Sci. USA, 93: 6025-6030. https://doi.org/10.1073/pnas.93.12.6025

Dulka, E.A., Anthony DeFazio, R. and Moenter, S.M., 2020. Chemogenetic suppression of GnRH neurons during pubertal development can alter adult GnRH neuron firing rate and reproductive parameters in female mice. eNeuro, 7: ENEURO.0223-20.2020. https://doi.org/10.1523/ENEURO.0223-20.2020

Engqvist, H., Parris, T.Z., Anikó, K., Nemes, S., Rönnerman, E.W., De Lara, S., Biermann, J., Sundfeldt, K., Karlsson, P. and Helou, K., 2019. Immunohistochemical validation of COL3A1, GPR158 and PITHD1 as prognostic biomarkers in early-stage ovarian carcinomas. BMC Cancer, 19: 928. https://doi.org/10.1186/s12885-019-6084-4

Foster, D.L. and Hileman, S.M., 2015. Puberty in the sheep. In: Knobil and Nll’s physiology of reproduction (Fourth Edition), pp. 1441-1485. https://doi.org/10.1016/B978-0-12-397175-3.00031-4

Ghaffari Novin, M., Allahveisi, A., Noruzinia, M., Farhadifar, F., Yousefian, E., Dehghani Fard, A., and Salimi, M., 2015. The relationship between transcript expression levels of nuclear encoded (TFAM, NRF1) and mitochondrial encoded (MT-CO1) genes in single human oocytes during oocyte maturation. Balkan J. med. Genet., 18: 39-46. https://doi.org/10.1515/bjmg-2015-0004

Gurskaya, N.G., Diatchenko, L., Chenchik, A., Siebert, P.D., Khaspekov, G.L., Lukyanov, K.A., Vagner, L.L., Ermolaeva, O.D., Lukyanov, S.A. and Sverdlov, E.D., 1996. Equalizing cDNA subtraction based on selective suppression of polymerase chain reaction: cloning of Jurkat cell transcripts induced by phytohemaglutinin and phorbol 12-myristate 13-acetate. Anal. Biochem., 240: 90-97. https://doi.org/10.1006/abio.1996.0334

Hirofumi, A., Tsutomu, M., Hiroyasu, K., Akiko, T., Keiko, I., Yasunari, F. and Tanri S., 2016. Usefulness of a management protocol for patients with cervical multicystic lesions: A retrospective analysis of 94 cases and the significance of GNAS mutation. J. Obstet. Gynaecol. Res., 42: 1588-1598. https://doi.org/10.1111/jog.13083

Hu, M.C., Hsu, N.C., El Hadj, N.B., Pai, C.I., Chu, H.P., Wang, C.K. and Chung, B.C., 2002. Steroid deficiency syndromes in mice with targeted disruption of Cyp11a1. Mol. Endocrinol., 16: 1943-1950. https://doi.org/10.1210/me.2002-0055

Jiang, M.X., Shi, Y., Sun, Z.G., Zhang, Z. and Zhu, Y., 2016. Inhibition of the binding between RGS2 and β-tubulin interferes with spindle formation and chromosome segregation during mouse oocyte maturation in vitro. PLoS One, 11: e0159535. https://doi.org/10.1371/journal.pone.0159535

Jiao, G.Z., Cao, X.Y., Cui, W., Lian, H.Y., Miao, Y.L., Wu, X.F., Han, D. and Tan, J.H., 2013. Developmental potential of prepubertal mouse oocytes is compromised due mainly to their impaired synthesis of glutathione. PLoS One, 8: e58018. https://doi.org/10.1371/journal.pone.0058018

John, B., Naczki, C., Patel, C., Ghoneum, A., Qasem, S., Salih, Z. and Said, N., 2019. Regulation of the bi-directional cross-talk between ovarian cancer cells and adipocytes by SPARC. Oncogene, 38: 4366-4383. https://doi.org/10.1038/s41388-019-0728-3

Kolkova, Z., Arakelyan, A., Bertil, C., Hansson, S. and Kriegova, E., 2013. Normalizing to GADPH jeopardises correct quantification of gene expression in ovarian tumours-IPO8 and RPL4 are reliable reference genes. J. Ovar. Res., 6: 60. https://doi.org/10.1186/1757-2215-6-60

Landis, C.A., Masters, S.B., Spada, A., Pace, A.M., Bourne, H.R. and Vallar, L., 1989. GTPase inhibiting mutations activate the alpha chain of Gs and stimulate adenylyl cyclase in human pituitary tumours. Nature, 340: 692-696. https://doi.org/10.1038/340692a0

Li, L., Mo, H., Zhang, J., Zhou, Y.X., Peng, X.H. and Luo, X.P., 2016. The role of heat shock protein 90B1 in patients with polycystic ovary syndrome. PLoS One, 11: e0152837. https://doi.org/10.1371/journal.pone.0152837

Liu, F., Zhao, H., Gong, L., Yao, L. and Li, Y.H., 2018. MicroRNA-129-3p functions as a tumor suppressor in serous ovarian cancer by targeting BZW1. Int. J. clin. exp. Pathol., 11: 5901-5908.

Liu, P., Zhong, Y., Cao, T., Sheng, X.J. and Huang, H., 2020. A frequent somatic mutation in the 3’UTR of GAPDH facilitates the development of ovarian cancer by creating a miR125b binding site. Oncol. Rep., 44: 887-896. https://doi.org/10.3892/or.2020.7663

Liu, Y., Xu, Y., Ding, L., Yu, L.L., Zhang, B.T. and Wei, D., 2020. LncRNA MEG3 suppressed the progression of ovarian cancer via sponging miR-30e-3p and regulating LAMA4 expression. Cancer Cell Int., 20: 181. https://doi.org/10.1186/s12935-020-01259-y

Lv, Y., 2020. Drosophila YBX1 homolog YPS promotes ovarian germ line stem cell development by preferentially recognizing 5-methylcytosine RNAs. Proc. natl. Acad. Sci. U.S.A., 117: 3603-3609. https://doi.org/10.1073/pnas.1910862117

Mariot, V., Wu, J.Y., Aydin, C., Mantovani, G., Mahon, M.J., Linglart, A., and Bastepe, M., 2011. Potent constitutive cyclic AMP-generating activity of XLαs implicates this imprinted GNAS product in the pathogenesis of McCune-Albright syndrome and fibrous dysplasia of bone. Bone, 48: 312-320. https://doi.org/10.1016/j.bone.2010.09.032

Matsuo, M., Coon, S.L. and Klein, D.C., 2013. RGS2 is a feedback inhibitor of melatonin production in the pineal gland. FEBS Lett., 587: 1392-1398. https://doi.org/10.1016/j.febslet.2013.03.016

Mehlmann, L.M., 2013. Losing Mom’s Message: Requirement for DCP1A and DCP2 in the degradation of maternal transcripts during oocyte maturation. Biol. Reprod., 88: 10. https://doi.org/10.1095/biolreprod.112.106591

Mehlmann, L.M., Kalinowski, R.R., Ross, L.F., Parlow, A.F., Hewlett, E.L. and Jaffe, L.A., 2006. Meiotic resumption in response to luteinizing hormone is independent of a Gi family G protein or calcium in the mouse oocyte. Dev. Biol., 299: 345-355. https://doi.org/10.1016/j.ydbio.2006.07.039

Meng-Chun, S., Nai-Chi, H., Chen-Che, H., Tzong-Shoon, W., Pao-Yen, L., Bon-chu Bon-chu., 2018. Mutation of mouse Cyp11a1 promoter caused tissue-specific reduction of gene expression and blunted stress response without affecting reproduction. Mol. Endocrinol., 4: 915-923.

Moenter, S.M., Caraty, A., Locatelli, A. and Karsch, F.J., 1991. Pattern of gonadotropin-releasing hormone (GnRH) secretion leading up to ovulation in the ewe: existence of a preovulatory GnRH surge. Endocrinology, 129: 1175-1182. https://doi.org/10.1210/endo-129-3-1175

Moravek, M.B., Shang, M., Menon, B. and Menon, K.M.J., 2016. HCG-mediated activation of mTORC1 signaling plays a crucial role in steroidogenesis in human granulosa lutein cells. Endocrine, 54: 217-224. https://doi.org/10.1007/s12020-016-1065-8

Nicolas, R., Arnaud, M., Nadia, C., Pauline, R.G., Harald, J. and Marie-Laure, K., 2013. Paternal GNAS mutations lead to severe intrauterine growth retardation (IUGR) and provide evidence for a role of XLαs in fetal development. J. clin. Endocrinol. Metab., 98: E1549-E1556. https://doi.org/10.1210/jc.2013-1667

Paley, P.J., Goff, B.A., Gown, A.M., Greer, B.E. and Sage, E.H., 2000. Alterations in SPARC and VEGF immunoreactivity in epithelial ovarian cancer. Gynecol. Oncol., 78: 336-341. https://doi.org/10.1006/gyno.2000.5894

Pawlak, P., Warzych, E., Hryciuk, M. and Lechniak, D., 2015. Transcript abundance, glutathione and apoptosis levels differ between porcine oocytes collected from prepubertal and cyclic gilts. Theriogenology, 84: 86-93. https://doi.org/10.1016/j.theriogenology.2015.02.016

Plant, T.M., 2015. Neuroendocrine control of the onset of puberty. Front. Neuroendocrinol., 38: 73-88. https://doi.org/10.1016/j.yfrne.2015.04.002

Preethi, S., Schomay, T.E., Aiello, K.A. and Alter, O., 2015. Tensor GSVD of patient- and platform-matched tumor and normal DNA copy-number profiles uncovers chromosome arm-wide patterns of tumor-exclusive platform-consistent alterations encoding for cell transformation and predicting ovarian cancer survival. PLoS One, 10: e0121396. https://doi.org/10.1371/journal.pone.0121396

Pusalkar, M., Meherji, P., Gokral, J., Chinnaraj, S. and Maitra, A., 2009. CYP11A1 and CYP17 promoter polymorphisms associate with hyperandrogenemia in polycystic ovary syndrome. Fertil. Steril., 92: 653-659. https://doi.org/10.1016/j.fertnstert.2008.07.016

Rafael, Cuesta, Adi, Y., Berman, Anya, Alayev, Marina, K. and Holz., 2018. Estrogen receptor α promotes protein synthesis by fine-tuning the expression of the eukaryotic translation initiation factor 3 subunit f (eIF3f). J. biol. Chem., 294: 2267-2278. https://doi.org/10.1074/jbc.RA118.004383

Rebrikov, D.V., Desai, S.M., Siebert, P.D. and Lukyanov, S.A., 2004. Suppression subtractive hybridization. Methods mol. Biol., 258: 107-134. https://doi.org/10.1385/1-59259-751-3:107

Said, N., Socha, M.J., Olearczyk, J.J., Elmarakby, A.A., Imig, J.D. and Motamed, K., 2007. Normalization of the ovarian cancer microenvironment by SPARC. Mol. Cancer Res., 5: 1015-1030. https://doi.org/10.1158/1541-7786.MCR-07-0001

Sasaki, T., Hohenester, E., Göhring, W. and Timpl, R., 1998. Crystal structure and mapping by site-directed mutagenesis of the collagen-binding epitope of an activated form of BM-40/SPARC/osteonectin. EMBO J., 17: 1625-1634. https://doi.org/10.1093/emboj/17.6.1625

Schutt, A.K., Blesson, C.S., Hsu, J.W., Valdes, C.T., Gibbons, W.E., Jahoor, F. and Yallampalli, C., 2019. Preovulatory exposure to a protein-restricted diet disrupts amino acid kinetics and alters mitochondrial structure and function in the rat oocyte and is partially rescued by folic acid. Reprod. Biol. Endocrinol., 17: 12. https://doi.org/10.1186/s12958-019-0458-y

Siu, M.K., Wong, O.G. and Cheung, A.N., 2009. TrkB as a therapeutic target for ovarian cancer. Exp. Opin. Therapeut. Targets, 13: 1169. https://doi.org/10.1517/14728220903196787

Smith, M., 2001. Neural signals that regulate GnRH neurones directly during the oestrous cycle. Reproduction, 122: 1. https://doi.org/10.1530/reprod/122.1.1

Song, N., Zhang, Y., Kong, F.F., Yang, H. and Ma, X.X., 2020. HOXA-AS2 promotes type I endometrial carcinoma via miRNA-302c-3p-mediated regulation of ZFX. Cancer Cell Int., 20: 359. https://doi.org/10.1186/s12935-020-01443-0

Terry, K., Mcgrath, M., Lee, I.M., Buring, J. and De Vivo, I., 2010. Genetic variation in CYP11A1 and StAR in relation to endometrial cancer risk. Gynecol. Oncol., 117: 255-259. https://doi.org/10.1016/j.ygyno.2010.02.002

Tumbarello, D.A., Andrews, M.R., Brenton, J.D. and Cukierman, E., 2016. SPARC regulates transforming growth factor beta induced (TGFBI) extracellular matrix deposition and paclitaxel response in ovarian cancer cells. PLoS One, 11: e0162698. https://doi.org/10.1371/journal.pone.0162698

Vali Pasha, K. and Vijayan E., 1989. Glutathione distribution in rat brain at different ages and the effect of intraventricular glutathione on gonadotropin levels in ovariectomized steroid primed rats. Brain Res. Bull., 22: 617-619. https://doi.org/10.1016/0361-9230(89)90079-8

Vali Pasha, K., 2007. Involvement of glutathione in puberty and FSH release. Neurosci. Lett., 423: 78-81. https://doi.org/10.1016/j.neulet.2007.06.024

Wallace, S.E., Regalado, E.S., Limin, G., Janda, A.L., Guo, D.C., Russo, C.F., Kulmacz, R.J., Hanna, N., Jondeau, G., Boileau, C., Arnaud, P., Lee, K., Leal, S.M., Hannuksela, M., Carlberg, B., Johnston, T., Antolik, C., Hostetler, E.M., Colombo, R. and Milewicz, D.M., 2019. MYLK pathogenic variants aortic disease presentation, pregnancy risk, and characterization of pathogenic missense variants. Genet. Med., 21: 144-151. https://doi.org/10.1038/s41436-018-0038-0

Wickenheisser, J.K., Biegler, J.M., Nelson-Degrave, V.L., Legro, R.S., Strauss, J.F., McAllister, J.M. and Franks, S., 2012. Cholesterol side-chain cleavage gene expression in theca cells: Augmented transcriptional regulation and mRNA stability in polycystic ovary syndrome. PLoS One, 7: e48963. https://doi.org/10.1371/journal.pone.0048963

Wu, Y.G., Liu, Y., Zhou, P., Lan, G.C., Han, D., Miao, D.Q. and Tan, J.H., 2007. Selection of oocytes for in vitro maturation by brilliant cresyl blue staining: a study using the mouse model. Cell Res., 17: 722-731. https://doi.org/10.1038/cr.2007.66

Xu, H., Sun, X., Huang, Y., Si, Q. and Li, M.K., 2020. Long noncoding RNA NEAT1 modifies cell proliferation, colony formation, apoptosis, migration and invasion via the miR4500/BZW1 axis in ovarian cancer. Mol. Med. Rep., 22: 3347–3357. https://doi.org/10.3892/mmr.2020.11408

Yang, K.T., Inoue, A., Lee, Y.J., Jiang, C.L. and Lin, F.J., 2019. Loss of Ikbkap/Elp1 in mouse oocytes causes spindle disorganization, developmental defects in preimplantation embryos and impaired female fertility. Sci. Rep., 9: 18875. https://doi.org/10.1038/s41598-019-55090-1

Zhen, X., Wu, B., Wang, J., Lu, C.L., Gao, H.F. and Qiao, J., 2015. Increased incidence of mitochondrial cytochrome C oxidase 1 gene mutations in patients with primary ovarian insufficiency. PLoS One, 10: e0132610. https://doi.org/10.1371/journal.pone.0132610