Abstract:Background Traditionally, Chinese indigenous sheep were classified geographically and morphologically into three groups: Mongolian, Kazakh and Tibetan. Herein, we aimed to evaluate the population structure and genome selection among 140 individuals from ten representative Chinese indigenous sheep breeds: Ujimqin, Hu, Tong, Large-Tailed Han and Lop breed (Mongolian group); Duolang and Kazakh (Kazakh group); and Diqing, Plateau-type Tibetan, and Valley-type Tibetan breed (Tibetan group). Results We analyzed the population using principal component analysis (PCA), STRUCTURE and a Neighbor-Joining (NJ)-tree. In PCA plot, the Tibetan and Mongolian groups were clustered as expected; however, Duolang and Kazakh (Kazakh group) were segregated. STRUCTURE analyses suggested two subpopulations: one from North China (Kazakh and Mongolian groups) and the other from the Southwest (Tibetan group). In the NJ-tree, the Tibetan group formed an independent branch and the Kazakh and Mongolian groups were mixed. We then used the d i statistic approach to reveal selection in Chinese indigenous sheep breeds. Among the 599 genome sequence windows analyzed, sixteen (2.7%) exhibited signatures of selection in four or more breeds. We detected three strong selection windows involving three functional genes: RXFP2 , PPP1CC and PDGFD. PDGFD , one of the four subfamilies of PDGF , which promotes proliferation and inhibits differentiation of preadipocytes, was significantly selected in fat type breeds by the Rsb (across pairs of populations) approach. Two consecutive selection regions in Duolang sheep were obviously different to other breeds. One region was in OAR2 including three genes ( NPR2 , SPAG8 and HINT2 ) the influence growth traits. The other region was in OAR 6 including four genes ( PKD2 , SPP1 , MEPE , and IBSP ) associated with a milk production quantitative trait locus. We also identified known candidate genes such as BMPR1B , MSRB3 , and three genes ( KIT , MC1R , and FRY ) that influence lambing percentage, ear size and coat phenotypes, respectively. Conclusions Based on the results presented here, we propose that Chinese native sheep can be divided into two genetic groups: the thin type (Tibetan group), and the fat type (Mongolian and Kazakh group). We also identified important genes that drive valuable phenotypes in Chinese indigenous sheep, especially PDGFD , which may influence fat deposition in fat type sheep.

Integration analysis of transcriptome and metabolome revealed the potential mechanism of spermatogenesis in Tibetan sheep (Ovis aries) at extreme high altitude

Transcriptome and metabolome revealed the effects of hypoxic environment on ovarian development of Tibetan sheep

Transcriptomic study of testicular hypoxia adaptation in Tibetan sheep

Genomic Insights into Tibetan Sheep Adaptation to Different Altitude Environments

Transcriptomic Study of Different Stages of Development in the Testis of Sheep

Expression and role of melatonin membrane receptors in the hypothalamic-pituitary-testicular axis of Tibetan sheep in a plateau pastoral area

Identification of Key Genes Related to Postnatal Testicular Development Based on Transcriptomic Data of Testis in Hu Sheep

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Transcriptomic Analysis of the HPG Axis-Related Tissues Reveals Potential Candidate Genes and Regulatory Pathways Associated with Testicular Size in Hu Sheep

Morphological and Scanning Electron Microscopic Study of the Gonadal Arterioles in the Tibetan Sheep.

iTRAQ-based proteomic analysis provides novel insight into the postnatal testicular development of Hu sheep

Multi-omic Analyses Shed Light on The Genetic Control of High-altitude Adaptation in Sheep

Single‐cell RNA sequencing of the Mongolia sheep testis reveals a conserved and divergent transcriptome landscape of mammalian spermatogenesis

Differential tissue expression of sex steroid-synthesizing enzyme CYP11A1 in male Tibetan sheep (Ovis aries)

Integrated Multi-Tissue Transcriptome Profiling Characterizes the Genetic Basis and Biomarkers Affecting Reproduction in Sheep (Ovis aries)

Genetic Signatures of High-Altitude Adaptation and Geographic Distribution in Tibetan Sheep

High-Altitude Stress Orchestrates Mrna Expression and Alternative Splicing of Ovarian Follicle Development Genes in Tibetan Sheep

Exploring the Molecular Characteristics and Role of PDGFB in Testis and Epididymis Development of Tibetan Sheep

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Genome-Wide Analysis Reveals Adaptation to High Altitudes in Tibetan Sheep

Response of Ruminal Microbiota-Host Gene Interaction to High-Altitude Environments in Tibetan Sheep