Hu sheep is a renowned prolific local sheep breed in China, widely distributed across the country due to its excellent reproductive performance. Deciphering the molecular mechanisms underlying the high fecundity of Hu sheep is crucial for improving the litter size of ewes. In this study, we genotyped 830 female Hu sheep using the Illumina OvineSNP50 BeadChip and performed genetic diversity analysis, selection signature detection, and a genome-wide association study (GWAS) for litter size. Our results revealed that the Hu sheep population exhibits relatively high genetic diversity. A total of 4927 runs of homozygosity (ROH) segments were detected, with the majority (74.73%) being short in length. Different genomic inbreeding coefficients (F_ROH, F_HOM, F_GRM, and F_UNI) ranged from −0.0060 to 0.0126, showing low levels of inbreeding in this population. Additionally, we identified 91 candidate genomic regions through three complementary selection signature methods, including ROH, composite likelihood ratio (CLR), and integrated haplotype score (iHS), and annotated 189 protein-coding genes. Moreover, we observed two significant SNPs related to the litter size of Hu sheep using GWAS analysis based on a repeatability model. Integrating the selection signatures and the GWAS results, we identified 15 candidate genes associated with litter size, among which BMPR1B and UNC5C were particularly noteworthy. These findings provide valuable insights for improving the reproductive performance and breeding of high-fecundity lines of Hu sheep.

Muscle development is a complex biological process involving an intricate network of multiple factor interactions. Through the analysis of transcriptome data and molecular biology confirmation, this study aims to reveal the molecular mechanism underlying sheep embryonic skeletal muscle development. The RNA sequencing of embryos was conducted, and microRNA (miRNA)-mediated competitive endogenous RNA (ceRNA) networks were constructed. qRT-PCR, siRNA knockdown, CCK-8 assay, scratch assay, and dual luciferase assay were used to carry out gene function identification. Through the analysis of the ceRNA networks, three miRNAs (miR-493-3p, miR-3959-3p, and miR-410-5p) and three genes (TEAD1, ZBTB34, and POGLUT1) were identified. The qRT-PCR of the DE-miRNAs and genes in the muscle tissues of sheep showed that the expression levels of the TEAD1 gene and miR-410-5p were correlated with the growth rate. The knockdown of the TEAD1 gene by siRNA could significantly inhibit the proliferation of sheep primary embryonic myoblasts, and the expression levels of SLC1A5, FoxO3, MyoD, and Pax7 were significantly downregulated. The targeting relationship between miR-410-5p and the TEAD1 gene was validated by a dual luciferase assay, and miR-410-5p can significantly downregulate the expression of TEAD1 in sheep primary embryonic myoblasts. We proved the regulatory relationship between miR-410-5p and the TEAD1 gene, which was related to the proliferation of sheep embryonic myoblasts. The results provide a reference and molecular basis for understanding the molecular mechanism of embryonic muscle development.

A genomic study was conducted to uncover the selection signatures in sheep that show extremely significant differences in growth traits under the same breed, age in months, nutrition level, and management practices. Hu sheep from Gansu Province and Gangba sheep from the Tibet Autonomous Region in China were selected. We collected whole-genome data from 40 sheep individuals (24 Hu sheep and 16 Gangba sheep), through whole-genome sequencing. Selection signals were analyzed using parameters such as F_ST, π _ratio, and Tajima’s D. We have identified several candidate genes that have undergone strong selection, particularly those associated with growth traits. Specifically, five growth-related genes were identified in both the Hu sheep group (HDAC1, MYH7B, LCK, ACVR1, GNAI2) and the Gangba sheep group (RBBP8, ACSL3, FBXW11, PLAT, CRB1). Additionally, in a genomic region strongly selected in both the Hu and Gangba sheep groups (Chr 22: 51,425,001-51,500,000), the growth-associated gene CYP2E1 was identified, further highlighting the genetic factors influencing growth characteristics in these breeds. This study analyzes the genetic basis for significant differences in sheep phenotypes, identifies candidate genes related to sheep growth traits, lays the foundation for molecular genetic breeding in sheep, and accelerates the genetic improvement in livestock.

As relatively new members of the non-coding RNA family, circRNAs play important roles in a variety of biological processes. However, the temporal expression pattern and the function of circRNAs during sheep skeletal muscle development remains unclear. This study aimed to identify circRNAs related to sheep skeletal muscle development and explore their roles in myoblast proliferation. The circRNA expression profiles of longissimus dorsi of sheep from F90, L30, and A3Y were obtained by the RNA-seq method. The function and mechanisms of the novel circCHRNG in muscle satellite cell proliferation were explored using CCK-8 assay, Western blot, qPCR, and dual-luciferase reporter assay. We identified 12,375 circRNAs, including 476, 133, and 233 DEcircRNAs found among three comparative groups. KEGG results showed that DEcircRNAs were enriched in muscle contraction, the regulation of cell proliferation, and the AMPK, insulin, and PI3K-Akt signaling pathways. Notably, a novel circRNA, termed circRNA CHRNG, acts as a miR-133 sponge to promote skeletal muscle satellite cell proliferation. Our study provides a systematic description of circRNAs of ovine skeletal muscle across fetal, lamb, and adult stages. GO and KEGG analyses showed that DEcircRNAs were enriched in multiple pathways associated with muscle development, such as the PI3K-Akt and AMPK signaling pathways. In addition, we propose that circCHRNG acts as a miR-133 sponge to upregulate the expression levels of SRF and MEF2A, thereby promoting myoblast proliferation.

Growth traits are influential factors that significantly affects the development of the sheep industry. A previous TMT proteomic analysis found that a key protein in the HIF signaling pathway, ARNT, may influence embryonic skeletal muscle growth and development in sheep. The purpose of this study was to better understand the association between the polymorphisms of ARNT and growth traits of sheep, and the potential function of ARNT. Real-time qPCR (qRT-PCR) of ARNT was carried out to compare its expression in different developmental stages of the muscle tissues and primary myoblasts in the Hu, Chinese merino, and Gangba sheep. The genetic variance of ARNT was detected using the Illumina Ovine SNP 50 K and 600 K BeadChip in the Hu and Ujimqin sheep populations, respectively. The CDS sequence of the ARNT gene was cloned in the Hu sheep using PCR technology. Finally, bioinformatic analytical methods were applied to characterize the genes and their hypothetical protein products. The qRT-PCR results showed that the ARNT gene was expressed significantly in the Chinese merino embryo after 85 gestation days (D85) (p < 0.05). Additionally, after the sheep were born, the expression of ARNT was significant at the weaning stage of the Hu sheep (p < 0.01). However, there was no difference in the Gangba sheep.In addition, six SNP loci were screened using 50 K and 600 K BeadChip. We found a significant association between rs413597480 A > G and the Hu sheep weight at weaning and backfat thickness in the 5-month-old sheep (p < 0.05), and four SNP loci (rs162298018 G > C, rs159644025 G > A, rs421351865 G > A, and rs401758103 A > G) were also associated with growth traits in the Ujimqin sheep (p < 0.05). Interestingly, we found that a G > C mutation at 1948 bp in the cloned ARNT CDS sequence of the Hu sheep was the same locus mutation as rs162298018 G > C identified using the 600 K BeadChip, which resulted in a nonconservative missense point mutation, leading to a change from proline to alanine and altering the number of DNA, protein-binding sites, and the α-helix of the ARNT protein. There was a strong linkage disequilibrium between rs162298018 G > C and rs159644025 G > A, and the ARNT protein was conserved among the goat, Hu sheep, and Texel sheep. And, we propose that a putative molecular marker for growth and development in sheep may be the G > C mutation at 1948 bp in the CDS region of the ARNT gene. Our study systematically analyzed the expression, structure, and function of the ARNT gene and its encoded proteins in sheep. This provides a basis for future studies of the regulatory mechanisms of the ARNT gene.