Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

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.

Expression, Polymorphism, and Potential Functional Sites of the BMPR1A Gene in the Sheep Horn

Sheep horns are composed of bone and sheaths, and the BMPR1A gene is required for cartilage and osteogenic differentiation. Therefore, the BMPR1A gene may have a function related to the sheep horn, but its relationship with the sheep horn remains unclear. In this study, we first utilized RNA sequencing (RNA-seq) data to investigate the expression of the BMPR1A gene in different tissues and breeds of sheep. Second, whole-genome sequencing (WGS) data were used to explore the functional sites of the BMPR1A gene. Lastly, the allele-specific expression of the BMPR1A gene was explored. Our results indicate that BMPR1A gene expression is significantly higher in the normal horn groups than in the scurred groups. Importantly, this trend is consistent across several sheep breeds. Therefore, this finding suggests that the BMPR1A gene may be related to horn type. A total of 43 Single-Nucleotide Polymorphisms (SNPs) (F-statistics > 0.15) and 10 allele-specific expressions (ASEs) exhibited difference between the large and small horn populations. It is probable that these sites significantly impact the size of sheep horns. Compared to other polled species, we discovered ten amino acid sites that could influence horn presence. By combining RNA-seq and WGS functional loci results, we identified a functional site at position 40574836 on chromosome 25 that is both an SNP and exhibits allele-specific expression. In conclusion, we demonstrated that the BMPR1A gene is associated with horn type and identified some important functional sites which can be used as molecular markers in the breeding of sheep horns.

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.

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 Resequencing of Ujimqin Sheep Identifies Genes Associated with Vertebral Number

The number of vertebrae is a crucial economic trait that can significantly impact the carcass length and meat production in animals. However, our understanding of the quantitative trait loci (QTLs) and candidate genes associated with the vertebral number in sheep (Ovis aries) remains limited. To identify these candidate genes and QTLs, we collected 73 Ujimqin sheep with increased numbers of vertebrae (T13L7, T14L6, and T14L7) and 23 sheep with normal numbers of vertebrae (T13L6). Through high-throughput genome resequencing, we obtained a total of 24,130,801 effective single-nucleotide polymorphisms (SNPs). By conducting a selective-sweep analysis, we discovered that the most significantly selective region was located on chromosome 7. Within this region, we identified several genes, including VRTN, SYNDIG1L, LTBP2, and ABCD4, known to regulate the spinal development and morphology. Further, a genome-wide association study (GWAS) performed on sheep with increased and normal vertebral numbers confirmed that ABCD4 is a candidate gene for determining the number of vertebrae in sheep. Additionally, the most significant SNP on chromosome 7 was identified as a candidate QTL. Moreover, we detected two missense mutations in the ABCD4 gene; one of these mutations (Chr7: 89393414, C > T) at position 22 leads to the conversion of arginine (Arg) to glutamine (Gln), which is expected to negatively affect the protein's function. Notably, a transcriptome expression profile in mouse embryonic development revealed that ABCD4 is highly expressed during the critical period of vertebral formation (4.5-7.5 days). Our study highlights ABCD4 as a potential major gene influencing the number of vertebrae in Ujimqin sheep, with promising prospects for future genome-assisted breeding improvements in sheep.

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.

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.

Whole-Genome Resequencing Reveals Selection Signatures of Abigar Cattle for Local Adaptation

Over time, indigenous cattle breeds have developed disease resistance, heat tolerance, and adaptability to harsh environments. Deciphering the genetic mechanisms underlying adaptive traits is crucial for their improvement and sustainable utilization. For the first time, we performed whole-genome sequencing to unveil the genomic diversity, population structure, and selection signatures of Abigar cattle living in a tropical environment. The population structure analysis revealed that Abigar cattle exhibit high nucleotide diversity and heterozygosity, with low runs of homozygosity and linkage disequilibrium, suggesting a genetic landscape less constrained by inbreeding and enriched by diversity. Using nucleotide diversity (Pi) and population differentiation (FST) selection scan methods, we identified 83 shared genes that are likely associated with tropical adaption. The functional annotation analysis revealed that some of these genes are potentially linked to heat tolerance (HOXC13, DNAJC18, and RXFP2), immune response (IRAK3, MZB1, and STING1), and oxidative stress response (SLC23A1). Given the wider spreading impacts of climate change on cattle production, understanding the genetic mechanisms of adaptation of local breeds becomes crucial to better respond to climate and environmental changes. In this context, our finding establishes a foundation for further research into the mechanisms underpinning cattle adaptation to tropical environments.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

Genetic diversity, tissue-specific expression, and functional analysis of the ATP7A gene in sheep

In humans, variation of the ATP7A gene may cause cranial exostosis, which is similar to "human horn," but the function of the ATP7A gene in sheep is still unknown. Tissue expression patterns and potential functional loci analysis of the ATP7A gene could help understand its function in sheep horn. In this study, we first identified tissue, sex, breed, and species-specific expression of the ATP7A gene in sheep based on the RNA-sequencing (RNA-seq) data. Second, the potential functional sites of the ATP7A gene were analyzed by using the whole genome sequencing (WGS) data of 99 sheep from 10 breeds. Last, the allele-specific expression of the ATP7A gene was explored. Our result showed the ATP7A gene has significantly higher expression in the big horn than in the small horn, and the ATP7A gene has high expression in the horn and skin, suggesting that this gene may be related to the horn. The PCA results show that the region around the ATP7A can distinguish horned and hornless groups to some extent, further indicating that the ATP7A may be related to horns. When compared with other species, we find seven ruminate specific amino acid sites of the ATP7A protein, which can be important to the ruminate horn. By analyzing WGS, we found 6 SNP sites with significant differences in frequency in horned and hornless populations, and most of these variants are present in the intron. But we still find some potential functional sites, including three missenses, three synonymous mutations, and four Indels. Finally, by combining the RNA-seq and WGS functional loci results, we find three mutations that showed allele-specific expression between big and small horns. This study shows that the ATP7A gene in sheep may be related to horn size, and several potential functional sites we identified here can be useful molecular markers for sheep horn breeding.

Whole genome sequencing reveals signals of adaptive admixture in Creole cattle

The Creole cattle from Guadeloupe (GUA) are well adapted to the tropical environment. Its admixed genome likely played an important role in such adaptation. Here, we sought to detect genomic signatures of selection in the GUA genome. For this purpose, we sequenced 23 GUA individuals and combined our data with sequenced genomes of 99 animals representative of European, African and indicine groups. We detect 17,228,983 single nucleotide polymorphisms (SNPs) in the GUA genome, providing the most detailed exploration, to date, of patterns of genetic variation in this breed. We confirm the higher level of African and indicine ancestries, compared to the European ancestry and we highlight the African origin of indicine ancestry in the GUA genome. We identify five strong candidate regions showing an excess of indicine ancestry and consistently supported across the different detection methods. These regions encompass genes with adaptive roles in relation to immunity, thermotolerance and physical activity. We confirmed a previously identified horn-related gene, RXFP2, as a gene under strong selective pressure in the GUA population likely owing to human-driven (socio-cultural) pressure. Findings from this study provide insight into the genetic mechanisms associated with resilience traits in livestock.

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.

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.

Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep

BACKGROUND

The detection of selective traits in different populations can not only reveal current mechanisms of artificial selection for breeding, but also provide new insights into phenotypic variation in new varieties and the search for genes associated with important traits. Panou sheep is a cultivated breed of Tibetan sheep in China with stable genetic performance, consistent appearance and fast growth and development after decades of artificial selection and cultivation. Due to long-term adaptation to the high altitude, cold and hypoxic environment in the plateau area, they may have formed a unique gene pool that is different from other Tibetan sheep breeds. To explore the genetic resources of Panou sheep, we used next-generation sequencing technology for the first time to investigate the genome-wide population structure, genetic diversity, and candidate signatures of positive selection in Panou sheep.

RESULTS

Comparative genomic analysis with the closely related species Oula sheep (a native breed of Tibetan sheep in China) was used to screen the population selection signal of Panou sheep. Principal component analysis and neighbor joining tree showed that Panou sheep and Oula sheep had differences in population differentiation. Furthermore, analyses of population structure, they came from the same ancestor, and when K = 2, the two populations could be distinguished. Panou sheep exhibit genetic diversity comparable to Oula sheep, as shown by observed heterozygosity, expected heterozygosity and runs of homozygosity. Genome-wide scanning using the Fst and π ratio methods revealed a list of potentially selected related genes in Panou sheep compared to Oula sheep, including histone deacetylase 9 (HDAC9), protein tyrosine kinase 2 (PTK2), microphthalmia-related transcription factor (MITF), vesicular amine transporter 1 (VAT1), trichohyalin-like 1 (TCHHL1), amine oxidase, copper containing 3 (AOC3), interferon-inducible protein 35 (IFI35).

CONCLUSIONS

The results suggest that traits related to growth and development and plateau adaptation may be selection targets for the domestication and breeding improvement of Tibetan sheep. This study provides the fundamental footprints for Panou sheep breeding and management.

Genomic adaptation of Ethiopian indigenous cattle to high altitude

The mountainous areas of Ethiopia represent one of the most extreme environmental challenges in Africa faced by humans and other inhabitants. Selection for high-altitude adaptation is expected to have imprinted the genomes of livestock living in these areas. Here we assess the genomic signatures of positive selection for high altitude adaptation in three cattle populations from the Ethiopian mountainous areas (Semien, Choke, and Bale mountains) compared to three Ethiopian lowland cattle populations (Afar, Ogaden, and Boran), using whole-genome resequencing and three genome scan approaches for signature of selection (iHS, XP-CLR, and PBS). We identified several candidate selection signature regions and several high-altitude adaptation genes. These include genes such as ITPR2, MB, and ARNT previously reported in the human population inhabiting the Ethiopian highlands. Furthermore, we present evidence of strong selection and high divergence between Ethiopian high- and low-altitude cattle populations at three new candidate genes (CLCA2, SLC26A2, and CBFA2T3), putatively linked to high-altitude adaptation in cattle. Our findings provide possible examples of convergent selection between cattle and humans as well as unique African cattle signature to the challenges of living in the Ethiopian mountainous regions.

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

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Whole-genome resequencing reveals domestication and signatures of selection in Ujimqin, Sunit, and Wu Ranke Mongolian sheep breeds

Objective: The current study aimed to perform whole-genome resequencing of Chinese indigenous Mongolian sheep breeds including Ujimqin, Sunit, and Wu Ranke sheep breeds (UJMQ, SNT, WRK) and deeply analyze genetic variation, population structure, domestication, and selection for domestication traits among these Mongolian sheep breeds.Methods: Blood samples were collected from a total of 60 individuals comprising 20 WRK, 20 UJMQ, and 20 SNT. For genome sequencing, about 1.5 μg of genomic DNA was used for library construction with an insert size of about 350 bp. Pair-end sequencing were performed on Illumina NovaSeq platform, with the read length of 150 bp at each end. We then investigated the domestication and signatures of selection in these sheep breeds.Results: According to the population and demographic analyses, WRK and SNT populations were very similar, which were different from UJMQ populations. Genome wide association study identified 468 and 779 significant loci from SNT vs UJMQ, and UJMQ vs WRK, respectively. However, only 3 loci were identified from SNT vs WRK. Genomic comparison and selective sweep analysis among these sheep breeds suggested that genes associated with regulation of secretion, metabolic pathways including estrogen metabolism and amino acid metabolism, and neuron development have undergone strong selection during domestication.Conclusion: Our findings will facilitate the understanding of Chinese indigenous Mongolian sheep breeds domestication and selection for complex traits and provide a valuable genomic resource for future studies of sheep and other domestic animal breeding.

Natural adaptation and human selection of northeast African sheep genomes

African sheep manifest diverse but distinct physio-anatomical traits, which are the outcomes of natural- and human-driven selection. Here, we generated 34.8 million variants from 150 indigenous northeast African sheep genomes sequenced at an average depth of ∼54× for 130 samples (Ethiopia, Libya) and ∼20× for 20 samples (Sudan). These represented sheep from diverse environments, tail morphology and post-Neolithic introductions to Africa. Phylogenetic and model-based admixture analysis provided evidence of four genetic groups corresponding to altitudinal geographic origins, tail morphotypes and possible historical introduction and dispersal of the species into and across the continent. Running admixture at higher levels of K (6 ≤ K ≤ 25), revealed cryptic levels of genome intermixing as well as distinct genetic backgrounds in some populations. Comparative genomic analysis identified targets of selection that spanned conserved haplotype structures overlapping clusters of genes and gene families. These were related to hypoxia responses, ear morphology, caudal vertebrae and tail skeleton length, and tail fat-depot structures. Our findings provide novel insights underpinning morphological variation and response to human-driven selection and environmental adaptation in African indigenous sheep.

Genomic analyses of wild argali, domestic sheep, and their hybrids provide insights into chromosome evolution, phenotypic variation, and germplasm innovation

Understanding the genetic mechanisms of phenotypic variation in hybrids between domestic animals and their wild relatives may aid germplasm innovation. Here, we report the high-quality genome assemblies of a male Pamir argali (O ammon polii, 2n = 56), a female Tibetan sheep (O aries, 2n = 54), and a male hybrid of Pamir argali and domestic sheep, and the high-throughput sequencing of 425 ovine animals, including the hybrids of argali and domestic sheep. We detected genomic synteny between Chromosome 2 of sheep and two acrocentric chromosomes of argali. We revealed consistent satellite repeats around the chromosome breakpoints, which could have resulted in chromosome fusion. We observed many more hybrids with karyotype 2n = 54 than with 2n = 55, which could be explained by the selfish centromeres, the possible decreased rate of normal/balanced sperm, and the increased incidence of early pregnancy loss in the aneuploid ewes or rams. We identified genes and variants associated with important morphological and production traits (e.g., body weight, cannon circumference, hip height, and tail length) that show significant variations. We revealed a strong selective signature at the mutation (c.334C > A, p.G112W) in TBXT and confirmed its association with tail length among sheep populations of wide geographic and genetic origins. We produced an intercross population of 110 F2 offspring with varied number of vertebrae and validated the causal mutation by whole-genome association analysis. We verified its function using CRISPR-Cas9 genome editing. Our results provide insights into chromosomal speciation and phenotypic evolution and a foundation of genetic variants for the breeding of sheep and other animals.

The Complex and Diverse Genetic Architecture of the Absence of Horns (Polledness) in Domestic Ruminants, including Goats and Sheep

Horns are the most obvious common feature of Bovidae. The naturally occurring absence of horns in these species, also known as polledness, is of surprisingly heterogeneous nature, although they are Mendelian traits. This review compares in detail the molecular differences among the causes of inherited polledness in the domestic ruminant species of cattle, yak, sheep, and goat based on the causal gene variants that have been discovered in recent years. The genetic causes for the lack of horns in small ruminants seem not only to be more complex, e.g., in sheep, breed-specific characteristics are still unexplained, but in goats, there is also the associated disorder of intersexuality—polled intersex syndrome (PIS). In connection with animal welfare and the associated discussion about a legal ban on the dehorning of all farm animals, naturally hornless animals and the causal genetic variants are of increasing research interest in the age of genome editing. However, the low acceptance of genetic engineering in livestock, especially in European societies, limits its use in food-producing animals. Therefore, genotype-based targeted selection of naturally occurring variants is still a widely used method for spreading this desired trait within and across populations, at least in cattle and sheep.

Whole genome sequencing of Luxi Black Head sheep for screening selection signatures associated with important traits

Objective

Luxi Black Head sheep (LBH) is the first crossbreed specialized for meat production and was developed by crossbreeding Black Head Dorper sheep (DP) and Small Tailed Han sheep (STH) in the farming areas of northern China. Research on the genomic variations and selection signatures of LBH caused by continuous artificial selection is of great significance for identifying the genetic mechanisms of important traits of sheep and for the continuous breeding of LBH.

Methods

We explored the genetic relationships of LBH, DP, and several Mongolian sheep breeds by constructing phylogenetic tree, principal component analysis and linkage disequilibrium analysis. In addition, we analysed 29 whole genomes of sheep. The genome-wide selection signatures have been scanned with four methods: heterozygosity (H_P), fixation index (F_ST), cross-population extended haplotype homozygosity (XP-EHH) and the nucleotide diversity (θ_π) ratio.

Results

The genetic relationships analysis showed that LBH appeared to be an independent cluster closer to DP. The candidate signatures of positive selection in sheep genome revealed candidate genes for developmental process (HoxA gene cluster, BCL2L11, TSHR), immunity (CXCL6, CXCL1, SKAP2, PTK6, MST1R), growth (PDGFD, FGF18, SRF, SOCS2), and reproduction (BCAS3, TRIM24, ASTL, FNDC3A). Moreover, two signalling pathways closely related to reproduction, the thyroid hormone signalling pathway and the oxytocin signalling pathway, were detected.

Conclusion

The selective sweep analysis of LBH genome revealed candidate genes and signalling pathways associated with developmental process, immunity, growth, and reproduction. Our findings provide a valuable resource for sheep breeding and insight into the mechanisms of artificial selection.

Whole-genome resequencing to investigate the determinants of the multi-lumbar vertebrae trait in sheep

Multi-lumbar vertebrae trait is a beneficial mutation that can significantly improve livestock meat production. However, the genetic basis of the multi-lumbar vertebrae in sheep is still unclear. Here, we analysed the number of lumbar vertebrae of Duolang sheep and found three different traits of lumbar vertebrae number. Compared with the normal sheep, the length and weight of animal carcass from the multi-lumbar vertebrae sheep increased by 2.21 cm and 0.78 kg, respectively. We performed high-throughput genome resequencing on multi-lumbar vertebrae (n = 18) and normal (n = 11) Duolang sheep and obtained a total of more than 528.87 GB data. We found that the most significantly selective region were located in the 49.68-49.74 MB of chromosome 4 by selective-sweep analysis. We annotated this region and found that it contains SFRP4 which is known to regulate bone development. We further used the PCR-SSCP technology to detect the single nucleotide polymorphism (SNP) of the putative candidate SFRP4 and found that the two SNPs (rs600370085:C > T and rs415133338: A > G) of this gene were significantly associated with the multi-lumbar vertebrae of Duolang sheep. Our study indicates that the SFRP4 may be a potential major gene that affects the number of lumbar vertebrae in Duolang sheep, and has the potential to be utilized for sheep breeding in the future.

Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci

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

Genome-wide association study reveals 14 new SNPs and confirms two structural variants highly associated with the horned/polled phenotype in goats

Background

There is a long-term interest in investigating the genetic basis of the horned/polled phenotype in domestic goats. Here, we report a genome-wide association study (GWAS) to detect the genetic loci affecting the polled phenotype in goats.

Results

We obtained a total of 13,980,209 biallelic SNPs, using the genotyping-by-sequencing data from 45 Jintang Black (JT) goats, which included 32 female and nine male goats, and four individuals with the polled intersex syndrome (PIS). Using a mixed-model based GWAS, we identified two association signals, which were located at 150,334,857–150,817,260 bp (P = 5.15 × 10^− 119) and 128,286,704–131,306,537 bp (P = 2.74 × 10^− 15) on chromosome 1. The genotype distributions of the 14 most significantly associated SNPs were completely correlated with horn status in goats, based on the whole-genome sequencing (WGS) data from JT and two other Chinese horned breeds. However, variant annotation suggested that none of the detected SNPs within the associated regions were plausible causal mutations. Via additional read-depth analyses and visual inspections of WGS data, we found a 10.1-kb deletion (CHI1:g. 129424781_129434939del) and a 480-kb duplication (CHI1:150,334,286–150,818,098 bp) encompassing two genes KCNJ15 and ERG in the associated regions of polled and PIS-affected goats. Notably, the 10.1-kb deletion also served as the insertion site for the 480-kb duplication, as validated by PCR and Sanger sequencing. Our WGS genotyping showed that all horned goats were homozygous for the reference alleles without either the structural variants (SVs), whereas the PIS-affected goats were homozygous for both the SVs. We also demonstrated that horned, polled, and PIS-affected individuals among 333 goats from JT and three other Chinese horned breeds can be accurately classified via PCR amplification and agarose gel electrophoresis of two fragments in both SVs.

Conclusion

Our results revealed that two genomic regions on chromosome 1 are major loci affecting the polled phenotypes in goats. We provided a diagnostic PCR to accurately classify horned, polled, and PIS-affected goats, which will enable a reliable genetic test for the early-in-life prediction of horn status in goats.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08089-w.

Combined approaches identify known and novel genes associated with sheep litter size and non-seasonal breeding

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.

Trends towards revealing the genetic architecture of sheep tail patterning: Promising genes and investigatory pathways

Different sheep breeds have evolved after initial domestication, generating various tail phenotypic patterns. The phenotypic diversity of sheep tail patterns offers ideal materials for comparative analysis of its genetic basis. Evolutionary biologists, animal geneticists, breeders, and producers have been curious to clearly understand the underlying genetics behind phenotypic differences in sheep tails. Understanding the causal gene(s) and mutation(s) underlying these differences will help probe an evolutionary riddle, improve animal production performance, promote animal welfare, and provide lessons that help comprehend human diseases related to fat deposition (i.e., obesity). Historically, fat tails have served as an adaptive response to aridification and climate change. However, the fat tail is currently associated with compromised mating and animal locomotion, fat distribution in the animal body, increased raising costs, reduced consumer preference, and other animal welfare issues such as tail docking. The developing genomic approaches provide unprecedented opportunities to determine causal variants underlying phenotypic differences among populations. In the last decade, researchers have performed several genomic investigations to assess the genomic causality underlying phenotypic variations in sheep tails. Various genes have been suggested with the prominence of several potentially significant causatives, including the BMP2 and PDGFD genes associated with the fat tail phenotype and the TBXT gene linked with the caudal vertebrae number and tail length. Although the potential genes related to sheep tail characteristics have been revealed, the causal variant(s) and mutation(s) of these high-ranking candidate genes are still elusive and need further investigation. The review discusses the potential genes, sheds light on a knowledge gap, and provides possible investigative approaches that could help determine the specific genomic causatives of sheep tail patterns. Besides, characterizing and revealing the genetic determinism of sheep tails will help solve issues compromising sheep breeding and welfare in the future.

De-Domestication: An Extension of Crop Evolution

De-domestication or feralization is an interesting phenomenon in crops and livestock. Previously, evidence for crop de-domestication was based mainly on studies using phenotypic and genotypic data from limited molecular markers or gene segments. Recent genomic studies in rice, barley, and wheat provide comprehensive landscapes of de-domestication on a whole-genome scale. Here, we summarize crop de-domestication processes, ecological roles of de-domesticates, mechanisms underlying crop de-domestication syndromes, and conditions potentially favoring de-domestication events. We further explain how recent de-domestication studies have expanded our understanding of the complexity of crop evolution, and highlight the genetic novelties of de-domesticates beneficial for modern crop breeding.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

The expression and mutation of BMPR1B and its association with litter size in small-tail Han sheep (Ovis aries)

Previous studies have shown that BMPR1B promotes follicular development and ovarian granulosa cell proliferation, thereby affecting ovulation in mammals. In this study, the expression and polymorphism of the BMPR1B gene associated with litter size in small-tail Han (STH) sheep were determined. The expression of BMPR1B was detected in 14 tissues of STH sheep during the follicular phase as well as in the hypothalamic–pituitary–gonadal (HPG) axis of monotocous and polytocous STH sheep during the follicular and luteal phases using quantitative polymerase chain reaction (qPCR). Sequenom MassARRAY^® single nucleotide polymorphism (SNP) technology was also used to detect the polymorphism of SNPs in seven sheep breeds. Here, BMPR1B was highly expressed in hypothalamus, ovary, uterus, and oviduct tissue during the follicular phase, and BMPR1B was expressed significantly more in the hypothalamus of polytocous ewes than in monotocous ewes during both the follicular and luteal phases (P<0.05). For genotyping, we found that genotype and allele frequencies of three loci of the BMPR1B gene were extremely significantly different (P<0.01) between the monotocous and polytocous groups. Association analysis results showed that the g.29380965A>G locus had significant negative effects on the litter size of STH sheep, and the combination of g.29380965A>G and FecB (Fec – fecundity and B – Booroola; A746G) at the BMPR1B gene showed that the litter size of AG–GG, AA–GG, and GG–GG genotypes was significantly higher compared with other genotypes (P<0.05). This is the first study to find a new molecular marker affecting litter size and to systematically analyze the expression of BMPR1B in different fecundity and physiological periods of STH sheep.

Polymorphism, expression and structure analysis of key genes in the ovarian steroidogenesis pathway in sheep (Ovis aries)

Background

Litter size is an important factor that significantly affects the development of the sheep industry. Our previous TMT proteomics analysis found that three key proteins in the ovarian steroidogenesis pathway, STAR, HSD3B1, and CYP11A1, may affect the litter size trait of Small Tail Han sheep.

Objective

The purpose of this study was to better understand the relationship between polymorphisms of these three genes and litter size.

Material and Method

Sequenom MassARRAY detected genetic variance of the three genes in 768 sheep. Real‐time qPCR of the three genes was used to compare their expression in monotocous and polytocous sheep in relevant tissues. Finally, bioinformatics analysis predicted the protein sequences of the different SNP variants.

Result

Association analysis showed that there was a significant difference in litter size among the genotypes at two loci of the CYP11A1 gene (p < 0.05), but no significant difference was observed in litter size among all genotypes at all loci of the STAR and HSD3B1 genes (p > 0.05). However, STAR expression was significantly different in polytocous and monotocous sheep in the pituitary (p < 0.01). Tissue‐specific expression in the ovary was observed for HSD3B1 (p < 0.05), but its expression was not different between polytocous and monotocous sheep. Bioinformatics analysis showed that the g.33217408C > T mutation of CYP11A1 resulted in major changes to the secondary and tertiary structures. In contrast, gene polymorphisms in STAR and HSD3B1 had minimal impacts on their protein structures.

Discussion

This may explain why the CYP11A1 variant impacted litter size while the others did not. The single nucleotide polymorphism of the CYP11A1 gene would serve as a good molecular marker when breeding to increase litter size in sheep. Our study provides a basis for further revealing the function of the ovarian steroidogenesis pathway in sheep reproduction and sheep breeding.

Shared Ancestry and Signatures of Recent Selection in Gotland Sheep

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

Expression, structure and function analysis of the sperm-oocyte fusion genes Juno and Izumo1 in sheep (Ovis aries)

Background

JUNO and IZUMO1 are the first receptor-ligand protein pairs discovered to be essential for sperm-oocyte fusion; their interaction is indispensable for fertilization.

Methods

PCR was used to clone the full-length DNA sequence of the Juno gene in sheep. The single nucleotide polymorphism (SNP) loci of Juno were genotyped by Sequenom MassARRAY®. PCR combined with rapid amplification of cDNA Ends were used to clone the full-length cDNA sequence of Juno and Izumo1. Reverse transcriptase-PCR (RT-PCR) and real time-quantitative-PCR (RT-qPCR) were used to analyze the genes’ expression in tissues of sheep, and single cell RNA-seq was used to analyze the genes’ expression in oocytes, granulosa cells and follicular theca of polytocous and monotocous Small Tail Han ewes. Bioinformatics was used to analyze advanced structure and phylogeny of JUNO and IZUMO1 proteins.

Results

The full-length DNA sequence of the Juno gene in sheep was cloned and nine SNPs were screened. We found a significant association between the g.848253 C > A locus of Juno and litter size of Small Tail Han sheep (P < 0.05). The full-length cDNA sequence of Juno and Izumo1 genes from Small Tail Han sheep were obtained. We found a new segment of the Izumo1 CDS consisting of 35 bp, and we confirmed the Izumo1 gene has 9 exons, not 8. RT-qPCR showed that Juno and Izumo1 genes were highly expressed in ovarian and testicular tissues, respectively (P < 0.01). Single cell RNA-seq showed Juno was specifically expressed in oocytes, but not in granulosa cells or follicular theca, while Izumo1 displayed little to no expression in all three cell types. There was no difference in expression of the Juno gene in oocyte and ovarian tissue in sheep with different litter sizes, indicating expression of Juno is not related to litter size traits. Bioinformatic analysis revealed the g.848253 C > A locus of Juno results in a nonconservative missense point mutation leading to a change from Phe to Leu at position 219 in the amino acid sequence.

Conclusions

For the first time, this study systematically analyzed the expression, structure and function of Juno and Izumo1 genes and their encoded proteins in Small Tail Han sheep, providing the basis for future studies of the regulatory mechanisms of Juno and Izumo1 genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00548-4.

Polymorphism Detection of GDF9 Gene and Its Association with Litter Size in Luzhong Mutton Sheep (Ovis aries)

Simple Summary

GDF9 and BMPR1B are two important reproduction genes. In this study, the whole coding region of GDF9 was sequenced, of which the mutations were detected in Luzhong mutton sheep. The results suggested that two single nucleotide polymorphisms (SNPs), g.41768501A > G and g.41768485 G > A in GDF9 gene were associated with litter size. The g.41768485 G > A is a missense mutation which is predicted to affect the tertiary structure of the protein. Thus, these two mutations may be potential effective genetic markers to improve the litter size in sheep.

Abstract

Litter size is one of the most important economic traits in sheep. GDF9 and BMPR1B are major genes affecting the litter size of sheep. In this study, the whole coding region of GDF9 was sequenced and all the SNPs (single nucleotide polymorphisms) were determined in Luzhong mutton ewes. The FecB mutation was genotyped using the Sequenom MassARRAY^®SNP assay technology. Then, the association analyses between polymorphic loci of GDF9 gene, FecB, and litter size were performed using a general linear model procedure. The results showed that eight SNPs were detected in GDF9 of Luzhong mutton sheep, including one novel mutation (g.41769606 T > G). The g.41768501A > G, g.41768485 G > A in GDF9 and FecB were significantly associated with litter size in Luzhong mutton ewes. The g.41768485 G > A is a missense mutation in the mature GDF9 protein region and is predicted to affect the tertiary structure of the protein. The results preliminarily demonstrated that GDF9 was a major gene affecting the fecundity of Luzhong mutton sheep and the two loci g.41768501A > G and g.41768485 G > A may be potential genetic markers for improving litter size.

Darwin's Fancy Revised: An Updated Understanding of the Genomic Constitution of Pigeon Breeds

Through its long history of artificial selection, the rock pigeon (Columba livia Gmelin 1789) was forged into a large number of domestic breeds. The incredible amount of phenotypic diversity exhibited in these breeds has long held the fascination of scholars, particularly those interested in biological inheritance and evolution. However, exploiting them as a model system is challenging, as unlike with many other domestic species, few reliable records exist about the origins of, and relationships between, each of the breeds. Therefore, in order to broaden our understanding of the complex evolutionary relationships among pigeon breeds, we generated genome-wide data by performing the genotyping-by-sequencing (GBS) method on close to 200 domestic individuals representing over 60 breeds. We analyzed these GBS data alongside previously published whole-genome sequencing data, and this combined analysis allowed us to conduct the most extensive phylogenetic analysis of the group, including two feral pigeons and one outgroup. We improve previous phylogenies, find considerable population structure across the different breeds, and identify unreported interbreed admixture events. Despite the reduced number of loci relative to whole-genome sequencing, we demonstrate that GBS data provide sufficient analytical power to investigate intertwined evolutionary relationships, such as those that are characteristic of animal domestic breeds. Thus, we argue that future studies should consider sequencing methods akin to the GBS approach as an optimal cost-effective approach for addressing complex phylogenies.

A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing

Background

Copy number variation (CNV) is an important source of genetic variation that has a significant influence on phenotypic diversity, economically important traits and the evolution of livestock species. In this study, the genome-wide CNV distribution characteristics of 32 fine-wool sheep from three breeds were analyzed using resequencing.

Results

A total of 1,747,604 CNVs were detected in this study, and 7228 CNV regions (CNVR) were obtained after merging overlapping CNVs; these regions accounted for 2.17% of the sheep reference genome. The average length of the CNVRs was 4307.17 bp. “Deletion” events took place more frequently than “duplication” or “both” events. The CNVRs obtained overlapped with previously reported sheep CNVRs to variable extents (4.39–55.46%). Functional enrichment analysis showed that the CNVR-harboring genes were mainly involved in sensory perception systems, nutrient metabolism processes, and growth and development processes. Furthermore, 1855 of the CNVRs were associated with 166 quantitative trait loci (QTL), including milk QTLs, carcass QTLs, and health-related QTLs, among others. In addition, the 32 fine-wool sheep were divided into horned and polled groups to analyze for the selective sweep of CNVRs, and it was found that the relaxin family peptide receptor 2 (RXFP2) gene was strongly influenced by selection.

Conclusions

In summary, we constructed a genomic CNV map for Chinese indigenous fine-wool sheep using resequencing, thereby providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07387-7.

Whole Genome Sequencing Reveals the Effects of Recent Artificial Selection on Litter Size of Bamei Mutton Sheep

Simple Summary

Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. There is large variation in the reproductive abilities of Bamei mutton sheep. After recent artificial selection, the average lambing rate of the Bamei mutton nucleus group was over 150%. We used the FST (Fixation Index) and XP-EHH (The Cross-Population Extended Haplotype Homozygosity) statistical approach to detect the selective sweeps between high- and low-fecundity Bamei mutton sheep groups. JUN (JUN proto-oncogene, AP-1 transcription factor subunit), ITPR3 (inositol 1,4,5-trisphosphate receptor type 3, PLCB2 (phospholipase C beta 2), HERC5 (HECT and RLD domain containing E3 ubiquitin protein ligase 5), and KDM4B (lysine demethylase 4B) were detected that are potential responsible for litter size. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection.

Abstract

Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on detecting candidate genes associated with the increasing of the litter size in this breeds under recent artificial selection to improve the efficiency of mutton production. We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH (gonadotropin-releasing hormone), oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 (forkhead box A1) and ESR1 (estrogen receptor 1). Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions that are responsible for litter size can in turn be used for further selection.

Genetic signatures of high-altitude adaptation and geographic distribution in Tibetan sheep

Most sheep breeding programs designed for the tropics and sub-tropics have to take into account the impacts of environmental adaptive traits. However, the genetic mechanism regulating the multiple biological processes driving adaptive responses remains unclear. In this study, we applied a selective sweep analysis by combing 1% top values of F_st and ZHp on both altitude and geographic subpopulations (APS) in 636 indigenous Tibetan sheep breeds. Results show that 37 genes were identified within overlapped genomic regions regarding F_st significantly associated with APS. Out of the 37 genes, we found that 8, 3 and 6 genes at chromosomes (chr.) 13, 23 and 27, respectively, were identified in the genomic regions with 1% top values of ZHp. We further analyzed the INDEL variation of 6 genes at chr.27 (X chromosome) in APS together with corresponding orthologs of 6 genes in Capra, Pantholops, and Bos Taurus. We found that an INDEL was located within 5′UTR region of HAG1 gene. This INDEL of HAG1 was strongly associated with the variation of APS, which was further confirmed by qPCR. Sheep breeds carrying “C-INDEL” of HAG1 have significantly greater body weight, shear amount, corpuscular hemoglobin and globulin levels, but lower body height, than those carrying “CA-INDEL” of HAG1. We concluded that “C-INDEL” variation of HAG1 gene confers better hypoxia tolerance in the highlands of Tibetan and explains well geographic distributions in this population. These results contribute to our understanding of adaptive responses to altitude and geographic adaptation in Tibetan sheep populations and will help to guide future conservation programs for Tibetan sheep native to Qinghai-Tibetan Plateau.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

Population structure, genetic diversity, and selective signature of Chaka sheep revealed by whole genome sequencing

Background

Chaka sheep, named after Chaka Salt Lake, are adapted to a harsh, highly saline environment. They are known for their high-grade meat quality and are a valuable genetic resource in China. Furthermore, the Chaka sheep breed has been designated a geographical symbol of agricultural products by the Chinese Ministry of Agriculture.

Results

The genomes of 10 Chaka sheep were sequenced using next-generation sequencing, and compared to that of additional Chinese sheep breeds (Mongolian: Bayinbuluke and Tan; Tibetan: Oula sheep) to explore its population structure, genetic diversity and positive selection signatures. Principle component analysis and a neighbor-joining tree indicated that Chaka sheep significantly diverged from Bayinbuluke, Tan, and Oula sheep. Moreover, they were found to have descended from unique ancestors (K = 2 and K = 3) according to the structure analysis. The Chaka sheep genome demonstrated comparable genetic diversity from the other three breeds, as indicated by observed heterozygosity (Ho), expected heterozygosity (He), runs of homozygosity (ROH), linkage disequilibrium (LD) decay. The enrichment analysis revealed that in contrast to Mongolian or Tibetan lineage groups, the genes annotated by specific missense mutations of Chaka sheep were enriched with muscle structure development (GO:0061061) factors including insulin-like growth factor 1 (IGF1), growth differentiation factor 3 (GDF3), histone deacetylase 9 (HDAC9), transforming growth factor beta receptor 2 (TGFBR2), and calpain 3 (CAPN3), among others. A genome-wide scan using Fst and XP-CLR revealed a list of muscle-related genes, including neurofibromin 1 (NF1) and myomesin 1 (MYOM1), under potential selection in Chaka sheep compared with other breeds.

Conclusions

The comprehensive genome-wide characterization provided the fundamental footprints for breeding and management of the Chaka sheep and confirmed that they harbor unique genetic resources.

Genetic diversity and genetic origin of Lanping black-boned sheep investigated by genome-wide single-nucleotide polymorphisms (SNPs)

Lanping black-boned sheep was first discovered in the 1950s in Lanping county of China and characterized by black pigmentation on skin and internal organs. Due to the novel and unique trait, the genetic background of Lanping black-boned sheep is of great interest. Here, we genotyped genome-wide SNPs (single nucleotide polymorphisms) of Lanping black-boned sheep and Lanping normal sheep using Illumina OvineSNP50 BeadChip to investigate the genetic diversity and genetic origin of Lanping black-boned sheep. We also downloaded a subset SNP dataset of two Tibet-lineage sheep breeds and four other sheep breeds from the International Sheep Genomics Consortium (ISGC) as a reference for interpreting. Lanping black-boned sheep had a lower genetic diversity level when compared to seven other sheep breeds. Principal component analysis (PCA) showed that Lanping black-boned sheep and Lanping normal sheep were clustered into the Asian group, but there was no clear separation between the two breeds. Structure analysis demonstrated a high ancestry coefficient in Lanping black-boned sheep and Lanping normal sheep. However, the two populations were separated into two distinct branches in a neighbor-joining (NJ) tree. We further evaluated the genetic divergence using population F ST , which showed that the genetic differentiation that existed between Lanping black-boned sheep and Lanping normal sheep was higher than that between Tibet sheep and Changthangi sheep, which revealed that Lanping black-boned sheep is a different breed from Lanping normal sheep on the genetic level. In addition, structure analysis and NJ tree showed that Lanping black-boned sheep had a relatively close relation with Tibet sheep. The results reported herein are a first step toward understanding the genetic background of Lanping black-boned sheep, and it will provide informative knowledge on the unique genetic resource conservation and mechanism of novel breed formation.

Integrating Genome-Wide CNVs Into QTLs and High Confidence GWAScore Regions Identified Positional Candidates for Sheep Economic Traits

Copy number variations (CNVs) are important source of genetic variation, which can affect diverse economic traits through a variety of mechanisms. In addition, genome scan can identify many quantitative trait loci (QTLs) for the economic traits, while genome-wide association studies (GWAS) can localize genetic variants associated with the phenotypic variations. Here, we developed a method called GWAScore which collected GWAS summary data to identify potential candidates, and integrated CNVs into QTLs and high confidence GWAScore regions to detect crucial CNV markers for sheep growth traits. We got 197 candidate genes which were overlapping with the candidate CNVs. Some crucial genes (MYLK3, TTC29, HERC6, ABCG2, RUNX1, etc.) showed significantly elevated GWAScore peaks than other candidate genes. In this study, we developed the GWAScore method to excavate the potential value of candidate genes as markers for the sheep molecular breeding.

The effect of SNP rs400827589 in exon 2 of the MTNR1B gene on reproductive seasonality and litter size in sheep

In mammals, the melatonin receptor gene has been widely studied since it has a great influence on reproductive traits. However, little is known about the association between polymorphism of the coding region of the MTNR1B gene and year‐round oestrus or the litter size in Small Tail Han sheep. To better understand the effects of single nucleotide polymorphism (SNP) rs400827589 in MTNR1B, a population polymorphism analysis was conducted using genotyping data in 45 sheep breeds around the world. The results indicated that TT was the dominant genotype in all sheep breeds. The associations of this SNP with reproductive seasonality and litter size in Small Tail Han sheep showed rs400827589 was correlated with fecundity as assessed by reproductive seasonality and litter size (p < .05). Bioinformatics analysis indicated the change in amino acid from Ile to Leu may affect the function of the MTNR1B protein by impacting the secondary and tertiary protein structures. The present results demonstrate that rs400827589 could be used in the marker‐assisted selection of the litter size in Small Tail Han sheep.

Genomic Diversity, Population Structure, and Signature of Selection in Five Chinese Native Sheep Breeds Adapted to Extreme Environments

Through long term natural and artificial selection, domestic sheep (Ovis aries) have become adapted to a diverse range of agro-ecological environments and display multiple phenotypic traits. Characterization of diversity and selection signature is essential for genetic improvement, understanding of environmental adaptation, as well as utilization and conservation of sheep genetic resources. Here, we aimed to assess genomic diversity, population structure, and genomic selection among five Chinese native sheep breeds using 600K high density SNP genotypes. A total of 96 animals of the five breeds were selected from different geographical locations with extremely dry or humid conditions. We found a high proportion of informative SNPs, ranging from 93.3% in Yabuyi to 95.5% in Wadi, Hu, and Hetian sheep. The average pairwise population differentiation (F_ST) between the breeds was 0.048%, ranging from 0.022% to 0.054%, indicating their low to moderate differentiation. PCA, ADMIXTURE, and phylogenetic tree analyses revealed a clustering pattern of the five Chinese sheep breeds according to their geographical distribution, tail type, coat color, body size, and breeding history. The genomic regions under putative selection identified by F_ST and XP-EHH approaches frequently overlapped across the breeds, and spanned genes associated with adaptation to extremely dry or humid environments, innate and adaptive immune responses, and growth, wool, milk, and reproduction traits. The present study offers novel insight into genomic adaptation to dry and humid climates in sheep among other domestic animals and provides a valuable resource for further investigation. Moreover, it contributes useful information to sustainable utilization and conservation of sheep genetic resources.

Genomic Scan for Selection Signature Reveals Fat Deposition in Chinese Indigenous Sheep with Extreme Tail Types

Simple Summary

According to the tail types, sheep can be briefly classified into three groups (fat-tailed, fat-rumped, and thin-tailed sheep). In this study, we used these three typical breeds from Chinese indigenous sheep breeds to perform a genome scan for selective sweeps using Ovine Infinium HD SNP BeadChip genotype data. Results showed that 25 genomic regions exhibited selection signals and harbored 73 positional candidate genes. These genes were documented not only to be associated with tail fat formation, but also be related to reproduction, body conformation, and appearance. Our findings contributed to understanding genetic basis of fat deposition in sheep tail and provide a reference for developing new sheep breeds with an ideal tail type.

Abstract

It is a unique feature that fat can be deposited in sheep tails and rumps. To elucidate the genetic mechanism underlying this trait, we collected 120 individuals from three Chinese indigenous sheep breeds with extreme tail types, namely large fat-tailed sheep (n = 40), Altay sheep (n = 40), and Tibetan sheep (n = 40), and genotyped them using the Ovine Infinium HD SNP BeadChip. Then genomic scan for selection signatures was performed using the hapFLK. In total, we identified 25 genomic regions exhibiting evidence of having been under selection. Bioinformatic analysis of the genomic regions showed that selection signatures related to multiple candidate genes had a demonstrated role in phenotypic variation. Nine genes have documented association with sheep tail types, including WDR92, TBX12, WARS2, BMP2, VEGFA, PDGFD, HOXA10, ALX4, and ETAA1. Moreover, a number of genes were of particular interest, including RXFP2 associated with the presence/absence and morphology of horns; MITF involved in coat color; LIN52 and SYNDIG1L related to the number of teats; MSRB3 gene associated with ear sizes; LTBP2 considered as a positional candidate genes for number of ribs; JAZF1 regulating lipid metabolism; PGRMC2, SPAG17, TSHR, GTF2A1, and LARP1B implicated with reproductive traits. Our findings provide insights into fat tail formation and a reference for carrying out molecular breeding and conservation in sheep.