Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits

Identification of AGO2 and PLEC genes polymorphisms in Hu sheep and their relationship with body size traits

The body size traits are major traits in livestock, which intuitively displays the development of the animal's bones and muscles. This study used PCR amplification, Sanger sequencing, KASPar genotyping, and quantitative real-time reverse transcription PCR (qRT-PCR) to analyze the Single-nucleotide polymorphism and expression characteristics of Argonaute RISC catalytic component 2 (AGO2) and Plectin (PLEC) genes in Hu sheep. Two intron mutations were found in Hu sheep, which were AGO2 g.51700 A > C and PLEC g.23157 C > T, respectively. Through association analysis of two mutation sites and body size traits, it was found that AGO2 g.51700 A > C mainly affects the chest and cannon circumference of Hu sheep of while PLEC g.23157 C mainly affects body height and body length. The combined genotypes of AGO2 and PLEC genes with body size traits showed SNPs at the AGO2 g.51700 A > C and PLEC g.23157 C > T loci significantly improved the body size traits of Hu sheep. In addition, the AGO2 gene has the highest expression levels in the heart, rumen, and tail fat, and the PLEC gene is highly expressed in the heart. These two loci can provide new research ideas for improving the body size traits of Hu sheep.

Genomic diversity of the locally developed Latvian Darkheaded sheep breed

The Latvian Darkheaded is the only locally developed sheep breed. The breed was formed at the beginning of the 20th century by crossing local coarse-wooled sheep with the British Shropshire and Oxfordshire breeds. The breed was later improved by adding Ile-de-France, Texel, German blackheads, and Finnsheep to achieve higher prolificacy and better meat quality. Previous studies have reported the Latvian Darkheaded sheep to be closely related to Estonian and Lithuanian Blackface breeds, according to microsatellite data. To expand our knowledge of the genetic resources of the Latvian Darkheaded breed, we conducted a whole-genome resequencing analysis on 40 native sheep. The investigation showed that local sheep harbor genetic diversity levels similar to those observed among other improved breeds of European origin, including Charollais and Suffolk. Genome-wide nucleotide diversity (π) in Latvian Darkheaded sheep was 3.91 × 10-3, whereas the average observed heterozygosity among the 40 animals was 0.267 and 0.438 within the subsample of unrelated individuals. The Ne has rapidly decreased to 200 ten generations ago with a recent drop to Ne 73 four generations ago. However, inbreeding levels based on runs of homozygosity were, on average, low, with FROH ranging between 0.016 and 0.059. The analysis of the genomic composition of the breed confirmed shared ancestry with sheep of British origin, reflecting the history of the breed. Nevertheless, Latvian Darkheaded sheep were genetically separable. The contemporary Latvian Darkheaded sheep population is genetically diverse with a low inbreeding rate. However, further development of breed management programs is necessary to prevent an increase in inbreeding, loss of genetic diversity, and depletion of breed-specific genetic resources, ensuring the preservation of the native Latvian Darkheaded sheep.

Relationship between Rumen Microbial Differences and Phenotype Traits among Hu Sheep and Crossbred Offspring Sheep

This experiment was conducted to investigate the effect of three-way hybrid sheep and Hu sheep on serum indicators, rumen fermentation, rumen enzyme activity, and microorganisms in sheep. Healthy and similar birth weights from three groups (Hu, n = 11; Charolais × Australian White × Hu, CAH, n = 11; Charolais × Dorper × Hu, CDH, n = 11) were selected to be fed by the ewes until 45 days of age. Subsequently, they were weaned intensively and underwent short-term fattening for 3 months along with selected male lambs fed intensively. During this period, they were fed and watered ad libitum. Blood and rumen fluid were collected and analyzed for serum indicators and rumen fluid microorganisms, enzyme activity, and VFA, respectively, at the end of the fattening period. Compared with Hu lamb, the offspring of the three-way hybrid lamb showed significant improvements in body weight, serum lactate dehydrogenase, and creatinine content. However, there was no significant effect on serum immunity and antioxidant indices. In addition, the rumen fluid volatile fatty acid (VFA) molar concentration and microcrystalline cellulose and lipase content were significantly lower in the three-way hybrid lamb compared to Hu lamb, but β-glucosidase, amylase, pepsin, and VFA molar ratio were not significantly affected. Subsequently, 16S rRNA sequencing diversity analysis revealed that three-way hybrid lamb significantly increased rumen microbial ACE and Chao1 indices compared to Hu lamb. Meanwhile, the abundance of Verrucomicrobiota and Synergistota significantly increased at the phylum level. Correlation analysis showed that Prevotella had the highest proportion, while Rikenellaceae_RC9_gut_group correlated most closely with others genus. The microbial communities isovaleric acid molar concentration and proportion were strongly correlated. In addition, there were significant differences in correlations between microbial communities and isobutyric acid, butyric acid and valeric acid content, and their molar proportion, but they were not significantly correlated with digestive enzymes. From the functional enrichment analysis, it was found that hybrid progeny were mainly enriched in the pyruvate metabolism, microbial metabolism in diverse environments, carbon metabolism, and quorum sensing pathways. In contrast, the Hu sheep were primarily enriched in the cysteine and methionine, amino sugar and nucleotide sugar, and biosynthesis of secondary metabolite pathways. These results suggest that hybridization can play a role in regulating organismal metabolism and improve animal production performance by influencing the structure and characteristics of microbial communities.

A time-resolved multi-omics atlas of transcriptional regulation in response to high-altitude hypoxia across whole-body tissues

High-altitude hypoxia acclimatization requires whole-body physiological regulation in highland immigrants, but the underlying genetic mechanism has not been clarified. Here we use sheep as an animal model for low-to-high altitude translocation. We generate multi-omics data including whole-genome sequences, time-resolved bulk RNA-Seq, ATAC-Seq and single-cell RNA-Seq from multiple tissues as well as phenotypic data from 20 bio-indicators. We characterize transcriptional changes of all genes in each tissue, and examine multi-tissue temporal dynamics and transcriptional interactions among genes. Particularly, we identify critical functional genes regulating the short response to hypoxia in each tissue (e.g., PARG in the cerebellum and HMOX1 in the colon). We further identify TAD-constrained cis-regulatory elements, which suppress the transcriptional activity of most genes under hypoxia. Phenotypic and transcriptional evidence indicate that antenatal hypoxia could improve hypoxia tolerance in offspring. Furthermore, we provide time-series expression data of candidate genes associated with human mountain sickness (e.g., BMPR2) and high-altitude adaptation (e.g., HIF1A). Our study provides valuable resources and insights for future hypoxia-related studies in mammals.

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.

Effects of nonsynonymous single nucleotide polymorphisms of the KIAA1217, SNTA1 and LTBP1 genes on the growth traits of Ujumqin sheep

Sheep body size can directly reflect the growth rates and fattening rates of sheep and is also an important index for measuring the growth performance of meat sheep. In this study, high-resolution resequencing data from four sheep breeds (Dorper sheep, Suffolk sheep, Ouessant sheep, and Shetland sheep) were analyzed. The nonsynonymous single nucleotide polymorphisms of three candidate genes (KIAA1217, SNTA1, and LTBP1) were also genotyped in 642 healthy Ujumqin sheep using MALDI-TOFMS and the genotyping results were associated with growth traits. The results showed that different genotypes of the KIAA1217 g.24429511T>C locus had significant effects on the chest circumferences of Ujumqin sheep. The SNTA1 g.62222626C>A locus had different effects on the chest depths, shoulder widths and rump widths of Ujumqin sheep. This study showed that these two sites can be used for marker-assisted selection, which will be beneficial for future precision molecular breeding.

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.

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 Aşıklı 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 Aşıklı 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 Aşıklı 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.

Are two beneficial mutations (p.Q249R and 90-bp Indel) within the ovine BMPRIB gene associated with growth traits?

Background

The problem of achieving economic efficiency in sheep breeding can be largely solved by increasing sheep productivity. Recently, the BMPRIB gene has been revealed by GWAS as a potential candidate gene for sheep body morphometric traits. Therefore, the present study aimed to investigate whether genetic polymorphisms (p.Q249R SNP and 90-bp deletion) in the BMPRIB gene are associated with sheep growth traits.

Methods

PCR-based genotyping was performed on 1,875 sheep, including 1,191 Guiqian semi-fine wool (GQSFW), 560 Luxi Blackhead (LXBH), 55 Lanzhou fat-tailed (LZFT), and 69 Weining (WN) sheep. Genotype-phenotype association was assessed using the independent samples t-test and ANOVA. The significance level was set at αoriginal < 0.05. The threshold p-value for significance was adjusted after correction for multiple comparisons using the Bonferroni correction.

Results

After the Bonferroni correction, it was found that individuals with FecB+/FecB+ genotypes of the p.Q249R had significantly better growth traits in LXBH ewe lambs, including the body length, chest width, paunch girth, cannon circumference, and hip width (P<0.0005). Meanwhile, associations were observed between 90-bp deletion polymorphism and several growth traits (body length, body height, chest depth, and canon circumference) in GQSFW ewe adults after the Bonferroni correction (P < 0.0002), and individuals with the "DD" genotypes had greater growth traits.

Conclusion

Our findings align with the experimental observations from GWAS, which identified the BMPRIB gene as a potential candidate gene for body measurement traits. These findings not only confirm the previous study's results but also expand on them. Therefore, further investigations regarding the impact of BMPRIB polymorphisms on growth traits are necessary in other sheep breeds.

Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits

The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveal those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.

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

BACKGROUND

Chinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations.

RESULTS

Based on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability.

CONCLUSIONS

Analysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.

Two mutations at KRT74 and EDAR synergistically drive the fine-wool production in Chinese sheep

INTRODUCTION

Fine-wool sheep are the most common breed used by the wool industry worldwide. Fine-wool sheep have over a three-fold higher follicle density and a 50% smaller fiber diameter than coarse-wool sheep.

OBJECTIVES

This study aims to clarify the underlying genetic basis for the denser and finer wool phenotype in fine-wool breeds.

METHOD

Whole-genome sequences of 140 samples, Ovine HD630K SNP array data of 385 samples, including fine, semi-fine, and coarse wool sheep, as well as skin transcriptomes of nine samples were integrated for genomic selection signature analysis.

RESULTS

Two loci at keratin 74 (KRT74) and ectodysplasin receptor (EDAR) were revealed. Fine-scale analysis in 250 fine/semi-fine and 198 coarse wool sheep narrowed this association to one C/A missense variant of KRT74 (OAR3:133,486,008, P = 1.02E-67) and one T/C SNP in the regulatory region upstream of EDAR (OAR3:61,927,840, P = 2.50E-43). Cellular over-expression and ovine skin section staining assays confirmed that C-KRT74 activated the KRT74 protein and specifically enlarged cell size at the Huxley's layer of the inner root sheath (P < 0.01). This structure enhancement shapes the growing hair shaft into the finer wool than the wild type. Luciferase assays validated that the C-to-T mutation upregulated EDAR mRNA expression via a newly created SOX2 binding site and potentially led to the formation of more hair placodes.

CONCLUSIONS

Two functional mutations driving finer and denser wool production were characterized and offered new targets for genetic breeding during wool sheep selection. This study not only provides a theoretical basis for future selection of fine wool sheep breeds but also contributes to improving the value of wool commodities.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Screening of Litter-Size-Associated SNPs in NOX4, PDE11A and GHR Genes of Sheep

In previous studies, NOX4, PDE11A and GHR genes have been screened as important candidate genes for litter size in sheep by using the GWAS method; however, neither their effects on litter size nor the loci associated with litter size have been identified. In this study, three candidate loci (c.1057-4C > T in NOX4, c.1983C > T in PDE11A and c.1618C > T in GHR) were first screened based on our previous resequencing data of 10 sheep breeds. After the three loci were genotyped using Sequenom MassARRAY technology, we carried out population genetics analysis on the three loci and performed association analysis between the polymorphism of the three loci and the litter size of sheep. The results of population genetics analysis suggested that c.1057-4C > T in NOX4 and c.1983C > T in PDE11A may be subject to natural or artificial selection. The results of association analysis indicated that litter size was significantly associated with c.1057-4C > T in NOX4 and c.1983C > T in PDE11A (p < 0.05) in Small Tail Han sheep, and there was no significant interaction effect between the two loci on the litter size. In summary, c.1057-4C > T in NOX4 and c.1983 C > T in PDE11A can be considered candidate molecular markers for improving litter size in sheep.

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

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

Natural selection and genetic diversity maintenance in a parasitic wasp during continuous biological control application

Aphidius gifuensis is a parasitoid wasp and primary endoparasitoid enemy of the peach potato aphid, Myzus persicae. Artificially reared, captive wasps of this species have been extensively and effectively used to control populations of aphids and limit crop loss. However, the consequences of large-scale releasing of captive A. gifuensis, such as genetic erosion and reduced fitness in wild populations of this species, remains unclear. Here, we sequence the genomes of 542 A. gifuensis individuals collected across China, including 265 wild and 277 human-intervened samples. Population genetic analyses on wild individuals recovered Yunnan populations as the ancestral group with the most complex genetic structure. We also find genetic signature of environmental adaptation during the dispersal of wild populations from Yunnan to other regions. While comparative genomic analyses of captive wasps revealed a decrease in genetic diversity during long-term rearing, population genomic analyses revealed signatures of natural selection by several biotic (host plants) or abiotic (climate) factors, which support maintenance of the gene pool of wild populations in spite of the introduction of captive wasps. Therefore, the impact of large-scale release is reduced. Our study suggests that A. gifuensis is a good system for exploring the genetic and evolutionary effects of mass rearing and release on species commonly used as biocontrol agents.

Hybrid de novo genome assembly of the sexually dimorphic Lady Amherst's pheasant

Pheasants are an important group of birds, valued for their economic benefit as poultry birds, game birds, and as ornamental species for their plumage. Lady Amherst's pheasant Chrysolophus amherstiae is an ornamental species, valued for its elaborate and beautiful plumage. In this study, we present a high-quality de novo hybrid genome assembly of C. amherstiae. Previous attempts to sequence the genome of this species resulted in draft-level assemblies, which are not available in the public domain. Using a combination of Illumina short reads and Oxford Nanopore's long-reads, we assembled a high-quality genome of N50 ~3.9 Mb and near complete BUSCO assessment. We observed a correlation between effective population size and past climatic conditions, with an increase in population size during the warm interglacial periods. We further observed significant fluctuations in genes involved with the immune system and visual perception. C. amherstiae is a highly dimorphic species, and significant fluctuations in gene families involved in immune response, visual perception, among others, suggesting a role of mate choice and sexual selection in the evolution and maintenance of exaggerated traits in the males.

Whole-genome resequencing reveals candidate genes associated with milk production trait in Guanzhong dairy goats

Milk production is one of the most important economic utility of goats. Guanzhong dairy goat is a local dairy goat in Shaanxi Province of China and has high milk yield and quality. However, there are relatively few studies on molecular markers of milk production traits in Guanzhong dairy goats. In this study, we sequenced the whole genomes of 20 Guanzhong dairy goats, 10 of which had high milk yield (HM) and 10 of which had low milk yield (LM). We detected candidate signatures of selection in HM goats using Fst and π-ratio statistics and identified several candidate genes including ANPEP, ADRA1A and PRKG1 associated with milk production. Our results provide the basis for molecular breeding of milk production traits in Guanzhong dairy goats.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

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.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

Identification of genetic associations and functional SNPs of bovine KLF6 gene on milk production traits in Chinese holstein

BACKGROUND

Our previous research identified the Kruppel like factor 6 (KLF6) gene as a prospective candidate for milk production traits in dairy cattle. The expression of KLF6 in the livers of Holstein cows during the peak of lactation was significantly higher than that during the dry and early lactation periods. Notably, it plays an essential role in activating peroxisome proliferator-activated receptor α (PPARα) signaling pathways. The primary aim of this study was to further substantiate whether the KLF6 gene has significant genetic effects on milk traits in dairy cattle.

RESULTS

Through direct sequencing of PCR products with pooled DNA, we totally identified 12 single nucleotide polymorphisms (SNPs) within the KLF6 gene. The set of SNPs encompasses 7 located in 5' flanking region, 2 located in exon 2 and 3 located in 3' untranslated region (UTR). Of these, the g.44601035G > A is a missense mutation that resulting in the replacement of arginine (CGG) with glutamine (CAG), consequently leading to alterations in the secondary structure of the KLF6 protein, as predicted by SOPMA. The remaining 7 regulatory SNPs significantly impacted the transcriptional activity of KLF6 following mutation (P < 0.005), manifesting as changes in transcription factor binding sites. Additionally, 4 SNPs located in both the UTR and exons were predicted to influence the secondary structure of KLF6 mRNA using the RNAfold web server. Furthermore, we performed the genotype-phenotype association analysis using SAS 9.2 which found all the 12 SNPs were significantly correlated to milk yield, fat yield, fat percentage, protein yield and protein percentage within both the first and second lactations (P < 0.0001 ~ 0.0441). Also, with Haploview 4.2 software, we found the 12 SNPs linked closely and formed a haplotype block, which was strongly associated with five milk traits (P < 0.0001 ~ 0.0203).

CONCLUSIONS

In summary, our study represented the KLF6 gene has significant impacts on milk yield and composition traits in dairy cattle. Among the identified SNPs, 7 were implicated in modulating milk traits by impacting transcriptional activity, 4 by altering mRNA secondary structure, and 1 by affecting the protein secondary structure of KLF6. These findings provided valuable molecular insights for genomic selection program of dairy cattle.

Genetic diversity, population structure, and selective signature of sheep in the northeastern Tarim Basin

Local sheep in the northeastern Tarim Basin can adapt to dry and low-rainfall regional environments. In this study, three local sheep breeds in the northeastern Tarim Basin, LOP (LOP) sheep, Bayinbuluke (BYK) sheep, and Kunlun (KUN, also known as the Qiemo sheep) sheep, and three introduced sheep breeds, Suffolk (SUF) sheep, Dorset (APD) sheep, and Texel (TEX) sheep, were analyzed for genetic diversity, population structure, and selective signature using the Illumina OvineSNP50K BeadChip. We found that LOP, BYK, and KUN had lower observed heterozygosity and expected heterozygosity than TEX, SUF, and ADP, which were differentiated based on geographic distribution. We performed fixation index (FST) analysis on three local sheep breeds in the northeastern Tarim Basin (LOP, BYK, and KUN) and introduced sheep breeds (TEX, SUF, and ADP) to measure genetic differentiation. Nucleotide diversity (PI) analysis was performed on single-nucleotide polymorphism (SNP) data of LOP, BYK, and KUN. A total of 493 candidate genes were obtained by taking the intersection at a threshold of 5%. Among them, SMAD2, ESR2, and HAS2 were related to reproductive traits. PCDH15, TLE4, and TFAP2B were related to growth traits. SOD1, TSHR, and DNAJB5 were related to desert environmental adaptation. Analyzing the genetic patterns of local sheep in the northeastern Tarim Basin can protect the germplasm resources of local sheep and promote the development and utilization of sheep genetic resources.

Distribution of Runs of Homozygosity and Their Relationship with Candidate Genes for Productivity in Kazakh Meat-Wool Sheep Breed

Increasing the fertility of sheep remains one of the crucial issues of modern sheep breeding. The Kazakh meat-wool sheep is an excellent breed with high meat and wool productivity and well adapted to harsh conditions. Nowadays, runs of homozygosity (ROHs) are considered a suitable approach for studying the genetic characteristics of farm animals. The aims of the study were to analyze the distribution of ROHs, describe autozygosity, and detect genomic regions with high ROH islands. In this study, we genotyped a total of 281 Kazakh meat-wool sheep using the Illumina iScan® system (EquipNet, Canton, MA, USA) via Ovine SNP50 BeadChip array. As a results, a total of 15,069 ROHs were found in the three Kazakh meat-wool sheep populations. The mean number of ROH per animal across populations varied from 40.3 (POP1) to 42.2 (POP2) in the category 1+ Mb. Furthermore, the number of ROH per animal in ROH1-2 Mb were much higher than ROH2-4 Mb and ROH8-16 Mb in the three sheep populations. Most of individuals had small number of ROH>16 Mb. The highest and lowest genomic inbreeding coefficient values were observed in POP2 and POP3, respectively. The estimated FROH presented the impact that recent inbreeding has had in all sheep populations. Furthermore, a set of interesting candidate genes (BMP2, BMPR2, BMPRIB, CLOCK, KDM2B, TIAM1, TASP1, MYBPC1, MYOM1, and CACNA2D1), which are related to the productive traits, were found. Collectively, these findings will contribute to the breeding and conservation strategies of the Kazakh meat-wool sheep breed.

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.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

Novel Heredity Basis of the Four-Horn Phenotype in Sheep Using Genome-Wide Sequence Data

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10-8) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10-9) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = -156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

A computational framework for improving genetic variants identification from 5,061 sheep sequencing data

BACKGROUND

Pan-genomics is a recently emerging strategy that can be utilized to provide a more comprehensive characterization of genetic variation. Joint calling is routinely used to combine identified variants across multiple related samples. However, the improvement of variants identification using the mutual support information from multiple samples remains quite limited for population-scale genotyping.

RESULTS

In this study, we developed a computational framework for joint calling genetic variants from 5,061 sheep by incorporating the sequencing error and optimizing mutual support information from multiple samples' data. The variants were accurately identified from multiple samples by using four steps: (1) Probabilities of variants from two widely used algorithms, GATK and Freebayes, were calculated by Poisson model incorporating base sequencing error potential; (2) The variants with high mapping quality or consistently identified from at least two samples by GATK and Freebayes were used to construct the raw high-confidence identification (rHID) variants database; (3) The high confidence variants identified in single sample were ordered by probability value and controlled by false discovery rate (FDR) using rHID database; (4) To avoid the elimination of potentially true variants from rHID database, the variants that failed FDR were reexamined to rescued potential true variants and ensured high accurate identification variants. The results indicated that the percent of concordant SNPs and Indels from Freebayes and GATK after our new method were significantly improved 12%-32% compared with raw variants and advantageously found low frequency variants of individual sheep involved several traits including nipples number (GPC5), scrapie pathology (PAPSS2), seasonal reproduction and litter size (GRM1), coat color (RAB27A), and lentivirus susceptibility (TMEM154).

CONCLUSION

The new method used the computational strategy to reduce the number of false positives, and simultaneously improve the identification of genetic variants. This strategy did not incur any extra cost by using any additional samples or sequencing data information and advantageously identified rare variants which can be important for practical applications of animal breeding.

Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds

Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BACKGROUND

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

RESULTS

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

CONCLUSIONS

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

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.

A Single-Cell Atlas of an Early Mongolian Sheep Embryo

Cell types have been established during organogenesis based on early mouse embryos. However, our understanding of cell types and molecular mechanisms in the early embryo development of Mongolian sheep has been hampered. This study presents the first comprehensive single-cell transcriptomic characterization at E16 in Ujumqin sheep and Hulunbuir short-tailed sheep. Thirteen major cell types were identified at E16 in Ujumqin sheep, and eight major cell types were identified at E16 in Hulunbuir short-tailed sheep. Function enrichment analysis showed that several pathways were significantly enriched in the TGF-beta signaling pathway, the Hippo signaling pathway, the platelet activation pathway, the riboflavin metabolism pathway, the Wnt signaling pathway, regulation of the actin cytoskeleton, and the insulin signaling pathway in the notochord cluster. Glutathione metabolism, glyoxylate, and dicarboxylate metabolism, the citrate cycle, thyroid hormone synthesis, pyruvate metabolism, cysteine and methionine metabolism, thermogenesis, and the VEGF signaling pathway were significantly enriched in the spinal cord cluster. Steroid biosynthesis, riboflavin metabolism, the cell cycle, the Hippo signaling pathway, the Hedgehog signaling pathway, the FoxO signaling pathway, the JAK-STAT signaling pathway, and the Wnt signaling pathway were significantly enriched in the paraxial mesoderm cluster. The notochord cluster, spinal cord cluster, and paraxial mesoderm cluster were found to be highly associated with tail development. Pseudo-time analysis demonstrated that the mesenchyme can translate to the notochord in Ujumqin sheep. Molecular assays revealed that the Hippo signaling pathway was enriched in Ujumqin sheep. This comprehensive single-cell map revealed previously unrecognized signaling pathways that will further our understanding of the mechanism of short-tailed sheep formation.

The Effects of DDI1 on Inducing Differentiation in Ovine Preadipocytes via Oar-miR-432

Reducing fat deposition in sheep (Ovis aries) tails is one of the most important ways to combat rising costs and control consumer preference. Our previous studies have shown that oar-miR-432 is differentially expressed in the tail adipose tissue of Hu (a fat-tailed sheep breed) and Tibetan (a thin-tailed sheep breed) sheep and is a key factor in the negative regulation of fat deposition through BMP2 in ovine preadipocytes. This study investigated the effect of oar-miR-432 and its target genes in ovine preadipocytes. A dual luciferase assay revealed that DDI1 is a direct target gene of oar-miR-432. We transfected an oar-miR-432 mimic and inhibitor into preadipocytes to analyze the expression of target genes. Overexpression of oar-miR-432 inhibits DDI1 expression, whereas inhibition showed the opposite results. Compared with thin-tailed sheep, DDI1 was highly expressed in the fat-tailed sheep at the mRNA and protein levels. Furthermore, we transfected the overexpression and knockdown target genes into preadipocytes to analyze their influence after inducing differentiation. Knockdown of DDI1 induced ovine preadipocyte differentiation into adipocytes but suppressed oar-miR-432 expression. Conversely, the overexpression of DDI1 significantly inhibited differentiation but promoted oar-miR-432 expression. DDI1 overexpression also decreased the content of triglycerides. Additionally, DDI1 is a nested gene in intron 1 of PDGFD. When DDI1 was overexpressed, the PDGFD expression also increased, whereas DDI1 knockdown showed the opposite results. This is the first study to reveal the biological mechanisms by which oar-miR-432 inhibits preadipocytes through DDI1 and provides insight into the molecular regulatory mechanisms of DDI1 in ovine preadipocytes. These results have important applications in animal breeding and obesity-related human diseases.

Whole-genome resequencing of the native sheep provides insights into the microevolution and identifies genes associated with reproduction traits

BACKGROUND

Sheep genomes undergo numerous genes losses, gains and mutation that generates genome variability among breeds of the same species after long time natural and artificial selection. However, the microevolution of native sheep in northwest China remains elusive. Our aim was to compare the genomes and relevant reproductive traits of four sheep breeds from different climatic environments, to unveil the selection challenges that this species cope with, and the microevolutionary differences in sheep genomes. Here, we resequenced the genomes of 4 representative sheep breeds in northwest China, including Kazakh sheep and Duolang sheep of native breeds, and Hu sheep and Suffolk sheep of exotic breeds with different reproductive characteristics.

RESULTS

We found that these four breeds had a similar expansion experience from ~ 10,000 to 1,000,000 years ago. In the past 10,000 years, the selection intensity of the four breeds was inconsistent, resulting in differences in reproductive traits. We explored the sheep variome and selection signatures by F_ST and θπ. The genomic regions containing genes associated with different reproductive traits that may be potential targets for breeding and selection were detected. Furthermore, non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds with different reproductive characteristics were found. We identified PAK1, CYP19A1 and PER1 as a likely causal gene for seasonal reproduction in native sheep through qPCR, Western blot and ELISA analyses. Also, the haplotype frequencies of 3 tested gene regions related to reproduction were significantly different among four sheep breeds.

CONCLUSIONS

Our results provide insights into the microevolution of native sheep and valuable genomic information for identifying genes associated with important reproductive traits in sheep.

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.

Genome-Wide Signal Selection Analysis Revealing Genes Potentially Related to Sheep-Milk-Production Traits

Natural selection and domestication have shaped modern sheep populations into a vast range of phenotypically diverse breeds. Among these breeds, dairy sheep have a smaller population than meat sheep and wool sheep, and less research is performed on them, but the lactation mechanism in dairy sheep is critically important for improving animal-production methods. In this study, whole-genome sequences were generated from 10 sheep breeds, including 57 high-milk-yield sheep and 44 low-milk-yield sheep, to investigate the genetic signatures of milk production in dairy sheep, and 59,864,820 valid SNPs (Single Nucleotide Polymorphisms) were kept after quality control to perform population-genetic-structure analyses, gene-detection analyses, and gene-function-validation analyses. For the population-genetic-structure analyses, we carried out PCA (Principal Component Analysis), as well as neighbor-joining tree and structure analyses to classify different sheep populations. The sheep used in our study were well distributed in ten groups, with the high-milk-yield-group populations close to each other and the low-milk-yield-group populations showing similar classifications. To perform an exact signal-selection analysis, we used three different methods to find SNPs to perform gene-annotation analyses within the 995 common regions derived from the fixation index (F_ST), nucleotide diversity (Ɵπ), and heterozygosity rate (ZHp) results. In total, we found 553 genes that were located in these regions. These genes mainly participate in the protein-binding pathway and the nucleoplasm-interaction pathway, as revealed by the GO- and KEGG-function-enrichment analyses. After the gene selection and function analyses, we found that FCGR3A, CTSK, CTSS, ARNT, GHR, SLC29A4, ROR1, and TNRC18 were potentially related to sheep-milk-production traits. We chose the strongly selected genes, FCGR3A, CTSK, CTSS, and ARNT during the signal-selection analysis to perform a RT-qPCR (Reale time Quantitative Polymerase Chain Reaction) experiment to validate their expression-level relationship with milk production, and the results showed that FCGR3A has a significant negative relationship with sheep-milk production, while other three genes did not show any positive or negative relations. In this study, it was discovered and proven that the candidate gene FCGR3A potentially contributes to the milk production of dairy sheep and a basis was laid for the further study of the genetic mechanism underlying the strong milk-production traits of sheep.

Insertion/Deletion (InDel) Variants within the Sheep Fat-Deposition-Related PDGFD Gene Strongly Affect Morphological Traits

Platelet-derived growth factor D (PDGFD) is a member of the PDGF gene family, and it plays an important role in the regulation of adipocyte development in mammals. Furthermore, genome-wide association studies (GWAS) have previously identified it as a candidate gene associated with fleece fiber variation, body size, and the fat-tail phenotype in domestic Chinese sheep. In this study, a total of 1919 indigenous Chinese sheep were genotyped to examine the association between nucleotide sequence variations in PDGFD and body morphology. Our results detected both a 14 bp insertion in intron 2 and a 13 bp deletion in intron 4 of PDGFD. Moreover, these two InDel loci had low to moderate polymorphism. Notably, the 13 bp deletion mutation of PDGFD was found to significantly affect sheep body size. Yearling rams in the Luxi black-headed sheep (LXBH) containing a heterozygous genotype (insertion/deletion, ID) were found to have larger body length, chest depth, and body weight than those with wild genotypes. Furthermore, adult ewes in the Guiqian semi-fine wool sheep (GSFW) containing a homozygous mutation (deletion/deletion, DD) were found to have smaller chest width than their peers. Moreover, yearling ewes in this group with the same homozygous mutation were found to have lower body weight, chest width, and cannon circumference compared to those of other individuals. This study demonstrates that PDGFD InDel polymorphisms have the potential to be effective molecular markers to improve morphological traits in domestic Chinese sheep.

Expression profiles of the CD274 and PLEKHH2 gene and association of its polymorphism with hematologic parameters in sheep

CD274 and PLEKHH2 genes have been identified as immune- and multiple diseases-related genes, and have recently garnered significant interest. However, their role in regulating immune functions in sheep remains largely unexplored. In this study, we aimed to investigate the effects of polymorphisms in CD274 and PLEKHH2 on hematologic parameters in 915 sheep. Our results showed that the CD274 and PLEKHH2 genes were most highly expressed in the spleen and tail fat, respectively, as determined by qRT-PCR. We also identified a G to A mutation (g 0.11858 G > A) in the exon 4 region of CD274, and a C to G mutation (g 0.38384 C > G) in the intron 8 region of PLEKH2. Association analysis revealed that CD274 g 0.11858 G > A was significantly associated with RBC, HCT, MCHC, and MCV (P < 0.05), while PLEKHH2 g 0.38384 C > G was significantly associated with HCT, MPV, MCHC, and MCV (P < 0.05). These results suggest that CD274 and PLEKHH2 genes may play a role in regulating blood physiological indicators and could be potential functional candidates for influencing immune traits in sheep breeding programs.

Identification of Critical Genes for Ovine Horn Development Based on Transcriptome during the Embryonic Period

Horns, also known as headgear, are a unique structure of ruminants. As ruminants are globally distributed, the study of horn formation is critical not only for increasing our understanding of natural and sexual selection but also for the breeding of polled sheep breeds to facilitate modern sheep farming. Despite this, a significant number of the underlying genetic pathways in sheep horn remain unclear. In this study, to clarify the gene expression profile of horn buds and investigate the key genes in horn bud formation, RNA-sequencing (RNA-seq) technology was utilized to investigate differential gene expression in the horn buds and adjacent forehead skin of Altay sheep fetuses. There were only 68 differentially expressed genes (DEGs) identified, consisting of 58 up-regulated genes and 10 down-regulated genes. RXFP2 was differentially up-regulated in the horn buds and had the highest significance (p-value = 7.42 × 10^-14). In addition, 32 DEGs were horn-related genes identified in previous studies, such as RXFP2, FOXL2, SFRP4, SFRP2, KRT1, KRT10, WNT7B, and WNT3. Further, Gene Ontology (GO) analysis showed that the DEGs were mainly enriched with regard to growth, development, and cell differentiation. Pathway analysis revealed that the Wnt signaling pathway may be responsible for horn development. Further, through combining the protein-protein interaction networks of the DEGs, it was found that the top five hub genes, namely, ACAN, SFRP2, SFRP4, WNT3, and WNT7B, were also associated with horn development. Our results suggest that only a few key genes, including RXFP2, are involved in bud formation. This study not only validates the expression of candidate genes identified at the transcriptome level in previous studies but also provides new possible marker genes for horn development, which may promote our understanding of the genetic mechanisms of horn formation.

Whole genome sequencing of simmental cattle for SNP and CNV discovery

BACKGROUD

The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity.

RESULTS

In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility.

CONCLUSIONS

In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

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.

Population Structure and Selection Signatures of Domestication in Geese

The goose is an economically important poultry species and was one of the first to be domesticated. However, studies on population genetic structures and domestication in goose are very limited. Here, we performed whole genome resequencing of geese from two wild ancestral populations, five Chinese domestic breeds, and four European domestic breeds. We found that Chinese domestic geese except Yili geese originated from a common ancestor and exhibited strong geographical distribution patterns and trait differentiation patterns, while the origin of European domestic geese was more complex, with two modern breeds having Chinese admixture. In both Chinese and European domestic geese, the identified selection signatures during domestication primarily involved the nervous system, immunity, and metabolism. Interestingly, genes related to vision, skeleton, and blood-O2 transport were also found to be under selection, indicating genetic adaptation to the captive environment. A forehead knob characterized by thickened skin and protruding bone is a unique trait of Chinese domestic geese. Interestingly, our population differentiation analysis followed by an extended genotype analysis in an additional population suggested that two intronic SNPs in EXT1, an osteochondroma-related gene, may plausibly be sites responsible for knob. Moreover, CSMD1 and LHCGR genes were found to be significantly associated with broodiness in Chinese domestic geese and European domestic geese, respectively. Our results have important implications for understanding the population structure and domestication of geese, and the selection signatures and variants identified in this study might be useful in genetic breeding for forehead knob and reproduction traits.

Identification of Genomic Signatures in Bullmastiff Dogs Using Composite Selection Signals Analysis of 23 Purebred Clades

Dog breeds represent canine sub-populations with distinctive phenotypic features and limited genetic diversity. We have established a resource to study breed-specific genetic diversity. Utilising genetic resources within our laboratory biobank, public domain genotype data and the phylogenetic framework of 23 breed clades, the primary objective for this study was to identify genomic regions that differentiate the Bullmastiff breed. Through application of a composite index analysis (CSS), genomic signatures were identified in Bullmastiffs when compared to the formative breeds, Mastiffs and Bulldogs, and to 22 other breed groups. Significant regions were identified on 15 chromosomes, with the most differentiated regions found on CFA1, CFA9, and CFA18. These regions may reflect genetic drift following establishment of the breed or the effects of selective breeding during development of the modern Bullmastiff. This was supported by analysis of genes from the identified genomic regions, including 458 genes from the multi-clade analysis, which revealed enriched pathways that may be related to characteristic traits and distinct morphology of the breed. The study demonstrates the utility of the CSS method in breed-specific genome analysis and advances our understanding of genetic diversity in Bullmastiff dogs.

Resequencing of a Pekin duck breeding population provides insights into the genomic response to short-term artificial selection

BACKGROUND

Short-term, intense artificial selection drives fast phenotypic changes in domestic animals and leaves imprints on their genomes. However, the genetic basis of this selection response is poorly understood. To better address this, we employed the Pekin duck Z2 pure line, in which the breast muscle weight was increased nearly 3-fold after 10 generations of breeding. We denovo assembled a high-quality reference genome of a female Pekin duck of this line (GCA_003850225.1) and identified 8.60 million genetic variants in 119 individuals among 10 generations of the breeding population.

RESULTS

We identified 53 selected regions between the first and tenth generations, and 93.8% of the identified variations were enriched in regulatory and noncoding regions. Integrating the selection signatures and genome-wide association approach, we found that 2 regions covering 0.36 Mb containing UTP25 and FBRSL1 were most likely to contribute to breast muscle weight improvement. The major allele frequencies of these 2 loci increased gradually with each generation following the same trend. Additionally, we found that a copy number variation region containing the entire EXOC4 gene could explain 1.9% of the variance in breast muscle weight, indicating that the nervous system may play a role in economic trait improvement.

CONCLUSIONS

Our study not only provides insights into genomic dynamics under intense artificial selection but also provides resources for genomics-enabled improvements in duck breeding.

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.

Selection signature analysis reveals RDH5 performed key function in vision during sheep domestication process

Abstract. As one of the most successful domesticated animals in the Neolithic age, sheep gradually migrated all over the world with human activities. During the domestication process, remarkable changes have taken place in morphology, physiology, and behavior, resulting in different breeds with different characters via artificial and natural selection. However, the genetic background responsible for these phenotypic variations remains largely unclear. Here, we used whole genome resequencing technology to compare and analyze the genome differences between Asiatic mouflon wild sheep (Ovis orientalis) and Hu sheep (Ovis aries). A total of 755 genes were positively selected in the process of domestication and selection, and the genes related to sensory perception had directional evolution in the autosomal region, such as OPRL1, LEF1, TAS1R3, ATF6, VSX2, MYO1A, RDH5, and some novel genes. A missense mutation of c.T722C/p.M241T in exon 4 of RDH5 existing in sheep were found, and the T allele was completely fixed in Hu sheep. In addition, the mutation with the C allele reduced the retinol dehydrogenase activity encoding by RDH5, which can impair retinoic acid metabolism and further influenced the visual cycle. Overall, our results showed significant enrichment for positively selected genes involved in sensory perception development during sheep domestication; RDH5 and its variants may be related to the retinal degeneration in sheep. We infer that the wild sheep ancestors with weaker visual sensitivity were weeded out by humans, and the mutation was selective, swept by the dual pressures of natural and artificial selection.

Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds

There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.