Journal cover Journal topic
Archives Animal Breeding Archiv Tierzucht
Journal topic
AAB | Articles | Volume 63, issue 1
Arch. Anim. Breed., 63, 193–201, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
Arch. Anim. Breed., 63, 193–201, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Original study 26 Jun 2020

Original study | 26 Jun 2020

Genetic diversity and genetic origin of Lanping black-boned sheep investigated by genome-wide single-nucleotide polymorphisms (SNPs)

Genetic diversity and genetic origin of Lanping black-boned sheep investigated by genome-wide single-nucleotide polymorphisms (SNPs)
Heli Xiong1,2, Xiaoming He1, Jing Li1,3, Xingneng Liu1, Chaochao Peng1, Dongmei Xi1, and Weidong Deng1 Heli Xiong et al.
  • 1Yunnan Provincial Key Laboratory of Animal Nutrition and Feed, Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, People's Republic of China
  • 2Yunnan Animal Science and Veterinary Institute, Kunming 650224, People's Republic of China
  • 3Yunnan Kunming Police Dog Base of Ministry of Public Security, Kunming 650201, People's Republic of China

Correspondence: Weidong Deng (

Back to toptop

Lanping black-boned sheep was first discovered in the 1950s in Lanping county of China and characterized by black pigmentation on skin and internal organs. Due to the novel and unique trait, the genetic background of Lanping black-boned sheep is of great interest. Here, we genotyped genome-wide SNPs (single nucleotide polymorphisms) of Lanping black-boned sheep and Lanping normal sheep using Illumina OvineSNP50 BeadChip to investigate the genetic diversity and genetic origin of Lanping black-boned sheep. We also downloaded a subset SNP dataset of two Tibet-lineage sheep breeds and four other sheep breeds from the International Sheep Genomics Consortium (ISGC) as a reference for interpreting. Lanping black-boned sheep had a lower genetic diversity level when compared to seven other sheep breeds. Principal component analysis (PCA) showed that Lanping black-boned sheep and Lanping normal sheep were clustered into the Asian group, but there was no clear separation between the two breeds. Structure analysis demonstrated a high ancestry coefficient in Lanping black-boned sheep and Lanping normal sheep. However, the two populations were separated into two distinct branches in a neighbor-joining (NJ) tree. We further evaluated the genetic divergence using population FST, which showed that the genetic differentiation that existed between Lanping black-boned sheep and Lanping normal sheep was higher than that between Tibet sheep and Changthangi sheep, which revealed that Lanping black-boned sheep is a different breed from Lanping normal sheep on the genetic level. In addition, structure analysis and NJ tree showed that Lanping black-boned sheep had a relatively close relation with Tibet sheep. The results reported herein are a first step toward understanding the genetic background of Lanping black-boned sheep, and it will provide informative knowledge on the unique genetic resource conservation and mechanism of novel breed formation.

1 Introduction
Back to toptop

Lanping black-boned sheep (LPBB) was first discovered in the 1950s and characterized by black pigmentation on skin and internal organs compared to the reddish coloration in Lanping normal sheep (LPN) (Deng et al., 2008; Li et al., 2018); the pigmentation pattern is similar to that in silky fowl of China (Li et al., 2018; Deng et al., 2006). Attributable to their tender and tasty meat quality and no special smell of mutton, LPBB stands out among Lanping local sheep breeds and has been designated as a novel genetic resource by the Chinese Ministry of Agriculture (Li et al., 2018).

LPBB has the same morphology with LPN such as coat color and horns and tail shape; there are slight phenotypic differences between the two breeds seen by careful visible inspection: LPBB has dark teeth and gums and visible mucosa. However, there are great differences between the two breeds: the skin, muscle, and inner organs – including kidney, heart, lung, and trachea of LPBB – were dark colored compared to the red coloration in LPN (Deng et al., 2006). The husbandry system of LPBB is similar to semi-feral sheep breeds of Tibet sheep, which graze freely with more natural selection and less human intervention such as selective breeding and nutrition complements (Deng et al., 2006; Pan et al., 2018). Due to the phenotypic similarity with LPN and less human intervention, it is supposed by scientists that the formation of LPBB was the result of advantageous genetic mutations from LPN in the process of adaptation to the local harsh environments, which include high altitude (approximately 3000 m), strong radiation, steep terrain, and low temperatures during winter (Li, 2009). Meanwhile, some scholars believe that gene flow is present in Lanping local sheep and Tibet sheep because Lanping is located on the Tea Horse Road leading to Tibet in Yunnan and is adjacent to Diqing Tibetan Autonomous Prefecture (Li, 2009). Furthermore, previous studies reported that Yunnan-Kweichow Plateau sheep has closer relation with Qinghai-Tibetan Plateau sheep compared to northern Chinese sheep (Wei et al., 2015; Yang et al., 2016; Hu et al., 2019). Thus, investigating the genetic relation of LPBB with these sheep breeds may provide informative knowledge on genetic origin of LPBB.

Genetic diversity is an important population genetics parameter that helps to explain the process of evolution (Notter, 1999) and is a major concern, considering the necessity of conservation in local breeds (Meloni et al., 2015). The ability of a population to respond adaptively to environmental changes depends on its level of genetic diversity, and a species without enough genetic diversity is thought to be unable to cope with changing environments or evolving competitors and parasites (Khodabakhshzadeh et al., 2016). Lanping black-boned sheep was found initially with about 200 individuals and the population number increased slowly to 2000 from 2001 to 2005 (Deng et al., 2006). Due to the small population and isolated habitat, LPBB was expected to have low genetic diversity and increased incidence of inbreeding, which was not good for this unique and rare sheep breed. Thus, investigating genetic diversity of LPBB will provide a guide for the strategies on this genetic resource conservation.

To date, genes related to the mechanism of black pigmentation traits on LPBB have been studied such as TYR (Deng et al., 2008), MC1R (Deng et al., 2009a), TYRP1, TYRP2 (Deng et al., 2009b), and EDN3 (Darwish et al., 2018). However, the genetic background of LPBB has not yet been reported. In this study, we aimed to investigate the genetic diversity and genetic origin of Lanping black-boned sheep by genotyping LPBB and LPN using Illumina OvineSNP50 BeadChip. We also downloaded a subset of publicly available SNP (single-nucleotide polymorphism) BeadChip data (Kijas et al., 2012a) from the Ovine HapMap project of International Sheep Genomics Consortium (ISGC) (Kijas et al., 2012a) to use this established resource as a reference for interpreting genetic diversity and phylogenetic relationship within Lanping local breeds and for developing assumption about their relationship with worldwide breeds (Beynon et al., 2015). The results obtained in this study will be informative for future research and conservation of this unique breed.

2 Materials and methods
Back to toptop

2.1 Animal materials

We collected 101 LPBB (48 males and 53 females) and 106 (19 males and 87 females) LPN sheep blood samples from Lanping County, Yunnan Province, China. Each individual was more than 1 year old and also carefully confirmed to match the breed standard. We then selected 15 LPBB (7 males and 8 females) and 15 LPN (7 males and 8 females) individuals randomly from the above collected blood samples to genotype; one LPN blood sample was lost for some reason. Genomic DNA was extracted from 500 µL of whole blood using Takara Blood Genome DNA Extraction Kit following the manufacture's guidelines.

All the experimental procedures were conducted under the International Animal Care and Use Committee of the Yunnan Agricultural University. The care and use of animals fully complied with local animal welfare laws, guidelines, and policies.

2.2 Genotyping and data quality control

In total, 29 genomic DNA samples were genotyped using Illumina OvineSNP50 BeadChip. The second dataset used in this study, including two Tibetan-lineage breeds (Tibetan and Changthangi sheep), another two Asian breeds (n=90), one southwest Asian (n=20) breed, and two southern and western European breeds (SW Europe) (n=43), was downloaded from ISGC and is available at (last access: 24 June 2020). Then, we merged and generated a 182-individual dataset (Table 1) with 49 034 overlapping SNPs for subsequent analysis. SNPs that cannot pass any of the following criteria were excluded from the analysis using PLINK 1.9 (Brito et al., 2017a; Johnston et al., 2016; Purcell et al., 2007): (1) with minor allele frequency (MAF) > 0.01; (2) with missing genotype data < 0.10; (3) individuals with missing genotype data  < 10 %; (4) with Hardy–Weinberg equilibrium P value > 0.00001; (5) included in the latest reference assembly of the sheep genome Oar_v4.1; and (6) located on autosomes. After filtering, a total of 45 754 autosomal SNPs were selected for further analyses. Individuals with potential kinships will also be excluded; we estimated the value of identity by state (IBS) between all samples using the “cluster – matrix” flag in PLINK to identify the genetic relatedness among sheep individuals (Gorkhali et al., 2016; Johnston et al., 2016).

Table 1Summary of animal resources and genetic diversity of eight sheep populations.

Pn represents the proportion of SNP displaying polymorphism, He represents average expected heterozygosity, Ho represents average observed heterozygosity, and F represents inbreeding coefficient.

Download Print Version | Download XLSX

2.3 Genetic diversity

Five metrics were used to evaluate levels of within-breed genetic diversity, and the values were estimated using PLINK 1.9 (Purcell et al., 2007). The proportion of polymorphic SNP (Pn) gives the fraction of total SNP that displayed both alleles within each population (Brito et al., 2017a). Pn was calculated as the proportion of SNPs with average minor allele frequency (MAF) greater than 1 % within each breed (Brito et al., 2017a). MAF is the frequency estimate of the least common allele per breed and estimated using the “freq” flag. Expected heterozygosity (He) and observed heterozygosity (Ho) were estimated using the “hardy” flag, and inbreeding coefficient (F) was estimated using the “het” flag (Beynon et al., 2015; Brito et al., 2017a).

2.4 Phylogenetic analysis

A pruned dataset of 182 sheep containing 36 711 SNPs, which excluded SNPs in LD (linkage disequilibrium) (PLINK – indep-pairwise 50 10 0.2; – extract, was used to investigate the genetic structure. Principal component analysis (PCA) was performed with GCTA software (version 1.26.0, genome-wide complex trait analysis) (Yang et al., 2011). Population structure was evaluated using Structure software (version 2.3.4) (Falush et al., 2007). The neighbor-joining (NJ) tree was constructed using MEGA software (version 7.0) on the basis of all SNPs (Hall, 2013). The pairwise r2 value within each population was calculated with parameter “– r2 – ld-window 99999 – ld-window-r2 0.2” in PLINK to compare LD patterns among breeds. Population divergence was calculated as FST using VCFtools (version 4.2) (Danecek et al., 2011).

3 Results
Back to toptop

3.1 Genetic diversity

Inbred individuals were not observed in the 182 studied sheep according to the IBS score (IBS < 0.9) (Yang et al., 2016). The level of Pn present in LPBB population was 0.8390, which was slightly lower than LPN population but greatly lower than the other six populations; this reflected that LPBB population had lower loci displaying polymorphism compared to the remaining populations (Table 1). The distribution of MAF per breed is given in Fig. 1 and shows that the LPBB population has the lowest level of highly variable SNPs (MAF > 0.3) with a proportion of 33.88 %, and the highest level is CME with a proportion of 47.41 %; LPBB population had an excess of low allele frequency SNPs (MAF = 0) compared with the other breeds. The value of He was close to Ho in all populations and the LPBB population had the lowest He (0.2861), followed by LPN population (0.2870); compared to pure breeds, the cross breed of CME had the highest level of He (0.3551). The inbreeding coefficients in all populations were detected with weak to moderate levels, ranging from 0.0267 to 0.2835 and the level in the LPBB population was moderate (0.1847).

Figure 1Minor allele frequency (MAF) of eight sheep breeds.


3.2 Phylogenetic analyses

We first performed principle component analysis (PCA) based on the pruned genotype data of 36 711 SNPs among all individuals to examine the genetic relationship between LPBB and its geographic neighbors and worldwide sheep breeds. The largest PC (6.931 % of total variation) separated eight breeds into three groups consistent with their genetic origin which were Asian breeds (PC1 < 0, LPBB, LPN, TIB, CHA and BGE), southwest Asian (0 < PC1 < 0.05, NDZ), and SW European breeds (PC1 > 0.05, CME, SAB). The second PC (4.467 % of total variation) separated Asian sheep breeds into three clusters consistent with their geographical distance: Lanping local breeds (PC2 <−0.05, LPBB and LPN), Tibet-lineage breeds (−0.05< PC2 < 0.05, TIB and CHA), and BGE, which is an Asian breed far away from the first two groups (PC2 > 0.05). Furthermore, LPBB and LPN populations showed relatively close relation with TIB and CHA populations. Based on PC1 and PC2, LPBB individuals clustered more tightly than LPN individuals; however, there was not a clear separation between LPBB and LPN populations. In addition, five individuals of CHA were outside of their expected population clusters (Fig. 2).

Figure 2Principal components 1 and 2 for the 182 sheep individuals.


Figure 3Population structure of 182 sheep individuals inferred from Structure 2.3.4 software. The length of each colored segment represents the proportion of the individual genome inferred from ancestral populations (K=2–8).


Figure 4The NJ tree was constructed using MEGA 7.0 and edited by iTOL.


To investigate ancestry and admixture proportion of LPBB breed, we performed population structure analysis based on pruned SNPs. Model-based clustering of the individuals was analyzed by assuming numbers of populations (K) between 1 and 8 (Beynon et al., 2015). Figure 3 shows the distribution of individuals into clusters for K=2–8. The clustering analysis for K=2 shows that 182 sheep were genetically divided into Asian and SW European breeds; NDZ population had an admixture component of Asian breeds and SW European breeds. When K=3, BGE breed was separated from Asian breeds, which was consistent with the geographical origin and the PCA analysis that PC2 separated BGE breeds from the remaining Asian breeds. When the K value became large, most breeds tended to be separated except LPBB and LPN. The ideal clustering of eight populations determined by the ΔK method (Nie et al., 2016) (Fig. 5) was K=7. LPBB and LPN breeds were clustered into one population at K=7; this revealed that LPBB and LPN breeds had a very similar genetic background, which was consistent with the PCA that LPBB and LPN were mixed. Figure 4 shows that at K=7 LPBB and LPN sheep are represented by one main cluster and have some features that are present in the TIB cluster. Although LPBB and LPN were clustered together, there are differences between the two populations: the admixture level of LPBB was slightly higher than that of LPN, LPBB had more features that presented in TIB clusters than LPN the representative clusters in CHA were present in LPN with a low ancestry coefficient. Meanwhile, a similar pattern was observed in TIB and CHA clusters; the representative clusters in LPBB and LPN were also present in TIB and CHA clusters with a low ancestry coefficient. This result demonstrated that LPBB and LPN breeds had shared ancestry with TIB breed or there was a gene flow between them.

Figure 5Optimal value of K determined by δK (ΔK).


To further resolve the phylogenetic relationship of LPBB with their geographical neighbors and worldwide sheep breeds, we constructed a neighbor-joining tree based on the whole genome-wide SNPs after filtering. Agreeing with PCA analysis, the phylogenetic tree split all breeds into three distinct branches which represent Asian, southwest Asian, and SW European sheep breeds (Fig. 3) and five individuals of CHA were outside of their clusters. The phylogenetic tree shows that LPBB and LPN breeds are separated into two clear clusters and positioned in one clade with TIB and CHA breed. Furthermore, LPBB and LPN breeds had the shortest branch length compared to the TIB breed, which revealed that LPBB and LPN had a relatively close relationship with TIB. The fact that they had common ancestry or gene flow between them was confirmed by the structure analysis.

Due to the morphology and genetic similarity of LPBB and LPN breeds, it was questionable that LPBB and LPN were distinct enough to be considered a different breed. To investigate their genetic difference, we used populations FST to evaluate the degree of genetic divergence between LPBB and LPN sheep breeds and compared them against the divergence that exists between populations recognized as separate breeds. The pairwise population FST of eight breeds was calculated using VCFtools. The FST value can range from zero (no genetic differentiation between the populations or from the ancestral population) to one (fully differentiated in each population). FST values ranging from 0.05 to 0.15 and 0.15 to 0.25 mean moderate and large differentiations, respectively (Grasso et al., 2014). Table 2 shows the FST that existed between each breed pair. The FST value of all pairwise was small to moderate, which was consistent with the report by Kijas et al. (2009). They studied 74 breeds from different parts of the world and showed that sheep breeds had maintained generally low genetic differentiation. The lowest level of divergence was found between TIB and CHA breeds (0.0285), which are both Tibet-lineage sheep breeds, while the greatest was observed between LPBB and SAB (0.1584), which are Asian and SW European breeds, respectively; this result was supported by the structure analysis result that they were from completely distinct clusters. Moderate genetic divergence was observed (0.0646) between LPBB and LPN breeds, which was higher than that between the two Tibet-lineage separate breeds. This strongly indicated that LPBB breed is sufficiently different with the LPN breed and can be considered a separate breed. The divergence level of LPBB-TIB and LPN-TIB was relatively lower with FST values of 0.0750 and 0.0759, respectively, which agreed with the NJ tree analysis that LPBB and LPN were more related to TIB than to the remaining breeds.

Table 2Population divergence measured as FST.

* FST is given above the diagonal and its standard deviation (SD) for each combination is given below.

Download Print Version | Download XLSX

To better understand the population variation, we performed linkage disequilibrium (LD) decay analysis, which can be informative for population demography and assessing the number of markers required to associate genetic variation with traits. Estimates of LD based on the r2 value were different between the eight populations. Compared to other breeds, LPBB population showed an overall slow decay rate and a high level of LD, which suggested that it was derived from a relatively small ancestral population (Fig. 6). Furthermore, the high LD level of LPBB population shows that there is strong linkage between SNP markers, and a lower marker density was required in the genome-wide association study.

Figure 6Linkage disequilibrium (LD) decay of eight sheep breeds, with one line per breed.


4 Discussion
Back to toptop

In this study, we evaluated for the first time the genetic diversity and population structure of LPBB with seven other populations. The mean values of the three metrics (Pn, He, F) of the downloaded six populations in this study were slightly different but had the same trend with the values reported in Kijas et al. (2012a); this might be the result of different filter criteria of SNPs.

Analysis of genetic diversity revealed that the LPBB population had the lowest level of diversity with an estimated genetic diversity (He) of 0.2861 and polymorphic SNPs (Pn) of 0.8390 (Table 1) when compared to other breeds. The lower level of genetic diversity of LPBB population was also detected in the distribution of MAF and confirmed by the lower rate of LD decay in LPBB. CME, a crossbreed of Australian Merino × Xinjiang sheep, had the highest level of diversity with an estimated He of 0.3551 and Pn of 0.9740; this was consistent with the report that the crossbreed has a higher genetic diversity level compared to the pure breed (Meadows et al., 2008; Brito et al., 2017b). Furthermore, Australian Merinos were reported as the most diverse sheep population (Meadows et al., 2008). The low genetic diversity in the LPBB population can be explained by their geographically separated habitat and small population size (Meloni et al., 2015). Lanping black-boned sheep have distributes in three towns of Yuping Mountain of Lanping. There was dozens of Lanping black-boned sheep individuals when first discovered in the 1950s, and until 2000 their number increased slowly to several hundred. In addition, due to the high elevation and steep road of Yuping Mountain, there is a geographically isolated area, and no sheep genes were introduced. The higher LD level and slow decay of LPBB population also demonstrated that the population size was relatively small. Genetic diversity is correlated with fitness (Silva et al., 2006): populations with low level of genetic diversity are expected to be less able to adapt to environmental changes (Meloni et al., 2015; Carrol and Fox, 2008) or evolving competitors and parasites (Khodabakhshzadeh et al., 2016); thus, it is suggested to conserve this unique breed.

Population structure was analyzed using PCA, structure, and NJ tree. The PCA analysis showed that the LPBB population was clustered into Asian group and mixed with LPN population, which was confirmed by the structure analysis that LPBB individuals were grouped into Asian population when K=2 and LPBB and LPN individuals were clustered into one population at the optimal K value (K=7). This result revealed that LPBB and LPN populations had great genetic similarity, which strongly explained the morphology similarity between LPBB and LPN individuals. This raises questions about whether the two populations are sufficiently different to be considered separate breeds. To search for evidence of genetic divergence, we used population FST to estimate genetic differentiation between the two populations and other breeds. The genetic differentiation observed between LPBB and LPN (FST=0.0646) was higher than that between the recognized separate breeds of TIB and CHA (FST=0.0285), which were both Tibet lineage. In addition, Kijas et al. (2012b) choose a subset of genotyped data from ISGC and grouped them into selection lines within breed, breed pairs of Mediterranean origin, and breed pairs of southern vs. northern European origin. FST values of these three groups were 0.017, 0.042, and 0.114, respectively, and this showed that the FST value that existed between LPBB and LPN was higher than selection lines within breed and breed pairs of Mediterranean origin. Therefore, on the basis of genetic data, we can conclude that LPBB can be considered a different breed from LPN.

LPBB population was clustered into Asian group and mixed with LPN population in both PCA and structure analysis, which demonstrated that LPBB has an Asian origin and had a very close relation with the LPN breed. Furthermore, the structure analysis revealed that LPBB and LPN breeds had shared ancestry with TIB breeds, which coincided with the NJ tree analysis that LPBB and LPN individuals were clustered into one clade with TIB individuals. This result agreed with the traditional classification that Lanping local sheep were Tibet lineage due to their phenotype traits such as thin tail (Hu et al., 2019) and carpet wool (Wei et al., 2015). To date, there are more than 42 native sheep breeds established in China, and they are classified into three major lineages by their geographic distribution known as northern China, Qinghai-Tibetan Plateau, and Yunnan-Kweichow Plateau sheep (Hu et al., 2019; Yang et al., 2016); Yunnan-Kweichow Plateau sheep had closer affinity with Qinghai-Tibet sheep compared with northern China sheep and could be explained by their thin-tail origin (Hu et al., 2019). Lanping is located on the Tea Horse Road leading to Tibet in Yunnan and is adjacent to Diqing Tibetan Autonomous Prefecture. This study provided the genetic evidence that LPBB sheep had close affinity with Tibet sheep. For a more accurate understanding of the origin of LPBB sheep breeds, the relation with other Yunnan-Kweichow Plateau breeds should be studied.

5 Conclusions
Back to toptop

In this study, we estimated for the first time the genetic diversity and genetic origin of Lanping black-boned sheep by using genome-wide SNP data. We observed that the LPBB population had the lowest genetic diversity compared to its geographical neighbors and four other world breeds. This finding improved our understanding of the genetic diversity in LPBB sheep breed and suggested that strategies should be implemented to maintain or increase genetic diversity in this breed. The PCA and structure analysis revealed that LPBB and LPN populations had great genetic similarity, whilst the NJ tree and genetic differentiation showed that there are differences between LPBB and LPN populations, and LPBB can be considered a specific separate breed. In addition, this study demonstrated that LPBB and LPN populations had common ancestry with the TIB population, which provided insight into the genetic origin of LPBB sheep and will be useful for the novel breed formation research.

Data availability
Back to toptop
Data availability. 

The genotyped data reported in this article are available upon request for research purposes.

Author contributions
Back to toptop
Author contributions. 

WD and HX designed and conceived the study; XL and CP performed sample collection and experiments; HX, XH, JL, and DX analyzed the data; and HX and WD wrote the article.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

The authors appreciated the help of James W. Kijas (CSIRO) for technical support, and the staff from the Animal Science and Veterinary Bureau in the Lanping County and Jianchuan County and local farmers Fuhai Yang and Sinanduojie Tang during sample collection.

Financial support
Back to toptop
Financial support. 

This research has been supported by the National Nature Science Foundation of China (grant no. 31660634) and the National Nature Science Foundation of Yunnan Province (grant no. 2016FA007).

Review statement
Back to toptop
Review statement. 

This paper was edited by Steffen Maak and reviewed by Xianyong Lan and Emiliano Lasagna.

Back to toptop

Beynon, S. E., Slavov, G. T., Farré, M., Sunduimijid, B., and Larkin, D. M.: Population structure and history of the Welsh sheep breeds determined by whole genome genotyping, BMC Genet., 16, 65,, 2015. 

Brito, L. F., Kijas, J. W., Ventura, R. V., Sargolzaei, M., Porto-Neto, L. R., Cánovas, A., Feng, Z., Jafarikia, M., and Schenkel, F. S.: Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers, BMC Genom., 18, 229,, 2017a. 

Brito, L. F., McEwan, J. C., Miller, S. P., Pickering, N. K., Bain, W. E., Dodds, K. G., Schenkel, F. S., and Clarke, S. M.: Genetic diversity of a New Zealand multi-breed sheep population and composite breeds' history revealed by a high-density SNP chip, BMC Genet., 18, 25,, 2017b. 

Carrol, S. P. and Fox, C. W.: Conservation biology: evolution in action, Oxford University Press, New York, 2008. 

Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., Handsaker, R. E., Lunter, G., Marth, G. T., Sherry, S. T., McVean, G., and Durbin, R.: The variant call format and VCFtools, Bioinformatics, 27, 2156–2158,, 2011. 

Darwish, H., Zhang, Y., Cui, K., Yang, Z., Han, D., Dong, X., Mao, H., Deng, W., and Deng, X.: Molecular cloning and characterization of the endothelin 3 gene in black bone sheep, J. Anim. Sci. Biotechnol., 9, 57,, 2018. 

Deng, W. D., Yang, S. L., Huo, Y. Q., Gou, X., Shi, X. W., and Mao, H. M.: Physiological and genetic characteristics of black-boned sheep (Ovis aries), Anim. Genet., 37, 586–588,, 2006. 

Deng, W. D., Xi, D. M., Gou, X., Yang, S. L., Shi, X. W., and Mao, H. M.: Pigmentation in Black-boned sheep (Ovis aries): association with polymorphism of the Tyrosinase gene, Molec. Biol. Rep., 35, 379–385,, 2008. 

Deng, W. D., Shu, W., Yang, S. L., Shi, X. W., and Mao, H. M.: Pigmentation in Black-boned sheep (Ovis aries): association with polymorphism of the MC1R gene, Molec. Biol. Rep., 36, 431–436,, 2009a. 

Deng, W. D., Tan, Y., Wang, X., Xi, D., He, Y., Yang, S., Mao, H., and Gao, S.: Molecular cloning, sequence characteristics, and polymorphism analyses of the tyrosinase-related protein 2/DOPAchrome tautomerase gene of black-boned sheep (Ovis aries), Genome, 52, 1001–1011,, 2009b. 

Falush, D., Stephens, M., and Pritchard, J. K.: Inference of population structure using multilocus genotype data: dominant markers and null alleles, Molec. Ecol. Notes, 7, 574–578,, 2007. 

Gorkhali, N. A., Dong, K., Yang, M., Song, S., Kader, A., Shrestha, B. S., He, X., Zhao, Q., Pu, Y., Li, X., Kijas, J., Guan, W., Han, J., Jiang, L., and Ma, Y.: Genomic analysis identified a potential novel molecular mechanism for high-altitude adaptation in sheep at the Himalayas, Scient. Rep., 6, 29963,, 2016. 

Grasso, A. N., Goldberg, V., Navajas, E. A., Iriarte, W., Gimeno, D., Aguilar, I., Medrano, J. F., Rincon, G., and Ciappesoni, G.: Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale, Merino and Creole sheep, Genet. Mol. Biol., 37, 389–395,, 2014. 

Hall, B. G.: Building Phylogenetic Trees from Molecular Data with MEGA, Molec. Biol. Evol., 30, 1229–1235,, 2013. 

Hu, X., Yang, J., Xie, X., Lv, F., Cao, Y., Li, W., Liu, M., Wang, Y., Li, J., Liu, Y., Ren, Y., Shen, Z., Wang, F., Hehua, E., Han, J., and Li, M.: The Genome Landscape of Tibetan Sheep Reveals Adaptive Introgression from Argali and the History of Early Human Settlements on the Qinghai-Tibetan Plateau, Molec. Biol. Evol., 36, 283–303,, 2019. 

Johnston, S. E., Bérénos, C., Slate, J., and Pemberton, J. M.: Conserved Genetic Architecture Underlying Individual Recombination Rate Variation in a Wild Population of Soay Sheep (Ovis aries), Genetics, 203, 583–598,, 2016. 

Khodabakhshzadeh, R., Mohammadabadi, M. R., Esmailizadeh, A. K., Moradi Shahrebabak, H., Bordbar, F., and Ansari Namin, S.: Identification of point mutations in exon 2 of GDF9 gene in Kermani sheep, Pol. J. Veterin. Sci., 19, 281–289,, 2016. 

Kijas, J. W., Townley, D., Dalrymple, B. P., Heaton, M. P., Maddox, J. F., McGrath, A., Wilson, P., Ingersoll, R. G., McCulloch, R., McWilliam, S., Tang, D., McEwan, J., Cockett, N., Oddy, V. H., Nicholas, F. W., and Raadsma, H.: A genome wide survey of SNP variation reveals the genetic structure of sheep breeds, PLoS One, 4, e4668,, 2009. 

Kijas, J. W., Lenstra, J. A., Hayes, B., Boitard, S., Porto Neto, L. R., San Cristobal, M., Servin, B., McCulloch, R., Whan, V., Gietzen, K., Paiva, S., Barendse, W., Ciani, E., Raadsma, H., McEwan, J., and Dalrymple, B.: Genome-wide analysis of the world's sheep breeds reveals high levels of historic mixture and strong recent selection, PLoS Biol., 10, e1001258,, 2012a. 

Kijas, J. W., Miller, J. E., Hadfield, T., McCulloch, R., Garcia-Gamez, E., Porto, N. L., and Cockett, N.: Tracking the emergence of a new breed using 49,034 SNP in sheep, PLoS One, 7, e41508,, 2012b. 

Li, G., Xiong, H., Xi, D., Memon, S., and Deng, W.: An examination of melanogenic traits and TYRP1 polymorphism in Nanping and Romney Marsh sheep breeds, Arch. Anim. Breed., 61, 131–141,, 2018.  

Li, W. B.: Research progress of Lanping black-boned sheep, Chin. Livest. Poult. Breed., 5, 58–60, 2009. 

Meadows, J. R., Chan, E. K., and Kijas, J. W.: Linkage disequilibrium compared between five populations of domestic sheep, BMC Genet., 9, 61,, 2008. 

Meloni, M., Reid, A., Caujapé-Castells, J., Soto, M., Fernández-Palacios, J. M., and Conti, E.: High genetic diversity and population structure in the endangered Canarian endemic Ruta oreojasme (Rutaceae), Genetica, 143, 571–580,, 2015. 

Nie, G., Huang, L., Zhang, X., Taylor, M., Jiang, Y., Yu, X., Liu, X., Wang, X., and Zhang, Y.: Marker-Trait Association for Biomass Yield of Potential Bio-fuel Feedstock Miscanthus sinensis from Southwest China, Front. Plant Sci., 7, 802,, 2016. 

Notter, D. R.: The importance of genetic diversity in livestock populations of the future, J. Anim, Sci., 77, 61–69,, 1999. 

Pan, Z., Li, S., Liu, Q., Wang, Z., Zhou, Z., Di, R., Miao, B., Hu, W., Wang, X., Hu, X., Xu, Z., Wei, D., He, X., Yuan, L., Guo, X., Liang, B., Wang, R., Li, X., Cao, X., Dong, X., Xia, Q., Shi, H., Hao, G., Yang, J., Luosang, C., Zhao, Y., Jin, M., Zhang, Y., Lv, S., Li, F., Ding, G., Chu, M., and Li, Y.: Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization, GigaScience, 7, 4,, 2018. 

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., Maller, J., Sklar, P., de Bakker, P. I. W., Daly, M. J., and Sham, P. C.: PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses, Am. J. Human Genet., 81, 559–575,, 2007. 

Silva, A. D., Luikart, G., Yoccoz, N. G., Cohas, A., and Allainé, D.: Genetic diversity-fitness correlation revealed by microsatellite analyses in European alpine marmots (Marmota marmota), Conserv. Genet., 7, 371–382,, 2006. 

Wei, C., Wang, H., Liu, G., Wu, M., Cao, J., Liu, Z., Liu, R., Zhao, F., Zhang, L., Lu, J., Liu, C., and Du, L.: Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds, BMC Genom., 16, 194,, 2015. 

Yang, J., Lee, S. H., Goddard, M. E., and Visscher, P. M.: GCTA: A Tool for Genome-wide Complex Trait Analysis, Am. J. Human Genet., 88, 76–82,, 2011. 

Yang, J., Li, W., Lv, F., He, S., Tian, S., Peng, W., Sun, Y., Zhao, Y., Tu, X., Zhang, M., Xie, X., Wang, Y., Li, J., Liu, Y., Shen, Z., Wang, F., Liu, G., Lu, H., Kantanen, J., Han, J., Li, M., and Liu, M.: Whole-Genome Sequencing of Native Sheep Provides Insights into Rapid Adaptations to Extreme Environments, Molec. Biol. Evol., 33, 2576–2592,, 2016. 

Short summary
Lanping black-boned sheep is a novel and rare sheep breed discovered in the 1950s in China and has the same hyperpigmentation trait as silk fowl. Lanping black-bones sheep has same morphology with Lanping native sheep such as coat color, and horns and tail shape, so it's interesting to investigate the genetic relationship of the two breeds and the genetic origin of Lanping black-boned sheep. Lanping black-boned sheep and Lanping native sheep have great genetic similarity.
Lanping black-boned sheep is a novel and rare sheep breed discovered in the 1950s in China and...