The dairy goat is one of the earliest dairy livestock species, which plays an important role in the economic development, especially for developing countries. With the development of agricultural civilization, dairy goats have been widely distributed across the world. However, few studies have been conducted on the specific characteristics of dairy goat. In this study, we collected the whole-genome data of 89 goat individuals by sequencing 48 goats and employing 41 publicly available goats, including five dairy goat breeds (Saanen, Nubian, Alpine, Toggenburg, and Guanzhong dairy goat; n = 24, 15, 11, 6, 6), and three goat breeds (Guishan goat, Longlin goat, Yunshang Black goat; n = 6, 15, 6). Through compared the genomes of dairy goat and non-dairy goat to analyze genetic diversity and selection characteristics of dairy goat. The results show that the eight goats could be divided into three subgroups of European, African, and Chinese indigenous goat populations, and we also found that Australian Nubian, Toggenburg, and Australian Alpine had the highest linkage disequilibrium, the lowest level of nucleotide diversity, and a higher inbreeding coefficient, indicating that they were strongly artificially selected. In addition, we identified several candidate genes related to the specificity of dairy goat, particularly genes associated with milk production traits (GHR, DGAT2, ELF5, GLYCAM1, ACSBG2, ACSS2), reproduction traits (TSHR, TSHB, PTGS2, ESR2), immunity traits (JAK1, POU2F2, LRRC66). Our results provide not only insights into the evolutionary history and breed characteristics of dairy goat, but also valuable information for the implementation and improvement of dairy goat cross breeding program.

2 SNPs were discovered in our previous genome-wide association study (GWAS): s58995.1 (rs420767326 A>G) in MEF2B gene and OAR3_115712045.1 (rs401775061 A>C) in TRHDE gene, which were significantly associated with post-weaning gain in sheep. Herein, we performed a replication experiment to investigate single nucleotide polymorphisms (SNPs) within the MEF2B and TRHDE gene exons, the 5′untranslated regions (within 1000bp), the 3′ untranslated regions (within 1000bp) and their associations with Ujumqin sheep growth traits in 4-month age and 6-month age, respectively. Finally,3 SNPs were selected to be investigated including 1 SNP in 3′untranslated regions in MEF2B gene (rs417014745 A>G) and 2 SNPs in TRHDE gene (rs426980328 T>C and rs430810656 G>A).The χ2 test showed all the 3 variations were in Hardy–Weinberg equilibrium (P>0.05) status. Association analysis suggested that rs426980328 T>C was significantly associated with body weight and chest girth in 4-month age (PA exhibited extremely significant association with body weight and chest girth in 4-month age (PG was extremely significantly associated with body weight and chest girth in 4-month age and chest girth in 6-month age (P<0.01), and it was also significantly associated with body weight in 6-month age (P<0.05). Combined effect analysis indicated significant associations between the combinations of rs426980328-rs417014745, rs430810656-rs417014745 and several growth traits (P<0.05). These results suggested MEF2B and TRHDE genes affected growth traits in Ujumqin sheep and the combination effect of the two genes also played a significant effective role. These SNPs might have potential value as genetic markers for growth traits and it could be used in Ujumqin sheep breeding in future. Further studies are necessary to confirm our findings.

Long non-coding RNA (lncRNA) has been extensively studied in many livestock. However, compared with other livestock breeds, there is less research regarding donkey lncRNA function. It has been reported that lncRNA plays an important role in the timing control of development, aging, and death of livestock. Therefore, the study of donkey skeletal muscle lncRNA is of great significance for exploring donkey meat production performance. In this study, RNA-Seq was used to perform high-throughput sequencing of skeletal muscle (longissimus dorsi and gluteus) of two Dezhou donkey strains (SanFen and WuTou). The differentially expressed lncRNAs were screened between different strains and tissues. Then candidate genes for conjoint analysis were screened based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, the accuracy of the sequencing data was verified by real-time quantitative polymerase chain reaction (RT-qPCR). Herein, we identified 3869 novel lncRNAs and 73 differentially expressed lncRNAs. Through the comparison between groups, the specific expression of lncRNAs were found in different strains and muscle tissues. Importantly, we constructed the lncRNA-miRNA-mRNA interaction network and found three important candidate lncRNAs (MSTRG.9787.1, MSTRG.3144.1, and MSTRG.9886.1) and four candidate genes (ACTN1, CDON, FMOD, and BMPR1B). Analysis of the KEGG pathway indicates that the TGF-β signaling pathway plays a pivotal role in the growth and development of skeletal muscle in Dezhou donkey strains.

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.

The growth and development of embryonic skeletal muscle plays a crucial role in sheep muscle mass. But proteomic analyses for embryonic skeletal development in sheep had been little involved in the past research. In this study, we explored differential abundance proteins during embryonic skeletal muscle development by the tandem mass tags (TMT) and performed a protein profile analyses in the longissimus dorsi of Chinese merino sheep at embryonic ages Day85 (D85N), Day105 (D105N) and Day135 (D135N). 5,520 proteins in sheep embryonic skeletal muscle were identified, and 1,316 of them were differential abundance (fold change ≥1.5 and p-value < 0.05). After the KEGG enrichment analyses, these differential abundance proteins were significant enriched in the protein binding, muscle contraction and energy metabolism pathways. After validation of the protein quantification with the parallel reaction monitoring (PRM), 41% (16/39) significant abundance proteins were validated, which was similar to the results of protein quantification with TMT. All results indicated that D85N to D105N was the stage of embryonic muscle fibers proliferation, while D105N to D135N was the stage of their hypertrophy. These findings provided a deeper understanding of the function and rules of proteins in different phases of sheep embryonic skeletal muscle growth and development.

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.