Chinese Indigenous Sheep Research Articles

Population structure and breed identification of Chinese indigenous sheep breeds using whole genome SNPs and InDels

BackgroundAccurate breed identification is essential for the conservation and sustainable use of indigenous farm animal genetic resources. In this study, we evaluated the phylogenetic relationships and genomic breed compositions of 13 sheep breeds using SNP and InDel data from whole genome sequencing. The breeds included 11 Chinese indigenous and 2 foreign commercial breeds. We compared different strategies for breed identification with respect to different marker types, i.e. SNPs, InDels, and a combination of SNPs and InDels (named SIs), different breed-informative marker detection methods, and different machine learning classification methods.ResultsUsing WGS-based SNPs and InDels, we revealed the phylogenetic relationships between 11 Chinese indigenous and two foreign sheep breeds and quantified their purities through estimated genomic breed compositions. We found that the optimal strategy for identifying these breeds was the combination of DFI_union for breed-informative marker detection, which integrated the methods of Delta, Pairwise Wright's FST, and Informativeness for Assignment (namely DFI) by merging the breed-informative markers derived from the three methods, and KSR for breed assignment, which integrated the methods of K-Nearest Neighbor, Support Vector Machine, and Random Forest (namely KSR) by intersecting their results. Using SI markers improved the identification accuracy compared to using SNPs or InDels alone. We achieved accuracies over 97.5% when using at least the 1000 most breed-informative (MBI) SI markers and even 100% when using 5000 SI markers.ConclusionsOur results provide not only an important foundation for conservation of these Chinese local sheep breeds, but also general approaches for breed identification of indigenous farm animal breeds.

Characterization of Flavor Compounds in Chinese Indigenous Sheep Breeds Using Gas Chromatography-Ion Mobility Spectrometry and Chemometrics.

This study analyzed the flavor compounds in the meat of four indigenous breeds of Chinese sheep through the use of gas chromatography-ion mobility spectrometry (GC-IMS). GC-IMS provided information on the characteristics and strength of 71 volatile flavor compounds (monomers and dimers), with aldehydes, alcohols and ketones being the most abundant in all types of sheep meat. The compounds with higher intensity peaks in the sheep meat were aldehydes (n-nonanal, octanal, heptanal, 3-methylbutanal, and hexanal), alcohols (1-octen-3-ol, hexanol, and pentanol), ketones (3-hydroxy-2-butanon, 2-butanone, and 2-propanone), esters (methyl benzoate), and thiazole (trimethylthiazole). The volatile flavor components in the meat of the different breeds of sheep obtained via GC-IMS were further differentiated using principal component analysis. In addition, orthogonal partial least squares discriminant analysis (OPLS-DA) and variable importance on projection (VIP) were used to determine the characteristic flavor compounds in the meats of different breeds of sheep, and 21 differentially volatile components were screened out based on having a VIP above 1. These results indicate that GC-IMS combined with multivariate analysis is a convenient and powerful method for characterizing and discriminating sheep meat.

Sheep litter size heredity basis using genome-wide selective analysis.

Sheep are important herbivorous domestic animal globally, and the Chinese indigenous sheep breed has a multitude of economically significant variations due to the diverse geographical and ecological conditions. In particular, certain native breeds exhibit a visible high litter size phenotype due to the selection pressure of natural and artificial for thousands of years, offering an ideal animal model for investigating sheep's fecundity. In this study, selective signal analysis was performed on public whole-genome sequencing data from 60 sheep across eight breeds to identify candidate genes related to litter size. Results revealed that a total of 34,065,017 single-nucleotide polymorphisms (SNPs) were identified from all sheep, and 65 candidate genes (CDGs) were pinpointed from the top 1% of interacted windows and SNPs between the pairwise fixation index (FST, >0.149543) and cross-population extended haplotype homozygosity (XP-EHH, >0.701551). A total of 41 CDGs (e.g. VRTN, EYA2 and MCPH1) were annotated to 576 GO terms, of which seven terms were directly linked to follicular and embryonic development (e.g. TBXT, BMPR1B, and BMP2). In addition, 73 KEGG pathways were enriched by 21 CDGs (e.g. ENTPD5, ABCD4 and RXFP2), mainly related to Hippo (TCF4, BMPR1B and BMP2), TGF-β (BMPR1B and BMP2), PI3K-Akt (ITGB4, IL4R and PPP2R5A) and Jak-STAT signalling pathways (IL20RA and IL4R). Notably, a series of CDGs was under strong selection in sheep with high litter size traits. These findings result could improve the comprehension of the genetic underpinnings of sheep litter size. Furthermore, it provides valuable CDGS for future molecular breeding.

Whole-genome resequencing of Chinese indigenous sheep provides insight into the genetic basis underlying climate adaptation

BackgroundChinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations.ResultsBased on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability.ConclusionsAnalysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep.

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Characterization and discrimination of the flavor profiles of Chinese indigenous sheep breeds via electronic sensory, smart instruments and chemometrics

This study was conducted to analyze the flavor profile of six commercially available Chinese indigenous breeds of sheep meat and to distinguish them based on their flavor characteristics. Various techniques such as gas chromatography-mass spectrometry (GC-MS), electronic nose (E-nose), high performance liquid chromatography (HPLC) and electronic tongue (E-tongue) were employed to analyze both volatile and non-volatile compounds in the meat samples. The results showed that, a total of 28 volatile compounds in sheep meat were detected, with 1-octen-3-ol, 1-pentanol, hexanal, octanal, and nonanal identified as the major contributors to its aroma. Tan Sheep and Oula Sheep had significantly higher concentrations of 1-pentanol, nonanal, 1-octen-3-ol, 1-octanol, and hexanoic acid than other breeds. Sweet amino acids were the most abundant free amino acids, comprising 72–85% of the total FAAs, and inosinic acid was the most abundant ATP-related compound. Small-Tailed Han sheep had significantly higher levels of SFAAs, total FAAs, and AMP, while Tan sheep had significantly higher levels of UFAAs, and Oula Sheep had significantly higher levels of GMP. The response values of differential electronic sensors W5S, W1W, NMS, PKS, and SCS were found to be useful in discriminating between breeds of sheep meat. The fusion electronic sensory signals were used to train support vector machine and random forest models, both achieving 100% accuracy in discriminating between sheep breeds.

Genome-wide DNA methylation analysis reveals different methylation patterns in Chinese indigenous sheep with different type of tail.

The study was aimed to analyze the difference of genome-wide DNA differential methylation in Lanzhou Large-tailed sheep, Altay sheep and Tibetan sheep, which the typical breeds with different type tails, as to screen the differentially methylated genes (DMGs) that affect the type of tails. In this study, three Lanzhou Large-tailed sheep, three Altay sheep and three Tibetan sheep were detected by whole genome bisulfite sequencing (WGBS). The degree of genome-wide DNA methylation, differentially methylated regions (DMRs) and DMGs were analyzed. The candidate genes affecting the tail type of sheep were identified by GO and KEGG pathway enrichment analysis of DMGs. we identified 68,603 different methylated regions (DMCs) and 75 differentially methylated genes (DMGs) associated with these DMCs. Functional analysis showed that these DMGs were mainly enriched in biological process, cellular component and molecular function, Some of the genes in these pathways are involved in fat metabolism: NFATC4, LPIN2, MGAT2 and MAT2B. Our results may help to further understand the epigenetic regulation mechanisms of deposition of fat in the tail of sheep and provide new basic data for the study of local sheep.

Genetic diversity and selection of Tibetan sheep breeds revealed by whole-genome resequencing.

This study aimed to elucidate the underlying gene regions responsible for productive, phenotypic or adaptive traits in different ecological types of Tibetan sheep and the discovery of important genes encoding valuable traits. We used whole-genome resequencing to explore the genetic relationships, phylogenetic tree, and population genetic structure analysis. In addition, we identified 28 representative Tibetan sheep single-nucleotide polymorphisms (SNPs) and genomic selective sweep regions with different traits in Tibetan sheep by fixation index (Fst) and the nucleotide diversity (θπ) ratio. The genetic relationships analysis showed that each breed partitioned into its own clades and had close genetic relationships. We also identified many potential breed-specific selective sweep regions, including genes associated with hypoxic adaptability (MTOR, TRHDE, PDK1, PTPN9, TMTC2, SOX9, EPAS1, PDGFD, SOCS3, TGFBR3), coat color (MITF, MC1R, ERCC2, TCF25, ITCH, TYR, RALY, KIT), wool traits (COL4A2, ERC2, NOTCH2, ROCK1, FGF5, SOX9), and horn phenotypes (RXFP2). In particular, a horn-related gene, RXFP2, showed the four most significantly associated SNP loci (g. 29481646 A>G, g. 29469024 T>C, g. 29462010 C>T, g. 29461968 C>T) and haplotypes. This finding demonstrates the potential for genetic markers in future molecular breeding programs to improve selection for horn phenotypes. The results will facilitate the understanding of the genetic basis of production and adaptive unique traits in Chinese indigenous Tibetan sheep taxa and offer a reference for the molecular breeding of Tibetan sheep.

Whole-genome resequencing reveals domestication and signatures of selection in Ujimqin, Sunit, and Wu Ranke Mongolian sheep breeds.

ObjectiveThe current study aimed to perform whole-genome resequencing of Chinese indigenous Mongolian sheep breeds including Ujimqin, Sunit, and Wu Ranke sheep breeds (UJMQ, SNT, WRK) and deeply analyze genetic variation, population structure, domestication, and selection for domestication traits among these Mongolian sheep breeds.MethodsBlood samples were collected from a total of 60 individuals comprising 20 WRK, 20 UJMQ, and 20 SNT. For genome sequencing, about 1.5 μg of genomic DNA was used for library construction with an insert size of about 350 bp. Pair-end sequencing were performed on Illumina NovaSeq platform, with the read length of 150 bp at each end. We then investigated the domestication and signatures of selection in these sheep breeds.ResultsAccording to the population and demographic analyses, WRK and SNT populations were very similar, which were different from UJMQ populations. Genome wide association study identified 468 and 779 significant loci from SNT vs UJMQ, and UJMQ vs WRK, respectively. However, only 3 loci were identified from SNT vs WRK. Genomic comparison and selective sweep analysis among these sheep breeds suggested that genes associated with regulation of secretion, metabolic pathways including estrogen metabolism and amino acid metabolism, and neuron development have undergone strong selection during domestication.ConclusionOur findings will facilitate the understanding of Chinese indigenous Mongolian sheep breeds domestication and selection for complex traits and provide a valuable genomic resource for future studies of sheep and other domestic animal breeding.

Whole-genome resequencing reveals molecular imprints of anthropogenic and natural selection in wild and domesticated sheep

The abundance of domesticated sheep varieties and phenotypes is largely the result of long-term natural and artificial selection. However, there is limited information regarding the genetic mechanisms underlying phenotypic variation induced by the domestication and improvement of sheep. In this study, to explore genomic diversity and selective regions at the genome level, we sequenced the genomes of 100 sheep across 10 breeds and combined these results with publicly available genomic data from 225 individuals, including improved breeds, Chinese indigenous breeds, African indigenous breeds, and their Asian mouflon ancestor. Based on population structure, the domesticated sheep formed a monophyletic group, while the Chinese indigenous sheep showed a clear geographical distribution trend. Comparative genomic analysis of domestication identified several selective signatures, including IFI44 and IFI44L genes and PANK2 and RNF24 genes, associated with immune response and visual function. Population genomic analysis of improvement demonstrated that candidate genes of selected regions were mainly associated with pigmentation, energy metabolism, and growth development. Furthermore, the IFI44 and IFI44L genes showed a common selection signature in the genomes of 30 domesticated sheep breeds. The IFI44 c. 54413058 C>G mutation was selected for genotyping and population genetic validation. Results showed that the IFI44 polymorphism was significantly associated with partial immune traits. Our findings identified the population genetic basis of domesticated sheep at the whole-genome level, providing theoretical insights into the molecular mechanism underlying breed characteristics and phenotypic changes during sheep domestication and improvement.

Isobaric Tags for Relative and Absolute Quantification-Based Proteomics Reveals Candidate Proteins of Fat Deposition in Chinese Indigenous Sheep With Morphologically Different Tails.

Background: Chinese indigenous sheep can be classified into two types according to their tail morphology: fat-rumped and thin-tailed sheep, of which the typical breeds are Altay sheep and Tibetan sheep, respectively.Methods: To identify the differentially expressed proteins (DEPs) underlying the phenotypic differences between tail types, we used isobaric tags for relative and absolute quantification (iTRAQ) combined with multi-dimensional liquid chromatography tandem-mass spectrometry (LC-MS/MS) technology to detect candidate proteins. We then subjected these to a database search and identified the DEPs. Finally, bioinformatics technology was used to carry out Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.Results: A total of 3,248 proteins were identified, of which 44 were up-regulated and 40 were down-regulated DEPs. Analyzing their GO function terms and KEGG pathways revealed that the functions of these DEPs are mainly binding, catalytic activity, structural molecule activity, molecular function regulator, and transporter activity. Among the genes encoding the DEPs, APOA2, GALK1, ADIPOQ, and NDUFS4 are associated with fat formation and metabolism.Conclusion: The APOA2, GALK1, ADIPOQ, and NDUFS4 genes may be involved in the deposition of fat in the tail of sheep. This study provides a scientific basis for the breeding of thin-tailed sheep.

Polymorphisms in BMPRIB gene affect litter size in Chinese indigenous sheep breed

The BMPRIB gene belongs to the TGF-β superfamily and is considered to be a regulator of sheep reproductive performance. Single nucleotide polymorphisms (SNPs) of BMPRIB gene in the Small Tail Han, Hu, Mongolian, Oula, Gansu Alpine Fine-wool, Dorper and Australian White sheep were detected by Sanger sequencing. Five SNPs (rs427897187 G > A, rs418841713 A > G, rs159952533 T > C, rs429416173 C > A and rs403555643 A > G) of BMPRIB gene were identified. For rs427897187 G > A, further analysis revealed that genotype GG and GA had 0.26 (p < 0.05) and 0.33 (p < 0.05) litter size less than those with genotype AA in Oula sheep. For rs403555643 A > G, further analysis revealed that genotype GG and AG had 0.65 (p < 0.05) and 0.38 (p < 0.05) litter size more than those with genotype AA in Oula sheep, and genotype GG had 0.56 (p < 0.05) litter size more than those with genotype AA in Mongolian sheep. The results showed that rs427897187 G > A and rs403555643 A > G are potential molecular markers wich could improve litter size of Chinese indigenous sheep and be used in Chinese indigenous sheep breeding.

Estimates of genomic inbreeding and identification of candidate regions that differ between Chinese indigenous sheep breeds

BackgroundA run of homozygosity (ROH) is a consecutive tract of homozygous genotypes in an individual that indicates it has inherited the same ancestral haplotype from both parents. Genomic inbreeding can be quantified based on ROH. Genomic regions enriched with ROH may be indicative of selection sweeps and are known as ROH islands. We carried out ROH analyses in five Chinese indigenous sheep breeds; Altay sheep (n = 50 individuals), Large-tailed Han sheep (n = 50), Hulun Buir sheep (n = 150), Short-tailed grassland sheep (n = 150), and Tibetan sheep (n = 50), using genotypes from an Ovine Infinium HD SNP BeadChip.ResultsA total of 18,288 ROH were identified. The average number of ROH per individual across the five sheep breeds ranged from 39 (Hulun Buir sheep) to 78 (Large-tailed Han sheep) and the average length of ROH ranged from 0.929 Mb (Hulun Buir sheep) to 2.544 Mb (Large-tailed Han sheep). The effective population size (Ne) of Altay sheep, Large-tailed Han sheep, Hulun Buir sheep, Short-tailed grassland sheep and Tibetan sheep were estimated to be 81, 78, 253, 238 and 70 five generations ago. The highest ROH-based inbreeding estimate (FROH) was 0.0808 in Large-tailed Han sheep, whereas the lowest FROH was 0.0148 in Hulun Buir sheep. Furthermore, the highest proportion of long ROH fragments (> 5 Mb) was observed in the Large-tailed Han sheep breed which indicated recent inbreeding. In total, 49 ROH islands (the top 0.1% of the SNPs most commonly observed in ROH) were identified in the five sheep breeds. Three ROH islands were common to all the five sheep breeds, and were located on OAR2: 12.2–12.3 Mb, OAR12: 78.4–79.1 Mb and OAR13: 53.0–53.6 Mb. Three breed-specific ROH islands were observed in Altay sheep (OAR15: 3.4–3.8 Mb), Large-tailed Han sheep (ORA17: 53.5–53.8 Mb) and Tibetan sheep (ORA5:19.8–20.2 Mb). Collectively, the ROH islands harbored 78 unique genes, including 19 genes that have been documented as having associations with tail types, adaptation, growth, body size, reproduction or immune response.ConclusionDifferent ROH patterns were observed in five Chinese indigenous sheep breeds, which reflected their different population histories. Large-tailed Han sheep had the highest genomic inbreeding coefficients and the highest proportion of long ROH fragments indicating recent inbreeding. Candidate genes in ROH islands could be used to illustrate the genetic characteristics of these five sheep breeds. Our findings contribute to the understanding of genetic diversity and population demography, and help design and implement breeding and conservation strategies for Chinese sheep.

Variations in the insulin receptor substrate 1 (IRS1) and its association with growth traits in Chinese black Tibetan sheep (Ovis aries)

This paper presented the results on the study of insulin receptor substrate 1 (IRS1) polymorphisms in Chinese black Tibetan sheep. Via DNA direct sequencing, four variations within 3′ untranslated region (UTR) of IRS1, including g.9382T > G, g.9569T > G, g.9572C > T and g.9695A > C were detected in the black Tibetan sheep population. Based on the χ 2 test, those four loci deviated from Hardy–Weinberg equilibrium (p < 0.05). In g.9569T > G locus, genotype of GG possessed advantage on body weight (p < 0.05). In g.9572C > T locus, individuals with genotype of TT homozygous mutation decreased significantly on body weight, withers height, body length and chest circumference (p < 0.05 or p < 0.01). In g.9695A > C locus, the body weight and chest circumference were also higher in AA carriers than in CC carriers (p < 0.05). Our results provided evidence that polymorphisms in IRS1 were associated with growth efficiency traits by quantitative genetic analysis, and may be used for marker-assisted selection in Chinese indigenous sheep.

Genome wide association study for the identification of genes associated with tail fat deposition in Chinese sheep breeds.

ABSTRACTChinese indigenous sheep can be classified into three types based on tail morphology: fat-tailed, fat-rumped, and thin-tailed sheep, of which the typical breeds are large-tailed Han sheep, Altay sheep, and Tibetan sheep, respectively. To unravel the molecular genetic basis underlying the phenotypic differences among Chinese indigenous sheep with these three different tail types, we used ovine high-density 600K single nucleotide polymorphism (SNP) arrays to detect genome-wide associations, and performed general linear model analysis to identify candidate genes, using genotyping technology to validate the candidate genes. Tail type is an important economic trait in sheep. However, the candidate genes associated with tail type are not known. The objective of this study was to identify SNP markers, genes, and chromosomal regions related to tail traits. We performed a genome-wide association study (GWAS) using data from 40 large-tailed Han sheep, 40 Altay sheep (cases) and 40 Tibetan sheep (controls). A total of 31 significant (P<0.05) SNPs associated with tail-type traits were detected. For significant SNPs' loci, we determined their physical location and performed a screening of candidate genes within each region. By combining information from previously reported and annotated biological functional genes, we identified SPAG17, Tbx15, VRTN, NPC2, BMP2 and PDGFD as the most promising candidate genes for tail-type traits. Based on the above identified candidate genes for tail-type traits, BMP2 and PDGFD genes were selected to investigate the relationship between SNPs within the tails in the Altay and Tibetan populations. rs119 T>C in exon1 of the BMP2 gene and one SNP in exon4 (rs69 C>A) of the PDGFD gene were detected. rs119 was of the TT genotype in Altay sheep, while it was of the CC genotype in Tibetan sheep. On rs69 of the PDGFD gene, Altay sheep presented with the CC genotype; however, Tibetan sheep presented with the AA genotype.

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

BackgroundCopy 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.ResultsA 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.ConclusionsIn 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.

Genome-Wide Analysis Revealed Homozygosity and Demographic History of Five Chinese Sheep Breeds Adapted to Different Environments.

Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1–5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits (5S_rRNA), the innate and adaptive immune response (HERC2 and CYFIP1), digestion and metabolism (CENPJ), growth (SPP1), body size and developments (GJB2 and GJA3). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.

Genetic variability analysis of 26 sheep breeds in the Czech Republic

Genetic variability analysis of 26 sheep breeds in the Czech Republic

Comparative transcriptome of reproductive axis in Chinese indigenous sheep with different FecB genotypes and prolificacies

In the present study, there were 128, 263, and 139 differentially expressed genes (DEGs) from the ovary, pituitary, and hypothalamus, respectively, with the FecB B+/FecB ++ genotype. Results indicated that there was a synchronized, and differential increase in mRNA transcript abundances for genes encoding for ovarian steroidogenesis-signaling proteins in the ovary, TGF-β-signaling in hypothalamus, cAMP-signaling-transduction, and dopaminergic-synapse pathways in the pituitary. The values for these variables were associated with mean prolificacy in indigenous sheep and results provided indications for candidate regulator proteins that coordinate processes related to prolificacy in ewes. Furthermore, the results when there was evaluation of the FecB-genotype-specific co-expression modules implied that CYP17 was the hub gene connecting the pathways and for modulation of folliculogenesis and ovulation in the ewe with the FecB B+ genotype.

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.