The heritability of litter size in sheep is low and controlled by multiple genes, but the research on its related genes is not sufficient. Here, to explore the expression pattern of multi-tissue genes in Chinese native sheep, we selected 10 tissues of the three adult ewes with the highest estimated breeding value in the early study of the prolific Xinggao sheep population. The global gene expression analysis showed that the ovary, uterus, and hypothalamus expressed the most genes. Using the Uniform Manifold Approximation and Projection (UMAP) cluster analysis, these samples were clustered into eight clusters. The functional enrichment analysis showed that the genes expressed in the spleen, uterus, and ovary were significantly enriched in the Ataxia Telangiectasia Mutated Protein (ATM) signaling pathway, and most genes in the liver, spleen, and ovary were enriched in the immune response pathway. Moreover, we focus on the expression genes of the hypothalamic–pituitary–ovarian axis (HPO) and found that 11,016 genes were co-expressed in the three tissues, and different tissues have different functions, but the oxytocin signaling pathway was widely enriched. To further explore the differences in the expression genes (DEGs) of HPO in different sheep breeds, we downloaded the transcriptome data in the public data, and the analysis of DEGs (Xinggao sheep vs. Sunite sheep in Hypothalamus, Xinggao sheep vs. Sunite sheep in Pituitary, and Xinggao sheep vs. Suffolk sheep in Ovary) revealed the neuroactive ligand–receptor interactions. In addition, the gene subsets of the transcription factors (TFs) of DEGs were identified. The results suggest that 51 TF genes and the homeobox TF may play an important role in transcriptional variation across the HPO. Altogether, our study provided the first fundamental resource to investigate the physiological functions and regulation mechanisms in sheep. This important data contributes to improving our understanding of the reproductive biology of sheep and isolating effecting molecular markers that can be used for genetic selection in sheep.

Spermatogenesis is a complicated course of several rigorous restrained steps that spermatogonial stem cells undergo to develop into highly specialized spermatozoa; however, specific genes and signal pathways, which regulate the amplification, differentiation and maturation of these cells, remain unclear. We performed bioinformatics analyses to investigate the dynamic changes of the gene expression patterns at three time points in the course of the first wave of rat spermatogenesis. Differently expressed genes (DEGs) were identified, and the features of DEGs were further analyzed with GO (Gene Ontology), KEGG (Kyoto Encyclopedia of Genes and Genomes) and Short Time-series Expression Miner (STEM). A total of 2954 differentially expressed genes were identified. By using STEM, the top 10 key genes were selected in the profile according to the enrichment results, and the distinguishable biological functions encoded by these DEGs were automatically divided into three parts. Genes from 6, 8 and 10 days were related to biosynthesis, immune response and cell junction, and genes from 14, 15 and 16 days were related to energy metabolic pathways. The results also suggest that genes from 29, 31 and 35 days may shift metabolic to sperm motility, sperm flagellum and cilium movement.

Xinggao sheep are a breed of Chinese domestic sheep that are adapted to the extremely cold climatic features of the Hinggan League in China. The economically vital reproductive trait of ewes (litter size, LS) and productive traits of lambs (birth weight, BWT; weaning weight, WWT; and average daily gain, ADG) are expressed in females and later in life after most of the selection decisions have been made. This study estimated the genetic parameters for four traits to explore the genetic mechanisms underlying the variation, and we performed genome-wide association study (GWAS) tests on a small sample size to identify novel marker trait associations (MTAs) associated with prolificacy and growth. We detected two suggestive significant single-nucleotide polymorphisms (SNPs) associated with LS and eight significant SNPs for BWT, WWT, and ADG. These candidate loci and genes also provide valuable information for further fine-mapping of QTLs and improvement of reproductive and productive traits in sheep.

Follicle-stimulating hormone (FSH) regulates ovarian follicle development through specific gene expression programs. Granulosa cells (GCs) are somatic cells surrounding the oocytes, secreting gonadotropins to regulate ovulation and promote follicular development. By analyzing the effects of different doses of FSH on the proliferation of GCs, we found that adding 10 ng/mL of FSH, as the optimal concentration, could promote the growth of GCs. Furthermore, we have successfully constructed the first CRISPR-Cas9 knockout library targeting the genes on chromosomes 2 and 3 and the X chromosomes of the sheep massively parallel coding gene, as well as an ovarian GCs knockout cell library. For the first time, we have exposed the knockout cell library to a concentration of 10 ng/mL FSH to explore the underlying mechanisms. Through this screening, we have identified 836 positive–negative screening genes that are responsive to FSH, thereby revealing the regulatory mechanisms and screening the functionality of candidate genes. Next, RNA-Seq of control (0 ng/mL), low (10 ng/mL), and high (100 ng/mL) doses of FSH revealed 1708 differentially expressed genes, and combined with 836 genes, we obtained 129 FSH dose-dependent genes with extremely significant differences. This enables us to delve deeper into investigating and identifying the mechanisms by which FSH regulates GCs. More generally, we have discovered new regulatory factors and identified reproductivity-associated major effectors. These findings provide novel research directions for further studies on sheep reproduction.