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.

Exome Sequencing Reveals Genetic Differentiation Due to High-Altitude Adaptation in the Tibetan Cashmere Goat ( Capra Hircus )

Sequence Characterization of DSG3 Gene to Know Its Role in High-Altitude Hypoxia Adaptation in the Chinese Cashmere Goat

Genomic Variation, Population History, and Long-Term Genetic Adaptation to High Altitudes in Tibetan Partridge (perdix Hodgsoniae).

Additional File 3: Table S2. of Exome Sequencing Reveals Genetic Differentiation Due to High-Altitude Adaptation in the Tibetan Cashmere Goat (capra Hircus)

Genomic Insights into Tibetan Sheep Adaptation to Different Altitude Environments

Markhor-derived introgression of a genomic region encompassing PAPSS2 confers high-altitude adaptability in Tibetan goats

Genome-wide comparative analyses reveal selection signatures underlying adaptation and production in Tibetan and Poll Dorset sheep

Genome-Wide Selection Signals Reveal Candidate Genes Associated with Plateau Adaptation in Tibetan Sheep

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

Genomic analysis identified a potential novel molecular mechanism for high-altitude adaptation in sheep at the Himalayas

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

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

Comparative Genome Analyses Reveal the Unique Genetic Composition and Selection Signals Underlying the Phenotypic Characteristics of Three Chinese Domestic Goat Breeds

Stress adaptation in Tibetan cashmere goats is governed by inherent metabolic differences and manifested through variable cashmere phenotypes

The Profiles and Functions of RNA Editing Sites Associated with High-Altitude Adaptation in Goats.

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

Variations in HIF-1α Contributed to High Altitude Hypoxia Adaptation via Affected Oxygen Metabolism in Tibetan Sheep

Expression and Variations in EPO Associated with Oxygen Metabolism in Tibetan Sheep.

LncRNA and Protein Expression Profiles Reveal Heart Adaptation to High-Altitude Hypoxia in Tibetan Sheep

Recent Selection and Introgression Facilitated High-Altitude Adaptation in Cattle

Physiology and Transcriptomics Analysis Reveal the Contribution of Lungs on High-Altitude Hypoxia Adaptation in Tibetan Sheep