High-throughput DNA sequence analysis elucidates novel insight into the genetic basis of adaptation in local sheep.

Understanding how evolutionary factors related to climate adaptation and human selection have influenced the genetic architecture of domesticated animals is of great interest in biology. In the current study, by using 304 whole genomes from different geographical regions (including Europe, north Africa, Southwest Asia, east Asia, west Africa, south Asia, east Africa, Australia and Turkey), We evaluate global sheep population dynamics in terms of genetic variation and population structure. We further conducted comparative population analysis to study the genetic underpinnings of climate adaption to local environments and also morphological traits. In order to identify genomic signals under selection, we applied fixation index (FST) and also nucleotide diversity (θπ) statistical measurements. Our results revealed several candidate genes on different chromosomes under selection for local climate adaptation (e.g. HOXC12, HOXC13, IRF1, FGD2 and GNAQ), body size (PDGFA, HMGA2, PDE3A) and also morphological related traits (RXFP2). The discovered candidate genes may offer newel insights into genetic underpinning of regional adaptation and commercially significant features in local sheep.

Geometric morphometric evaluation of mandibles of four sheep breeds: Bardoka, Ivesi, Polish Mountain sheep and Turcana.

The enduring relationship between humans and domestic sheep has evolved over millennia, showcasing diverse uses such as meat, milk, wool, leather and fur, shaped by geographical, historical, cultural and social factors. The sheep breeds discussed include the Ivesi from Southeastern Anatolia, known for its varied animal products; the resilient Turcana breed of Romania; Kosovo's Bardoka, valued for its triple-purpose characteristics; and Poland's Polish Mountain Sheep, uniquely utilized for milk production in cheese making. Sheep, with their enduring relationship with humans and significant economic importance, have attracted scientific interest in morphometric studies of their mandibles, yielding valuable data applicable across various fields including basic anatomy, veterinary clinical anatomy, zooarchaeology and veterinary forensic medicine. Traditional morphometric studies rely on statistical methods to compare length, depth and angular ratios between anatomical formations, often highlighting differences between specific points but not fully revealing shape variations between distinct groups. Geometric morphometric analysis has emerged as a preferred method in recent years, enabling shape analyses using coordinate data from various imaging techniques, facilitating a comprehensive examination of mandibular morphometrics among sheep breeds across different countries. This study involved four sheep breeds from different countries, namely Ivesi from Turkey, Bardoka from Kosovo, Polish Mountain Sheep from Poland and Turcana from Romania, with a total of 70 mandibles sourced from various veterinary faculties. Mandibular photographs were meticulously captured, focusing on the right side of mandible pairs and placing landmarks and semi-landmarks along the entire edge, enabling geometric morphometric analysis using tpsUtil, tpsDig2 and MorphoJ software. The analysis included principal component analysis, canonical variate analysis and discriminant function analysis for pairwise comparisons, facilitating a comprehensive examination of mandibular shape variations among the different sheep breeds. Using geometric morphometric methods, this study analysed mandibles from four distinct sheep breeds sourced from different countries, revealing notable variations in regions such as the ramus mandibula, angulus mandibula and incisive areas, attributed to genetic, geographical and dietary influences, highlighting the importance of continued research to better comprehend these shape differences.

Assessing the performance of Moghani crossbred lambs derived from different mating systems with Texel and Booroola sheep.

In our ongoing project, which focuses on the introgression of Booroola/FecB gene and the myostatin (MSTN) gene into purebred Moghani sheep, we assessed the performance of second-generation Moghani crossbreds such as second crossbreds (F2) and initial backcross generation (BC1). These crossbreds were generated through different mating systems, including in-breeding, outcrossing, first paternal backcrossing (PBC1), and first maternal backcrossing (MBC1). Notably, F2 strains exhibited lean tail, woolly fleece and a higher percentage of white coat color compared to BC1. The impact of mating systems and birth types on pre-weaning survival rates was found to be statistically significant (P < 0.0001), with singleton offspring resulting from paternal backcross showing a particularly substantial effect. The F2 crossbred lambs carrying the Booroola gene did not show a statistically significant difference in survivability compared to those carrying the MSTN gene, implying the Booroola prolificacy gene had no significant impact on survival outcomes. However, the occurrence of multiple births had a significant negative impact on lamb survival (P < 0.0001). The PBC1 sheep strains, specifically Texel Tamlet ram strains carrying the MSTN mutation, exhibited superior growth rates compared to others (P < 0.05). Interestingly, the MSTN mutation in the homozygous variant genotype significantly impacts growth rate before weaning compared to other genotypes and pure Moghani sheep (P < 0.05). In conclusion, this study objectively underscores the pivotal role of genetic factors, specifically through strategic mating systems like paternal backcrossing, in enhancing desired traits and growth rates in Moghani sheep, thereby contributing valuable insights to the field of sheep breeding programs.

Identification of the FGB gene polymorphism and analysis of its association with fat deposition traits in Hu sheep.

As a crucial economic trait, fat deposition is directly related to carcass quality and feed efficiency in sheep. The purpose of this study was to investigate the polymorphisms of the FGB gene related to fat deposition and detect the expression features of the FGB gene in different adipose tissues of sheep by using Sanger sequencing, MassARRAY® SNP technique, and quantitative real-time PCR (qRT-PCR). Results showed that in the intron region of the FGB gene, a SNP g. 3378953 A > T has been identified, and significant association was found between perirenal fat weight, perirenal fat relative weight, mesenteric fat weight, and mesenteric fat relative weight (P < 0.05). Moreover, qRT-PCR analysis showed that FGB was expressed in all three adipose tissues, and FGB gene expression level in the AA genotype was significantly lower than that in the AT or TT genotypes (P < 0.05). Therefore, the FGB gene can be used as a candidate gene to reduce fat deposition in Hu sheep breeding, and the selection of the AA genotype in Hu sheep in production practice is more conducive to improving production efficiency.

Population Genomic History of the Endangered Anatolian and Cyprian Mouflons in Relation to Worldwide Wild, Feral, and Domestic Sheep Lineages.

Once widespread in their homelands, the Anatolian mouflon (Ovis gmelini anatolica) and the Cyprian mouflon (Ovis gmelini ophion) were driven to near extinction during the 20th century and are currently listed as endangered populations by the International Union for Conservation of Nature. While the exact origins of these lineages remain unclear, they have been suggested to be close relatives of domestic sheep or remnants of proto-domestic sheep. Here, we study whole genome sequences of n = 5 Anatolian mouflons and n = 10 Cyprian mouflons in terms of population history and diversity, comparing them with eight other extant sheep lineages. We find reciprocal genetic affinity between Anatolian and Cyprian mouflons and domestic sheep, higher than all other studied wild sheep genomes, including the Iranian mouflon (O. gmelini). Studying diversity indices, we detect a considerable load of short runs of homozygosity blocks (<2 Mb) in both Anatolian and Cyprian mouflons, reflecting small effective population size (Ne). Meanwhile, Ne and mutation load estimates are lower in Cyprian compared with Anatolian mouflons, suggesting the purging of recessive deleterious variants in Cyprian sheep under a small long-term Ne, possibly attributable to founder effects, island isolation, introgression from domestic lineages, or differences in their bottleneck dynamics. Expanding our analyses to worldwide wild and feral Ovis genomes, we observe varying viability metrics among different lineages and a limited consistency between viability metrics and International Union for Conservation of Nature conservation status. Factors such as recent inbreeding, introgression, and unique population dynamics may have contributed to the observed disparities.

Genome wide association studies for carcass traits measured by video image analysis in crossbred lambs.

This is the first UK genome wide association study investigating potential links between Video Image Analysis (VIA) carcass traits and molecular polymorphisms in crossbred sheep. Phenotypic and genotypic data were collected from two crossbred lamb populations: Texel x Scotch Mule (TxSM, n = 2330) and Texel x Lleyn (TxL, n = 3816). Traits measured included live weights at birth, eight weeks and weaning (∼15 weeks). VIA-predicted traits included total weights and weights of fat, muscle and bone in the whole carcass and primal (hind leg, saddle, shoulder) regions. Within-breed heritabilities estimated for the VIA traits ranged from 0.01 to 0.70, indicating potential for inclusion of some traits in breeding programmes. The two crossbred populations differed in SNPs associated with different traits. Two SNPs on chromosomes two (s74618.1) and eight (s68536.1), respectively, reached genome-wise significance for TxSM, explaining <1% of trait variance, for whole carcass fat and muscle weights, hind leg and saddle fat weights and shoulder bone weights. For TxL, four SNPs reached genome-wise significance, on chromosome two for hind leg muscle weight (OAR2_117,959,202 and OAR2_11804335), on chromosome 10 for whole carcass bone weight (OAR19_8,995,957.1), and on chromosome 19 for weaning weight (s40847.1), each explaining <1% of trait genetic variation. Differences in apparent genetic control of carcass traits may be influenced by the lambs' cross-breed, but also by management decisions affecting environmental variance and trait definitions, which should be understood in order to define protocols for incorporation of carcass traits into (cross)breeding programmes. IMPLICATIONS: Combining VIA-measured carcass traits with conventional production traits in a breeding programme could potentially improve the production and product quality of meat sheep. Phenotypes for VIA traits could be collected relatively easily if VIA machines were present at all abattoir sites. The current study and future Genome Wide Association Studies may help to identify potentially informative molecular markers, that explain large proportions of the genetic variance observed in VIA-measured carcass traits. Including this information in the estimation of breeding values could increase the accuracy of prediction, increasing the potential rate of genetic improvement for product quality. This study confirms the polygenic architecture of the investigated carcass traits, with a small number of molecular markers that each explain a small amount of genetic variation. Further studies across breed types are recommended to further test and validate molecular markers for traits related to lamb carcass quality, as measured by video image analysis.

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries).

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Asikli Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Asikli Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Asikli Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network.

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

Identification and screening of circular RNAs during adipogenic differentiation of ovine preadipocyte by RNA-seq.

Circular RNAs (circRNAs) are a class of non-coding RNAs that play important roles in preadipocyte differentiation and adipogenesis. However, little is known about genome-wide identification, expression profile, and function of circRNAs in sheep. To investigate the role of circRNAs during ovine adipogenic differentiation, the subcutaneous adipose tissue of Tibetan rams was collected in June 2022. Subsequently, the preadipocytes were immediately isolated from collected adipose tissue and then induced to begin differentiation. The adipocytes samples cultured on days 0, 2, and 8 of preadipocytes differentiation were used to perform RNA sequencing (RNA-seq) analysis to construct the expression profiles of circRNAs. Subsequently, the function of differentially expressed circRNAs was investigated by performing the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of their parent genes. Finally, a circRNAs-miRNAs-mRNAs network involved in adipogenic differentiation was been analyzed. As a result, a total of 6,449 candidate circRNAs were identified in ovine preadipocytes. Of these circRNAs identified, 63 candidate circRNAs were differentially expressed among the three differentiation stages and their parent genes were mainly enriched in acetyl-CoA metabolic process, positive regulation of lipid biosynthetic process, positive regulation of steroid biosynthetic process, and focal adhesion pathway (P < 0.05). Based on a circRNAs-miRNAs-mRNAs regulatory network constructed, circ_004977, circ_006132 and circ_003788 were found to function as competing endogenous RNAs (ceRNAs) to regulate ovine preadipocyte differentiation and lipid metabolism. The results provide an improved understanding of functions and molecular mechanisms of circRNAs underlying ovine adipogenesis in sheep.

The moderate fat deposition contributes to improve mutton quality, which is associated with the differentiation of preadipocytes. To investigate roles of circular RNAs (circRNAs) in preadipocyte differentiation, we identified circRNAs on days 0, 2, and 8 of preadipocytes differentiation and compared the expression profile of circRNAs at different adipogenic differentiation stages. A total of 6,449 candidate circRNAs were identified, among which 63 candidate circRNAs were differentially expressed among the three differentiation stages. The parent genes of differentially expressed circRNAs were enriched in several biological process and pathways related to lipid metabolism and synthesis. In addition, several circRNAs may regulate ovine preadipocyte differentiation by interacting with microRNAs (miRNAs). The results reveal the potential roles of circRNAs in adipogenic differentiation of sheep.

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep.

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

Increasing the meat productivity of young sheep based on the use of the gene pool of the Dorper and Hissar breeds.

In the pursuit of enhanced mutton production, improving the genetic reservoir of sheep with early maturation and high meat productivity is imperative. This study aims to assess the efficacy of integrating Dorper and Hissar rams into the breeding program of Kazakh fat-tailed coarse-haired ewes for generating young mutton. The research involved forming three groups, each comprising 40 ewes of the Kazakh fat-tailed coarse-haired breed, based on analog pairs. Ewes in Group I were inseminated with Dorper ram semen, those in Group II were inseminated with Hissar ram semen, and Group III served as a control with purebred Kazakh fat-tailed coarse-haired sheep breeding. Results revealed that crossbred rams in Group II achieved a significantly higher live weight of 45.2 kg at 120 days of age, surpassing the other groups by 9.7 kg and 10.6 kg. Crossbred gimmers in Group II reached a live weight of 42.0 kg by 4 months, outpacing the other groups by 12.2 kg. The crossbred lambs exhibited an expansive, deep, and sturdy physique, indicative of elevated meat productivity. Physique index analysis displayed that crossbred rams exhibited elongated limbs, bulkiness, and massiveness compared to purebred Kazakh fat-tailed coarse-haired lambs. In the 4.0-4.5-month age range, crossbred rams demonstrated a higher carcass muscle yield than their purebred counterparts, albeit the latter exhibited a 0.18% greater bone yield. Moreover, the meat of groups I and II sheep contained 19.6% and 20.1% protein content, respectively, surpassing the local Kazakh fat-tailed sheep population by 0.7% and 1.2% in absolute terms.

Inheritance of breeding traits in Karakul sheep under different selection approaches.

In the previous century, each intrabreed type of the Karakul sheep breed was characterized by significant numbers, representing a super population with rich genetic diversity. However, over time, the genetic diversity within the breed's gene pool has undergone significant depletion. At present, the Karakul breed is predominantly composed of only two small populations, distinguished by their fur colors: black and gray. Consequently, under such circumstances, genetic advancements in breeding endeavors are likely to be relatively limited, especially given the potential risk of these populations disappearing altogether in the future. Hence, the preservation and judicious utilization of the available genetic resources within the black and gray Karakul sheep populations hold paramount importance in breeding efforts. The primary objective of our research was to investigate the heritability of breeding traits among gray lambs through various selection options. The study was conducted at the "Kumkent" base farm in the Sozak district of the Kyzylorda region. Our findings revealed that the inheritance of gray and black fur colors across the different selection options occurred in a consistent ratio. In the first selection variant (a gray ram with even silver marking ♂ x a black jacket fur type ewe with intense pigmentation ♀), the proportion of gray offspring was 50.6%. Similarly, in the second selection variant (a black jacket fur type ram with intense pigmentation ♂ x a gray ewe with even silver marking ♀), the proportion of gray offspring was 49.6%. The percentage of black lambs obtained in both selection options was nearly equivalent, with 49.4% and 50.4% in the first and second variants, respectively.

Wool production and quality traits of pure- and crossbred Merino-type sheep.

Very little information is available on the quality of wool produced by terminal crosses out of wool producing dam lines. This study was therefore undertaken to elucidate the wool quality traits of four crossbred lines relative to Merinos and Dohne Merinos, which served as dam lines. Each dam line was mated to rams of their own breed as well as to Dormer or Ile de France rams to create four crossbred lines, namely, Dohne x Dormer, Dohne x Ile de France, Merino x Dormer, and Merino x Ile de France, in addition to the purebred Dohne Merino and Merino lines. Four rams and four ewes from each of these lines were reared up to one year of age under optimal growth conditions and shorn at the end of that time. Fleeces were weighed and samples collected for analysis. Neither sex nor genotype influenced clean fleece weight or clean yield percentage, but fibre diameter differed between genotypes. Purebred Merino had the finest wool (18.26 µm) and the Merino x Dormer cross the coarsest (26.01 µm). However, all lines still produced wool that could be used in manufacturing clothing, while fleeces showed good uniformity. The two purebred lines had the highest comfort factor (> 98%). The crossbred lines produced a similar quantity of wool as the purebreds, although of poorer quality. All genotypes except Merino x Dormer appear to produce wool that could be processed into garments, although the crossbred wool would only be suitable for outerwear.

Ruminal microbiota-host crosstalks promote ruminal epithelial development in neonatal lambs with alfalfa hay introduction.

Ruminal microbiota is gradually established after birth, while microbiota maturation could be highly diverse because of varied solid dietary accessibility. However, how the ruminal microbiota accreted from postnatal hay diets alters rumen epithelial development, and how this affects animal health remains largely unknown. Here, neonatal lambs were introduced to starchy corn-soybean starter or corn-soybean starter + alfalfa hay (AH) to investigate the influences of early life ruminal microbiome on rumen epithelial development using integrated 16s rRNA sequencing-metagenome-transcriptome approaches. The results showed that AH introduction elevated average daily weight gain, rumen weight and volume, rumen epithelial papillae length, and rumen muscle layer thickness. Meanwhile, the relative abundance of fibrolytic bacteria (Christensenellaceae R-7 group, Prevotellaceae UCG-001, and Succinivibrio), acetate producer (Acetitomaculum and Mitsuokella), and propionate producer Succiniclasticum was increased in the rumen content by AH supplementation (P < 0.05). Moreover, AH introduction decreased the relative abundance of total CAZymes, CBM, and GH and increased the abundance of KO genes related to volatile fatty acid (VFA) generation in the rumen content. AH lambs had a higher relative abundance of Succiniclasticum, Megasphaera, Succinivibrio, and Suttonella (P < 0.05), while a lower relative abundance of Cloacibacillus, Desulfovibrio, Dialister, Intestinimonas, Parabacteroides, and Pseudoscardovia (P < 0.05) in the rumen epithelial samples. Furthermore, these alterations in ruminal microbial structure and function resulted in ruminal epithelial cell proliferation and development pathways activation. In summary, AH introduction benefited ruminal fiber degradation and VFA generation bacteria colonization and promoted ruminal epithelial development. These findings provide new insights into ruminal microbial-host interactions in the early life.IMPORTANCEWhile it is established that a fiber-rich diet promotes rumen development in lambs, further research is needed to investigate the precise response of rumen microbiota and epithelium to high-quality alfalfa hay. Here, we observed that the inclusion of alfalfa hay led to a discernible alteration in the developmental trajectory of the rumen. Notably, there was a favorable shift in the rumen's volume, morphology, and the development of rumen papillae. Furthermore, ruminal microbial structure and function resulted in ruminal epithelial cell proliferation and development pathways activation, collectively provide compelling evidence supporting the capacity of alfalfa hay to enhance rumen development and health through ruminal micrbiota-host crosstalks. Our findings elucidate the functional response of the rumen to alfalfa hay introduction, providing new insights into strategies for promoting healthy development of the rumen in young ruminants.

Distinguishing host responses, extensive viral dissemination and long-term viral RNA persistence in domestic sheep experimentally infected with Crimean-Congo haemorrhagic fever virus Kosovo Hoti.

Crimean-Congo haemorrhagic fever orthonairovirus (CCHFV) is a tick-borne, risk group 4 pathogen that often causes a severe haemorrhagic disease in humans (CCHF) with high case fatality rates. The virus is believed to be maintained in a tick-vertebrate-tick ecological cycle involving numerous wild and domestic animal species; however the biology of CCHFV infection in these animals remains poorly understood. Here, we experimentally infect domestic sheep with CCHFV Kosovo Hoti, a clinical isolate representing high pathogenicity to humans and increasingly utilized in current research. In the absence of prominent clinical signs, the infection leads to an acute viremia and coinciding viral shedding, fever and markers for potential impairment in liver and kidney functions. A number of host responses distinguish the subclinical infection in sheep versus fatal infection in humans. These include an early reduction of neutrophil recruitment and its chemoattractant, IL-8, in the blood stream of infected sheep, whereas neutrophil infiltration and elevated IL-8 are features of fatal CCHFV infections reported in immunodeficient mice and humans. Several inflammatory cytokines that correlate with poor disease outcomes in humans and have potential to cause vascular dysfunction, a primary hallmark of severe CCHF, are down-regulated or restricted from increasing in sheep. Of particular interest, the detection of CCHFV RNA (including full-length genome) in a variety of sheep tissues long after the acute phase of infection indicates a widespread viral dissemination in the host and suggests a potentially long-term persisting impact of CCHFV infection. These findings reveal previously unrecognized aspects of CCHFV biology in animals.

Phenotype transition from wild mouflon to domestic sheep.

The domestication of animals started around 12,000 years ago in the Near East region. This "endless process" is characterized by the gradual accumulation of changes that progressively marked the genetic, phenotypic and physiological differences between wild and domesticated species. The main distinctive phenotypic characteristics are not all directly attributable to the human-mediated selection of more productive traits. In the last decades, two main hypotheses have been proposed to clarify the emergence of such a set of phenotypic traits across a variety of domestic species. The first hypothesis relates the phenotype of the domesticated species to an altered thyroid hormone-based signaling, whereas the second one relates it to changes in the neural crest cells induced by selection of animals for tameness. These two hypotheses are not necessarily mutually exclusive since they may have contributed differently to the process over time and space. The adaptation model induced by domestication can be adopted to clarify some aspects (that are still controversial and debated) of the long-term evolutionary process leading from the wild Neolithic mouflon to the current domestic sheep. Indeed, sheep are among the earliest animals to have been domesticated by humans, around 12,000 years ago, and since then, they have represented a crucial resource in human history. The aim of this review is to shed light on the molecular mechanisms and the specific genomic variants that underlie the phenotypic variability between sheep and mouflon. In this regard, we carried out a critical review of the most recent studies on the molecular mechanisms that are most accredited to be responsible for coat color and phenotype, tail size and presence of horns. We also highlight that, in such a complicate context, sheep/mouflon hybrids represent a powerful and innovative model for studying the mechanism by which the phenotypic traits related to the phenotypic responses to domestication are inherited. Knowledge of these mechanisms could have a significant impact on the selection of more productive breeds. In fact, as in a journey back in time of animal domestication, the genetic traits of today's domestic species are being progressively and deliberately shaped according to human needs, in a direction opposite to that followed during domestication.

Genetic profile of adaptive immune traits and relationships with parasite resistance and productivity in Scottish Blackface sheep.

Gastrointestinal (GI) parasites cause significant production losses in grazing ruminants which can be mitigated by breeding animals resistant to disease. Lymphocyte cytokine production and parasite-specific Immunoglobulin A (IgA) are adaptive immune traits associated with immunity to GI parasites. To explore the utility of these traits for selective breeding purposes, this study estimated the genetic parameters of the immune traits in sheep and assessed their relationship with disease and productivity traits. Whole blood stimulation assays were performed on 1 040 Scottish Blackface lambs at two months of age in 2016-2017. Blood was stimulated with either pokeweed mitogen (PWM), a non-specific activator of lymphocytes, and Teladorsagia circumcincta (T-ci) larval antigen to activate parasite-specific T lymphocytes. The type of adaptive immune response was determined by quantifying production of cytokines interferon-gamma (IFN-γ), interleukin (IL)-4, and IL-10, which relate to T-helper type (Th) 1, Th2 and regulatory T cell responses, respectively. Serum T-ci specific IgA was also quantified. Heritabilities were estimated for each immune trait by univariate analyses. Genetic and phenotypic correlations were estimated between different immune traits, and between immune traits vs. disease and productivity traits that were recorded at three months of age. Disease phenotypes were expressed as faecal egg counts (FEC) of nematode parasites (Strongyles and Nematodirus), faecal oocyst counts (FOC) of coccidian parasites, and faecal soiling score; production was measured as lamb live weight. Significant genetic variation was observed in all immune response traits. Heritabilities of cytokine production varied from low (0.14 ± 0.06) to very high (0.77 ± 0.09) and were always significantly greater than zero (P < 0.05). IgA heritability was found to be moderate (0.41 ± 0.09). Negative associations previously identified between IFN-γ production and FOC, and IL-4 production and strongyle FEC, were not evident in this study, potentially due to the time-lag between immune and parasitology measures. Instead, a positive genetic correlation was found between FOC and PWM-induced IFN-γ production, while a negative genetic correlation was found between FOC and T-ci induced IL-10. Live weight was negatively genetically correlated with IFN-γ responses. Overall, IFN-γ and IL-4 responses were positively correlated, providing little evidence of cross-regulation of Th1 and Th2 immunity within individual sheep. Furthermore, T-ci specific IgA was highly positively correlated with PWM-induced IL-10, indicating a possible role for this cytokine in IgA production. Our results suggest that while genetic selection for adaptive immune response traits is possible and may be beneficial for parasite control, selection of high IFN-γ responsiveness may negatively affect productivity.

Exploring genetic variation in potential indicators of resilience in sheep using fibre diameter measured along the wool staple.

Production animals are increasingly exposed to a wide variety of disturbances that can compromise their productivity, health and well-being. As a result, there is a growing need to be able to select animals that are more resilient to environmental disturbances. Fibre diameter variation measured along a wool staple is expected to contain information about how resilient sheep are to the disturbances of their internal and external environment. This study aimed to develop potential resilience indicators from fibre diameter variation, estimate their genetic parameters and assess whether these traits are genetically correlated across three age stages. The study used 6 140 Merino sheep from the Sheep Cooperative Research Centre Information Nucleus Flocks recorded at yearling, 2 years old, and adult ages. Eight potential traits were defined based on theory, literature and exploratory analysis, which were suggested to capture the animal's ability to resist, respond and recover from potential disturbances. Genetic evaluation of the traits was conducted using pedigree-based animal models. The traits were shown to be low to moderately heritable (0.01-0.33) when examined at each of the three age stages. The potential indicators were generally well correlated with one another within age stages. Further, the genetic correlation between the same trait measured at different age stages was moderate to high between yearling and 2 years old (0.35-0.94) and between 2 years old and adults (0.18-0.70), while slightly lower between yearling and adult estimates (0.09-0.62). These results suggest that selection for resilience indicators from fibre diameter is possible; however, further studies are warranted to refine the trait definitions and validate these indicators against other measures of health, fitness and productive performance.

Genome-wide association studies of parasite resistance, productivity and immunology traits in Scottish Blackface sheep.

Gastrointestinal parasitism represents a global problem for grazing ruminants, which can be addressed sustainably by breeding animals to be more resistant against infection by parasites. The aim of this study was to assess the genetic architecture underlying traits associated with gastrointestinal parasite resistance, immunological profile and production in meat sheep, and identify and characterise candidate genes affecting these traits. Data on gastrointestinal parasite infection (faecal egg counts for Strongyles (FECS) and Nematodirus (FECN) and faecal oocyst counts for Coccidia, along with faecal soiling scores (DAG), characterised by the accumulation of faeces around the perineum) and production (live weight (LWT)) were gathered from a flock Scottish Blackface lambs at three and four months of age. Data on the immune profile were also collected from a subset of these lambs at two and five months of age. Immune traits included the production of Interferon-γ (IFN-γ), Interleukin (IL)-4 and IL-10 following stimulation of whole blood with pokeweed mitogen (PWM) or antigen from the gastric parasite Teladorsagia circumcincta (T-ci), and serum levels of T. circumcincta-specific immunoglobulin A (IgA). Animals were genotyped with genome-wide DNA arrays, and a total of 1 766 animals and 45 827 Single Nucleotide Polymorphisms (SNPs) were retained following quality control and imputation. Genome-wide association studies were performed for 24 traits. The effects of individual markers with significant effects were estimated, and the genotypic effect solutions were used to estimate additive and dominance effects, and the proportion of additive genetic variance attributed to each SNP locus. A total of 15 SNPs were associated at least at a suggestive level with FECS, FECN, DAG, IgA, PWM-induced IFN-γ and IL-4, and T-ci-induced IL-10. This study uncovered 52 genes closely related to immune function in proximity to these SNPs. A number of genes encoding C-type lectins and killer cell lectin-like family members were close to a SNP associated with FECN, while several genes encoding IL-1 cytokine family members were found to be associated with IgA. Potential candidate genes belonging to or in close proximity with the Major Histocompatibility Complex (MHC) were revealed, including Homeostatic Iron Regulator and butyrophilin coding genes associated with IFN-γ(PWM), and IL-17 coding genes associated with IgA. Due to the importance of the MHC in the control of immune responses, these genes may play an important role in resistance to parasitic infections. Our results reveal a largely complex and polygenic genetic profile of the studied traits in this Scottish Blackface sheep population.

Amplification of Ovis aries papillomavirus type 2 DNA from an ovine cutaneous fibropapilloma.

Seven of 60 Perendale sheep within a flock developed single or multiple exophytic masses on their distal hind limbs. A mass was excised from one sheep and histological evaluation revealed epidermal and mesenchymal proliferation, papillomavirus-induced keratinocyte changes and marked keratohyalin clumping. Ovis aries papillomavirus type 2 DNA sequences were amplified using PCR.

Sept des 60 moutons Perendale d'un troupeau ont développé des masses exophytiques uniques ou multiples sur leurs membres postérieurs distaux. Une masse a été excisée sur un mouton et l'évaluation histologique a révélé une prolifération épidermique et mésenchymateuse, des modifications kératinocytaires induites par le papillomavirus et une agglutination marquée de kératohyaline. Les séquences d'ADN du papillomavirus Ovis aries de type 2 ont été amplifiées par PCR.

Sete de 60 ovelhas Perendale de um rebanho desenvolveram massas exofíticas na porção distal dos seus membros posteriores. Uma massa foi removida de uma ovelha e a avaliação histopatológica revelou proliferação mesenquimal e epidérmica, alterações queratinocíticas induzidas por papilomavírus e aglomeração queratohialina. Sequências de papilomavírus Ovis aries tipo 2 foram amplificadas utilizando PCR.

Siete de 60 ovejas Perendale dentro de un rebaño desarrollaron masas exofíticas únicas o múltiples en sus extremidades traseras distales. Se extirpó una masa de una oveja y la evaluación histológica reveló proliferación epidérmica y mesenquimal, cambios de queratinocitos inducidos por el virus del papiloma y marcada acumulación de queratohialina. Mediante PCR se amplificaron secuencias de DNA del virus del papiloma Ovis aries tipo 2.