Long divergent haplotypes introgressed from wild sheep are associated with distinct morphological and adaptive characteristics in domestic sheep.

The worldwide sheep population comprises more than 1000 breeds. Together, these exhibit a considerable morphological diversity, which has not been extensively investigated at the molecular level. Here, we analyze whole-genome sequencing individuals of 1,098 domestic sheep from 154 breeds, and 69 wild sheep from seven Ovis species. On average, we detected 6.8%, 1.0% and 0.2% introgressed sequence in domestic sheep originating from Iranian mouflon, urial and argali, respectively, with rare introgressions from other wild species. Interestingly, several introgressed haplotypes contributed to the morphological differentiations across sheep breeds, such as a RXFP2 haplotype from Iranian mouflon conferring the spiral horn trait, a MSRB3 haplotype from argali strongly associated with ear morphology, and a VPS13B haplotype probably originating from urial and mouflon possibly associated with facial traits. Our results reveal that introgression events from wild Ovis species contributed to the high rate of morphological differentiation in sheep breeds, but also to individual variation within breeds. We propose that long divergent haplotypes are a ubiquitous source of phenotypic variation that allows adaptation to a variable environment, and that these remain intact in the receiving population probably due to reduced recombination.

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.

Construction of a high-density genetic map and QTL localization of body weight and wool production related traits in Alpine Merino sheep based on WGR.

BACKGROUND: The Alpine Merino is a new breed of fine-wool sheep adapted to the cold and arid climate of the plateau in the world. It has been popularized in Northwest China due to its superior adaptability as well as excellent production performance. Those traits related to body weight, wool yield, and wool fiber characteristics, which are economically essential traits in Alpine Merino sheep, are controlled by QTL (Quantitative Trait Loci). Therefore, the identification of QTL and genetic markers for these key economic traits is a critical step in establishing a MAS (Marker-Assisted Selection) breeding program. RESULTS: In this study, we constructed the high-density genetic linkage map of Alpine Merino sheep by sequencing 110 F1 generation individuals using WGR (Whole Genome Resequencing) technology. 14,942 SNPs (Single Nucleotide Polymorphism) were identified and genotyped. The map spanned 2,697.86 cM, with an average genetic marker interval of 1.44 cM. A total of 1,871 high-quality SNP markers were distributed across 27 linkage groups, with an average of 69 markers per LG (Linkage Group). Among them, the smallest genetic distance is 19.62 cM for LG2, while the largest is 237.19 cM for LG19. The average genetic distance between markers in LGs ranged from 0.24 cM (LG2) to 3.57 cM (LG17). The marker density in the LGs ranged from LG14 (39 markers) to LG1 (150 markers). CONCLUSIONS: The first genetic map of Alpine Merino sheep we constructed included 14,942 SNPs, while 46 QTLs associated with body weight, wool yield and wool fiber traits were identified, laying the foundation for genetic studies and molecular marker-assisted breeding. Notably, there were QTL intervals for overlapping traits on LG4 and LG8, providing potential opportunities for multi-trait co-breeding and further theoretical support for selection and breeding of ultra-fine and meaty Alpine Merino 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.

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 diversity of United States Rambouillet, Katahdin and Dorper sheep.

BACKGROUND: Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs). RESULTS: Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms. CONCLUSIONS: In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

A single base pair duplication in the SLC33A1 gene is associated with fetal losses and neonatal lethality in Manech Tête Rousse dairy sheep.

We recently discovered that the Manech Tête Rousse (MTR) deficient homozygous haplotype 2 (MTRDHH2) probably carries a recessive lethal mutation in sheep. In this study, we fine-mapped this region through whole-genome sequencing of five MTRDHH2 heterozygous carriers and 95 non-carriers from various ovine breeds. We identified a single base pair duplication within the SLC33A1 gene, leading to a frameshift mutation and a premature stop codon (p.Arg246Alafs*3). SLC33A1 encodes a transmembrane transporter of acetyl-coenzyme A that is crucial for cellular metabolism. To investigate the lethality of this mutation in homozygous MTR sheep, we performed at-risk matings using artificial insemination (AI) between heterozygous SLC33A1 variant carriers (SLC33A1_dupG). Pregnancy was confirmed 15 days post-AI using a blood test measuring interferon Tau-stimulated MX1 gene expression. Ultrasonography between 45 and 60 days post-AI revealed a 12% reduction in AI success compared with safe matings, indicating embryonic/fetal loss. This was supported by the MX1 differential expression test suggesting fetal losses between 15 and 60 days of gestation. We also observed a 34.7% pre-weaning mortality rate in 49 lambs born from at-risk matings. Homozygous SLC33A1_dupG lambs accounted for 47% of this mortality, with deaths occurring mostly within the first 5 days without visible clinical signs. Therefore, appropriate management of SLC33A1_dupG with an allele frequency of 0.04 in the MTR selection scheme would help increase overall fertility and lamb survival.

MFSD2A frameshift variant in Kerry Hill sheep with microcephaly.

Microcephaly is a rare neurodevelopmental disorder characterized by reduced skull circumference and brain volume that occurs sporadically in farm animals. We investigated an early-onset neurodegenerative disorder observed in seven lambs of purebred Kerry Hill sheep. Clinical signs included inability to stand or severe ataxia, convulsions, and early death. Diagnostic imaging and brain necropsy confirmed microcephaly. The pedigree of the lambs suggested monogenic autosomal recessive inheritance. We sequenced the genome of one affected lamb, and comparison with 115 control genomes revealed a single private protein-changing variant. This frameshift variant, MFSD2A: c.285dupA, p.(Asp96fs*9), represents a 1-bp duplication predicted to truncate 80% of the open reading frame. MFSD2A is a transmembrane protein that is essential for maintaining blood-brain barrier homeostasis and plays a key role in regulating brain lipogenesis. Human MFSD2A pathogenic variants are associated with a neurodevelopmental disorder with progressive microcephaly, spasticity, and brain imaging abnormalities (NEDMISBA, OMIM 616486). Here we present evidence for the occurrence of a recessively inherited form of microcephaly in sheep due to a loss-of-function variant in MFSD2A (OMIA 002371-9940). To the best of our knowledge, this is the first report of a spontaneous MFSD2A variant in domestic animals.

ASReml-based estimation of early genetic parameters in hybrid mutton sheep populations.

The aim of this study was to analyze the effects of non-genetic factors on the estimation of genetic parameters of early growth traits in hybrid mutton sheep using ASReml software, in order to provide theoretical basis for screening the optimal hybriding combinations and accelerating the breeding process of new breeds of specialized housed-feeding mutton sheep. We selected the wellgrown hybrid Southhu (Southdown × Hu sheep) and Dorhu (Dorset × Hu sheep) sheep as the research objects, constructed weight correction formulae for SH and DH sheep at 60 and 180 days; and used ASReml software to investigate the effects of non-genetic factors on the estimation of genetic parameters of early growth traits in hybrid sheep. The results showed that the birth month and birth type were found significant for all traits (p < 0.001); the heritability of maternal effects ranged from 0.0709 to 0.1859. It was found that both SH and DH sheep emerged the potential for rapid early growth and development, early growth traits are significantly affected by maternal genetic effects, thereby the maternal effect should be taken into consideration for the purpose of improving accuracy in parameter estimations and therefore increasing the success of 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.

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.

Genome-wide detection of structural variation in some sheep breeds using whole-genome long-read sequencing data.

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep.

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.

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.

Comparative assessment of growth curves in Brazilian white Morada Nova hair sheep breed.

The Morada Nova sheep breed is essential for the economy of the semi-arid region of Northeast Brazil, standing out for its adaptability, resistance to parasites and reproductive ability. However, the white variant is endangered, highlighting the importance of studies on its productivity to support conservation efforts. This study focuses on the growth curve of the Morada Nova sheep breed, using nonlinear models and analyzing flock profiles. Total of 764 observations of 165 animals from four farms in Ceará and Rio Grande do Norte, Brazil, were analyzed. Canonical discriminant analysis (CDA) was used for the exploratory analysis and four nonlinear models were used to study the growth curve. Weight from birth to 270 days of age, absolute growth rate (AGR), and the impact of sex on growth curves were assessed. Sex and farm are significant discriminating variables (P < 0.05) for the studied effects (weight and age). Weight was the primary phenotypic biomarker that discriminated between the two indicators, while age was a discriminating indicator only for the core effect. The Gompertz model was the most efficient, presenting the lowest residuals and greatest convergence. The study reveals new information about the growth of Morada Nova sheep, the white variety, including weight differences between the sexes at all analyzed ages and an inflection point before 90 days of age. These discoveries contribute to the understanding of the breed's growth and help in the formulation of conservation strategies.

Gene expression in the Longissimus dorsi muscle related to meat quality from tropical hair lambs.

The present study describes the expression of genes in the Longissimus dorsi muscle related to meat quality of hair lambs finished in an Integration Crop-Livestock system. Twenty-eight non-castrated lambs of two breeds, Somalis Brasileira and Santa Inês, at 120 ± 15 days of age, with an average initial live weight of 18 ± 3.1 kg, were kept in a pasture-based finishing system with supplementation. Upon reaching 28 kg body weight, animals were sent for slaughter. Samples of the Longissimus dorsi and Biceps femoris muscle were harvested for analyses of gene expression and physicochemical properties. Significant differences were detected between the breeds for tissue and chemical composition, whereas the physical aspects did not differ. We observed the expression of six genes related to lipid synthesis (acetyl-CoA carboxylase [ACACA], fatty acid synthase [FAS], stearoyl-CoA desaturase [SCD], lipoprotein lipase [LPL], cell death-inducing DFFA-like effector A [CIDEA], and thyroid hormone responsive [THRSP]) and six genes related to molecular synthesis (myostatin [MSTN], growth differentiation factor 8 [GDF8], insulin-like growth factor 1 [IGF1], insulin-like growth factor 2 [IGF2], delta-like 1 homolog [DLK1], and growth hormone receptor [GHr]) in both breeds. The Santa Inês breed and the Somalis Brasileira showed similar expression patterns of genes related to lipogenesis and myogenesis of the Longissimus dorsi muscle, with the exception of the THRSP gene, in which the Somalis Brasileira have more receptors for the action of thyroid hormones, which resulted in greater thickness of fat in the carcass (subcutaneous fat) and higher lipid content in the chemical composition of the meat.

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.

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity, Introgression, and Agronomically Important Loci.

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3'-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Evidence of novel Treponema phylotypes implicated in contagious ovine digital dermatitis and association of treponemes with major lameness causing foot pathogens.

In this study 269 swabs collected from 254 ovine foot lesions and 15 apparently healthy ovine feet were screened by PCR for the presence of major lameness causing foot pathogens viz. Treponema species, D. nodosus, F. necrophorum and T. pyogenes with the presumption that ovine foot lesion positive for Treponema species alone or in association with other three pathogens were categorized as contagious ovine digital dermatitis (CODD). While samples positive for D. nodosus alone or its combination with F. necrophorum and T. pyogenes were considered as footrot (FR) and samples in which F. necrophorum or T. pyogenes was found either alone or in combination were considered as interdigital dermatitis (ID). The overall occurrence of Treponema sp. in ovine foot lesions was 48.0%, and ranged from 33 to 58%. In Treponema positive samples D. nodosus, F. necrophorum and T. pyogenes were present in 34 (27.4%), 66 (54.4%) and 84 (68.5%) in contrast to Treponema negative samples in which these were present in 15 (11.1%), 20 (14.12%) and 17 (12.6%) samples, respectively. The data signifies that Treponema sp. are significantly associated with these foot pathogens and their different combinations with Treponema sp. influence the severity of CODD lesion. The identification of Treponema phylotypes was done by sequencing the 16S rRNA gene fragment of ten representative samples. Out of ten sequences, four (Trep-2, Trep-4, Trep-7 and Trep-10) were identical to Treponema sp. phylotype 1 (PT1) that belongs to phylogroup T. refringens-like, one sequence (Trep-1) was genetically close (90% sequence homology) to Treponema brennaborense while five sequences (Trep-3, Trep-5, Trep-6, Trep-8 and Trep-9) matched with uncultured bacterium clones of treponemes forming separate monophyletic group in phylogenetic tree and could represent new digital dermatitis phylogroup presently containing five ovine specific phylotypes. This is the first report on the presence of Treponema phylotypes other than three digital dermatitis (DD) Treponema phylogroups viz. T. phagedenis-like, T. medium/T. vincentii-like, and T. pedis-like that are frequently detected in CODD lesions. Metagenomic analysis of two representative samples revealed the abundance of genus Treponema in CODD lesion while this genus was absent in swab collected from clinically healthy foot suggesting that it might play primary role in producing CODD. These findings may further aid in understanding the etiopathogenesis of CODD and could help to develop appropriate treatment and mitigation strategies to combat the disease.

A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds.

BACKGROUND: To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. RESULTS: The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. CONCLUSIONS: To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.