Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

Analysing the genetic diversity of three sheep breeds in Turkey and nearby countries using 50 K SNPs data

This study analysed the genetic diversity and population structure of eight sheep breeds in Turkey and nearby countries. Moderate genetic diversity was observed, with the Sakiz (SKZ) exhibiting the highest diversity based on heterozygosity and allelic richness (AR) values. Genetic distances revealed differentiation between the populations, with the most significant divergence between the Cyprus Fat Tail (CFT) and SKZ breeds. PCA demonstrated SKZ and Chios (CHI) clustering together, indicating genetic similarity. Karakas (KRS), Norduz (NDZ), Afshari (AFS), Moghani (MOG) and others showed overlap, reflecting genetic relationships. Ancestry analysis found that KRS was predominantly inherited from the second ancestral population, while SKZ and NDZ were primarily derived from the first and second ancestral lineages. This illustrated the populations' diverse origins. Most genetic variation (96.84%) was within, not between, populations. The phi-statistic (PhiPT) indicated moderate differentiation overall. Phylogenetic analysis further demonstrated the genetic distinctiveness of the SKZ breed. ROH and FROH analyses showed that SKZ exhibited the highest homozygosity and inbreeding, while KRS displayed the lowest. This study elucidates these breeds' genetic diversity, structure and relationships. Key findings include moderate diversity, evidence of differentiation between breeds, diverse ancestral origins and distinct ROH patterns. This provides insights into the population's genetic characteristics and conservation requirements.

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.

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.

Melanophilin Polymorphism in Ferrets of Different Color

In mammals, the main contribution to the variability of pigmentation is made by two groups of genes directly related to the metabolic pathways of pigment synthesis and controlling the transport of melanosomes in melanocytes to keratinocytes. In order to identify the genetic basis of pigmentation variants, the nucleotide sequences of the melanophilin gene were compared in two groups of ferrets-silver-colored and wild-type animals-using sequencing of 16 exons. In carriers of silver color, a single nucleotide deletion was detected in the 9th exon, leading to a shift in the reading frame and the formation of a stop codon downstream. The protein encoded by the mutant allele is almost completely devoid of the C terminal domain of the protein responsible for the contact of melanosomes with actin during their moving to the periphery of melanocytes, but it retains the leading domain involved in the formation of melanosomes. The combination of the preservation of the N domain and the defect of the C domain of the mutant protein for the first time makes it possible to explain the incomplete dominance of the wild-type protein in heterozygotes.

Assessment of morphological features of skin and hair cover of Karakul lambs (South of Kazakhstan)

The relevance of the present research is related to the great importance of Karakul sheep breeding in Kazakhstan and the insufficient investigation of morphological features of the skin and hair cover of Karakul lambs in modern literature. The purpose of the research was to explore the specific features of hair pigmentation and histological structure of the skin of Karakul lambs of different colour groups: black, sur, grey, and brown. The work involved the standard methods of assessing the pigment content in the hair of lambs and methods of histological studies. The study utilised statistical methods including the Kolmogorov-Smirnov criterion, one-factor analysis of variance, Fisher's F-criterion, and the Bonferroni method, with data processed through the "PASW Statistics 17" and "Statistica 6.0" software packages. The results of the research demonstrate that the pigment content in the hair of new-born Karakul lambs depends on their colouring. Differences in pigment precipitation between different colour groups of lambs were established: pigment grains in black hairs and suras were larger, in comparison with brown and grey hairs, in which pigment was represented by the smallest grains. At the age of 19 days, no significant changes in the colouring of the lambs examined were found (p > 0.05), whereas on reaching 4 months of age, depigmentation (p < 0.05) was noted in all colour groups, which was slower in black lambs.

Single-cell sequencing reveals the reproductive variations between primiparous and multiparous Hu ewes

BACKGROUND

In the modern sheep production systems, the reproductive performance of ewes determines the economic profitability of farming. Revealing the genetic mechanisms underlying differences in the litter size is important for the selection and breeding of highly prolific ewes. Hu sheep, a high-quality Chinese sheep breed, is known for its high fecundity and is often used as a model to study prolificacy traits. In the current study, animals were divided into two groups according to their delivery rates in three consecutive lambing seasons (namely, the high and low reproductive groups with ≥ 3 lambs and one lamb per season, n = 3, respectively). The ewes were slaughtered within 12 h of estrus, and unilateral ovarian tissues were collected and analyzed by 10× Genomics single-cell RNA sequencing.

RESULTS

A total of 5 types of somatic cells were identified and corresponding expression profiles were mapped in the ovaries of each group. Noticeably, the differences in the ovary somatic cell expression profiles between the high and low reproductive groups were mainly clustered in the granulosa cells. Furthermore, four granulosa cell subtypes were identified. GeneSwitches analysis revealed that the abundance of JPH1 expression and the reduction of LOC101112291 expression could lead to different evolutionary directions of the granulosa cells. Additionally, the expression levels of FTH1 and FTL in mural granulosa cells of the highly reproductive group were significantly higher. These genes inhibit necroptosis and ferroptosis of mural granulosa cells, which helps prevent follicular atresia.

CONCLUSIONS

This study provides insights into the molecular mechanisms underlying the high fecundity of Hu sheep. The differences in gene expression profiles, particularly in the granulosa cells, suggest that these cells play a critical role in female prolificacy. The findings also highlight the importance of genes such as JPH1, LOC101112291, FTH1, and FTL in regulating granulosa cell function and follicular development.

Impacts of Population Size and Domestication Process on Genetic Diversity and Genetic Load in Genus Ovis

In theoretical biology, a prevailing hypothesis posits a profound interconnection between effective population size (Ne), genetic diversity, inbreeding, and genetic load. The domestication and improvement processes are believed to be pivotal in diminishing genetic diversity while elevating levels of inbreeding and increasing genetic load. In this study, we performed a whole genome analysis to quantity genetic diversity, inbreeding, and genetic load across seven wild Ovis species and five domesticated sheep breeds. Our research demonstrates that the genetic load and diversity of species in the genus Ovis have no discernible impact on recent Ne, and three species within the subgenus Pachyceros tend to carry a higher genetic load and lower genetic diversity patterns. The results coincide with these species' dramatic decline in population sizes within the subgenus Pachyceros ~80-250 thousand years ago. European mouflon presented with the lowest Ne, lower genetic diversity, and higher individual inbreeding coefficient but a lower genetic load (missense and LoF). This suggests that the small Ne of European mouflon could reduce harmful mutations compared to other species within the genus Ovis. We showed lower genetic diversity in domesticated sheep than in Asiatic mouflon, but counterintuitive patterns of genetic load, i.e., lower weak genetic load (missense mutation) and no significant difference in strong genetic load (LoF mutation) between domestic sheep and Asiatic mouflon. These findings reveal that the "cost of domestication" during domestication and improvement processes reduced genetic diversity and purified weak genetic load more efficiently than wild species.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Whole-Genome Resequencing Reveals Selection Signal Related to Sheep Wool Fineness

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

Heritability and recursive influence of host genetics on the rumen microbiota drive body weight variance in male Hu sheep lambs

BACKGROUND

Heritable rumen microbiota is an important modulator of ruminant growth performance. However, no information exists to date on host genetics-rumen microbiota interactions and their association with phenotype in sheep. To solve this, we curated and analyzed whole-genome resequencing genotypes, 16S rumen-microbiota data, and longitudinal body weight (BW) phenotypes from 1150 sheep.

RESULTS

A variance component model indicated significant heritability of rumen microbial community diversity. Genome-wide association studies (GWAS) using microbial features as traits identified 411 loci-taxon significant associations (P < 10-8). We found a heritability of 39% for 180-day-old BW, while also the rumen microbiota likely played a significant role, explaining that 20% of the phenotypic variation. Microbiota-wide association studies (MWAS) and GWAS identified four marker genera (Bonferroni corrected P < 0.05) and five novel genetic variants (P < 10-8) that were significantly associated with BW. Integrative analysis identified the mediating role of marker genera in genotype influencing phenotype and unravelled that the same genetic markers have direct and indirect effects on sheep weight.

CONCLUSIONS

This study reveals a reciprocal interplay among host genetic variations, the rumen microbiota and the body weight traits of sheep. The information obtained provide insights into the diverse microbiota characteristics of rumen and may help in designing precision microbiota management strategies for controlling and manipulating sheep rumen microbiota to increase productivity. Video Abstract.

Genomic Landscape of Copy Number Variations and Their Associations with Climatic Variables in the World's Sheep

Sheep show characteristics of phenotypic diversity and adaptation to diverse climatic regions. Previous studies indicated associations between copy number variations (CNVs) and climate-driven adaptive evolution in humans and other domestic animals. Here, we constructed a genomic landscape of CNVs (n = 39,145) in 47 old autochthonous populations genotyped at a set of high-density (600 K) SNPs to detect environment-driven signatures of CNVs using a multivariate regression model. We found 136 deletions and 52 duplications that were significantly (Padj. < 0.05) associated with climatic variables. These climate-mediated selective CNVs are involved in functional candidate genes for heat stress and cold climate adaptation (e.g., B3GNTL1, UBE2L3, and TRAF2), coat and wool-related traits (e.g., TMEM9, STRA6, RASGRP2, and PLA2G3), repairing damaged DNA (e.g., HTT), GTPase activity (e.g., COPG), fast metabolism (e.g., LMF2 and LPIN3), fertility and reproduction (e.g., SLC19A1 and CCDC155), growth-related traits (e.g., ADRM1 and IGFALS), and immune response (e.g., BEGAIN and RNF121) in sheep. In particular, we identified significant (Padj. < 0.05) associations between probes in deleted/duplicated CNVs and solar radiation. Enrichment analysis of the gene sets among all the CNVs revealed significant (Padj. < 0.05) enriched gene ontology terms and pathways related to functions such as nucleotide, protein complex, and GTPase activity. Additionally, we observed overlapping between the CNVs and 140 known sheep QTLs. Our findings imply that CNVs can serve as genomic markers for the selection of sheep adapted to specific climatic conditions.

The identification of a novel SNP in the resistin (RETN) gene associated with growth traits in Karakul and Awassi sheep

Resistin is one of the most important adipocytokines in mammalian cells due to its involvement in insulin resistance, obesity, and autoimmune diseases. Resistin is encoded by RETN gene that is primarily expressed in adipose tissues. Mutations in this gene have been associated with several productive traits in animals. This study was conducted to assess the possible biomarker capacity of RETN by evaluating its association with growth traits in two economically important sheep in the Middle East. Genotyping was conducted using PCR-single strand conformation polymorphism (SSCP), and the polymorphism of RETN was associated with several growth traits for three months intervals starting from birth until one year of age. In a total of 190 Karakul sheep and 245 Awassi sheep, only one SNP (233A > C) was detected in the RETN gene. The identified novel SNP showed significant associations with all growth traits at the ages of six, nine, and twelve months. At the age of six months onward, lambs with AC and CC genotypes showed respectively lower body weight and length, chest and abdominal circumferences, and wither and rump heights than those with AA genotype. Due to the remarkable association between RETN;233A > C and lower growth traits, this genotype is suggested as a promising marker to assess growth traits in Karakul and Awassi sheep. This is the first study that demonstrated the importance of RETN as a possible tool for evaluating growth traits in two breeds of sheep with a possibility to be applied to other breeds via large-scale association analysis.

Whole-body adipose tissue multi-omic analyses in sheep reveal molecular mechanisms underlying local adaptation to extreme environments

The fat tail of sheep is an important organ that has evolved to adapt to extreme environments. However, the genetic mechanisms underlying the fat tail phenotype remain poorly understood. Here, we characterize transcriptome and lipidome profiles and morphological changes in 250 adipose tissues from two thin-tailed and three fat-tailed sheep populations in summer and winter. We implement whole-genome selective sweep tests to identify genetic variants related to fat-tails. We identify a set of functional genes that show differential expression in the tail fat of fat-tailed and thin-tailed sheep in summer and winter. These genes are significantly enriched in pathways, such as lipid metabolism, extracellular matrix (ECM) remodeling, molecular transport, and inflammatory response. In contrast to thin-tailed sheep, tail fat from fat-tailed sheep show slighter changes in adipocyte size, ECM remodeling, and lipid metabolism, and had less inflammation in response to seasonal changes, indicating improved homeostasis. Whole-genome selective sweep tests identify genes involved in preadipocyte commitment (e.g., BMP2, PDGFD) and terminal adipogenic differentiation (e.g., VEGFA), which could contribute to enhanced adipocyte hyperplasia. Altogether, we establish a model of regulatory networks regulating adipose homeostasis in sheep tails. These findings improve our understanding of how adipose homeostasis is maintained, in response to extreme environments in animals.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise F_ST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and F_ST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.