Search

The Booroola fecundity gene (FecB) has a mutation that was found to increase the ovulation rate and litter size in Booroola Merino sheep. This mutation is also associated with the fecundity of small-tail han (STH) sheep, an important maternal breed used to produce hybrid offspring for mutton production in China. Previous research showed that the FecB gene affects reproduction in STH sheep, based on litter size records. However, the effects of this gene on estrus, ovulation, and endocrine characteristics in these sheep remain unclear. Here, we analyzed the traits mentioned earlier and compared them among the three FecB genotypes of STH ewes using estrus synchronization. Overall, 53 pluriparous ewes were selected from among 890 STH ewes and subjected to FecB genotyping for experiments to characterize estrous and ovulation rates. FecB heterozygous (+B) ewes presented an earlier onset of estrus (42.9 ± 2.2 h) and a shorter estrous cycle (17.2 ± 0.2 days) (P ≤ 0.05). The ovulation rates increased with the increasing copy number of the B allele (P ≤ 0.01). Ovulation time showed no significant differences among the three FecB genotypes. The serum concentrations of follicle-stimulating hormone (FSH), luteinizing hormone, estrogen (E2), and progesterone (P4) were measured in 19 of the ewes. Serum concentrations of E2 and FSH dramatically varied around the time of behavioral estrus. In FecB mutant homozygous (BB) ewes, E2 concentration had two peaks, which were higher (P ≤ 0.05) than those of ++ genotypes. FSH concentration of BB ewes was higher (P ≤ 0.05) than that of the ++ ewes just after estrus. The expression of the estrogen receptor 1 (ESR1) gene in the +B genotype was higher than in the other genotypes. Based on the data for the reproductive performance of STH ewes with the three FecB genotypes, our study suggests that the development of follicles in ewes with the B allele is dependent on the response to FSH regulated by E2 in the early stage. +B ewes, exhibiting moderate ovulation and litter size and a shorter estrous cycle, can be highly recommended in sheep crossbreeding systems for commercial mutton production. Moreover, this study provides useful information to conserve better and use the genetic resources of STH sheep in China.

MicroRNA (miRNA) is a kind of noncoding RNA whose function involved in various biological processes in neuronal maturation and adipocyte cells, such as differentiation, proliferation, development, apoptosis, and metabolism. Herein, miRNA-Seq was used to identify miRNAs in the tail fat tissue of Hu sheep (short-fat-tailed) and Tibetan sheep (short-thin-tailed). In this study, 155 differentially expression miRNAs (DE miRNAs) were identified, including 78 up-regulated and 77 down-regulated. Among these DE miRNAs, 17 miRNAs were reported and related with lipid metabolism. MiRanda and RNAhybrid software were used to predict the target genes of DE miRNAs, obtaining the number of targeting relationships is 38553. Target genes were enriched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). 742 terms and 302 single pathways are enriched, including lipid metabolic process, response to lipid, cellular lipid catabolic process, lipid catabolic process, cellular lipid metabolic process, inositol lipid-mediated signaling, calcium channel activity, PI3K-Akt signaling pathway, MAPK signaling pathway, ECM-receptor interaction, AMPK signaling pathway, Wnt signaling pathway and TGF-beta signaling pathway. Notably, miR-379-5p was associated with tail fat deposition of sheep. Dual-Luciferase reporter assays showed miR-379-5p and HOXC9 had targeted relationship. The result of RT-qPCR showed that the expression trend of miR-379-5p and HOXC9 was opposite. miR-379-5p was down-regulated and highly expressed in tail adipose tissue of Tibetan sheep. HOXC9 was highly expressed in adipose tissue of Hu sheep. These results could provide a meaningful theoretical basis for studying the molecular mechanisms of sheep tail adipogenesis.

Most sheep are seasonal estrus, and they breed in autumn when the days get shorter. Seasonal estrus is an important factor that affects the productivity and fertility of sheep. The key point to solve this problem is to explore the regulation mechanism of estrus in sheep. Therefore, in this study, transcriptomic sequencing technology was used to identify differentially expressed mRNAs in the hypothalamus, pituitary and ovary of Small Tail Han sheep (year-round estrus) and tan sheep (seasonal estrus) among luteal, proestrus and estrus stages. There were 256,923,304,156 mRNAs being identified in the hypothalamus, pituitary and ovary, respectively. Functional analysis showed that the photosensor, leucine and isoleucine biosynthesis pathways were enriched significantly. It is speculated that photoperiod may initiate estrus by stimulating the corresponding pathways in hypothalamus. ODC1, PRLH, CRYBB2, SMAD5, OPN1SW, TPH1 are believed to be key genes involved in the estrogen process. In conclusion, this study expanded the database of indigenous sheep breeds, and also provided new candidate genes for future genetic and molecular studies on the seasonal estrus trait in sheep.

The pituitary is a remarkably dynamic organ with roles in hormone (FSH and LH) synthesis and secretion. In animals with the

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.

The pituitary pars tuberalis (PT) is the regulating center of seasonal reproduction, which can sense the melatonin signal and eventually cause downstream changes of GnRH secretion through TSHβ. Recently, lncRNAs have been identified in animal reproductive-related tissues, and they play important roles in reproductive regulation. Therefore, in this study, we expect to identify photoperiod-induced lncRNAs and genes in pituitary PT of sheep by comparison of expression profiles between short photoperiod (SP) and long photoperiod (LP). Through RNA-Seq, a total of 55,472 lncRNAs were identified in pituitary PT of Sunite ewes. The number of differentially expressed (DE) genes and lncRNAs between SP and LP increased gradually with the extension of LP (from LP7 to LP42). The notable LP-induced candidate genes included EYA3, TSHB, SIX1, DCT, VMO1, AREG, SUV39H2, and EZH2, and SP-induced genes involved ENSOARG00000012585, CHGA, FOS, SOCS3, and TH. In enriched pathways for DE genes and lncRNA target genes between SP and LP, the reproduction- and circadian-related pathways were highlighted. In addition, the interactome analysis of lncRNAs and their targets implied that MSTRG.209166 and its trans-target TSHB, MSTRG.288068 and its cis-target SIX1, and ENSOARG00000026131 and its cis-target TH might participate in regulation of seasonal reproduction. Together, these results will help to determine important photoperiod-induced lncRNAs and genes and give us some new insights into the epigenetic regulation of seasonal reproduction in sheep.

Bone is continuously remodeled in a dynamic process maintained by osteoclasts and osteoblasts, and imbalances in the relative activities of these cell types can cause various pathological conditions, including rheumatoid arthritis and osteoporosis. Osteoclasts are multinucleated cells that serve an important role in regulating the development of osteoporosis. Furthermore, P2X7 receptor activation has a vital role in physiological and pathological reactions in bone, including bone disease. Therefore, the present study aimed to investigate the effect of P2X7 receptor on osteoclast differentiation and to explore the underlying molecular mechanism by western blotting and tartrate‑resistant acid phosphatase staining. The results indicated that the expression levels of P2X7 receptor and intracellular Ca2+ concentration levels were very high in mature osteoclasts. Furthermore, P2X7 receptor overexpression increased the number of multinucleated osteoclasts and the expression of osteoclastogenesis‑related proteins. P2X7 receptor overexpression was also associated with downstream activation of Ca

Although it has been investigated for many years, the physiological processes regulating prolificacy in sheep remains unclear because of regulation by many genes. To better understand the effects of three single nucleotide polymorphisms (SNPs) comprising g.48462350CT in BMP2, g.58171856CG and g.58171886AC in BMP7, a population genetic analysis was conducted using data obtained from genotyping in 768 sheep from six breeds (three polytocous and three monotocous). The results indicate that all the sheep breeds were considered to conform to the Hardy-Weinberg equilibrium (P 0.05). The associations of these three SNPs with litter size in 384 Small Tail Han sheep were analyzed, therefore, and found to be correlated with fecundity as assessed by mean litter size (P 0.05). Bioinformatic analysis indicated there was a transmembrane domain change that occurred after a mutation in BMP2 at g.48462350CT, and changes involving transcription factors such as USF1, USF2 and INMS1 in the BMP7 promoter region might be involved in greater sheep prolificacy.

Qira black sheep is a famous indigenous sheep breed in China. The objectives of this study are to identify candidate genes related to body size, and to estimate the level of inbreeding depression on body size based on runs of homozygosity in Qira black sheep. Here, 188 adult Qira black sheep were genotyped with a high density (630 K) SNP chip and genome-wide association study for body weight and body size traits (including withers height, body slanting length, tail length, chest girth, chest width, and chest depth) were performed using an additive linear model. In consequence, 12 genome- and chromosome-wide significant SNPs and, accordingly, six candidate genes involved in muscle differentiation, metabolism and cell processes were identified. Of them, ZNF704 (zinc finger protein 704) was identified for body weight; AK2 (adenylate kinase 2) and PARK2 (parkin RBR E3 ubiquitin protein ligase) for tail length; MOCOS (molybdenum cofactor sulfurase) and ELP2 (elongator acetyltransferase complex subunit 2) for chest width; and MFAP1 (microfibril associated protein 1) for chest girth. Additionally, inbreeding depressions on body size were observed in the current herd. These results will provide insightful understandings into the genetic mechanisms of adult body size, and into the conservation and utilization of Qira black sheep.

Simple Summary Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. There is large variation in the reproductive abilities of Bamei mutton sheep. After recent artificial selection, the average lambing rate of the Bamei mutton nucleus group was over 150%. We used the FST (Fixation Index) and XP-EHH (The Cross-Population Extended Haplotype Homozygosity) statistical approach to detect the selective sweeps between high- and low-fecundity Bamei mutton sheep groups. JUN (JUN proto-oncogene, AP-1 transcription factor subunit), ITPR3 (inositol 1,4,5-trisphosphate receptor type 3, PLCB2 (phospholipase C beta 2), HERC5 (HECT and RLD domain containing E3 ubiquitin protein ligase 5), and KDM4B (lysine demethylase 4B) were detected that are potential responsible for litter size. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. Abstract Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on detecting candidate genes associated with the increasing of the litter size in this breeds under recent artificial selection to improve the efficiency of mutton production. We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH (gonadotropin-releasing hormone), oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 (forkhead box A1) and ESR1 (estrogen receptor 1). Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions that are responsible for litter size can in turn be used for further selection.

Ovine litter size (LS) is an important trait showing variability within breeds. It remains largely unknown whether inbreeding depression on LS exists based on genomic homozygous regions, and whether the homozygous regions resulted from inbreeding are significantly associated with LS in sheep. We here reanalyze a set of single nucleotide polymorphism (SNP) chip of six breeds to characterize the patterns of runs of homozygosity (ROH), to evaluate inbreeding levels and inbreeding depressions on LS, and to identify candidate homozygous regions responsible for LS. Consequently, unique ROH patterns were observed among six sheep populations. Inbreeding depression on LS was only found in Hu sheep, where a significant reduction of 0.016, 0.02, and 0.02 per 1% elevated inbreeding FROH4&ndash, 8, FROH&gt, 8 and the total inbreeding measure was observed, respectively. Nine significantly homozygous regions were found for LS in Hu sheep, where some promising genes for LS possibly via regulation of the development of oocytes (NGF, AKT1, and SYCP1), fertilization (SPAG17, MORC1, TDRD9, ZFYVE21, ADGRB3, and CKB), embryo implantation (PPP1R13B, INF2, and VANGL1) and development (DPPA2, DPPA4, CDCA4, CSDE1, and ADSSL1), and reproductive health (NRG3, BAG5, CKB, and XRCC3) were identified. These results from the present study would provide insights into the genetic management and complementary understandings of LS in sheep.

Long non-coding RNA (lncRNA) accounts for a large proportion of RNA in animals. The thyroid gland has been established as an important gland involved in animal reproduction, however, little is known of its gene expression patterns and potential roles in the sheep. Herein, RNA-Seq was used to detect reproduction-related differentially expressed lncRNAs (DELs) and mRNAs (DEGs) in the follicular phase (FT) FecBBB (MM) and FecB++ (ww) genotypes of Small Tail Han (STH) sheep thyroids. Overall, 29 DELs and 448 DEGs in thyroid between MM and ww sheep were screened. Moreover, GO and KEGG enrichment analysis showed that targets of DELs and DEGs were annotated in biological transitions, such as cell cycle, oocyte meiosis and methylation, which in turn affect reproductive performance in sheep. In addition, we constructed co-expression and networks of lncRNAs-mRNAs. Specifically, XLOC_075176 targeted MYB, XLOC_014695 targeted VCAN, 106991527 targeted CASR, XLOC_075176 targeted KIFC1, XLOC_360232 targeted BRCA2. All these differential lncRNAs and mRNAs expression profiles in the thyroid provide a new resource for elucidating the regulatory mechanism underlying STH sheep prolificacy.

The thyroid gland is a vital endocrine organ involved in the reproduction of animals via the regulation of hormone synthesis and secretion. LncRNAs have been proven to play important roles in reproductive regulation; however, the associated mechanism in the thyroid gland has not been clarified. In this study, we investigated to identify photoperiod-induced lncRNAs and mRNAs in the thyroid gland in Sunite ewes by comparing the expression profiles of short photoperiod (SP) and long photoperiods (LP). A total of 41,088 lncRNAs were identified in the thyroid gland through RNA-Seq. Functional analysis of differentially expressed lncRNAs using the R package revealed that reproductive hormone- and photoperiod response-related pathways, including the prolactin signaling, cAMP signaling, and circadian rhythm pathways, were significantly enriched. An mRNA-lncRNA interaction analysis suggested that the lncRNA LOC1056153S88 trans targets ARG2 and CCNB3, and the lncRNA LOC105607004 trans targets DMXL2, both of these might be involved in seasonal sheep breeding reproduction. Together, these results will provide resources for further studies on seasonal reproduction in sheep.

The objective of the current study was to analyze expression levels of synapse differentiation inducing 1-like (SYNDIG1L) and unc-13 homolog C (UNC13C) genes in different tissues, while single-nucleotide polymorphisms (SNPs) of two genes were associated with multiple thoracic vertebrae traits in both Small-tailed Han sheep (STH) and Sunite sheep (SNT). The expression levels of SYNDIG1L and UNC13C were analyzed in the brain, cerebellum, heart, liver, spleen, lung, kidney, adrenal gland, uterine horn, longissimus muscle, and abdominal adipose tissues of two sheep breeds with different thoracic vertebral number (TVN) sheep (T13 groups and T14 groups) by real-time quantitative polymerase chain reaction (RT-qPCR). Meanwhile, the polymorphisms of UNC13C gene g.52919279C>T and SYNDIG1L gene g.82573325C>A in T14 and T13 were genotyped by the Sequenom MassARRAY® SNP assay, and association analysis was performed with the TVN. The results demonstrated that UNC13C gene was extensively expressed in 11 tissues. The expression of UNC13C gene in longissimus muscle of T14 groups of STH was significantly higher than that of T13 groups (PA was significantly correlated with the TVN in both STH (P

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction, however, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. Therefore, in this study, we described the hypothalamic AS events and variants in differentially expressed genes (DEGs) in Small-tailed Han sheep without the FecB mutation at polytocous sheep in the follicular phase vs. monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase vs. monotocous sheep in the luteal phase (PL vs. ML) via an RNA-seq study for the first time. We found 39 DEGs with AS events (AS DEGs) in PF vs. MF, while 42 AS DEGs were identified in PL vs. ML. No DEGs with single nucleotide polymorphisms (SNPs) were observed in PF vs. MF, but five were identified in PL vs. ML. We also performed a correlation analysis of transcriptomics and proteomics, and the results suggested several key DEGs/differentially expressed proteins (DEPs), such as LGALS3 in PF vs. MF and aspartoacylase (ASPA) and transthyretin (TTR) in PL vs. ML, could be candidate genes influencing ovine litter size. In addition, further analyses suggested that AS events, SNPs and miRNA-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. All in all, this study provides a new insight into ovine and even other mammalian reproduction.

Follistatin‐like 3 (FSTL3) is a regulator of cellular apoptosis and was previously identified via RNA‐Seq to be associated with follicular development in mammalian ovaries. However, the mechanism underlying the FSTL3 regulation of oestrus in sheep remained poorly understood. In this study, the oestrogen (E2) and progesterone (P4) concentrations in blood were detected, and the expression level and functional analysis of FSTL3 in the ovary were studied during the different reproductive stage in Aohan fine wool sheep (seasonal breeding breed in China). The concentrations of E2 and P4 at the anestrus were significantly lower compared to dioestrus, proestrus and oestrus stages. Higher expression levels of FSTL3 were observed in the sheep ovary, hypothalamus, and thyroid. During different reproductive stages, higher expression levels were found during the stages of dioestrus and proestrus, while lower levels were found during the oestrus and anestrus stages. Functional analysis of FSTL3 was performed in primary granulosa cells (GCs) of sheep. The concentration of E2 increased significantly after RNAi interference of FSTL3, while the P4 level decreased. FSTL3 can decrease P4 levels, which might be involved in mediating oestrous cycle in sheep.

Simple Summary Cashmere goats are a unique husbandry resource in China. These goats are well known for producing the highest cashmere yield and best fiber quality in the world. Although cashmere is highly valued and also known as “fiber gem” and “soft gold”, few studies have examined the genetic basis of cashmere traits in cashmere goats. Here, we identified selection signals by comparing Fst and XP-EHH (the cross population extend haplotype homozygosity test) of a non-cashmere breed (Huanghuai goat) with those of two cashmere breeds (Inner Mongolia and Liaoning cashmere goats). Two genes (WNT10A and CSN3) were potentially associated with cashmere traits. This information may be valuable for studying the genetic uniqueness of cashmere goats and elucidating the mechanisms underlying cashmere traits in cashmere goats. Abstract Inner Mongolia and Liaoning cashmere goats in China are well-known for their cashmere quality and yield. Thus, they are great models for identifying genomic regions associated with cashmere traits. Herein, 53 Inner Mongolia cashmere goats, Liaoning cashmere goats and Huanghuai goats were genotyped, and 53,347 single-nucleotide polymorphisms (SNPs) were produced using the Illumina Caprine 50K SNP chip. Additionally, we identified some positively selected SNPs by analyzing Fst and XP-EHH. The top 5% of SNPs had selection signatures. After gene annotation, 222 and 173 candidate genes were identified in Inner Mongolia and Liaoning cashmere goats, respectively. Several genes were related to hair follicle development, such as TRPS1, WDR74, LRRC14, SPTLC3, IGF1R, PADI2, FOXP1, WNT10A and CSN3. Gene enrichment analysis of these cashmere trait-associated genes related 67 enriched signaling pathways that mainly participate in hair follicle development and stem cell pluripotency regulation. Furthermore, we identified 20 overlapping genes that were selected in both cashmere goat breeds. Among these overlapping genes, WNT10A and CSN3, which are associated with hair follicle development, are potentially involved in cashmere production. These findings may improve molecular breeding of cashmere goats in the future.

Tibetan goat is an ancient breed, which inhabits the adverse conditions of the plateaus in China. To investigate the role of selection in shaping its genomes, we genotyped Tibetan goats (Nagqu Prefecture, above 4500 m) and three lowland populations (Xinjiang goats, Taihang goats and Huanghuai goats). The result of PCA, neighbor-joining (N-J) tree and model-based clustering showed that the genetic structure between the Tibetan goat and the three lowland populations has significant difference. As demonstrated by the di statistic, we found that some genes were related to the high-altitude adaptation of Tibetan goats. Functional analysis revealed that these genes were enriched in the VEGF (vascular endothelial growth factor) signaling pathway and melanoma, suggesting that nine genes (FGF2, EGFR, AKT1, PTEN, MITF, ENPEP, SIRT6, KDR, and CDC42) might have important roles in the high-altitude adaptation of Nagqu Tibetan goats. We also found that the LEPR gene was under the strongest selection (di value = 16.70), and it could induce upregulation of the hypoxic ventilatory response. In addition, five genes (LEPR, LDB1, EGFR, NOX4 and FGF2) with high di values were analyzed using q-PCR. Among them, we found that LEPR, LDB1 and FGF2 exhibited higher expression in the lungs of the Tibetan goats, LEPR, EGFR and LDB1 exhibited higher expression in the hearts of the Huanghuai goat. Our results suggest that LEPR, LDB1, EGFR and FGF2 genes may be related to the high-altitude adaptation of the goats. These findings improve our understanding of the selection of the high-altitude adaptability of the Nagqu Tibetan goats and provide new theoretical knowledge for the conservation and utilization of germplasm resources.

Background: Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on increasing the litter size of this breed to improve the efficiency of mutton production.Results: We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH, oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 and ESR1. Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B.Conclusion: These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions of KDM4B that are responsible for litter size can in turn be used for further selection.

The objective of this study was to assess the genetic diversity and population structure of nine sheep populations (four famous mutton breeds worldwide, three popular mutton populations and two most famous high prolific breeds in China) and predict the heterosis among these populations based on microsatellite markers. Results showed the relatively rich genetic diversity in Chinese sheep populations. The results of UPGMA phylogenetic tree and Structure analysis indicated that the nine investigated populations can be divided into two groups. Suffolk and Dorset were clustered in one group, and the other group can be further divided into three clusters: GMM-BAS-BAM, HUS-STH and DOS-DOP. This is consistent with sheep breeding history. TreeMix analysis also hinted the possible gene flow from GMM to SUF. Additionally, based on the calculation of genetic distance, the optimal cross combinations are recommended for two popular mutton sheep populations (BAM and BAS) in the southern Mongolian Plateau and the most famous high prolific breeds (Small Tail Han sheep and Hu sheep) in China, respectively. Together, these results bring us an in-depth view of genetic diversity within sheep population and the genetic relationship among populations.

The pituitary is a vital endocrine organ that regulates animal seasonal reproduction by controlling the synthesis and secretion of the hormone. The change of photoperiod is the key factor affecting the function of the pituitary in animals, but the mechanism is unclear. Here, we studied the transcriptomic variation in pars distalis (PD) of the pituitary between short photoperiod (SP) and long photoperiod (LP) using RNA sequencing based on the OVX+E2 sheep. 346 differentially expressed (DE) lncRNAs and 186 DE-mRNA were found in the PD. Moreover, function annotation analysis indicated that the reproductive hormones and photoperiod response-related pathways including aldosterone synthesis and secretion, insulin secretion, thyroid hormone synthesis, and circadian entrainment were enriched. The interaction analysis of mRNA-lncRNA suggested that MSTRG.240648, MSTRG.85500, MSTRG.32448, and MSTRG.304959 targeted CREB3L1 and DUSP6, which may be involved in the photoperiodic regulation of the PD. These findings provide resources for further study on the seasonal reproductive in ewes.

Previous studies showed that the NLR family pyrin domain-containing 5 (NLRP5) and NLRP9 genes are two important reproductive genes, however, their effects on sheep litter size are unknown. Therefore, in this study, we first genotyped seven sheep breeds via the MassARRAY&reg, SNP system at the loci g.60495375A &gt, G, g.60495363G &gt, A, and g.60499690C &gt, A in NLRP5, and g.59030623T &gt, C and g.59043397A &gt, C in NLRP9. Our results revealed that each locus in most sheep breeds contained three genotypes. Then, we conducted population genetic analysis of single nucleotide polymorphisms in NLRP5 and NLRP9, and we found that the polymorphism information content value in all sheep breeds ranged from 0 to 0.36, and most sheep breeds were under Hardy&ndash, Weinberg equilibrium (p &gt, 0.05). Furthermore, association analysis in Small Tail Han sheep indicated that two loci, g.60495363G &gt, A in NLRP5 and g.59030623T &gt, C in NLRP9, were highly associated with litter size. The mutation in g.60495363G &gt, A may decrease interactions of NLRP5 with proteins, such as GDF9, whereas the mutation in g.59030623T &gt, C may enhance the combining capacity of NLRP9 with these proteins, consequently, these mutations may influence the ovulation rate and even litter size. The findings of our study provide valuable genetic markers that can be used to improve the breeding of sheep and even other mammals.

Photoperiod is one of the important factors leading to seasonal reproduction of sheep. However, the molecular mechanisms underlying the photoperiod regulation of seasonal reproduction remain poorly understood. In this study, we compared the expression profiles of mRNAs, lncRNAs, and circRNAs in uterine tissues from Sunite sheep during three different photoperiods, namely, the short photoperiod (SP), short transfer to long photoperiod (SLP), and long photoperiod (LP). The results showed that 298, 403, and 378 differentially expressed (DE) mRNAs, 171, 491, and 499 DE lncRNAs, and 124, 270, and 400 DE circRNAs were identified between SP and LP, between SP and SLP, and between LP and SLP, respectively. Furthermore, functional enrichment analysis showed that the differentially expressed RNAs were mainly involved in the GnRH signaling pathway, thyroid hormone synthesis, and thyroid hormone signaling pathway. In addition, co-expression networks of lncRNA&ndash, mRNA were constructed based on the correlation analysis between the differentially expressed RNAs. Our study provides new insights into the expression changes of RNAs in different photoperiods, which might contribute to understanding the molecular mechanisms of seasonal reproduction in sheep.

Simple Summary The NLR family pyrin domain-containing 5 (NLRP5) and NLRP9 genes are two important reproductive genes; however, their effects on litter size in sheep are unknown. In this study, we conducted population genetic and association analyses on five NLRP5 and NLRP9 loci of sheep. Our results suggested that a mutation in g.60495363G > A may decrease interactions of NLRP5 with proteins, such as the growth differentiation factor 9 (GDF9), whereas a mutation in g. g.59030623T > C may enhance the NLRP9-combining capacity with these proteins. Consequently, these mutations may lead to differences in ovulation rate and even litter size. Overall, this study provided useful genetic markers that can be used to improve sheep breeding. Abstract Previous studies showed that the NLR family pyrin domain-containing 5 (NLRP5) and NLRP9 genes are two important reproductive genes; however, their effects on sheep litter size are unknown. Therefore, in this study, we first genotyped seven sheep breeds via the MassARRAY® SNP system at the loci g.60495375A > G, g.60495363G > A, and g.60499690C > A in NLRP5, and g.59030623T > C and g.59043397A > C in NLRP9. Our results revealed that each locus in most sheep breeds contained three genotypes. Then, we conducted population genetic analysis of single nucleotide polymorphisms in NLRP5 and NLRP9, and we found that the polymorphism information content value in all sheep breeds ranged from 0 to 0.36, and most sheep breeds were under Hardy–Weinberg equilibrium (p > 0.05). Furthermore, association analysis in Small Tail Han sheep indicated that two loci, g.60495363G > A in NLRP5 and g.59030623T > C in NLRP9, were highly associated with litter size. The mutation in g.60495363G > A may decrease interactions of NLRP5 with proteins, such as GDF9, whereas the mutation in g.59030623T > C may enhance the combining capacity of NLRP9 with these proteins; consequently, these mutations may influence the ovulation rate and even litter size. The findings of our study provide valuable genetic markers that can be used to improve the breeding of sheep and even other mammals.

MicroRNA (miRNA) is a sort of endogenous ~20-25 nt non-coding RNAs, and it can regulate a variety of biological events. We found the miR-378 may involve in regulating the muscle development of sheep during our previous research. However, the molecular mechanism of miR-378 regulating myoblast proliferation is still unclear. In this research, we predicted that BMP2 (Bone morphogenetic protein 2) was the target gene of miR-378 and the BMP-Smad signal pathway that BMP2 participated in playing an important role in the muscle development. Therefore, we tried to determine whether miR-378 influence myoblast proliferation of sheep through the BMP-Smad signal pathway. The results indicated that inhibit BMP-Smad signal pathway by interfering Smad4 to promote proliferation of sheep myoblasts; promote BMP-Smad signal pathway by interfering Smad7 to inhibit proliferation of sheep myoblasts; over-expression miR-378 promotes BMP-Smad signal pathway and myoblast proliferation in sheep; interfering miR-378 inhibits BMP-Smad signal pathway and myoblast proliferation in sheep. However, when both of which functioned at the myoblast, miR-378 could not fully depend on BMP-Smad signal pathway to regulate myoblast proliferation. In sum, both miR-378 and BMP-Smad can influence the proliferation of myoblast, but miR-378 does not target the 3' UTR of sheep BMP2.

Inner Mongolia and Liaoning cashmere goats are two outstanding Chinese multipurpose breeds that adapt well to the semi-arid temperate grassland. These two breeds are characterized by their soft cashmere fibers, thus making them great models to identify genomic regions that are associated with cashmere fiber traits. Whole-genome sequencing of 70 cashmere goats produced more than 5.52 million single-nucleotide polymorphisms and 710,600 short insertions and deletions. Further analysis of these genetic variants showed some population-specific molecular markers for the two cashmere goat breeds that are otherwise phenotypically similar. By analyzing FST and θπ outlier values, we identified 135 genomic regions that were associated with cashmere fiber traits within the cashmere goat populations. These selected genomic regions contained genes, which are potential involved in the production of cashmere fiber, such as FGF5, SGK3, IGFBP7, OXTR, and ROCK1. Gene ontology enrichment analysis of identified short insertions and deletions also showed enrichment in keratinocyte differentiation and epidermal cell differentiation. These findings demonstrate that this genomic resource will facilitate the breeding of cashmere goat and other Capra species in future.

CircRNA plays important roles in cell proliferation, differentiation, autophagy and apoptosis during development. However, there are few reports on circRNAs related to livestock reproduction. In this study, we identified circRNAs by deep sequencing and analyzed their expression in the uteri of polytocous and monotocous sheep (FecB++) during follicular and luteal phases. There were 147 and 364 circRNAs with differential expression in the follicular and luteal phases, respectively. GO and KEGG enrichment analysis was performed for the host genes of the circRNAs to predict the functions of differentially expressed circRNAs. These source genes were mainly involved in the estrogen signaling pathway, TGF&beta, signaling pathway, GnRH signaling pathway, oxytocin signaling pathway, pentose phosphate pathway, and starch and sucrose metabolism related to reproduction and energy metabolism. CircRNA expression patterns were validated by RT-qPCR. Our findings provide a solid foundation for the identification and characterization of key important circRNAs involved in reproduction.