Polymorphism Analysis of NOTCH2 and CD1A Genes and Their Association with Wool Traits in Subo Merino Sheep

TL;DR

This study identifies genetic variations in NOTCH2 and CD1A genes linked to wool traits in Subo Merino sheep, offering insights for breeding better wool quality.

Contribution

The study discovers specific SNPs in NOTCH2 and CD1A genes associated with wool traits and their potential impact on protein structures.

Findings

Specific SNPs in NOTCH2 are linked to fiber diameter variation and greasy fleece weight.

CD1A SNPs are significantly associated with fiber diameter standard deviation and crimp number.

CD1A gene is highly expressed in fine wool, suggesting a role in regulating wool quality.

Abstract

To identify the genetic factors influencing wool quality in Subo Merino sheep, this study focused on the NOTCH2 and CD1A genes. Through genetic analysis, we detected six single nucleotide polymorphisms (SNPs) within these genes and examined their associations with key wool traits. Our results indicated that specific SNPs in NOTCH2 are closely linked to the coefficient of variation of fibre diameter and greasy fleece weight, while SNPs in CD1A are significantly associated with the standard deviation of fibre diameter and crimp number. We also found that these genetic variations may alter the secondary and tertiary structures of the proteins encoded by NOTCH2 and CD1A, potentially affecting their biological functions. Furthermore, qPCR results demonstrated that CD1A is higher expressed in fine wool, reinforcing its potential role in regulating wool quality. Overall, this research identifies promising candidate SNPs that could serve as molecular markers, providing practical guidance for breeding Subo Merino sheep with improved wool traits and enhancing our understanding of the genetic mechanisms underlying wool quality.

To identify molecular markers associated with wool traits in fine-wool sheep, we examined genetic polymorphisms in the NOTCH2 and CD1A genes in 944 Subo Merino sheep in this study. Subsequently, we performed association analyses between mutation sites in the NOTCH2 and CD1A genes and wool traits using SAS 9.4 software, followed by linkage disequilibrium (LD) analysis of different mutation sites using Haploview 4.2 software. Additionally, bioinformatics tools were employed to predict the potential impacts of missense mutations on protein secondary and tertiary structures. Finally, quantitative PCR (qPCR) was used to assess the expression levels of the NOTCH2 and CD1A genes. Genetic analysis revealed six polymorphic sites in NOTCH2 and CD1A, all of which were missense mutations. Two SNPs in NOTCH2 (SNP1 and SNP2) showed significant associations with the coefficient of variation of fibre diameter, and SNP1 was also associated with greasy fleece weight. Four SNPs in CD1A (SNP3–SNP6) were significantly associated with fibre diameter standard deviation, and SNP3, SNP4, and SNP5 were additionally associated with crimp number. LD analysis revealed that SNP3, SNP4, and SNP5 were closely linked. Bioinformatics analysis indicated that the mutations caused alterations in the secondary and tertiary structures of the NOTCH2 and CD1A proteins. qPCR results showed that the CD1A gene was highly expressed in the fine wool fibre group compared with the ultra-fine wool fibre group. In conclusion, this study revealed a genetic association between NOTCH2 and CD1A and wool traits. The results are expected to provide a theoretical foundation for breeding wool traits in Subo Merino sheep, thereby enhancing the economic value of fine wool.

Figures

Full text

1. Introduction

Fine-wool sheep are predominantly found in the pastoral and semi-pastoral regions of central and western China, as well as in ethnic border areas. They play a vital role in Xinjiang’s animal husbandry and have become a key industry for improving the livelihoods of local farmers and herders [1]. However, in recent years, the fine-wool sheep industry in Xinjiang has faced significant challenges due to the increasing dominance of meat sheep production. The widening price gap between meat and wool has reduced profitability in fine-wool breeding, leading to declining wool output, low breeding efficiency, and inferior wool quality. Although domestic wool is influenced by many factors, as a high-grade raw material in the textile industry, the long-term viability of wool is an inevitable trend [2]. Currently, in the wool market, the price of ultrafine wool is two to three times that of ordinary wool. Therefore, high-quality fine wool is the focus of market attention and a key target for scientific research breakthroughs [3]. However, traditional breeding methods can no longer meet the demands of wool production. To accelerate the development of the wool industry, molecular breeding is particularly important in the breeding and selection of fine-wool sheep [4].

Single nucleotide polymorphisms (SNPs) are considered the most promising molecular markers due to their abundance, wide distribution, and rich informational content. With the continual reduction in the costs of chip and sequencing technologies, SNP marker technology has been widely adopted in animal breeding research. The polymorphisms of numerous genes and their associations with wool traits have gradually been elucidated. Studies have identified significant associations between specific SNP variants in the PTPN3, KRT83, KRT39, CCSER1, RPS6KC1, KCNRG, KCNK9, and CLYBL genes and wool or cashmere traits. These genes may serve as valuable candidates for improving wool and cashmere quality, providing important references for breeding programmes and research into genetic mechanisms [5,6,7,8].

In preliminary studies, our research group identified several candidate genes (FZD3, ARPP21, LMNB1, RASA1, PAK1, IFNAR2, FAT3, CD1A, and NOTCH2) through genome-wide association analysis [9]. Notably, members of the Notch family play a crucial role in hair follicle morphogenesis, regulating the proliferation and differentiation of follicular cells across multiple developmental stages. Specifically, NOTCH1, NOTCH2, and NOTCH3 exhibit distinct spatial expression patterns within different hair follicle cell layers [10,11,12]. He et al. [13] demonstrated that NOTCH1 expression peaks at embryonic day 105 (E105) during ovine skin development, suggesting its pivotal role in folliculogenesis. Complementary findings by Vauclair et al. [14] confirmed NOTCH1’s essential functions in late-stage embryonic follicle development and postnatal hair cycle regulation. Similarly, NOTCH2 has been implicated in cutaneous differentiation and follicular development [15]. Regarding CD1A, this transmembrane protein, a member of the CD1 family, plays dual roles in immune regulation and antigen presentation. As a dominant skin antigen-presenting molecule, CD1A-targeted neutralising antibodies show therapeutic potential for dermatological and systemic disorders [16]. Recently, it has been found that changes in CD1A-dependent T cell responses are associated with the pathogenesis of various inflammatory skin diseases [17]. Therefore, we conclude that the distribution of CD1A in skin immune cells is closely related to the pathological state of the skin, which may affect the immune environment surrounding hair follicles. In summary, the NOTCH2 and CD1A genes play important roles in hair growth and can be used as candidate genes affecting fine wool traits for SNP mining.

Subo Merino sheep, a novel superfine-wool breed developed in 2014, exhibit remarkable stress resistance, high reproductive survival rates, superior wool characteristics (17–19 μm fibre diameter), and high wool yield. This breed has filled China’s gap in fine-wool sheep production (80 s wool grade), enriched the nation’s fine-wool germplasm resources, and enhanced both domestic wool quality and international market competitiveness [18,19,20]. Therefore, in this study, the SNPs of NOTCH2 and CD1A genes were identified in Subo Merino sheep, and the genetic effects of these SNPs on wool traits were analysed in depth, aiming to provide new markers for the molecular breeding of wool traits in fine-wool sheep.

2. Materials and Methods

2.1. Phenotype Measurement and Sample Collection

In this study, total of 944 one-year-old Subo Merino ewes were obtained from two regions in Xinjiang: Yili Gongnaisi Sheep Farm (n = 473) and Aksu Baicheng Sheep Farm (n = 471). Sample collection involved the following four steps:

• (1)Blood samples were collected from the jugular vein of these 944 Subo Merino sheep using heparin anticoagulant tubes. After heparin is mixed with blood, the blood samples were immediately placed in an icebox, further transported to the laboratory, and finally stored at −20 °C refrigerator for DNA extraction;
• (2)While taking blood samples, the greasy fleece weight (GFW), live weight before shearing (LWBS), and live weight after shearing (LWAS) of these sheep were measured and recorded. Additionally, the staple length (SL), fineness count (FC), crimp, hair length (HL), and crimp number (CN) of these sheep were also measured and recorded [21];
• (3)Wool samples were collected from a site 10 cm posterior to the left scapular edge (midline region). Wool samples were further washed using the conventional washing process and allowed to dry naturally. Measurements of the mean fibre diameter (MFD), coefficient of variation of fibre diameter (CVFD), and fibre diameter standard deviation (FDSD) were taken in a laboratory maintained at a constant temperature (20 ± 2 °C) and humidity (65 ± 4%) using a fibre diameter optical analyser (OFDA2000, Ningbo Jiangnan Instrument Factory, Ningbo, China) [22]. The parameters measured included. Excel 2019 was used to compile the data on wool traits. SPSS 27.0 software [23] was employed to perform descriptive statistical analyses on the relevant wool trait data;
• (4)Based on the results of MFD measurements, 10 sheep with the smallest MFD were designated as the ultra-fine wool fibre group (UFW, 16.21 ± 0.46 μm), and 10 with the largest MFD as the fine wool fibre group (FW, 20.68 ± 0.93 μm). For these two groups, additional wool samples were collected from the left forelimb and 5 cm posterior to the scapula and were further used to measure the MFD. 20 skin tissue samples (approximately 2 cm × 2 cm) from these two groups were collected using a skin sampler, immediately frozen in liquid nitrogen, and were further stored at −80 °C refrigerator for RNA extraction.

2.2. DNA Extraction and SNP Typing

Total of 944 blood genomic DNA was extracted using the Blood/Cell/Tissue Genomic DNA Extraction Kit (DP304) from Tiangen Biochemical Technology (Beijing) Co., Ltd. (Beijing, China). The quality of the DNA was assessed by 1.0% agarose gel electrophoresis, and its concentration was measured using an NanoDrop™ 2000 devices (Thermo Fisher Scientific, Waltham, MA, USA). Samples with a concentration greater than 20 ng/μL and an OD260/OD280 ratio between 1.7 and 1.9 were deemed suitable for the experiment. Subsequently, SNP typing was performed using Fluidigm’s Biomark™ HD system (Biomark™ HD, San Francisco, CA, USA).

2.3. Genetic Diversity Analysis

Excel 2019 was utilised to organise the data and to calculate genotype frequencies, allele frequencies, and genetic diversity parameters, including homozygosity (Ho), heterozygosity (He), effective number of alleles (Ne), polymorphism information content (PIC), and Hardy–Weinberg equilibrium for each genotype within the population [24,25]. Additionally, when the Hardy–Weinberg equilibrium test yields a p > 0.05, it indicates that the population is in Hardy–Weinberg equilibrium.

The calculation formulas for Ho, He, Ne, PIC and $[eqn]$ are as follows:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where m is the number of alleles, $[eqn]$ and $[eqn]$ are the frequencies of the i-th and j-th alleles in the population, $[eqn]$ is the expected frequency, and $[eqn]$ is the observed frequency.

Subsequently, linkage disequilibrium (LD) analysis of the mutation sites was conducted using Haploview 4.2 software [26]. Finally, haplotype analysis was carried out using the geneHapR package [27] in the R language 4.4.2.

2.4. Correlation Analysis

According to the results of wool fibre diameter measurements, the correlation between different genotypes and wool traits was analysed using the GLM procedure in SAS 9.4 [28]. The least-squares variance analysis method was employed, with genotype and field effects treated as fixed effects. During the GLM analysis, t-tests were conducted between levels. Concurrently, Tukey’s method was employed for multiple comparison correction. The results are presented as least-squares means ± standard error. p < 0.05 indicates a significant difference, and p < 0.01 indicates a highly significant difference.

The linear model is as follows:

$$[eqn]$$

In the form, $[eqn]$ : sheep individual phenotypic value; $[eqn]$ : group mean $[eqn]$ : genotypic SNP effect; $[eqn]$ : field effect; $[eqn]$ : random error.

2.5. Biological Function Prediction

Non-synonymous SNPs result in changes to the protein sequence before and after mutation. Using information on different genotypes at non-synonymous mutation sites, the wild-type and mutant gene sequences at these sites in the NOTCH2 and CD1A genes were retrieved from the Ensembl nucleic acid database. SOPMA [29] was employed to predict the secondary structure of the protein before and after mutation, while SWISS-MODEL [30] was used to predict the tertiary structure, with all parameters set to default.

2.6. RNA Extraction

Total RNA was extracted using the TRIzol method. The OD260/280 and OD260/230 ratios of the total RNA were measured with a NanoDrop™ 2000 devices (Thermo Fisher Scientific, MA, USA), and the integrity of the RNA (RIN) was assessed using an Agilent 2100 Bioanalyser (Agilent Technologies, Santa Clara, CA, USA). The OD260/280 ratio ranged from 1.8 to 2.0, the OD260/230 ratio was greater than 2.0, and the RIN value ranged from 7.0 to 8.5, indicating suitability for subsequent experiments. Next, reverse transcription was performed using the Takara reverse transcription kit (DRR036A) from Bao Bioengineering (Dalian) Co., Ltd. (Dalian, China).

2.7. Primer Design and qPCR

According to the NCBI database, the sheep NOTCH2 gene sequence (accession number: XM_060401816.1) and the CD1A gene sequence (accession number: XM_042256151.1) were obtained. Primer pairs were designed using Primer Premier 5.0 and validated via NCBI Primer-BLAST (https://blast.ncbi.nlm.nih.gov, accessed on 30 May 2025) to select suitable candidates for the assay. The primers were then commercially synthesised by Bioengineering (Shanghai) Co., Ltd. (Shanghai, China). Simultaneously, the GAPDH gene was used as the internal reference gene [31]. The list of primers is shown in Table 1. Subsequently, the reaction conditions for real-time PCR were established by referring to the system and parameters recommended in the Talent qPCR PreMix (SYBR Green) kit manual (FP209) from Tiangen Biochemical Technology (Beijing) Co., Ltd. (Beijing, China). qPCR detection was performed using the CFX96™ Real-Time System instrument (Bio-Rad, Hercules, CA, USA). The amplification procedure was as follows: pre-denaturation at 95 °C for 3 min, followed by denaturation at 95 °C for 5 s, annealing at 60 °C for 10 s, and extension at 72 °C for 15 s, with fluorescence signal collection. The last three steps were repeated 40 times, followed by a dissociation curve stage. The qPCR results were calculated using the 2^−ΔΔCt^ algorithm [32].

3. Results

3.1. Descriptive Statistics of Wool Traits

Descriptive statistics for wool traits in Subo Merino sheep were calculated, and the results are presented in Table 2. The MFD was 17.71 μm, ranging from 12.80 to 24.00 μm, with FDSD of 1.84 μm and CVFD of 10.39%. The average values and ranges of each trait correspond with objective data, and the small standard deviation indicates that the variables are relatively concentrated around the mean.

3.2. Analysis of Mutation Sites in the NOTCH2 and CD1A Genes

Six mutation sites were identified in the NOTCH2 and CD1A genes of Subo Merino sheep, all of which were missense mutations. SNP1 is located in exon 30 of the NOTCH2 gene, and SNP2 is located in exon 24 of the NOTCH2 gene, with A → C and G → A mutations, respectively. SNP3 is located in the second exon of the CD1A gene, with a G → A mutation. SNP4 and SNP5 are located in the third exon of the CD1A gene, with A → C and G → A mutations, respectively. SNP6 is located in the fourth exon of the CD1A gene, with an A → T mutation. Detailed information on the specific mutation sites is provided in Table 3.

3.3. Analysis of Genotype Frequencies and Allele Frequencies of NOTCH2 and CD1A Genes

As shown in Table 4, for SNP1 of the NOTCH2 gene, the GG genotype and G allele were predominant, whereas for SNP2, the CC genotype and C allele were predominant. For the CD1A gene, SNP3 and SNP5 shared the same predominant genotype (GG) and allele (G); SNP4 was characterised by the AA genotype and A allele; and SNP6 was characterised by the TT genotype and T allele. The chi-square test indicated that SNP1, SNP2, SNP3, SNP4, and SNP5 were in Hardy–Weinberg equilibrium (p > 0.05), whereas SNP6 deviated from Hardy–Weinberg equilibrium (p < 0.05).

3.4. Population Genetic Analysis of NOTCH2 and CD1A Genes

The genetic diversity of SNPs in Subo Merino sheep was analysed, and the results are presented in Table 5. As shown, the Ho for the SNPs in the NOTCH2 and CD1A genes were 0.739, 0.783, 0.559, 0.556, 0.554, and 0.753, respectively. The Ne for these SNPs were 1.36, 1.285, 1.737, 1.741, 1.764, and 1.797, respectively. With regard to PIC, the values for SNP3, SNP4, SNP5, and SNP6 ranged between 0.25 and 0.5, indicating moderate polymorphism. In contrast, the PIC values for SNP1 and SNP2 were less than 0.25, corresponding to low polymorphism at these loci. These results suggest that the CD1A gene loci exhibit relatively high polymorphism.

3.5. Genetic Effects of NOTCH2 and CD1A Genes on Wool Traits

3.5.1. Variance Analysis of Different Genotypes of NOTCH2 and CD1A Genes on Wool Traits

The effects of SNPs in the NOTCH2 and CD1A genes on the wool traits of Subo Merino sheep were analysed using least-squares analysis of variance. The results are presented in Table 6. SNP1 had a significant effect on MFD (p < 0.05), SNP2 significantly affected the CVFD (p < 0.05), and SNP4 significantly influenced SL (p < 0.05). SNP6 had a significant effect on FDSD and CVFD in this group (p < 0.05), while the remaining SNPs showed no significant effects (p > 0.05).

3.5.2. Association Analysis of NOTCH2 and CD1A Genes with Wool Traits

Table 7 presents the results of association analyses between SNPs in the NOTCH2 and CD1A genes and wool traits. For FDSD, significant associations were observed with SNP3, SNP4, and SNP5 (p < 0.05), a highly significant association with SNP6 (p < 0.01), and no significant associations with SNP1 or SNP2 (p > 0.05). MFD, SL, FC, LWAS, and LWBS showed no significant associations with any SNPs (p > 0.05). Regarding SNP3, individuals with GA or GG genotypes exhibited significant associations with HL and crimp compared to those with the AA genotype (p < 0.05), while AA or GA genotypes were significantly associated with CN compared to the GG genotype (p < 0.05). For SNP4, AA and AC genotypes were significantly associated with crimp relative to the CC genotype (p < 0.05), and CN differed significantly between AC and CC genotypes (p < 0.05). For SNP5, AA and GA genotypes showed significant associations with CN compared to the GG genotype (p < 0.05). In SNP6, AT and TT genotypes were significantly associated with CVFD relative to the AA genotype (p < 0.05). Regarding SNP1, CVFD differed significantly among GG, TG, and TT genotypes (p < 0.05), and GFW in the TT genotype showed significant differences compared to the GG and TG genotypes (p < 0.05). In SNP2, CVFD in the CC genotype differed significantly from that in the TT and TC genotypes (p < 0.05), while no significant associations were detected between other wool traits and genotypes (p > 0.05).

3.6. LD and Haplotype Analysis

LD analysis was conducted on six SNPs of the NOTCH2 and CD1A genes using Haploview 4.2 software. The linkage between SNPs is illustrated in Figure 1A. The colour of each block ranges from light to dark (white to red), indicating the degree of linkage from low to high; a darker colour signifies stronger linkage. The D′ and r^2^ values range from 0 to 1, with higher values indicating a greater degree of linkage. An r^2^ value of 1 denotes complete LD, while r^2^ > 0.6 indicates strong LD. As shown in Table 8, the D’ values for SNP1 and SNP2, SNP3 and SNP4, SNP3 and SNP5, and SNP4 and SNP5 are all 1, indicating a high correlation and close linkage between these SNP pairs.

A haplotype refers to the combination of alleles at multiple loci that are inherited together on the same chromosome. In the Subo Merino sheep population, three haplotypes of the NOTCH2 gene were identified (Figure 1B). Among these, the H001 haplotype, primarily GC, had a frequency of 69.44%; H002, primarily TC, had a frequency of 28.74%; and H003 was mainly GT. Additionally, four haplotypes were observed in the CD1A gene (Figure 1C). The H001 haplotype, predominantly GAGA, had a frequency of 45.76%; H002, mainly GAGT, had a frequency of 44.55%; H003 was primarily ACAT; and H004 was mainly GAAT.

3.7. Analysis of Protein Structure Changes

An analysis of the amino acid substitutions at the mutation sites and their impact on protein secondary structure revealed that the proteins encoded by the NOTCH2 and CD1A genes comprise three types of secondary structures, primarily random coils, α-helices, and extended strands (Table 9). The mutation in SNP1 resulted in the amino acid at position 1805 changing from methionine to leucine, which led to a decrease in the proportion of α-helices and an increase in the proportion of extended strands. The mutation in SNP2 caused the amino acid at position 1307 to change from valine to methionine, resulting in a reduction in the proportion of random coils and an increase in extended strands. The mutation in SNP4 led to the substitution of histidine by proline at position 123, causing an increase in the proportion of α-helices and a decrease in random coils. The mutation in SNP6 resulted in the amino acid at position 270 changing from glutamic acid to valine, which increased the proportion of random coils and decreased the proportion of extended strands. The mutation in SNP3 replaced aspartic acid at position 78 with asparagine, while the mutation in SNP5 substituted alanine at position 127 with threonine. However, these two mutations did not alter the proportions of the protein’s secondary structures. Homology modelling of the protein tertiary structures before and after the missense mutations, performed using SWISS-MODEL, showed that the predicted tertiary structures were consistent with the secondary structure predictions, with random coils predominating in the protein structure (Figure 2).

3.8. qPCR Results of NOTCH2 and CD1A Genes

Based on wool fibre diameter measurements, the ten individuals with the lowest MFD were classified as the FW (mean: 20.68 ± 0.93 μm), while the ten with the highest MFD were designated as the UFW (mean: 16.21 ± 0.46 μm). As shown in Figure 3, the expression levels associated with different MFD in Subo Merino sheep exhibited highly significant differences between the UFW and FW (p < 0.01, Figure 3A). Furthermore, the expression levels of the NOTCH2 and CD1A genes were analysed in both groups. The results demonstrated that both genes were expressed at higher levels in the FW compared to the UFW. However, while NOTCH2 showed no statistically significant difference (p > 0.05, Figure 3B), CD1A exhibited a higher significant differential expression (p < 0.01, Figure 3C). In conclusion, the NOTCH2 and CD1A genes may contribute to the observed differences in MFD between the two groups.

4. Discussion

The growth and development of wool are regulated by the expression of related genes at the molecular level, with variations in wool traits resulting from mutations in key genes [33]. Studies by Ma et al. [34,35], Wang et al. [36], Zhao et al. [37], and Yue et al. [38] have demonstrated, through genome-wide association studies and gene polymorphism validation, that SNPs in genes such as ALX4, keratin-associated protein 2-1, KIF16B, SLIT3, and ZNF280B are associated with traits including SL, crimp, MFD, and wool production. In summary, analysing the genetic mechanisms underlying wool quality at the molecular level is of great significance for the breeding of fine-wool sheep.

Studies have shown that Notch signalling pathway is one of the top ten gene families implicated in hair follicle development and villus growth [39]. Bai et al. [40] analysed the correlation between NOTCH2 gene polymorphisms and villus and growth traits in Shaanbei white cashmere goats. Their results demonstrated that two mutation sites, rs665021370 and rs653705114, in the NOTCH2 gene had significant effects on villus and growth traits. Dilinare [41] investigated the genetic effects of the KRT85, NOTCH2, and ADAM9 genes on wool traits in Chinese Merino sheep (Xinjiang type). The findings revealed that the g.96438799 C>T locus of the NOTCH2 gene had an extremely significant effect on SL (p < 0.01), and the g.96432471 T>G locus had an extremely significant effect on CN (p < 0.01). The results of this study strongly support the significant role of the NOTCH2 gene in regulating wool traits. Specifically, NOTCH2 gene SNPs 1 and 2 were significantly associated with CVFD, while SNP1 was significantly associated with GFW. However, no significant association was found between the NOTCH2 gene and CN or SL, which is inconsistent with the findings of Dilinare [41]. We therefore speculate that this discrepancy may be attributable to variations in breed, sample size, and rearing conditions. To explore this further, qPCR analysis was conducted. Although NOTCH2 gene expression tended to be higher in the FW group, no significant difference was observed between the FW and UFW groups. This suggests that NOTCH2 gene may not play a major role in regulating MFD at the transcriptional level in this context. Instead, its potential influence on CVFD or GFW is likely mediated through post-transcriptional mechanisms or gene interactions. Since the CD1A gene was identified through GWAS screening and functional prediction as influencing the fineness of fine-wool sheep in our earlier work [9], further exploration of the CD1A gene is warranted. The qPCR results of this study also support the association of the CD1A gene with wool MFD. Mitchell et al. [42] identified that CD1A T cells may be involved in the pathology of various diseases, including cancer and autoimmune disorders. The results of this study showed that SNP3, SNP4, and SNP5 of the CD1A gene were significantly correlated with FDSD and CN, while SNP6 was significantly correlated with FDSD. SNP3 was significantly associated with HL (p < 0.05). Additionally, SNP3 and SNP4 were significantly associated with crimp (p < 0.05), and SNP6 was significantly associated with CVFD. These findings further support the important role of the CD1A gene in regulating wool traits.

PIC and heterozygosity are indicators used to assess genetic polymorphism within a population [43]. In the Subo Merino sheep population, the PIC values for the two mutation sites in the NOTCH2 gene ranged from 0.197 to 0.230, indicating low polymorphism. In contrast, the four mutation sites in the CD1A gene exhibited PIC values between 0.334 and 0.345, reflecting moderate polymorphism and suggesting a moderate level of genetic variation with relatively abundant polymorphism. Notably, SNP6 was found to deviate from Hardy–Weinberg equilibrium, which may indicate that factors such as natural selection, genetic drift, non-random mating, or mutation are influencing the population’s gene frequencies. This warrants further investigation into the causal factors and their impact on the population’s adaptability and viability. LD analysis revealed strong LD between SNP1 and SNP2, SNP3 and SNP4, SNP3 and SNP5, and SNP4 and SNP5, indicating non-random allelic associations and close linkage due to their physical proximity on the same chromosome, making them less likely to separate during genetic transmission. Additionally, this study found that SNP6 was in LD with other loci. This phenomenon may result from the unique genetic effect of this locus or the presence of other linked loci acting together on wool traits. Therefore, it is necessary to evaluate these SNPs separately, considering their different linkage relationships, to determine their breeding value and specific applications in breeding practices.

The spatial structure of proteins is closely related to their function and provides valuable biological information [44]. Studies have shown that SNPs located in gene exons can directly cause changes in the encoded amino acids, thereby affecting protein structure, function, and expression levels [45]. Existing research has confirmed that missense mutations can influence phenotypes by altering protein structure [46]. For example, missense mutations in the TLR2 gene can affect immune function [47]; missense mutations in the MYO7B, KIF13B, and LOC101121854 genes are associated with sheep body weight [48]; and missense mutations in the KRTAP36-2 and FST genes regulate wool yield and wool fibre diameter, respectively [49,50]. Therefore, it is necessary to investigate the relationship between changes in the secondary and tertiary structures of NOTCH2 and CD1A gene proteins and wool-related traits. In this study, mutations in the NOTCH2 gene at SNP1 and SNP2 resulted in changes such as a decrease in the proportion of α-helices, an increase in the proportion of extended strands, and a decrease in the proportion of random coil. As NOTCH2 is involved in signalling pathways related to hair follicle development, changes in its protein structure may disrupt the accuracy and efficiency of signal transmission, ultimately affecting hair follicle morphogenesis, cyclic growth, and wool traits such as hair thickness and density. The secondary structure of the CD1A gene remained unchanged before and after mutations at the SNP3 and SNP5 sites. However, following mutations at the SNP4 and SNP6 sites, the secondary structure underwent changes including an increase in the proportion of α-helices, a decrease in the proportion of random coil, and a decrease in the proportion of extended strands. We speculate that these changes may affect the recognition and binding of the CD1A protein with other proteins, which is crucial for maintaining the stability of the hair follicle microenvironment.

5. Conclusions

In this study, six mutation sites were identified in the NOTCH2 and CD1A genes of Subo Merino sheep. SNP1 and SNP2 of the NOTCH2 gene were significantly correlated with CVFD and GFW, while SNP3, SNP4, SNP5, and SNP6 of the CD1A gene were significantly correlated with FDSD and CN. The results indicate that the NOTCH2 and CD1A genes can influence their functions by altering the secondary and tertiary structures of the proteins. Furthermore, qPCR analysis demonstrated that the CD1A gene is higher expressed in the FW. In summary, the findings of this study provide a foundation for further investigation into the genetic mechanisms underlying the wool traits of Subo Merino sheep. The identified mutation sites are expected to serve as potential molecular markers affecting the wool performance of this breed and offer a reference for molecular breeding aimed at improving wool quality.