Two mutations at KRT74 and EDAR synergistically drive the fine-wool production in Chinese sheep.

INTRODUCTION: Fine-wool sheep are the most common breed used by the wool industry worldwide. Fine-wool sheep have over a three-fold higher follicle density and a 50% smaller fiber diameter than coarse-wool sheep. OBJECTIVES: This study aims to clarify the underlying genetic basis for the denser and finer wool phenotype in fine-wool breeds. METHOD: Whole-genome sequences of 140 samples, Ovine HD630K SNP array data of 385 samples, including fine, semi-fine, and coarse wool sheep, as well as skin transcriptomes of nine samples were integrated for genomic selection signature analysis. RESULTS: Two loci at keratin 74 (KRT74) and ectodysplasin receptor (EDAR) were revealed. Fine-scale analysis in 250 fine/semi-fine and 198 coarse wool sheep narrowed this association to one C/A missense variant of KRT74 (OAR3:133,486,008, P = 1.02E-67) and one T/C SNP in the regulatory region upstream of EDAR (OAR3:61,927,840, P = 2.50E-43). Cellular over-expression and ovine skin section staining assays confirmed that C-KRT74 activated the KRT74 protein and specifically enlarged cell size at the Huxley's layer of the inner root sheath (P < 0.01). This structure enhancement shapes the growing hair shaft into the finer wool than the wild type. Luciferase assays validated that the C-to-T mutation upregulated EDAR mRNA expression via a newly created SOX2 binding site and potentially led to the formation of more hair placodes. CONCLUSIONS: Two functional mutations driving finer and denser wool production were characterized and offered new targets for genetic breeding during wool sheep selection. This study not only provides a theoretical basis for future selection of fine wool sheep breeds but also contributes to improving the value of wool commodities.

Transcriptome profiling of differentiating adipose-derived stem cells across species reveals new genes regulating adipogenesis.

Adipose-derived stem cells (ADSCs) that are enriched in adipose tissue with multilineage differentiation potential have become an important tool in therapeutic research and tissue engineering. Certain breeds of sheep exhibit a unique fat tail trait such that tail tissue accounts for approximately 10 % of body weight and can provide an excellent source of ADSCs. Here, we describe isolation of primary ADSCs from ovine embryonic fat tail tissues that displayed high self-renewal capacity, multilineage differentiation and excellent adipogenic ability. Through transcriptome analysis covering ADSCs differentiating into adipocytes, 37 transcription factors were involved in early transcriptional events that initiate a regulatory cascade of adipogenesis; the entire adipogenic activity consists of a reduction in proliferation ability and upregulation of genes related to lipid generation and energy metabolism, as well as several genes associated with myogenesis. Furthermore, Comparative transcriptome analysis across species (sheep, human, and mouse) revealed enhanced basal metabolic ability in differentiating ovine ADSCs, which may relate to the excellent adipogenic capability of these cells. We also identified a small evolutionarily conserved gene set, consisting of 21 and 22 genes exhibiting increased and decreased expression, respectively. Almost half (20) of these genes have not previously been reported to regulate adipogenesis in mammals. In this study, we identified important regulators that trigger ovine adipocyte differentiation, main biological pathways involved in adipogenesis as well as the evolutionarily conserved genes governing adipogenic process across species. Our study provides a novel excellent biomaterial and novel genes regulating adipogenesis for cellular transplantation therapy and investigations of fat metabolism.

Genome-Wide Signal Selection Analysis Revealing Genes Potentially Related to Sheep-Milk-Production Traits.

Natural selection and domestication have shaped modern sheep populations into a vast range of phenotypically diverse breeds. Among these breeds, dairy sheep have a smaller population than meat sheep and wool sheep, and less research is performed on them, but the lactation mechanism in dairy sheep is critically important for improving animal-production methods. In this study, whole-genome sequences were generated from 10 sheep breeds, including 57 high-milk-yield sheep and 44 low-milk-yield sheep, to investigate the genetic signatures of milk production in dairy sheep, and 59,864,820 valid SNPs (Single Nucleotide Polymorphisms) were kept after quality control to perform population-genetic-structure analyses, gene-detection analyses, and gene-function-validation analyses. For the population-genetic-structure analyses, we carried out PCA (Principal Component Analysis), as well as neighbor-joining tree and structure analyses to classify different sheep populations. The sheep used in our study were well distributed in ten groups, with the high-milk-yield-group populations close to each other and the low-milk-yield-group populations showing similar classifications. To perform an exact signal-selection analysis, we used three different methods to find SNPs to perform gene-annotation analyses within the 995 common regions derived from the fixation index (FST), nucleotide diversity (Ɵπ), and heterozygosity rate (ZHp) results. In total, we found 553 genes that were located in these regions. These genes mainly participate in the protein-binding pathway and the nucleoplasm-interaction pathway, as revealed by the GO- and KEGG-function-enrichment analyses. After the gene selection and function analyses, we found that FCGR3A, CTSK, CTSS, ARNT, GHR, SLC29A4, ROR1, and TNRC18 were potentially related to sheep-milk-production traits. We chose the strongly selected genes, FCGR3A, CTSK, CTSS, and ARNT during the signal-selection analysis to perform a RT-qPCR (Reale time Quantitative Polymerase Chain Reaction) experiment to validate their expression-level relationship with milk production, and the results showed that FCGR3A has a significant negative relationship with sheep-milk production, while other three genes did not show any positive or negative relations. In this study, it was discovered and proven that the candidate gene FCGR3A potentially contributes to the milk production of dairy sheep and a basis was laid for the further study of the genetic mechanism underlying the strong milk-production traits of sheep.

Transcriptome Profiling of Developing Ovine Fat Tail Tissue Reveals an Important Role for MTFP1 in Regulation of Adipogenesis.

Fat-tail sheep exhibit a unique trait whereby substantial adipose tissue accumulates in the tail, a phenotype that is advantageous in many agroecological environments. In this study, we conducted histological assays, transcriptome analysis and functional assays to examine morphogenesis, characterize gene expression, and elucidate mechanisms that regulate fat tail development. We obtained the microstructure of tail before and after fat deposition, and demonstrated that measurable fat deposition occurred by the 80-day embryo (E80) stage, earlier than other tissues. Transcriptome profiling revealed 1,058 differentially expressed genes (DEGs) with six markedly different expression trends. GSEA enrichment and other downstream analyses showed important roles for genes and pathways involving in metabolism and that mitochondrial components were specifically overexpressed in the fat tail tissue of the 70-day embryo (E70). One hundred and eighty-three genes were further identified by leading edge gene analysis, among which, 17 genes have been reported in previous studies, including EEF1D, MTFP1, PPP1CA, PDGFD. Notably, the MTFP1 gene was highly correlated with the expression of other genes and with the highest enrichment score and gene expression change. Knockdown of MTFP1 in isolated adipose derived stem cells (ADSCs) inhibited cell proliferation and migration ability, besides, promoted the process of adipogenesis in vitro.

A study using single-locus and multi-locus genome-wide association study to identify genes associated with teat number in Hu sheep.

The multiple teats trait is common in many species of mammals and is considered related to lactation ability in swine. However, in Hu sheep, related gene research is still relatively limited. In this study, a genome-wide association study was used to identify genetic markers and genes related to the number of teats in the Hu sheep population, a native Chinese sheep breed. A single marker method and several multi-locus methods were utilized. A total of 61 SNPs were found to be related to the number of teats. Among these, 11 SNPs and one SNP were consistently detected by two and three multi-locus models respectively. Four SNPs were concordantly identified between the single marker and multi-locus methods. We also performed quantitative real-time PCR testing of these identified candidate genes, identifying three genes with significantly different expression. Our study suggested that the LHFP, DPYSL2, and TDP-43 genes may be related to the number of teats in sheep. The combination of single and multi-locus GWAS detected additional SNPs not found with only one model. Our results provide new and important insights into the genetic mechanisms of the mammalian multiparous teat phenotype. These findings may be useful for future breeding and understanding the genetics of sheep and other livestock.

Expression, structure and function analysis of the sperm-oocyte fusion genes Juno and Izumo1 in sheep (Ovis aries).

BACKGROUND: JUNO and IZUMO1 are the first receptor-ligand protein pairs discovered to be essential for sperm-oocyte fusion; their interaction is indispensable for fertilization. METHODS: PCR was used to clone the full-length DNA sequence of the Juno gene in sheep. The single nucleotide polymorphism (SNP) loci of Juno were genotyped by Sequenom MassARRAY®. PCR combined with rapid amplification of cDNA Ends were used to clone the full-length cDNA sequence of Juno and Izumo1. Reverse transcriptase-PCR (RT-PCR) and real time-quantitative-PCR (RT-qPCR) were used to analyze the genes' expression in tissues of sheep, and single cell RNA-seq was used to analyze the genes' expression in oocytes, granulosa cells and follicular theca of polytocous and monotocous Small Tail Han ewes. Bioinformatics was used to analyze advanced structure and phylogeny of JUNO and IZUMO1 proteins. RESULTS: The full-length DNA sequence of the Juno gene in sheep was cloned and nine SNPs were screened. We found a significant association between the g.848253 C > A locus of Juno and litter size of Small Tail Han sheep (P < 0.05). The full-length cDNA sequence of Juno and Izumo1 genes from Small Tail Han sheep were obtained. We found a new segment of the Izumo1 CDS consisting of 35 bp, and we confirmed the Izumo1 gene has 9 exons, not 8. RT-qPCR showed that Juno and Izumo1 genes were highly expressed in ovarian and testicular tissues, respectively (P < 0.01). Single cell RNA-seq showed Juno was specifically expressed in oocytes, but not in granulosa cells or follicular theca, while Izumo1 displayed little to no expression in all three cell types. There was no difference in expression of the Juno gene in oocyte and ovarian tissue in sheep with different litter sizes, indicating expression of Juno is not related to litter size traits. Bioinformatic analysis revealed the g.848253 C > A locus of Juno results in a nonconservative missense point mutation leading to a change from Phe to Leu at position 219 in the amino acid sequence. CONCLUSIONS: For the first time, this study systematically analyzed the expression, structure and function of Juno and Izumo1 genes and their encoded proteins in Small Tail Han sheep, providing the basis for future studies of the regulatory mechanisms of Juno and Izumo1 genes.

Correction: Genomic mapping identifies two genetic variants in the MC1R gene for coat colour variation in Chinese Tan sheep.

[This corrects the article DOI: 10.1371/journal.pone.0235426.].

Genomic mapping identifies two genetic variants in the MC1R gene for coat colour variation in Chinese Tan sheep.

Coat colour is one of the most important economic traits of sheep and is mainly used for breed identification and characterization. This trait is determined by the biochemical function, availability and distribution of phaeomelanin and eumelanin pigments. In our study, we conducted a genome-wide association study to identify candidate genes and genetic variants associated with coat colour in 75 Chinese Tan sheep using the ovine 600K SNP BeadChip. Accordingly, we identified two significant SNPs (rs409651063 at 14.232 Mb and rs408511664 at 14.228 Mb) associated with coat colour in the MC1R gene on chromosome 14 with -log10(P) = 2.47E-14 and 1.00E-13, respectively. The consequence of rs409651063 was a missense variant (g.14231948 G>A) that caused an amino acid change (Asp105Asn); however, the second SNP (rs408511664) was a synonymous substitution and is an upstream variant (g.14228343G>A). Moreover, our PCR analysis revealed that the genotype of white sheep was exclusively homozygous (GG), whereas the genotypes of black-head sheep were mainly heterozygous (GA). Interestingly, allele-specific expression analysis (using the missense variant for the skin cDNA samples from black-head sheep) revealed that only the G allele was expressed in the skin covered with white hair, while both the G and A alleles were expressed in the skin covered with black hair. This finding indicated that the missense mutation that we identified is probably responsible for white coat colour in Tan sheep. Furthermore, qPCR analysis of MC1R mRNA level in the skin samples was significantly higher in black-head than white sheep and very significantly higher in GA than GG individuals. Taken together, these results help to elucidate the genetic mechanism underlying coat colour variation in Chinese indigenous sheep.

Molecular Cloning of the B4GALNT2 Gene and Its Single Nucleotide Polymorphisms Association with Litter Size in Small Tail Han Sheep.

A new fecundity gene named the FecL (mutation), which regulates the ovulation rate, was discovered in French Lacaune sheep. The B4GALNT2 (beta-1, 4-N-acetyl-galactosaminyl transferase 2) gene was considered as the potential FecL mutation gene. This study explores whether the effect of the FecL mutation exists in other sheep breeds, and the features of the B4GALNT2 gene in terms of the molecular structure and its expression profile. Using Sanger sequencing, we found that high and low fecundity breeds from among 11 measured sheep breeds all had no variation in the three specific mutation sites, which were linked with the FecL mutation. However, two mutations of g.36946470C > T and g.36933082C > T in the exon of B4GALNT2 had a significant effect on litter size in the first parity for Small Tail Han (STH) Sheep (p < 0.05). Two transcription start sites (TSS) of B4GALNT2 in its 5'-flanking region were discovered in ovine granule cells in vitro, through the RACE (Rapid amplification of cDNA ends) method. Except for in the kidney and oviduct, no significant difference in expression levels had been found between STH sheep and Tan sheep breeds. The B4GALNT2 gene, as a candidate for FecL, may have a relationship with the differences in litter size in STH sheep. B4GALNT2 is mainly expressed in the ovine ovary, which also suggests that B4GALNT2 plays an important role in sheep reproduction.

Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization.

Background: Animal domestication has been extensively studied, but the process of feralization remains poorly understood. Results: Here, we performed whole-genome sequencing of 99 sheep and identified a primary genetic divergence between 2 heterogeneous populations in the Tibetan Plateau, including 1 semi-feral lineage. Selective sweep and candidate gene analysis revealed local adaptations of these sheep associated with sensory perception, muscle strength, eating habit, mating process, and aggressive behavior. In particular, a horn-related gene, RXFP2, showed signs of rapid evolution specifically in the semi-feral breeds. A unique haplotype and repressed horn-related tissue expression of RXFP2 were correlated with higher horn length, as well as spiral and horizontally extended horn shape. Conclusions: Semi-feralization has an extensive impact on diverse phenotypic traits of sheep. By acquiring features like those of their wild ancestors, semi-feral sheep were able to regain fitness while in frequent contact with wild surroundings and rare human interventions. This study provides a new insight into the evolution of domestic animals when human interventions are no longer dominant.