The Transcriptional Cell Atlas of Testis Development in Sheep at Pre-Sexual Maturity.

Sheep testes undergo a dramatic rate of development with structural changes during pre-sexual maturity, including the proliferation and maturation of somatic niche cells and the initiation of spermatogenesis. To explore this complex process, 12,843 testicular cells from three males at pre-sexual maturity (three-month-old) were sequenced using the 10× Genomics ChromiumTM single-cell RNA-seq (scRNA-seq) technology. Nine testicular somatic cell types (Sertoli cells, myoid cells, monocytes, macrophages, Leydig cells, dendritic cells, endothelial cells, smooth muscle cells, and leukocytes) and an unknown cell cluster were observed. In particular, five male germ cell types (including two types of undifferentiated spermatogonia (Apale and Adark), primary spermatocytes, secondary spermatocytes, and sperm cells) were identified. Interestingly, Apale and Adark were found to be two distinct states of undifferentiated spermatogonia. Further analysis identified specific marker genes, including UCHL1, DDX4, SOHLH1, KITLG, and PCNA, in the germ cells at different states of differentiation. The study revealed significant changes in germline stem cells at pre-sexual maturation, paving the way to explore the candidate factors and pathways for the regulation of germ and somatic cells, and to provide us with opportunities for the establishment of livestock stem cell breeding programs.

Selective Sweeps Uncovering the Genetic Basis of Horn and Adaptability Traits on Fine-Wool Sheep in China.

Long-term natural and artificial selection leads to change in certain regions of the genome, resulting in selection signatures that can reveal genes associated with selected traits, such as horns (i.e., polled/horned), high-quality wool traits, and high-altitude hypoxia adaptability. These are complex traits determined by multiple genes, regulatory pathways, and environmental factors. A list of genes with considerable effects on horn and adaptability traits has not been found, although multiple quantitative trait loci (QTL) have been identified. Selection signatures could be identified using genetic differentiation (F ST ), polymorphism levels θπ, and Tajima's D. This study aimed to identify selection signatures in fine-wool sheep and to investigate the genes annotated in these regions, as well as the biological pathways involved in horn and adaptability traits. For this purpose, the whole-genome sequence of 120 individuals from four breeds, which come from different elevations and habitats in China, was used to analyze selection signatures for horn and adaptability traits. Annotation of the consensus regions of F ST and θπ ratios revealed a list of identified genes associated with polled/horned and high-altitude hypoxia adaptability traits, such as RXPF2, EERFC4, MSH6, PP1R12A, THBS1, ATP1B2, RYR2, and PLA2G2E. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified genes related primarily to mismatch repair, metabolism, vascular smooth muscle contraction, and cardiac muscle contraction. This is the first study to demonstrate that selection signatures play an important role in the polled/horned and high-altitude hypoxia adaptability traits of fine-wool sheep breeds that have undergone high-intensity selection and adapted to different ecological environments in China. Changes observed in the genome of fine-wool sheep may have acted on genomic regions that affect performance traits and provide a reference for genome design and breeding.

Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs.

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs' expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17ß4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing.

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.

Comprehensive Analysis of Long Non-coding RNA and mRNA Transcriptomes Related to Hypoxia Adaptation in Tibetan Sheep.

Tibetan sheep have lived on the Qinghai-Tibet Plateau for a long time, and after long-term natural selection, they have shown stable genetic adaptability to high-altitude environments. However, little is known about the molecular mechanisms of the long non-coding (lnc)RNAs involved in the adaptation of Tibetan sheep to hypoxia. Here, we collected lung tissues from high-altitude Tibetan sheep and low-altitude Hu sheep for RNA sequencing to study the regulatory mechanisms of the lncRNAs and mRNAs in the adaptation of Tibetan sheep to hypoxia. We identified 254 differentially expressed lncRNAs and 1,502 differentially expressed mRNAs. We found 20 pairs of cis-regulatory relationships between 15 differentially expressed lncRNAs and 14 protein-coding genes and two pairs of trans-regulatory relationships between two differentially expressed lncRNAs and two protein-coding genes. These differentially expressed mRNAs and lncRNA target genes were mainly enriched in pathways related to lipid metabolism and immune function. Interaction network analysis showed that 17 differentially expressed lncRNAs and 15 differentially expressed mRNAs had an interactive relationship. Additionally, we used six differentially expressed lncRNAs and mRNAs to verify the accuracy of the sequencing data via qRT-PCR. Our results provide a comprehensive overview of the expression patterns of the lncRNAs and mRNAs involved in the adaptation of Tibetan sheep to hypoxia, laying a foundation for further analysis of the adaptations of plateau animals.

Liquid Chromatography Quadrupole Time-of-Flight Mass Spectrometry and Rapid Evaporative Ionization Mass Spectrometry Were Used to Develop a Lamb Authentication Method: A Preliminary Study.

A untargeted metabolomics approach was proposed in this study based on ultra-high performance liquid chromatography quadrupole time-of-flight (UHPLC-QTOF) and rapid evaporative ionization mass spectrometry (REIMS) to discriminate lamb and mutton meat and investigate their subtle metabolic differences, considering the higher popularity of lamb meat than mutton in the market. Multivariate statistical analysis was performed for data processing in order to distinguish between the two sample types. A total of 42 potential metabolites (20 in positive and 22 in negative ion mode) were defined for UHPLC-QTOF analysis, which provided references for discriminating the two kinds of meat. Furthermore, three potential markers were tentatively identified using LC/MS data against chemical databases. In addition, 14 potential metabolites were putatively identified in negative ion mode using the LipidMaps database. Meanwhile, the data-driven soft independent modeling of class analogy (DD-SIMCA) model was established, which could rapidly differentiate non-pretreated lamb meat and mutton with 92% specificity, rendering REIMS a promising technique for meat identification.

Analysis of dynamic and widespread lncRNA and miRNA expression in fetal sheep skeletal muscle.

The sheep is an economically important animal, and there is currently a major focus on improving its meat quality through breeding. There are variations in the growth regulation mechanisms of different sheep breeds, making fundamental research on skeletal muscle growth essential in understanding the regulation of (thus far) unknown genes. Skeletal muscle development is a complex biological process regulated by numerous genes and non-coding RNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). In this study, we used deep sequencing data from sheep longissimus dorsi (LD) muscles sampled at day 60, 90, and 120 of gestation, as well as at day 0 and 360 following birth, to identify and examine the lncRNA and miRNA temporal expression profiles that regulate sheep skeletal myogenesis. We stained LD muscles using histological sections to analyse the area and circumference of muscle fibers from the embryonic to postnatal development stages. Our results showed that embryonic skeletal muscle growth can be characterized by time. We obtained a total of 694 different lncRNAs and compared the differential expression between the E60 vs. E90, E90 vs. E120, E120 vs. D0, and D0 vs. D360 lncRNA and gene samples. Of the total 701 known sheep miRNAs we detected, the following showed a wide range of expression during the embryonic stage: miR-2387, miR-105, miR-767, miR-432, and miR-433. We propose that the detected lncRNA expression was time-specific during the gestational and postnatal stages. GO and KEGG analyses of the genes targeted by different miRNAs and lncRNAs revealed that these significantly enriched processes and pathways were consistent with skeletal muscle development over time across all sampled stages. We found four visual lncRNA-gene regulatory networks that can be used to explore the function of lncRNAs in sheep and may be valuable in helping improve muscle growth. This study also describes the function of several lncRNAs that interact with miRNAs to regulate myogenic differentiation.

Genetic signatures of high-altitude adaptation and geographic distribution in Tibetan sheep.

Most sheep breeding programs designed for the tropics and sub-tropics have to take into account the impacts of environmental adaptive traits. However, the genetic mechanism regulating the multiple biological processes driving adaptive responses remains unclear. In this study, we applied a selective sweep analysis by combing 1% top values of Fst and ZHp on both altitude and geographic subpopulations (APS) in 636 indigenous Tibetan sheep breeds. Results show that 37 genes were identified within overlapped genomic regions regarding Fst significantly associated with APS. Out of the 37 genes, we found that 8, 3 and 6 genes at chromosomes (chr.) 13, 23 and 27, respectively, were identified in the genomic regions with 1% top values of ZHp. We further analyzed the INDEL variation of 6 genes at chr.27 (X chromosome) in APS together with corresponding orthologs of 6 genes in Capra, Pantholops, and Bos Taurus. We found that an INDEL was located within 5'UTR region of HAG1 gene. This INDEL of HAG1 was strongly associated with the variation of APS, which was further confirmed by qPCR. Sheep breeds carrying "C-INDEL" of HAG1 have significantly greater body weight, shear amount, corpuscular hemoglobin and globulin levels, but lower body height, than those carrying "CA-INDEL" of HAG1. We concluded that "C-INDEL" variation of HAG1 gene confers better hypoxia tolerance in the highlands of Tibetan and explains well geographic distributions in this population. These results contribute to our understanding of adaptive responses to altitude and geographic adaptation in Tibetan sheep populations and will help to guide future conservation programs for Tibetan sheep native to Qinghai-Tibetan Plateau.

Comparative proteomics reveals genetic mechanisms underlying secondary hair follicle development in fine wool sheep during the fetal stage.

The aim of this study was to investigate the genetic mechanisms underlying wool production by characterizing the skin protein profile and determining the proteomic changes that occur as a consequence of development in wool-producing sheep using a label-free proteomics approach. Samples were collected at four stages during gestation (87, 96, 102, and 138 days), and every two consecutive stages were statistically compared (87 versus 96, 96 versus 102, and 102 versus 138 days). We identified 227 specific proteins in the sheep proteome that were present in all four stages, and 123 differentially abundant proteins (DAPs). We also observed that the microstructure of the secondary follicles changed significantly during the development of the fetal skin hair follicle. The screened DAPs were strictly related to metabolic and skin development pathways, and were associated with pathways such as the glycolysis/gluconeogenesis. These analyses indicated that the wool production of fine wool sheep is regulated via a variety of pathways. These findings provide an important resource that can be used in future studies of the genetic mechanisms underlying wool traits in fine wool sheep, and the identified DAPs should be further investigated as candidate markers for predicting wool traits in sheep. SIGNIFICANCE: Wool quality (fiber diameter, length, etc.) is an important economic trait of fine wool sheep that is determined by secondary follicle differentiation and re-differentiation. Secondary follicles of fine wool sheep developed from a bud (87 days), and underwent differentiation (96 days) and rapid growth (102 days) until maturity (138 days) during gestation. Comparative analysis based on differential proteomics of these four periods could provide a better understanding of the wool growth mechanism of fine wool sheep and offer novel strategies for improving fine wool quality by breeding.

Genome-Wide Association Study of Body Weight Traits in Chinese Fine-Wool Sheep.

Body weight is an important economic trait for sheep and it is vital for their successful production and breeding. Therefore, identifying the genomic regions and biological pathways that contribute to understanding variability in body weight traits is significant for selection purposes. In this study, the genome-wide associations of birth, weaning, yearling, and adult weights of 460 fine-wool sheep were determined using resequencing technology. The results showed that 113 single nucleotide polymorphisms (SNPs) reached the genome-wide significance levels for the four body weight traits and 30 genes were annotated effectively, including AADACL3, VGF, NPC1, and SERPINA12. The genes annotated by these SNPs significantly enriched 78 gene ontology terms and 25 signaling pathways, and were found to mainly participate in skeletal muscle development and lipid metabolism. These genes can be used as candidate genes for body weight in sheep, and provide useful information for the production and genomic selection of Chinese fine-wool sheep.

Evaluation of Crossbreeding of Australian Superfine Merinos with Gansu Alpine Finewool Sheep to Improve Wool Characteristics.

Crossbreeding of Australian Superfine Merinos (ASMs) with Gansu Alpine Finewool (GAF) sheep and an evaluation of the potential benefits of this genetic cross has not been previously conducted. 13 ASMs were crossbred with GAF sheep over a five year period with backcrossing designed to assess heterosis. Data from 11,178 lambs sired by 189 rams were used in the study. Genotype, birth year, birth type, dam age, sex and/or management group, and record age were fitted as fixed effects and within-genotype sire fitted as a random effect. Crossbreeds of 1/2 ASM expressed the most desirable effects for improving average fiber diameter (AFD), clean fleece weight (CFW), yield, coefficient of variation of AFD (CVAFD), yearling staple length (YSL) to AFD ratio (YSL/AFD), and CFW to metabolic yearling bodyweight (YWT0.75) ratio (CFW/YWT0.75) but showed the least post-weaning average daily gain (powADG) and YWT. Genotype of backcrossing with 1/4 ASM obtained moderate improvements in AFD, CFW, CVAFD, and YSL/AFD but the highest YSL, WWT, and prwADG. Except for yield (-1.42%) and CFW/YWT0.75 (-1%), heterosis estimates were generally low and positive, and ranged from 0.1% for CVAFD to 4% for powADG, which indicates the potential to improve relevant traits through exploiting heterosis to a varying extent. The ASMs sampled in this study were found to be superior to GAFs for AFD, CFW, yield, and CVAFD by 19.82%, 11.68%, 14.47%, and 6.99%, respectively, but inferior for YSL, PowADG, and YWT by 4.36%, 50.97%, and 16.93%, respectively. ASMs also appeared to be more efficient than GAFs in clean wool production (25.34%) and staple length growth (16.17%). The results of our study strongly suggest that an infusion of ASM genes via crossbreeding is an effective and appropriate approach to improve wool microns and wool production from GAF sheep, and we make recommendations to tackle the undesirable traits of YWT and YSL from ASM introduction.

Genetic Diversity and Phylogenetic Evolution of Tibetan Sheep Based on mtDNA D-Loop Sequences.

The molecular and population genetic evidence of the phylogenetic status of the Tibetan sheep (Ovis aries) is not well understood, and little is known about this species' genetic diversity. This knowledge gap is partly due to the difficulty of sample collection. This is the first work to address this question. Here, the genetic diversity and phylogenetic relationship of 636 individual Tibetan sheep from fifteen populations were assessed using 642 complete sequences of the mitochondrial DNA D-loop. Samples were collected from the Qinghai-Tibetan Plateau area in China, and reference data were obtained from the six reference breed sequences available in GenBank. The length of the sequences varied considerably, between 1031 and 1259 bp. The haplotype diversity and nucleotide diversity were 0.992±0.010 and 0.019±0.001, respectively. The average number of nucleotide differences was 19.635. The mean nucleotide composition of the 350 haplotypes was 32.961% A, 29.708% T, 22.892% C, 14.439% G, 62.669% A+T, and 37.331% G+C. Phylogenetic analysis showed that all four previously defined haplogroups (A, B, C, and D) were found in the 636 individuals of the fifteen Tibetan sheep populations but that only the D haplogroup was found in Linzhou sheep. Further, the clustering analysis divided the fifteen Tibetan sheep populations into at least two clusters. The estimation of the demographic parameters from the mismatch analyses showed that haplogroups A, B, and C had at least one demographic expansion in Tibetan sheep. These results contribute to the knowledge of Tibetan sheep populations and will help inform future conservation programs about the Tibetan sheep native to the Qinghai-Tibetan Plateau.

Integrated Analysis of the Roles of Long Noncoding RNA and Coding RNA Expression in Sheep (Ovis aries) Skin during Initiation of Secondary Hair Follicle.

Initiation of hair follicle (HF) is the first and most important stage of HF morphogenesis. However the precise molecular mechanism of initiation of hair follicle remains elusive. Meanwhile, in previous study, the more attentions had been paid to the function of genes, while the roles of non-coding RNAs (such as long noncoding RNA and microRNA) had not been described. Therefore, the roles of long noncoding RNA(LncRNA) and coding RNA in sheep skin during the initiation of sheep secondary HF were integrated and analyzed, by using strand-specific RNA sequencing (ssRNA-seq).A total of 192 significant differentially expressed genes were detected, including 67 up-regulated genes and 125 down-regulated genes between stage 0 and stage 1 of HF morphogenesis during HF initiation. Only Wnt2, FGF20 were just significant differentially expressed among Wnt, Shh, Notch and BMP signaling pathways. Further expression profile analysis of lncRNAs showed that 884 novel lncRNAs were discovered in sheep skin expression profiles. A total of 15 lncRNAs with significant differential expression were detected, 6 up-regulated and 9 down-regulated. Among of differentially expressed genes and LncRNA, XLOC002437 lncRNA and potential target gene COL6A6 were all significantly down-regulated in stage 1. Furthermore, by using RNAhybrid, XLOC005698 may be as a competing endogenous RNA ''sponges" oar-miR-3955-5p activity. Gene Ontology and KEGG pathway analyses indicated that the significantly enriched pathway was peroxisome proliferator-activated receptors (PPARs) pathway (corrected P-value < 0.05), indicating that PPAR pathway is likely to play significant roles during the initiation of secondary HF.Results suggest that the key differentially expressed genes and LncRNAs may be considered as potential candidate genes for further study on the molecular mechanisms of HF initiation, as well as supplying some potential values for understanding human hair disorders.

Exploring Differentially Expressed Genes and Natural Antisense Transcripts in Sheep (Ovis aries) Skin with Different Wool Fiber Diameters by Digital Gene Expression Profiling.

Wool fiber diameter (WFD) is the most important economic trait of wool. However, the genes specifically controlling WFD remain elusive. In this study, the expression profiles of skin from two groups of Gansu Alpine merino sheep with different WFD (a super-fine wool group [FD = 18.0 ± 0.5 µm, n=3] and a fine wool group [FD=23.0 ± 0.5 µm, n=3]) were analyzed using next-generation sequencing-based digital gene expression profiling. A total of 40 significant differentially expressed genes (DEGs) were detected, including 9 up-regulated genes and 31 down-regulated genes. Further expression profile analysis of natural antisense transcripts (NATs) showed that more than 30% of the genes presented in sheep skin expression profiles had NATs. A total of 7 NATs with significant differential expression were detected, and all were down-regulated. Among of 40 DEGs, 3 DEGs (AQP8, Bos d2, and SPRR) had significant NATs which were all significantly down-regulated in the super-fine wool group. In total of DEGs and NATs were summarized as 3 main GO categories and 38 subcategories. Among the molecular functions, cellular components and biological processes categories, binding, cell part and metabolic process were the most dominant subcategories, respectively. However, no significant enrichment of GO terms was found (corrected P-value >0.05). The pathways that were significantly enriched with significant DEGs and NATs were mainly the lipoic acid metabolism, bile secretion, salivary secretion and ribosome and phenylalanine metabolism pathways (P < 0.05). The results indicated that expression of NATs and gene transcripts were correlated, suggesting a role in gene regulation. The discovery of these DEGs and NATs could facilitate enhanced selection for super-fine wool sheep through gene-assisted selection or targeted gene manipulation in the future.