Genetic characterization and implications for conservation of the last autochthonous Mouflon population in Europe.

Population genetic studies provide accurate information on population structure, connectivity, and hybridization. These are key elements to identify units for conservation and define wildlife management strategies aimed to maintain and restore biodiversity. The Mediterranean island of Sardinia hosts one of the last autochthonous mouflon populations, descending from the wild Neolithic ancestor. The first mouflon arrived in Sardinia ~ 7000 years ago and thrived across the island until the twentieth century, when anthropogenic factors led to population fragmentation. We analysed the three main allopatric Sardinian mouflon sub-populations, namely: the native sub-populations of Montes Forest and Mount Tonneri, and the reintroduced sub-population of Mount Lerno. We investigated the spatial genetic structure of the Sardinian mouflon based on the parallel analysis of 14 highly polymorphic microsatellite loci and mitochondrial D-loop sequences. The Montes Forest sub-population was found to harbour the ancestral haplotype in the phylogeny of European mouflon. We detected high levels of relatedness in all the sub-populations and a mitochondrial signature of hybridization between the Mount Lerno sub-population and domestic sheep. Our findings provide useful insights to protect such an invaluable genetic heritage from the risk of genetic depletion by promoting controlled inter-population exchange and drawing informed repopulation plans sourcing from genetically pure mouflon stocks.

An epidemiological study of gastrointestinal nematode and Eimeria coccidia infections in different populations of Kazakh sheep.

This is an epidemiological study on the gastrointestinal nematode (GIN) and Eimeria coccidia infections in Kazakh sheep and the F1 and F2 generations of Kazakh × Texel sheep crosses. A total of 7599 sheep fecal samples were collected from the Zhaosu County and Nilka County in Ili Kazakh Autonomous Prefecture in the four seasons-spring, summer, autumn, and winter of 2019. The parasite causing the infection was identified by the saturated saline floating method, and the infection intensity was calculated by the modified McMaster method. SPSS19.0 was used to evaluate the differences in the fecal egg count (FEC) of for GIN and the fecal oocyst count (FOC) value of for coccidia per sample. The results showed that there were nine types of sheep GIN infections and Eimeria coccidia in these two counties of Ililocations, with the dominant parasite species of Haemonchus contortus, Trichostrongylus spp., and Ostertagia spp as the predominant parasites in the sheep. Most of the GIN and coccidia infections in these two regions were mild and moderate. The mean log (FEC) of GIN infection in the Zhaosu area was significantly higher than that in the Nilka area, whereas the mean log (FOC) of coccidia infection in Zhaosu was significantly lower than that of Nilka. The mean log (FEC) of GIN infection in the four seasons was the highest in spring, followed by in summer, then in autumn, and the lowest in winter. The mean log (FOC) of coccidia infection was the highest in spring, followed by in autumn, and was the lowest in summer and winter. The mean log (FEC) of GIN infection and log (FOC) of coccidia infection of Kazakh sheep was significantly higher than the F1 generation, which was then significantly higher than the F2 generation of summer. A positive correlation was found between the EPG and OPG levels in the sheep. These results showed that the GIN and coccidia infection intensities of the F1 generation sheep of Kazakh ×Texel crosses were significantly lower than that of Kazakh sheep paving the way for marker-based resistance selection.

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

Chinese 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.

Shared Ancestry and Signatures of Recent Selection in Gotland Sheep.

Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation FST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the FST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.

Paternal Origins and Migratory Episodes of Domestic Sheep.

The domestication and subsequent global dispersal of livestock are crucial events in human history, but the migratory episodes during the history of livestock remain poorly documented [1-3]. Here, we first developed a set of 493 novel ovine SNPs of the male-specific region of Y chromosome (MSY) by genome mapping. We then conducted a comprehensive genomic analysis of Y chromosome, mitochondrial DNA, and whole-genome sequence variations in a large number of 595 rams representing 118 domestic populations across the world. We detected four different paternal lineages of domestic sheep and resolved, at the global level, their paternal origins and differentiation. In Northern European breeds, several of which have retained primitive traits (e.g., a small body size and short or thin tails), and fat-tailed sheep, we found an overrepresentation of MSY lineages y-HC and y-HB, respectively. Using an approximate Bayesian computation approach, we reconstruct the demographic expansions associated with the segregation of primitive and fat-tailed phenotypes. These results together with archaeological evidence and historical data suggested the first expansion of early domestic hair sheep and the later expansion of fat-tailed sheep occurred ∼11,800-9,000 years BP and ∼5,300-1,700 years BP, respectively. These findings provide important insights into the history of migration and pastoralism of sheep across the Old World, which was associated with different breeding goals during the Neolithic agricultural revolution.

Vascular morphology and expression of angiogenic factors in sheep ovaries and their relation with the lambing rate.

The study aims to present the vascular anatomy of ovarian arterial system in sheep and to explore its relation with the angiogenic factors vascular endothelial growth factor (VEGF) and angiopoietin-2 (Ang-2), environmental factors and the lambing rate. Small-tailed Han sheep, Hu sheep, Aohan finewool sheep and Tibetan sheep who are from different habitats in China and with different lambing rates were investigated. The vascular corrosion casting technique and scanning electron microscopy were used to understand arterial architecture, and real-time quantitative PCR (qRT-PCR) and Western blotting were used for assessing VEGF and Ang-2 expression. Results showed that high-prolificity breeds had significantly larger vessel diameters in ovaries than low-prolificity ones. VEGF gene transcript and protein expression were related with the lambing rate, while Ang-2 gene transcript was negatively related with the lambing rate. Results indicated that high-prolificity breeds have more abundant arterial blood supply, which is associated with VEGF expression.

Understanding the action of muscle iron concentration on dark cutting: An important aspect affecting consumer confidence of purchasing meat.

We investigated the association of muscle iron concentration, in addition to ultimate pH (pHU), on dark meat formation in sheep of different breeds fed forage-based diets. At 1 h simulated display, redness of meat (a*-value) increased (P < .0001) by about 3 units as the iron concentration increased from 10 to 22 mg/kg of meat, whereas the a*-value decreased by 2 units as pHU increased from 5.5 to 6.2 in fresh meat (P < .0001). After 90 days storage the corresponding responses were about 2 units increase for iron concentration and about 1 unit decrease for pHU, respectively. The results clearly show that increased muscle iron concentration was strongly associated with reduced dark cutting in fresh and stored meat evaluated at 1 h simulated display. We conclude that it may be desirable to measure iron concentration, along with pHU, for evaluation of the potential for carcasses to produce dark cutting meat, and for the meat to turn brown during display.

Prevalence and multilocus genotyping of Giardia duodenalis in Tan sheep (Ovis aries) in northwestern China.

Giardia duodenalis is a common intestinal protozoa, which can cause the occurrence of diarrhea, weight loss, and even death in animals or human, this threatens the husbandry industry and public health. It can infect virtually humans and all domestic animals including sheep. Tan sheep is one of the most important sheep breeds, which is short-tailed indigenous sheep breed used for production of high quality meat and pelts in China. However, there are no report regarding the occurrence and multilocus genotyping of G. duodenalis in Tan sheep in northwestern China. Thus, the objective of the present study was to investigate the prevalence and multilocus genotypes of G. duodenalis in Tan sheep. 1014 fecal samples were collected from Tan sheep from Ningxia Hui Autonomous Region, and three loci (ß-giardin (bg), glutamate dehydrogenase (gdh) and triosephosphate isomerase (tpi) genes) were amplified by nested PCR. The prevalence of G. duodenalis in Tan sheep was 14.5% (147/1014), two assemblages (assemblage A, n = 43; and E, n = 90) were detected, including one novel assemblage A at bg locus, one novel assemblage A at tpi locus, and 10 and 11 novel subtypes of assemblage E were detected at the bg and gdh loci, respectively. One MLGs was formed based on sequence variation among the three loci. Moreover, 9 Tan sheep were infected with two assemblages (A and E) based on the three loci. These findings expand the host range of G. duodenalis and revealed genetic diversity of G. duodenalis assemblages in Tan sheep.

Genetic and genomic analyses underpin the feasibility of concomitant genetic improvement of milk yield and mastitis resistance in dairy sheep.

Milk yield is the most important dairy sheep trait and constitutes the key genetic improvement goal via selective breeding. Mastitis is one of the most prevalent diseases, significantly impacting on animal welfare, milk yield and quality, while incurring substantial costs. Our objectives were to determine the feasibility of a concomitant genetic improvement programme for enhanced milk production and resistance to mastitis. Individual records for milk yield, and four mastitis-related traits (milk somatic cell count, California Mastitis Test score, total viable bacterial count in milk and clinical mastitis presence) were collected monthly throughout lactation for 609 ewes of the Chios breed. All ewes were genotyped with a mastitis specific custom-made 960 single nucleotide polymorphism (SNP) array. We performed targeted genomic association studies, (co)variance component estimation and pathway enrichment analysis, and characterised gene expression levels and the extent of allelic expression imbalance. Presence of heritable variation for milk yield was confirmed. There was no significant genetic correlation between milk yield and mastitis traits. Environmental factors appeared to favour both milk production and udder health. There were no overlapping of SNPs associated with mastitis resistance and milk yield in Chios sheep. Furthermore, four distinct Quantitative Trait Loci (QTLs) affecting milk yield were detected on chromosomes 2, 12, 16 and 19, in locations other than those previously identified to affect mastitis resistance. Five genes (DNAJA1, GHR, LYPLA1, NUP35 and OXCT1) located within the QTL regions were highly expressed in both the mammary gland and milk transcriptome, suggesting involvement in milk synthesis and production. Furthermore, the expression of two of these genes (NUP35 and OXCT1) was enriched in immune tissues implying a potentially pleiotropic effect or likely role in milk production during udder infection, which needs to be further elucidated in future studies. In conclusion, the absence of genetic antagonism between milk yield and mastitis resistance suggests that simultaneous genetic improvement of both traits be achievable.

Identification of two TYRP1 loss-of-function alleles in Valais Red sheep.

The Valais Red sheep breed is a local breed of the Swiss canton Valais. Although the breed is characterised by its brown colour, black animals occasionally occur and the objective of this study was to identify the causative genetic variants responsible for the obvious difference. A GWAS using high-density SNP data to compare 51 brown and 38 black sheep showed a strong signal on chromosome 2 at the TYRP1 locus. Haplotype analyses revealed three different brown-associated alleles. The WGS of three sheep revealed four protein-changing variants within the TYRP1 gene. Three of these variants were associated with the recessively inherited brown coat colour. This includes the known missense variant TYRP1:c.869G>T designated as bS oay and two novel loss-of-function variants. We propose to designate the frame-shift variant TYRP1:c.86_87delGA as bVS 1 and the nonsense variant TYRP1:c.1066C>T as bVS 2 . Interestingly, the bVS 1 allele occurs only in local breeds of Switzerland whereas the bVS 2 allele seems to be more widespread across Europe.

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations.

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short-tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D5 ) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.

Two linked TBXT (brachyury) gene polymorphisms are associated with the tailless phenotype in fat-rumped sheep.

T-box transcription factor T (TBXT), encoding the brachyury protein, is an embryonic nuclear transcription factor involved in mesoderm formation and differentiation. Previous studies indicate that TBXT mutations are responsible for the tailless or short-tailed phenotype of many vertebrates. To verify whether the tailless phenotype in fat-rumped sheep is associated with TBXT mutations, exon 2 of the TBXT gene for 301 individuals belonging to 13 Chinese and Iranian sheep breeds was directly sequenced. Meanwhile, 380 samples were used to detect the genotypes of the candidate variations by mapping to their reads databases in the Sequence Read Archive repository of GenBank. The results showed that one missense mutation, c.334G>T (GGG>TGG) with a completely linked synonymous variant c.333G>C (CCG>CCC) was found to be associated with the 'tailless' characteristic in typical fat-rumped sheep breeds. The c.334G>T transversion led to the conversion of glycine to tryptophan at the 112th amino acid in the T-box domain of the brachyury protein. In addition, crossbreeding experiments for long-tailed and tailless sheep showed that CT/CT allele of nucleotides (nt) 333 and 334, a recessive mutation, would cause sheep tails to be shorter, suggesting that these two linked variants at nucleotides 333 and 334 in TBXT are probably causative mutations responsible for the tailless phenotype in sheep.

Low genome-wide homozygosity in 11 Spanish ovine breeds.

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (FROH ) are generally small (FROH  < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.

Inferring the population structure of the Maghreb sheep breeds using a medium-density SNP chip.

North Africa has a great diversity of indigenous sheep breeds whose origin is linked to its environmental characteristics and to certain historical events that took place in the region. To date, few genome-wide studies have been conducted to investigate the population structure of North African indigenous sheep. The objective of the present study was to provide a detailed assessment of the genetic structure and admixture patterns of six Maghreb sheep populations using the Illumina 50K Ovine BeadChip and comparisons with 22 global populations of sheep and mouflon. Regardless of the method of analysis used, patterns of multiple hybridization events were observed within all North African populations, leading to a heterogeneous genetic architecture that varies according to the breed. The Barbarine population showed the lowest genetic heterogeneity and major southwest Asian ancestry, providing additional support to the Asian origin of the North African fat-tailed sheep. All other breeds presented substantial Merino introgression ranging from 15% for D'man to 31% for Black Thibar. We highlighted several signals of ancestral introgression between North African and southern European sheep. In addition, we identified two opposite gradients of ancestry, southwest Asian and central European, occurring between North Africa and central Europe. Our results provide further evidence of the weak global population structure of sheep resulting from high levels of gene flow among breeds occurring worldwide. At the regional level, signs of recent admixture among North African populations, resulting in a change of the original genomic architecture of minority breeds, were also detected.

Relationships between body reserve dynamics and rearing performances in meat ewes1.

The main objective of this work was to study the relationships between body reserve (BR) dynamics and rearing performance (PERF) traits in ewes from a Romane meat sheep flock managed extensively on "Causse" rangelands in the south of France. Flock records were used to generate data sets covering 14 lambing years (YR). The data set included 1,146 ewes with 2 ages of first lambing (AGE), 3 parities (PAR), and 4 litter sizes (LS). Repeated measurements of the BW and BCS were used as indicators of BR. The ewe PERF traits recorded were indirect measurements for maternal abilities and included prolificacy, litter weight and lamb BW at lambing and weaning, ADG at 1, 2, and 3 mo after lambing, and litter survival from lambing to weaning. The effects of different BW and BCS trajectories (e.g., changes in BW and BCS across the production cycle), previously been characterized in the same animals, on PERF traits were investigated. Such trajectories reflected different profiles at the intraflock level in the dynamics of BR mobilization-accretion cycles. Genetic relationships between BR and PERF traits were assessed. All the fixed variables considered (i.e., YR, AGE, PAR, LS, and SEX ratio of the litter) have significant effects on the PERF traits. Similarly, BW trajectories had an effect on the PERF traits across the 3 PARs studied, particularly during the first cycle (PAR 1). The BCS trajectories only affected prolificacy, lamb BW at birth, and litter survival. Most of the PERF traits considered here showed moderate heritabilities (0.17-0.23) except for prolificacy, the lamb growth rate during the third month and litter survival which showed very low heritabilities. With exception of litter survival and prolificacy, ewe PERF traits were genetically, strongly, and positively correlated with BW whatever the physiological stage. A few weak genetic correlations were found between BCS and PERF traits. As illustrated by BW and BCS changes over time, favorable genetic correlations were found, even if few and moderate, between BR accretion or mobilization and PERF traits, particularly for prolificacy and litter weight at birth. In conclusion, our results show significant relationships between BR dynamics and PERF traits in ewes, which could be considered in future sheep selection programs aiming to improve robustness.

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds.

A dataset consisting of 787 animals with high-density SNP chip genotypes (346 774 SNPs) and 939 animals with medium-density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between FPED and FROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean FROH estimates from medium-density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make-up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.

Genome-wide association studies for milk production traits in Valle del Belice sheep using repeated measures.

Genome-wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test-day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.

Detecting selection signatures in three Iranian sheep breeds.

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype-phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP-EHH (cross-population extend haplotype homozygosity) and FST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori-Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success.

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA-miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm-egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

Genetic differentiation and phylogeny of 27 sheep populations based on structural gene loci.

To explore the genetic divergence and phylogeny of Chinese indigenous sheep breeds, in the current study, we analyzed the polymorphisms of 5 structural loci in ten sheep populations, including Sishui Fur, Sunite, Wurank, Bayinbuluke, Altay, Small-Tailed Han, Wadi, Tan, Tong and Hu sheep. The data were then compared with those from an additional 13 Asian and 4 European sheep populations acquired by the same experimental method. Based on the genetic distance and the results of a cluster analysis, we constructed the phylogenetic relationship of 27 populations. The results showed that the sheep populations in this study could be classified into four genetic groups: "Mongolian", "Tibetan", "South-Southeast Asian" and "European" sheep groups. All 10 Chinese sheep breeds belonged to the "Mongolian sheep" lineage; however, Finnish Landrace sheep and Yunnan sheep could not be classified into any of the four groups. These results could provide a good reference for the protection and utilization of primary breed resources in China and phylogenic research on Asian sheep populations.