Next Article in Journal
Analysis of Runs of Homozygosity in Aberdeen Angus Cattle
Next Article in Special Issue
Characterisation of Ovine KRTAP19-3 and Its Impact on Wool Traits in Chinese Tan Sheep
Previous Article in Journal
The Interaction between Canine Semen Bacteria and Semen Quality Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep

1
Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810001, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Qinghai Sheep Breeding and Promotion Service Center, Gangcha 812300, China
*
Author to whom correspondence should be addressed.
Animals 2024, 14(15), 2152; https://doi.org/10.3390/ani14152152
Submission received: 21 June 2024 / Revised: 15 July 2024 / Accepted: 22 July 2024 / Published: 24 July 2024
(This article belongs to the Special Issue Genetics and Breeding in Ruminants)

Abstract

:

Simple Summary

The horned type of sheep, including the polled phenotype and scurs phenotype, emerged during the process of artificial domestication. Sheep serve as an optimal animal model for investigating the correlation between quality and quantity traits of horns, as well as the multi-gene regulatory mechanism underlying these traits. While significant progress has been made in understanding the effects of natural selection and sex selection on horn phenotypes, there is limited research on polymorphisms resulting from artificial selection. Through genome-wide association analysis, this study reveals the genetic mapping of horn phenotypes, polled genes, and superior haplotypes to provide valuable insights into studying the genetic mechanisms governing horn traits and other specific characteristics in ruminants while also laying a foundation for breeding.

Abstract

The types and morphology of sheep horns have been extensively researched, yet the genetic foundation underlying the emergence of diverse horn characteristics during the breeding of polled Tibetan sheep has remained elusive. Genome-wide association analysis (GWAS) was performed on 103 subtypes (normal large horn, scurs, and polled) differentiated from G2 (offspring (G2) of parent (G1) of polled) of the polled core herd. Six single nucleotide polymorphisms (SNPs) located on chromosome 10 of the relaxin family peptide receptor 2 (RXFP2) gene exhibited positive correlations with horn length, horn base circumference, and horn base interval. Furthermore, in genotyping 382 G2 individuals, significant variations were observed for each specific horn type. Three additional mutations were identified near the target SNP upstream of the amplification product. Finally, the RXFP2-specific haplotype associated with the horned trait effectively maintained horn length, horn base circumference, and horn base interval in Tibetan sheep, as confirmed by population validation of nine loci in a sample size of 1125 individuals. The present study offers novel insights into the genetic differentiation of the horned type during improvement breeding and evolution, thereby establishing a robust theoretical foundation for polled Tibetan sheep breeding and providing valuable guidance for practical production.

Graphical Abstract

1. Introduction

In cases where strong directional selection is evident, genetic traits often fail to exhibit the expected evolutionary response [1]. Therefore, elucidating the genes and genomic regions that contribute to trait variation provides an invaluable opportunity for comprehending the intricacies of genetic diversity processes within populations. The presence of adaptive variation at the genotype level may not be readily apparent when focusing on phenotypes for analysis, potentially due to genetic correlations between the phenotype of interest and other fitness-related traits [2]. Genetic resources obtained through directional selection [3] were utilized in this study to identify the genetic polymorphisms associated with horn traits within the population. Additionally, the effects of genetic variation on specific traits were determined, thereby establishing a crucial foundation for understanding the relationship between genotype and phenotype traits following artificial intervention.
The cranial appendages of ruminant animals, commonly referred to as headgear, consist of a pneumatized bony core fused with the frontal bone of the skull and are covered by a continuously growing keratin sheath [4]. Horn injuries to cattle caused by too-frequent and severe bumps are listed as one of the top 10 challenges to the U.S. beef supply [5]. To reduce injury risk to handlers and other animals as well as economic losses, dehorning procedures are performed [6]. Therefore, these appendages have become unpopular in modern captive breeding programs. The market-oriented sheep industry values its economic value because of animal welfare concerns [6,7,8]. More and more artificial interventions have been employed in recent years to genetically breed sheep without horns; the occurrence of this has been documented in historical records dating back to ancient Egypt [9]. However, there are still a few native breeds that naturally possess this trait. Prolonged human-mediated artificial selection has resulted in phenotypic variation within the genetic structure, which can be identified through patterns of genomic variation such as reduced nucleotide diversity and regions with very low heterozygosity. Importantly, the detection of numerous single nucleotide polymorphisms (SNPs) provides crucial information for identifying potential sites responsible for phenotypic variation. Therefore, the comprehensive identification and strategic deployment of specific SNPs via whole-genome sequencing (WGS) analysis promises substantial advancements in the domains of animal husbandry and molecular biology, particularly for the enhancement of production traits. A growing number of researchers are conducting extensive studies of various patterns of the presence and absence of horns, multi-horns, and deformed horns from around the world. With human domestication, domestic sheep breeds have also developed a hornless phenotype; breeders have been breeding polled sheep since the 16th century [10]. Earlier studies proposed that the autosomal 10 locus HoP, Ho+, and HohL alleles inheritance model controls the phenotype of horns, with the HoP allele determining polledness [10,11,12]. The horn allele has been mapped to OAR 10 such that there are at least two loci affecting the presence or absence of horns in Merino × Romney crosses [13]. The horn locus has been fine mapped to a 50 kb region in the sheep genome and predictive haplotypes for the polled HP allele have been identified using a segregating (Merino × Romney) × Merino resource [13]. RXFP2 is the principle candidate locus which supports a new model of horn-type inheritance in Soay sheep [14]. The analysis of the sequence revealed a 1833-base-pair genomic insertion located in the 3′-UTR region of the RXFP2 gene in polled animals [15,16]. The main candidate gene associated with the type and size of sheep horns, RXFP2, was found to have several new mutation sites. Among these mutations, the synonymous mutation p.P375 (c.1125A>G) was identified as a potential indicator for determining the presence or absence of horns in Tan sheep [17]. So far, the impact of two causal mutations on the process of horn formation remains unknown and, despite their influence on polled traits [18], no genetic variation has been identified that is causative of horn phenotype. Although the genetic loci associated with polledness have been mapped successfully via WGS, the observed variations in hornless sheep populations across different breeds are notable. Additionally, horn type in certain breeds is associated with variations in the RXFP2 gene region, indicating a highly complex genetic mechanism underlying the polled trait. These studies have only explained a small portion of the overall genetic variation; therefore, further experimentation is necessary to identify any causal sites related to the polled trait in ruminants.
Tibetan sheep are one of the prominent animal breeding resources unique to the Qinghai–Tibet Plateau. In order to elucidate the genetic characteristics of horn development, we established polled male and female groups (G1) based on phenotype in the adult phase and performed targeted mating, with the resulting offspring (G2) serving as our research subjects. Until now, there have been no reports on the genetic basis underlying horned differentiation enhanced breeding.
In this study, we conducted GWAS to analyze three horn phenotypes in G2 and subsequently identified specific haplotypes through population verification. By investigating candidate genes associated with polled traits at the genomic level, we were able to assess the molecular genetic characteristics of polled traits and provide fundamental data for localizing functional genes related to horn traits. This will serve as a crucial reference for the future molecular breeding of polled sheep.

2. Materials and Methods

2.1. Samples Collection and Whole-Genome Sequencing

The core group exhibited phenotypic variations in the second generation (G2), where polled rams and ewes (G1) were mated resulting in individuals with either normal large horns, scurs, or no horns. After continuously observing the horn status of a core population with a well-defined pedigree, we meticulously selected a total of 103 sheep for sequencing analysis, including 48 polled sheep, 29 sheep with normal large horns, and 26 sheep with scurs. The genotyping of 382 G2 individuals was conducted within a non-horned core group. The phenotypic values of the compiled horns were collected, and their whole blood was obtained and stored in EDTA vacuum tubes. The husbandry and management of these animals was consistently maintained at the same level. The QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) was utilized to extract total genomic DNA from 200 µL of sheep blood, followed by assessment of DNA quality. Illumina’s paired-end DNA Sample Prep kit standard library building procedure was employed to prepare sequencing libraries using a minimum of 3 µg genomic DNA. The Agilent 2100 bioanalyzer (Agilent, Santa Clara, CA, USA) was used for examination and real-time PCR quantification of the sequencing libraries prior to sequencing on the Illumina NovaSeq6000 NGS platform (Guangzhou, China).

2.2. Quality Control and Annotation

The Illumina data underwent filtration using FASTP (version 0.18.0) with the following criteria: (1) reads containing ≥10% unknown nucleotides (N) were removed; (2) reads with ≥50% of bases having Phred quality scores ≤ 20 were removed; and (3) reads with adapters were deleted. The resulting filtered clean reads were utilized for assembly analysis. The filtered reads were aligned to the reference genome using the Burrows–Wheeler Aligner (BWA) (version 0.7.12) alignment software (ncbi_GCF_002742125.1) with the settings ‘mem 4 - k 32 - M’, - k as the minimum seed length, and -M as an option used to mark shorter split alignment hits as secondary alignments [19]. Variant calling was performed for all samples using GATK’s Unified Genotyper (version 3.4-46) [20]. SNPs and InDels were filtered using GATK’s Variant Filtration with proper standards (-Window 4, -filter QD < 2.0 || FS > 60.0 || MQ < 40.0, -G_filter “GQ < 20”), and those exhibiting segregation distortion or sequencing errors were discarded. To determine the physical positions of each SNP, the variants were annotated using the ANNOVAR software (version 2) tool [21] to align and annotate SNPs.

2.3. Genome-Wide Association Analysis

The implementation of REGENIE [22] and SAIGE [23] was chosen due to the presence of rare variants and highly imbalanced case–control ratios. In an effort to enhance the power of the association test, we employed a gene-region-based approach (or other methods for generating variant sets) and utilized GMMAT (version 0.9.3) [24] software for conducting gene/set-based association analysis. We included three tests in the framework: the burden test, SKAT (sequence kernel association test) [25], and SKAT-O (optimal sequence kernel association test) [26]. The SKAT method employed a multiple regression model to directly assess the relationship between the phenotypes and genetic variants in a variant set, as well as covariates. This approach allows for variations in both direction and magnitude of effects among different variants. In the burden test, SKATO estimated the weight and SKAT statistics that maximize power using the formula TSKAT−0 = (1 − ρ)TSKAT + ρTburden. Additionally, SKAT-O required a grid search over ρ to determine the minimum p-value. The Bonferroni correction threshold (p-value = 0.01/marker number or 0.05/marker number) was employed to identify statistically significant associations. Candidate genes (CAGs) located within a 50 kb region upstream or downstream of the significantly associated markers were identified.

2.4. Enrichment Analysis of Candidate Genes

The candidate genes were mapped to GO terms in the Gene Ontology database [27] (http://www.geneontology.org/, accessed on 20 April 2024). Gene numbers were calculated for each term, and significantly enriched GO terms were determined using a hypergeometric test. The p-value calculation formula is as follows:
P = 1 i = 0 m 1 M i N M n i N n
The number of all genes with GO annotation is represented by N, while n represents the number of genes/CAGs in N. Similarly, M denotes the number of all genes annotated to specific pathways, and m represents the number of genes/CAGs in M. The calculated p-value underwent false discovery rate (FDR) correction, with a threshold set at FDR ≤ 0.05. Pathways that meet this condition are defined as significantly enriched GO terms in genes/CAGs.

2.5. Group Validation of Significantly Correlated SNPs

The blood samples of 1125 sheep from a large group were collected in a random manner, and the horn phenotype data (Figure 1, Supplementary Table S11), including measurements of horn length, horn base circumference, and horn base interval, were recorded using the following methods:
Horn length: the arc length along the longitudinal edge of the back of the horn, from its base to its end;
Horn base circumference: the circumference of the longitudinal base encircling the horn base;
Horn base interval: the linear distance between the vertical bases on the posterior side of the horn base.
The nine SNPs detected in exons of genes identified by screening analysis were validated through Sanger sequencing and SNaPshot mini-sequencing, and primers were designed using Primer3 v0.4.0 (1) [28] while single-base extension primers were designed using the online tool Primer 3 plus (Supplemental Tables S1 and S2).
The amplification and screening of SNPs were conducted in 1125 sheep to analyze the correlation between candidate mutations and horn traits. Initially, the amplification of exons from genomic DNA was achieved using primers. PCR reactions were performed in 10 μL volumes containing 5 μL Taq PCR Mix, 1 μL Primer Mix, 1 μL DNA, and 3 μL ddH2O. A PCR program consisting of an initial denaturation at 94 °C for 5 min, followed by 35 cycles at 94 °C for 20 s, 60 °C for 30 s, and finally at 72 °C for 30 s was carried out using a Veriti™96-Well Thermal Cycler (Applied Bio-systems, Foster City, CA, USA). All PCR products were purified with ExoSAPIT® (USB Corporation, Cleveland, OH, USA) at 37 °C for 40 min followed by incubation at 85 °C for 15 min. Subsequently, all SNaPshot® reactions were performed in a volume of 5 μL including 0.5 μL SNaPshotMix, 3 μL Pooled PCR Products, 1 μL Pooled Primers, and 0.5 μL ddH2O with a program consisting of an initial denaturation at 95 °C for 30 s and 35 cycles comprising 95 °C for 5 s, 52 °C for 5 s, and 60 °C for 5 s. All SNaPshot®-specific products were purified with Shrimp Alkaline Phosphatase-SAP (Thermo Fisher Scientific Inc., Waltham, MA USA) and incubated at 95 °C for 5 min. Finally, the Genemarker Software (v. 1.90) was used to analyze all the SNaPshot®-specific results. Allele frequency and genotype frequency were assessed using web-based software to determine their adherence to the Hardy–Weinberg equilibrium (http://scienceprimer.com/hardy-weinberg-equilibrium-calculator, accessed on 20 April 2024).
The association between various genotypes or haplotypes and horn traits was assessed using a sampling approach based on general linear models:
YijTl = ξ + αi + βj + γT + δl + θ
Here:
YijTl is the horned trait measured observation;
ξ is the overall mean;
αi is the genetic effect;
βj is the individual age effect;
γT is the year effect;
δl is the permanent environment effect;
θ is the random residual effect and is assumed to be independent, N (0, σ2) distribution.
The data were analyzed using the General Linear Model (GLM) in SPSS Statistics (version 29, IBM, Armonk, NY, USA), and statistical significance was determined at a p-value of less than 0.05.

3. Results

3.1. Summary Statistics of Phenotype Data and Quality Control

We established a polled core population and observed polymorphisms in the horn phenotypes of G2 generation sheep through directional selection, including degenerated, deformed horns known as scurs, polled individuals, and individuals with normal large horns. The descriptive statistics of phenotypic values associated with horn-type traits in a population of 103 Tibetan sheep are provided in Supplementary Table S3.
After implementing quality control and management measures for 36,023,978 SNP loci, the following conditions were excluded: (1) 14,513,444 SNPs with a secondary gene frequency less than 0.05; (2) 200,718 SNPs with a deletion rate exceeding 0.5; and (3) 30,361 SNPs with a heterozygosity ratio greater than 0.8. Subsequently, a total of 21,279,455 independent SNPs remained for the subsequent analysis. The final count revealed a total of 33,065,497 single nucleotide polymorphisms (SNPs) and 2,958,481 insertions/deletions (indels). The majority of the high-quality single nucleotide polymorphisms (SNPs) (61.61%) were observed in intergenic regions, characterized by T/C and A/G substitutions, while only 0.70% were located within exonic regions. The remaining SNPs were situated upstream (0.61%) and downstream (0.63%) of the open reading frame, primarily found within intronic regions (33.97%) (Supplementary Table S4).

3.2. Population Structures and Genome-Wide Association Analysis

According to the population structure, admixture v1.3 was employed with a range of K values from 2 to 5 (Figure 2A). Specifically, a K value of 2 was utilized and the corresponding outcome is depicted in Figure 2B. Kinship estimation and principal component analysis (PCA) were conducted on all individuals, confirming the validity of the sampling (Figure 2C,D). The bivariables in the Kinship plot are perfectly aligned along a straight line, indicating a strong linear relationship. Moreover, the sample exhibits high repeatability, suggesting excellent data reproducibility.
According to the number of independent effective SNPs, the genome-wide significance and suggestive significance values were determined as a p-value of 4.06 × 10−10 corresponding to a 1% significance level and a p-value of 2.34 × 10−9 corresponding to a 5% significance level, respectively. SNPs with p-values lower than these thresholds were identified as being associated with Tibetan sheep horn traits. The p-value Manhattan map shows significant correlations between SNPs and horn length (Figure 3A), horn base circumference (Figure 3C), and horn base interval (Figure 3E), with statistically significant regions covering SNPs sites (Figure 3B,D,F). These aforementioned SNPs exhibit a strong positive correlation with the three horn-related traits (Table 1, Supplementary Tables S5–S7), surpassing the predetermined threshold. Further verification will be conducted subsequently. Concurrently, we conducted a comparison of the reported sites in the literature pertinent to RXFP2 and observed that only a single site (OAR10_29461968) aligned with our findings (Table 2).
The genetic polymorphism of horn type was observed in both sexes, with three distinct morphological characteristics: large normal horns, scurs, and polled (polled accounted for 31.12% and 54.55% of all recorded cases in male and female sheep, respectively) (Table 3). Each horn type exhibited a significant number of quantitative variation sites. Differences in the dominance and expression of specific alleles were evident between the sexes, posing challenges to inferring an individual’s genotype.

3.3. Enrichment Analysis

To further elucidate candidate genes associated with significant SNPs, we conducted additional gene enrichment analysis. The pathways enriched in the RXFP2 gene were further subjected to additional recordings. The GO/KEGG pathway analysis revealed the presence of certain pathways associated with horn characteristics. Specifically, the relaxin signaling pathway and neuroactive ligand–receptor interaction were found to be significantly KEGG co-enriched in terms of horn length, horn base circumference, and horn base interval. Furthermore, these three traits collectively exhibited enrichment in 33 GO term pathways. The candidate genes associated with horn length (Figure 4A, Supplementary Table S8), horn base circumference (Figure 4B, Supplementary Table S9), and horn base interval (Figure 4C, Supplementary Table S10) exhibit co-enrichment during gonad development, oocyte maturation, sex differentiation, protein–hormone receptor activity, cell maturation and differentiation processes, as well as positive regulation of CAMP-mediated signaling pathways. These findings suggest that these three traits are likely regulated by shared biological pathways. The enrichment analysis revealed that the RXFP2 gene exhibited the highest level of enrichment among the Gene Ontology terms associated with biological processes.

3.4. Population Verification of SNPs Significantly Correlated with Tibetan Sheep Horn Traits

The SNPs identified by the GWAS mentioned above were subjected to significant testing between horned and polled Tibetan sheep, revealing a significant correlation between the RXFP2 gene and polled Tibetan sheep, which has been extensively reported in the literature. Furthermore, the descriptive statistical information of phenotypic values obtained from a large population through the random selection of 1125 ewes (Supplementary Table S11) provided population verification for these polymorphisms at the aforementioned RXFP2 sites. Genotyping, population genetic analysis, and association determination between SNPs and horn traits were conducted for all loci.
GWAS predicted six individually detected SNPs, while three additional mutations were identified in close proximity to the target SNPs upstream of amplification products (Supplementary Table S12). The value of PIC at SNP001 was <0.25, which was indicative of low polymorphism. The other SNPs we identified exhibited moderate levels of polymorphism (0.25 < PIC < 0.5) (Supplementary Table S13).
Nine SNP verification results revealed significant correlations between three mutation sites and the horn length, horn base circumference, and horn base interval of Tibetan sheep (Table 4). The genotype of horns had a significant effect on horn length, horn base circumference, and horn base interval, while there might be a polygenic influence associated with the genotypic horn.

3.5. Linkage Disequilibrium and Haplotype Block Association Analysis

The analysis results of linkage disequilibrium (LD) indicate a strong correlation between the SNP sites of the two haplotype modules (D’ > 0.85), with a coefficient of 0.81 observed for the correlation between the two haplotype modules (Figure 5A). Subsequently, a correlation analysis was conducted to examine the haplotypes of the two modules, revealing significant effects of different haplotypes on the horn type of Tibetan sheep (Table 5 and Table 6). In the 19 kb LD block, twelve haplotypes containing three SNPs were identified, while in the 51 kb LD block, twenty-two haplotypes containing six SNPs were found (Figure 5B,C). In Module 1, three haplotypes, namely AAGGCT (n = 11, p-value < 0.01), GAGACT (n = 9, p-value < 0.01), and GAGGTT (n = 8, p-value < 0.01) (Table 5), exhibited significant correlations with polled traits. In Module 2, the haplotypes AACCAAGGCTCT (n = 2, p-value < 0.01), AATTGAAATTTT (n = 18, p-value < 0.01), AATTGGAATTTT (n = 3, p-value < 0.01) GACTGAGATTCT (n = 4, p-value < 0.01), and GATTGGAATTTT (n = 15, p-value < 0.01) (Table 6) were found to be significantly associated with polled traits. The maintenance mechanism of horn characteristic variation is determined by the co-breeding of haplotypes associated with these three traits.

4. Discussion

The vestigialization of horns occurs in most domestic lines due to the diminishing relevance of traits that ensure fitness in natural environments under artificial breeding conditions [30]. Domestic animals have undergone domestication and improvement to produce similar phenotypic transformations [34], a phenomenon known as “domestication syndrome” [35,36,37]. However, our findings indicate that artificial intervention has led to the differentiation of three distinct horn phenotypes (normal large horns, scur, and polled) in polled-orientated progeny. Furthermore, we have determined that all types of horns, including horn length, horn base circumference, and horn base interval, are regulated by the RXFP2 gene located on chromosome 10 [33,38]. We can distinguish phenotypes that have the same horn but have different genotypes, suggesting that the genotype of the horn may be responsible not only for discrete variation in the diagonal type, but also for many quantitative genetic variations in the diagonal size [39], and that the quantified variation differences in the size of the horn appear to be caused by the genotype of the horn. This study establishes a foundation for investigating genotypic-level selection of horn types and morphology and enhances our understanding of the mechanisms underlying the maintenance and differentiation of horn polymorphisms in this population.
Interestingly, horn polymorphisms also occur in natural populations, such as Soay sheep (males exhibit both normal and scurred horns, while females can have normal, scurred, or polled horns) [39]. The same phenomenon is observed in the natural population of Tibetan sheep. It remains unclear whether this phenotype is associated with genetic introgressions or a free combination of RXFP2 superior alleles linked to horn types within the population. This hypothesis can be tested in future studies by generating genotypic data for additional sheep breeds.
The development of horns in sheep is influenced not only by genetic and sexual factors, but also by hormonal and environmental influences. Seasonal horn growth typically occurs during spring, when prolactin secretion reaches its peak. There are notable variations in the timing and pace of horn development across different life stages [40]. Sustained low levels of testosterone during the inactive period are essential for promoting male-type horn growth, while high levels of testosterone prior to estrus can actually impede growth—a mechanism that also regulates antler development in deer [41,42]. Sheep exhibited cryptorchidism after the RXFP2 gene was knocked out using CRISPR/Cas9 technology [43]. The growth rate of sheep horns is regulated by thyroid activity, with the peak concentration of thyroid hormone (T3, T4) coinciding with the rapid growth period [44]. Our findings demonstrated a significant correlation between the concentration of RXFP2 and gonadal development and maturation, as well as the regulatory impact of reproductive hormones on sheep horn development. Further investigation into the hormonal regulatory mechanisms will be conducted.
The inheritance of horns in sheep is closely associated with sex, as indicated by Wood’s findings [45]. Specifically, the horned and hornless genes exhibit distinct effects on males and females. In this regard, it can be observed that large horns are dominant in males while recessive in females [45]. A study conducted by Dominik et al. revealed the presence of a single nucleotide polymorphism located in a specific genomic region, which exhibited a prominent maternal imprinting effect in both males and females [12]. QTL analysis of diagonal morphology revealed a significant interaction between horn shape traits and sex, indicating that the same chromosomal region accounted for genetic variation in both horn length and base circumference angle. Furthermore, it was determined that these two traits were subject to genotypic (synergistic) selection [38]. The fixation of the ideal individual’s characteristics can be rapidly achieved by mating a homozygous hornless male with a homozygous hornless female [45], resulting in all offspring being hornless.
Prediction within a single family will incorporate both linkage information and LD information, assuming that non-genetic factors have been accurately accounted for. However, in the case of cross-family prediction, reliance can only be placed on LD information [46]. Haplotypes encompass a greater amount of information (heterozygosity) compared to a single SNP, thereby significantly enhancing accuracy and robustness in gene localization. Conversely, in small sample sizes, a single SNP fails to establish an association between a specific phenotype and gene, whereas haplotypes possess this capability [47]. Horn phenotypes can be accurately predicted with an accuracy of approximately 0.7 when utilizing a SNP for prediction [33]. Consequently, investigations into phenotypic correlation with genetic variation are more likely to yield successful outcomes if the SNPs employed have been demonstrated to exhibit linkage imbalance through methodologies such as haplotype analysis [48,49,50]. The alleles located on OAR10 at position 29,461,968 were found to be genetically linked with the neighboring loci in the gene map [14,30]. LD was observed to be lower between the 1.78 kb insert fragment of the RXFP2 gene and the SNP located at OAR10_29511510.1 (position 29,476,678) [16], whereas a higher LD was observed between the SNP inserted with OAR10_29546872.1) [33]. In our study, the haplotype combinations of AAGGCT, GAGACT, and GAGGTT as well as AACCAAGGCTCT, AATTGAAATTTT, AATTGGAATTTT, GACTGAGATTCT, and GATTGGAATTTT are the predominant ones. The presence of fitness disparities may also occur among sheep exhibiting similar phenotypes but possessing different underlying genotypes due to their association with other alleles at adjacent sites in a gene-linkage imbalance. We have successfully identified the haplotype associated with the horn trait, indicating that the absence of horns in artificially bred sheep is influenced by the interaction of multiple SNPs. This finding also highlights the synergistic genetic basis underlying these three horn traits. It is important to note that single SNP loci alone cannot provide sufficient evidence for establishing a correlation between variations in the RXFP2 gene and the expression of the horn trait.
Polled rams exhibit relatively low breeding values and were predominantly heterozygous, the proportion of polled individuals has greater potential in the G2 population, where ewes outnumber rams. It is evident that there exists a common genetic basis for both horn type and length, with at least 76% of the genetic variation attributed to the “horn” genotype [14]. Consequently, genome selection can significantly accelerate the rate of genetic improvement and positively impact both the breeding value and quantity of polled rams. Quantitative trait genes in animals were directly selected using molecular marker technology, and whole-genome association analysis based on the genome variation map was employed to enhance breed improvement more effectively. A breeding program thus may well produce polled sire with breeding values comparable to the top sire in the near future [51]. Based on these findings, we anticipate that a majority of lambs born in generations G3–G4 will be polled, thus establishing a new prototype for a polled strain.

5. Conclusions

In summary, our GWAS identified six single nucleotide polymorphisms (SNPs) on the RXFP2 gene located on chromosome 10 that exhibited significant associations with horned traits. Each horn type displayed a substantial number of variation sites. Furthermore, three additional mutations were subsequently detected near the target SNP upstream of the amplification product. Population validation of the specific haplotype is considered indicative of distinguishing between horn and polled phenotypes while maintaining horn length, horn base circumference, and horn base interval. This study aims to provide novel insights into the molecular characteristics of polled sheep under human intervention and can offer practical guidance.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani14152152/s1.

Author Contributions

K.Z. designed the experiment, while D.T. conducted data analysis and wrote the paper. X.L. and B.H. (Buying Han) were responsible for DNA sample preparation and data analysis. Z.Z. and B.H. (Bin Huang) collected the samples. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Natural Science Foundation of Qinghai Province (2022-ZJ-901). National Breeding Joint Research Project.

Institutional Review Board Statement

The present study adhered to the guidelines outlined in the “Regulations for the Management of Experimental Animals” (Ministry of Science and Technology, China, revised in June 2004). Ethical approval was obtained from the Animal Care and Use Committees of the Northwest Institute of Plateau Biology, Chinese Academy of Sciences (NWIPB2021311, 11 March 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data utilized in this study are accessible within this article as well as the supplementary information.

Acknowledgments

We greatly appreciate our collaborators for their assistance in sample collection.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Merilä, J.; Sheldon, B.C.; Kruuk, L.E.B. Explaining stasis: Microevolutionary studies in natural populations. Genetica 2001, 112, 199–222. [Google Scholar] [CrossRef] [PubMed]
  2. Lande, R.; Arnold, S.J. The measurement of selection on correlated characters. Evolution 1983, 37, 1210–1226. [Google Scholar] [CrossRef] [PubMed]
  3. Coltman, D.W.; O’Donoghue, P.; Hogg, J.T.; Festa-Bianchet, M. Selection and genetic (co) variance in bighorn sheep. Evolution 2005, 59, 1372–1382. [Google Scholar] [PubMed]
  4. Allais-Bonnet, A.; Grohs, C.; Medugorac, I.; Krebs, S.; Djari, A.; Graf, A.; Fritz, S.; Seichter, D.; Baur, A.; Russ, I.; et al. Novel Insights into the Bovine Polled Phenotype and Horn Ontogenesis in Bovidae. PLoS ONE 2013, 8, e63512. [Google Scholar] [CrossRef] [PubMed]
  5. Koene, P. Breeding and Feeding For Animal Health and Welfare In Organic Livestock Systems-Animal Welfare and Genetics In Organic Farming of Layers: The Example Of Cannibalism. In Proceedings of the fourth NAHWOA Workshop, Wageningen, The Netherlands, 24–27 March 2001; pp. 62–85. [Google Scholar]
  6. Misch, L.J.; Duffield, T.F.; Millman, S.T.; Lissemore, K.D. An investigation into the practices of dairy producers and veterinarians in dehorning dairy calves in Ontario. Can. Vet. J. 2007, 48, 1249–1254. [Google Scholar] [PubMed]
  7. Graf, B.; Senn, M. Behavioural and physiological responses of calves to dehorning by heat cauterization with or without local anaesthesia. Appl. Anim. Behav. Sci. 1999, 62, 153–171. [Google Scholar] [CrossRef]
  8. Prayaga, K.C. Genetic options to replace dehorning in beef cattle—A review. Aust. J. Agric. Res. 2007, 58, 1–8. [Google Scholar] [CrossRef]
  9. Roman, A.V. L’élevage Bovin en Égypte Antique. Ph.D. Thesis, École Nationale Vétérinaire d’Alfort, Maisons-Alfort, France, 2004. [Google Scholar]
  10. Dolling, C. Hornedness and polledness in sheep. I. The inheritance of polledness in the Merino. Aust. J. Agric. Res. 1960, 11, 427–438. [Google Scholar] [CrossRef]
  11. Coltman, D.W.; Pemberton, J.M. Inheritance of Coat Colour and Horn Type in Hirta Soay Sheep; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  12. Dominik, S.; Henshall, J.M.; Hayes, B.J. A single nucleotide polymorphism on chromosome 10 is highly predictive for the polled phenotype in Australian Merino sheep. Anim. Genet. 2012, 43, 468–470. [Google Scholar] [CrossRef]
  13. Montgomery, G.W.; Henry, H.M.; Dodds, K.G.; Beattie, A.E.; Wuliji, T.; Crawford, A.M. Mapping the Horns (Ho) locus in sheep: A fur-ther locus controlling horn development in domestic animals. J Hered. 1996, 87, 358–363. [Google Scholar] [CrossRef]
  14. Johnston, S.E.; Mcewan, J.C.; Pickering, N.K.; Kijas, J.W.; Beraldi, D.; Pilkington, J.G.; Pemberton, J.M.; Slate, J. Genome–wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population. Mol. Ecol. 2011, 20, 2555–2566. [Google Scholar] [CrossRef]
  15. Wiedemar, N.; Drögemüller, C. A 1.8-kb insertion in the 3′-UTR of RXFP2 is associated with polledness in sheep. Anim. Genet. 2015, 46, 457–461. [Google Scholar] [CrossRef] [PubMed]
  16. Lühken, G.; Krebs, S.; Rothammer, S.; Küpper, J.; Mioč, B.; Russ, I.; Medugorac, I. The 1.78-kb insertion in the 3′-untranslated region of RXFP2 does not segregate with horn status in sheep breeds with variable horn status. Genet. Sel. Evol. 2016, 48, 1–14. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, X.; Zhou, G.; Li, Q.; Zhao, D.; Chen, Y. Discovery of SNPs in RXFP2 related to horn types in sheep. Small Rumin. Res. 2014, 116, 133–136. [Google Scholar] [CrossRef]
  18. Medugorac, I.; Seichter, D.; Graf, A.; Russ, I.; Blum, H.; Göpel, K.H.; Rothammer, S.; Förster, M.; Krebs, S. Bovine Polledness—An Autosomal Dominant Trait with Allelic Heterogeneity. PLoS ONE 2012, 7, e39477. [Google Scholar] [CrossRef]
  19. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows—Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  20. Van Der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  21. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  22. Mbatchou, J.; Barnard, L.; Backman, J.; Marcketta, A.; Kosmicki, J.A.; Ziyatdinov, A.; Benner, C.; O’dushlaine, C.; Barber, M.; Boutkov, B.; et al. Computationally efficient whole-genome regression for quantitative and binary traits. Nat. Genet. 2021, 53, 1097–1103. [Google Scholar] [CrossRef]
  23. Dey, R.; Schmidt, E.M.; Abecasis, G.R.; Lee, S. A Fast and Accurate Algorithm to Test for Binary Phenotypes and Its Application to PheWAS. Am. J. Hum. Genet. 2017, 101, 37–49. [Google Scholar] [CrossRef]
  24. Chen, H.; Huffman, J.E.; Brody, J.A.; Wang, C.; Lee, S.; Li, Z.; Gogarten, S.M.; Sofer, T.; Bielak, L.F.; Bis, J.C.; et al. Efficient Variant Set Mixed Model Association Tests for Continuous and Binary Traits in Large-Scale Whole-Genome Sequencing Studies. Am. J. Hum. Genet. 2019, 104, 260–274. [Google Scholar] [CrossRef]
  25. Wu, M.C.; Lee, S.; Cai, T.; Li, Y.; Boehnke, M.; Lin, X. Rare-Variant Association Testing for Sequencing Data with the Sequence Kernel Association Test. Am. J. Hum. Genet. 2011, 89, 82–93. [Google Scholar] [CrossRef] [PubMed]
  26. Lee, S.; Wu, M.C.; Lin, X. Optimal tests for rare variant effects in sequencing association studies. Biostatistics 2012, 13, 762–775. [Google Scholar] [CrossRef] [PubMed]
  27. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  28. Koressaar, T.; Remm, M. Enhancements and modifications of primer design program Primer. Bioinformatics 2007, 23, 1289–1291. [Google Scholar] [CrossRef]
  29. Sim, Z.; Coltman, D.W. Heritability of Horn Size in Thinhorn Sheep. Front. Genet. 2019, 10, 959. [Google Scholar] [CrossRef]
  30. Pan, Z.; Li, S.; Liu, Q.; Wang, Z.; Zhou, Z.; Di, R.; Miao, B.; Hu, W.; Wang, X.; Hu, X.; et al. Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization. GigaScience 2018, 7, giy019. [Google Scholar] [CrossRef]
  31. Tian, D.; Han, B.; Li, X.; Liu, D.; Zhou, B.; Zhao, C.; Zhang, N.; Wang, L.; Pei, Q.; Zhao, K. Genetic diversity and selection of Tibetan sheep breeds revealed by whole-genome resequencing. Anim. Biosci. 2023, 36, 991–1002. [Google Scholar] [CrossRef]
  32. Kardos, M.; Luikart, G.; Bunch, R.; Dewey, S.; Edwards, W.; McWilliam, S.; Stephenson, J.; Allendorf, F.W.; Hogg, J.T.; Kijas, J. Whole-genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep. Mol. Ecol. 2015, 24, 5616–5632. [Google Scholar] [CrossRef] [PubMed]
  33. Duijvesteijn, N.; Bolormaa, S.; Daetwyler, H.D.; Van Der Werf, J.H.J. Genomic prediction of the polled and horned phenotypes in Merino sheep. Genet. Sel. Evol. 2018, 50, 28. [Google Scholar] [CrossRef]
  34. Yang, J.; Wang, D.F.; Huang, J.H.; Zhu, Q.H.; Luo, L.Y.; Lu, R.; Xie, X.L.; Salehian-Dehkordi, H.; Esmailizadeh, A.; Liu, G.E.; et al. Structural variant landscapes reveal convergent signatures of evolution in sheep and goats. Genome Biol. 2024, 25, 1–46. [Google Scholar] [CrossRef] [PubMed]
  35. Sakuma, S.; Salomon, B.; Komatsuda, T. The Domestication Syndrome Genes Responsible for the Major Changes in Plant Form in the Triticeae Crops. Plant Cell Physiol. 2011, 52, 738–749. [Google Scholar] [CrossRef] [PubMed]
  36. Wilkins, A.S.; Wrangham, R.W.; Fitch, W.T. The Domestication Syndrome in Mammals: A Unified Explanation Based on Neural Crest Cell Behavior and Genetics. Genetics 2014, 197, 795–808. [Google Scholar] [CrossRef] [PubMed]
  37. Lord, K.A.; Larson, G.; Coppinger, R.P.; Karlsson, E.K. The History of Farm Foxes Undermines the Animal Domestication Syndrome. Trends Ecol. Evol. 2020, 35, 125–136. [Google Scholar] [CrossRef] [PubMed]
  38. Poissant, J.; Davis, C.S.; Malenfant, R.M.; Hogg, J.T.; Coltman, D.W. QTL mapping for sexually dimorphic fitness-related traits in wild bighorn sheep. Heredity 2012, 108, 256–263. [Google Scholar] [CrossRef] [PubMed]
  39. Johnston, S.E.; Beraldi, D.; McRae, A.F.; Pemberton, J.M.; Slate, J. Horn type and horn length genes map to the same chromosomal region in Soay sheep. Heredity 2010, 104, 196–205. [Google Scholar] [CrossRef] [PubMed]
  40. Lincoln, G. Reproductive seasonality and maturation throughout the complete life-cycle in the mouflon ram (Ovis musimon). Anim. Reprod. Sci. 1998, 53, 87–105. [Google Scholar] [CrossRef] [PubMed]
  41. Goss, R.J. Deer Antlers: Regeneration, Function and Evolution; Academic Press: Cambridge, MA, USA, 2012. [Google Scholar]
  42. Lincoln, G.A. Teeth, horns and antlers: The weapons of sex. In The Differences between the Sexes; Cambridge University Press: Cambridge, UK, 1994; pp. 131–158. [Google Scholar]
  43. Gao, Y.; Xi, S.; Cai, B.; Wu, T.; Wang, Q.; Kalds, P.; Huang, S.; Wang, Y.; Han, S.; Pan, M.; et al. Sheep with Partial RXFP2 Knockout Exhibit Normal Horn Phenotype but Unilateral Cryptorchidism. J. Integr. Agric. 2023. [Google Scholar] [CrossRef]
  44. Lincoln, G.A.; Klandorf, H.; Anderson, N. Photoperiodic Control of Thyroid Function and Wool and Horn Growth in Rams and the Effect of Cranial Sympathectomy. Endocrinology 1980, 107, 1543–1548. [Google Scholar] [CrossRef]
  45. Wood, T.B. The Inheritance of Horns and Face Colour in Sheep. J. Agric. Sci. 1909, 3, 145–154. [Google Scholar] [CrossRef]
  46. Lee, S.H.; Van Der Werf, J.H.J.; Hayes, B.J.; Goddard, M.E.; Visscher, P.M. Predicting Unobserved Phenotypes for Complex Traits from Whole-Genome SNP Data. PLoS Genet. 2008, 4, e1000231. [Google Scholar] [CrossRef] [PubMed]
  47. Drysdale, C.M.; McGraw, D.W.; Stack, C.B.; Stephens, J.C.; Judson, R.S.; Nandabalan, K.; Arnold, K.; Ruano, G.; Liggett, S.B. Complex promoter and coding region β2 -adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness. Proc. Natl. Acad. Sci. USA 2000, 97, 10483–10488. [Google Scholar] [CrossRef]
  48. Stephens, J.C.; Schneider, J.A.; Tanguay, D.A.; Choi, J.; Acharya, T.; Stanley, S.E.; Jiang, R.; Messer, C.J.; Chew, A.; Han, J.H.; et al. Haplotype variation and linkage disequilibrium in 313 human genes. Science 2001, 293, 489–493. [Google Scholar] [CrossRef] [PubMed]
  49. Akey, J.; Jin, L.; Xiong, M. Haplotypes vs single marker linkage disequilibrium tests: What do we gain? Eur. J. Hum. Genet. Nat. 2001, 9, 291–300. [Google Scholar] [CrossRef] [PubMed]
  50. Martin, E.R.; Lai, E.H.; Gilbert, J.R.; Rogala, A.R.; Afshari, A.J.; Riley, J.; Finch, K.L.; Stevens, J.F.; Livak, K.J.; Slotterbeck, B.D.; et al. SNPing away at complex diseases: Analysis of sin-gle-nucleotide polymorphisms around APOE in Alzheimer disease. Am. J. Hum. Genet. 2000, 67, 383–394. [Google Scholar] [CrossRef] [PubMed]
  51. Windig, J.J.; Hoving-Bolink, R.A.; Veerkamp, R.F. Breeding for polledness in Holstein cattle. Livest. Sci. 2015, 179, 96–101. [Google Scholar] [CrossRef]
Figure 1. Measurement of horn length, horn base circumference, and horn base interval. The X axis represents the measured position, while the Y axis represents the measured value, which is depicted through a combination of a box plot and normal curve.
Figure 1. Measurement of horn length, horn base circumference, and horn base interval. The X axis represents the measured position, while the Y axis represents the measured value, which is depicted through a combination of a box plot and normal curve.
Animals 14 02152 g001
Figure 2. Population structure analysis in the nucleus herd of polled Tibetan sheep. (A) Population structure with K from 2 to 5; (B) cross-validation plot for determining the best K; (C) kinship plot of 103 Tibetan sheep; (D) principal component analysis.
Figure 2. Population structure analysis in the nucleus herd of polled Tibetan sheep. (A) Population structure with K from 2 to 5; (B) cross-validation plot for determining the best K; (C) kinship plot of 103 Tibetan sheep; (D) principal component analysis.
Animals 14 02152 g002
Figure 3. Manhattan plots and QQ plots of (A) horn length, (C) horn base circumference, and (E) horn base interval traits. (B) QQ plots of the horn length, (D) QQ plots of the horn base circumference, (F) QQ plots of the horn base interval.
Figure 3. Manhattan plots and QQ plots of (A) horn length, (C) horn base circumference, and (E) horn base interval traits. (B) QQ plots of the horn length, (D) QQ plots of the horn base circumference, (F) QQ plots of the horn base interval.
Animals 14 02152 g003
Figure 4. Gene Ontology (GO) and KEGG pathway analysis of the RXFP2 gene in terms of (A) horn length, (B) horn base circumference, and (C) horn base interval.
Figure 4. Gene Ontology (GO) and KEGG pathway analysis of the RXFP2 gene in terms of (A) horn length, (B) horn base circumference, and (C) horn base interval.
Animals 14 02152 g004
Figure 5. Linkage disequilibrium plot (r2) and haplotype blocks for SNPs of the RXFP2 gene. (A) A 19-kb LD block and 51-kb LD block; (B) four haplotypes within the sampled sheep population; (C) five haplotypes within the sampled sheep population.
Figure 5. Linkage disequilibrium plot (r2) and haplotype blocks for SNPs of the RXFP2 gene. (A) A 19-kb LD block and 51-kb LD block; (B) four haplotypes within the sampled sheep population; (C) five haplotypes within the sampled sheep population.
Animals 14 02152 g005
Table 1. Analysis of related SNPs associated with horn traits.
Table 1. Analysis of related SNPs associated with horn traits.
TraitChrRelated SNPsPosition (bp)p-ValueCandidate Gene
horn length10NC_040261.130,953,9721.19 × 10−18RXFP2 within
30,953,9975.98 × 10−18
30,957,8523.76 × 10−17
30,986,8798.96 × 10−17
30,953,0701.72 × 10−16
horn base circumference10NC_040261.130,953,9722.42 × 10−19
30,953,9973.24 × 10−18
30,953,0703.36 × 10−17
30,986,8798.84 × 10−17
horn base interval10NC_040261.130,935,8076.24 × 10−10
Table 2. Analysis of RXFP2 site characteristics.
Table 2. Analysis of RXFP2 site characteristics.
TraitsLoicCurrent Name (Breed)References
horn typec.789C>TTan sheep[17]
c.1117A>G
c.1125A>G
c.2059C>T
horn sizeOAR10_29685536Thinhorn sheep[29]
horn shapeOAR10 2946196889 Chinese sheep[30,31]
horn size29.45–29.55 MBighorn sheep[32]
horn phenotypeOAR10_29415140Soay sheep[14]
OAR10_29455959
OAR10_29511510
polled phenotypeOAR10_29389966Merino sheep[12]
horn phenotypeOAR10_29448537Merino sheep[33]
OAR10_29458450
OAR10_29546872
horn phenotypeOAR10_29481646Tibetan sheepa, [31]
OAR10_29469024
OAR10_30935807Tibetan sheepa
OAR10_30953070
OAR10_30953972
OAR10_30953997
OAR10_30957852
OAR10_30986879
a—this work.
Table 3. Phenotypic distribution and potential genotypes of horn type in genetically matched polled offspring.
Table 3. Phenotypic distribution and potential genotypes of horn type in genetically matched polled offspring.
SexHorn TypeSNPGenotypePolled Frequency (%)
WildMutantHomozygous
Male
(n = 193)
Normalg.30935807851043.88
g.3095307012362
g.3095397202462
g.3095399712362
g.3095785212560
g.3098687944339
Scurredg.309358073612125.00
g.3095307033412
g.3095397213612
g.3095399733412
g.3095785233412
g.3098687913279
Polledg.309358074119131.12
g.3095307021400
g.309539720610
g.3095399721400
g.3095785222390
g.3098687927340
Female
(n = 189)
Normalg.30935807572031.55
g.3095307001346
g.3095397246130
g.3095399701346
g.3095785201544
g.3098687912830
Scurredg.30935807232113.90
g.309530702816
g.3095397201016
g.309539972816
g.309578522816
g.309868793158
Polledg.309358077725054.55
g.30953070177411
g.3095397209111
g.30953997177411
g.30957852187410
g.3098687939594
Table 4. Association analyses of SNPs and genotypes in RXFP2 with horn traits.
Table 4. Association analyses of SNPs and genotypes in RXFP2 with horn traits.
SNPGenotypeHorn Length (cm)Horn Base Circumference (cm)Horn Base Interval (cm)
g.30935807AA (962)16.51 ± 14.77 Aa7.72 ± 5.57 Aa3.36 ± 2.09 Aa
GA (157)2.85 ± 7.13 Bb1.61 ± 3.92 Bb0.93 ± 1.90 Bb
GG (6)3.50 ± 4.28 ABb3.00 ± 4.69 ABb1.50 ± 2.35 ABb
p = 5.945 × 10−21p = 2.3087 × 10−17p = 0.008457
g.30953070TT (72)1.98 ± 5.57 C1.10 ± 3.29 C0.97 ± 1.77 C
CT (419)8.65 ± 13.34 B4.14 ± 5.65 B2.13 ± 2.35 B
CC (634)19.85 ± 13.95 A9.28 ± 4.73 A3.83 ± 1.79 A
p = 2.3448 × 10−9p = 2.0722 × 10−23p = 1.0313 × 10−45
g.30953972GG (28)2.88 ± 7.51 c1.55 ± 3.93 c0.57 ± 1.44 C
AG (465)8.12 ± 13.10 Bb3.86 ± 5.56 Bb2.06 ± 2.32 B
AA (632)19.78 ± 13.91 Aa9.26 ± 4.73 Aa3.83 ± 1.79 A
p = 0.012412p = 6.8297 × 10−15p = 7.1798 × 10−49
g.30953997AA (73)2.25 ± 5.99 C1.22 ± 3.43 C1.01 ± 1.80 C
GA (422)8.86 ± 13.54 B4.21 ± 5.68 B2.15 ± 2.35 B
GG (630)19.76 ± 13.92 A9.26 ± 4.73 A3.82 ± 1.80 A
p = 1.0693 × 10−8p = 4.3225 × 10−23p = 1.0569 × 10−45
g.30957852TT (77)1.85 ± 5.41 C1.03 ± 3.19 C0.96 ± 1.76 C
CT (424)8.81 ± 13.30 B4.26 ± 5.67 B2.17 ± 2.35 B
CC (624)20.00 ± 13.97 A9.31 ± 4.71 A3.84 ± 1.79 A
p = 1.7251 × 10−10p = 2.2845 × 10−26p = 1.267 × 10−47
g.30986879TT (176)3.33 ± 8.47 C1.77 ± 4.15 C1.08 ± 1.99 C
CT (490)11.63 ± 13.25 B5.79 ± 5.73 B2.80 ± 2.32 B
CC (459)21.94 ± 14.31 A9.90 ± 4.44 A3.98 ± 1.60 A
p = 5.5701 × 10−14p = 1.2202 × 10−39p = 5.9519 × 10−57
g.29469024AA (131)2.90 ± 8.96 C1.36 ± 3.78 C0.60 ± 1.57 C
GA (463)12.31 ± 14.61 B5.80 ± 5.82 B2.54 ± 2.34 B
GG (531)19.36 ± 13.88 A9.10 ± 4.92 A4.02 ± 1.560 A
p = 9.9762 × 10−15p = 51.2842 × 10−35p = 1.1715 × 10−73
g.29481646TT (833)16.56 ± 14.47 A7.82 ± 5.47 A3.45 ± 2.04 A
CT (27)5.85 ± 11.51 B2.65 ± 5.12 B1.13 ± 1.98 B
CC (265)9.05 ± 14.10 B4.20 ± 5.77 B1.85 ± 2.31 B
p = 1.4899 × 10−14p = 0.01p = 5.599 × 10−11
g.29461968GG (306)2.72 ± 8.85 C1.25 ± 3.64 C0.56 ± 1.53 C
AG (434)16.70 ± 14.65 B7.88 ± 5.54 B3.45 ± 2.07 B
AA (385)21.49 ± 12.84 A10.11 ± 3.97 A4.47 ± 0.87 A
p = 6.7264 × 10−34p = 2.828 × 10−86p = 75.0538 × 10−47
Within the same line, the mean values of different superscript lowercase letters were significantly different (p < 0.05). Within the same line, the mean values of different superscript uppercase letters differ significantly (p < 0.01).
Table 5. Association analysis of haplotypes in RXFP2 with horn traits in Module 1.
Table 5. Association analysis of haplotypes in RXFP2 with horn traits in Module 1.
LociHaplotypesNumbersHorn Length (cm)Horn Base Circumference (cm)Horn Base Interval (cm)
g.29481646
g.29469024
g.29461968
AAGGCT110.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
AAGGTT1203.16 ± 9.32 B1.48 ± 3.92 B0.65 ± 1.63 B
GAGACT90.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
GAGATT32416.46 ± 14.55 A7.79 ± 5.46 A3.41 ± 2.08 Aa
GAGGCC1222.99 ± 9.73 B1.34 ± 3.77 B0.59 ± 1.60 B
GAGGTT80.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
GGAACC323.83 ± 20.52 A10.33 ± 4.75 A4.50 ± 1.00 Aa
GGAACT722.57 ± 11.58 A10.21 ± 4.84 A4.36 ± 0.75 Aa
GGAATT37521.45 ± 12.83 A10.11 ± 3.96 A4.48 ± 0.88 Aa
GGGACC9619.53 ± 14.65 A9.12 ± 5.32 A4.00 ± 1.72 Aa
GGGATT58.30 ± 11.85 A4.40 ± 6.12 A1.80 ± 2.49 Ab
GGGGCC442.00 ± 6.68 B0.98 ± 3.20 B0.48 ± 1.38 B
Within the same line, the mean values of different superscript lowercase letters were significantly different (p < 0.05). Within the same line, the mean values of different superscript uppercase letters differ significantly (p < 0.01).
Table 6. Association analysis of haplotypes in RXFP2 with horn traits in Module 2.
Table 6. Association analysis of haplotypes in RXFP2 with horn traits in Module 2.
LociHaplotypesNumbersHorn Length (cm)Horn Base Circumference (cm)Horn Base Interval (cm)
g.30935807
g.30953070
g.30953972
g.30953997
g.30957852
g.30986879
AACCAAGGCCCC36322.17 ± 13.79 Aa10.04 ± 4.14 Aa4.07 ± 1.51 Aa
AACCAAGGCCCT22217.49 ± 13.23 B8.68 ± 4.95 B3.73 ± 1.95 Ab
AACCAAGGCCTT2512.18 ± 15.80 B5.50 ± 6.68 B1.96 ± 2.34 B
AACCAAGGCTCC1110.68 ± 8.21 B6.91 ± 4.89 Ab3.41 ± 2.22 Aa
AACCAAGGCTCT20.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
AACTGAGACCCC33.33 ± 5.77 B2.33 ± 4.04 B2.00 ± 3.46 Aa
AACTGAGACTCC7023.78 ± 16.71 Aa10.16 ± 5.28 Aa3.86 ± 1.70 Aa
AACTGAGACTCT1778.19 ± 11.83 B4.06 ± 5.42 B2.43 ± 2.38 B
AACTGAGACTTT483.57 ± 8.35 B2.02 ± 4.40 B1.38 ± 2.09 B
AATTGAAATTCT60.33 ± 0.82 B0.42 ± 1.02 B2.17 ± 2.40 Aa
AATTGAAATTTT180.00 ± 0.00 B0.00 ± 0.00 B1.11 ± 1.84 B
AATTGGAATTCT512.60 ± 12.80 Aa6.70 ± 6.32 Aa2.40 ± 2.27 Aa
AATTGGAATTTT30.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
GACCAAGGCCCC616.33 ± 13.93 Aa8.08 ± 6.48 Aa2.83 ± 2.32 Aa
GACTGAGACTCT663.85 ± 8.63 B2.07 ± 4.64 B1.07 ± 1.95 B
GACTGAGACTTT471.48 ± 3.29 B1.20 ± 2.70 B0.90 ± 2.07 B
GACTGAGATTCT40.00 ± 0.00 B0.00 ± 0.00 B2.00 ± 2.31 Ab
GATTGAAATTTT161.63 ± 3.34 B0.69 ± 2.75 B0.50 ± 1.37 B
GATTGGAATTCT20.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
GATTGGAATTTT150.00 ± 0.00 B0.00 ± 0.00 B0.00 ± 0.00 B
GGCTGAGACTCT23.50 ± 4.95 Ab4.00 ± 5.66 Aa2.00 ± 2.83 Aa
GGTTGGAATTTT43.50 ± 4.73 B2.50 ± 5.00 B1.25 ± 2.50 B
Within the same line, the mean values of different superscript lowercase letters were significantly different (p < 0.05). Within the same line, the mean values of different superscript uppercase letters differ significantly (p < 0.01).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, D.; Zhang, Z.; Huang, B.; Han, B.; Li, X.; Zhao, K. Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep. Animals 2024, 14, 2152. https://doi.org/10.3390/ani14152152

AMA Style

Tian D, Zhang Z, Huang B, Han B, Li X, Zhao K. Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep. Animals. 2024; 14(15):2152. https://doi.org/10.3390/ani14152152

Chicago/Turabian Style

Tian, Dehong, Zian Zhang, Bin Huang, Buying Han, Xue Li, and Kai Zhao. 2024. "Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep" Animals 14, no. 15: 2152. https://doi.org/10.3390/ani14152152

APA Style

Tian, D., Zhang, Z., Huang, B., Han, B., Li, X., & Zhao, K. (2024). Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep. Animals, 14(15), 2152. https://doi.org/10.3390/ani14152152

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop