Background: The quality and yield of wool determine the economic value of the fine-wool sheep. Therefore, discovering markers or genes relevant to wool traits is the cornerstone for the breeding of fine-wool sheep. In this study, we used the Illumina HiSeq X Ten platform to re-sequence 460 sheep belonging to four different fine-wool sheep breeds, namely, Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Aohan fine-wool sheep (AHS) and Qinghai fine-wool sheep (QHS). Eight wool traits, including fiber diameter (FD), fiber diameter coefficient of variance (FDCV), fiber diameter standard deviation (FDSD), staple length (SL), greasy fleece weight (GFW), clean wool rate (CWR), staple strength (SS) and staple elongation (SE) were examined. A genome-wide association study (GWAS) was performed to detect the candidate genes for the eight wool traits. Results: A total of 8.222 Tb of raw data was generated, with an average of approximately 8.59X sequencing depth. After quality control, 12,561,225 SNPs were available for analysis. And a total of 57 genome-wide significant SNPs and 30 candidate genes were detected for the desired wool traits. Among them, 7 SNPs and 6 genes are related to wool fineness indicators (FD, FDCV and FDSD), 10 SNPs and 7 genes are related to staple length, 13 SNPs and 7 genes are related to wool production indicators (GFW and CWR), 27 SNPs and 10 genes associated with staple elongation. Among these candidate genes, UBE2E3 and RHPN2 associated with fiber diameter, were found to play an important role in keratinocyte differentiation and cell proliferation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results, revealed that multitude significant pathways are related to keratin and cell proliferation and differentiation, such as positive regulation of canonical Wnt signaling pathway (GO:0090263). Conclusion: This is the first GWAS on the wool traits by using re-sequencing data in Chinese fine-wool sheep. The newly detected significant SNPs in this study can be used in genome-selective breeding for the fine-wool sheep. And the new candidate genes would provide a good theoretical basis for the fine-wool sheep breeding.

Background: Copy number variation (CNV) is an important source of genetic variation that has a significant influence on phenotypic diversity, economically important traits and the evolution of livestock species. In this study, the genome-wide CNV distribution characteristics of 32 fine-wool sheep from three breeds were analyzed using resequencing. Results: A total of 1,747,604 CNVs were detected in this study, and 7228 CNV regions (CNVR) were obtained after merging overlapping CNVs; these regions accounted for 2.17% of the sheep reference genome. The average length of the CNVRs was 4307.17 bp. “Deletion” events took place more frequently than “duplication” or “both” events. The CNVRs obtained overlapped with previously reported sheep CNVRs to variable extents (4.39–55.46%). Functional enrichment analysis showed that the CNVR-harboring genes were mainly involved in sensory perception systems, nutrient metabolism processes, and growth and development processes. Furthermore, 1855 of the CNVRs were associated with 166 quantitative trait loci (QTL), including milk QTLs, carcass QTLs, and health-related QTLs, among others. In addition, the 32 fine-wool sheep were divided into horned and polled groups to analyze for the selective sweep of CNVRs, and it was found that the relaxin family peptide receptor 2 (RXFP2) gene was strongly influenced by selection. Conclusions: In summary, we constructed a genomic CNV map for Chinese indigenous fine-wool sheep using resequencing, thereby providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep.

Dorper sheep (Ovis aries) (DPS), developed in the 1930s by crossing Dorset Horn and Blackhead Persian sheep in South Africa, is a world-famous composite breed for mutton production. The genetic basis underlying this breed is yet to be elucidated. Here, we report the sequencing and assembly of a highly contiguous Dorper sheep genome via integration of Oxford Nanopore Technology (ONT) sequencing and Hi-C (chromatin conformation capture) approaches. The assembled genome was around 2.64 Gb with a contig N50 of 73.33 Mb and 140 contigs in total. More than 99.5% of the assembled sequences could be anchored to 27 chromosomes and they were annotated with 20,450 protein-coding genes. Allele-specific expression (ASE) genes of Dorper sheep were revealed through ASE analysis and they were involved in the immune system, lipid metabolism, and environmental adaptation. A total of 5,701 and 456 allelic sites were observed in the SNP and indels loci identified from relevant whole-genome resequencing data. These allelic SNP and INDEL sites were annotated in 1,002 and 294 genes, respectively. Moreover, we calculated the number of variant sites and related genes derived from the maternal and paternal ancestors, revealing the genetic basis of outstanding phenotypic performance of Dorper sheep. In conclusion, this study reports the first reference genome of Dorper sheep and reveals its genetic basis through ASE. This study also provides a pipeline for mining genetic information of composite breeds, which has an implication for future hybrid-breeding practices.

Hair follicle development and wool shedding in sheep are poorly understood. This study investigated the population structures and genetic differences between sheep with different wool types to identify candidate genes related to these traits. We used Illumina ovine SNP 50K chip genotyping data of 795 sheep populations comprising 27 breeds with two wool types, measuring the population differentiation index (Fst), nucleotide diversity (θπ ratio), and extended haplotype homozygosity among populations (XP-EHH) to detect the selective signatures of hair sheep and fine-wool sheep. The top 5% of the Fst and θπ ratio values, and values of XP-EHH < −2 were considered strongly selected SNP sites. Annotation showed that the PRX, SOX18, TGM3, and TCF3 genes related to hair follicle development and wool shedding were strongly selected. Our results indicated that these methods identified important genes related to hair follicle formation, epidermal differentiation, and hair follicle stem cell development, and provide a meaningful reference for further study on the molecular mechanisms of economically important traits in sheep.

To investigate single nucleotide polymorphism (SNP) loci associated with yearling wool traits of fine-wool sheep for optimizing marker-assisted selection and dissection of the genetic architecture of wool traits, we conducted a genome-wide association study (GWAS) based on the fixed and random model circulating probability unification (FarmCPU) for yearling staple length (YSL), yearling mean fiber diameter (YFD), yearling greasy fleece weight (YGFW), and yearling clean fleece rate (YCFR) by using the whole-genome re-sequenced data (totaling 577 sheep) from the following four fine-wool sheep breeds in China: Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Qinghai fine-wool sheep (QHS), and Aohan fine-wool sheep (AHS). A total of 16 SNPs were detected above the genome-wise significant threshold (P = 5.45E-09), and 79 SNPs were located above the suggestive significance threshold (P = 5.00E-07) from the GWAS results. For YFD and YGFW traits, 7 and 9 SNPs reached the genome-wise significance thresholds, whereas 10 and 12 SNPs reached the suggestive significance threshold, respectively. For YSL and YCFR traits, none of the SNPs reached the genome-wise significance thresholds, whereas 57 SNPs exceeded the suggestive significance threshold. We recorded 14 genes located at the region of ±50-kb near the genome-wise significant SNPs and 59 genes located at the region of ±50-kb near the suggestive significant SNPs. Meanwhile, we used the Average Information Restricted Maximum likelihood algorithm (AI-REML) in the “HIBLUP” package to estimate the heritability and variance components of the four desired yearling wool traits. The estimated heritability values (h²) of YSL, YFD, YGFW, and YCFR were 0.6208, 0.7460, 0.6758, and 0.5559, respectively. We noted that the genetic parameters in this study can be used for fine-wool sheep breeding. The newly detected significant SNPs and the newly identified candidate genes in this study would enhance our understanding of yearling wool formation, and significant SNPs can be applied to genome selection in fine-wool sheep breeding.

Wool is produced via synthetic processes of wool follicles, which are embedded in the skin of sheep. The development of new-generation sequencing and RNA sequencing provides new approaches that may elucidate the molecular regulation mechanism of wool follicle development and facilitate enhanced selection for wool traits through gene-assisted selection or targeted gene manipulation. We performed de novo transcriptome sequencing of skin using the Illumina Hiseq 2000 sequencing system in sheep (Ovis aries). Transcriptome de novo assembly was carried out via short-read assembly programs, including SOAPdenovo and ESTScan. The protein function, clusters of orthologous group function, gene ontology function, metabolic pathway analysis, and protein coding region prediction of unigenes were annotated by BLASTx, BLAST2GO, and ESTScan. More than 26,266,670 clean reads were collected and assembled into 79,741 unigene sequences, with a final assembly length of 35,447,962 nucleotides. A total of 22,164 unigenes were annotated, accounting for 36.27% of the total number of unigenes, which were divided into 25 classes belonging to 218 signaling pathways. Among them, there were 17 signal paths related to hair follicle development. Based on mass sequencing data of sheepskin obtained by RNA-Seq, many unigenes were identified and annotated, which provides an excellent platform for future sheep genetic and functional genomic research. The data could be used for improving wool quality and as a model for human hair follicle development or disease prevention.