Expression Localization of the KRT32 Gene and Its Association of Genetic Variation with Wool Traits

Changes in keratin gene expression and spatiotemporal regulation determine the compositional content and cellular localization of wool keratin, thereby affecting wool traits. Therefore, keratin gene family member 32 (KRT32) was selected for a study using RT-qPCR, immunofluorescence, and penta-primer amplification refractory mutation system (PARMS) techniques. The results showed that KRT32 mRNA was highly expressed in the skin and localized to the inner root sheath (IRS), outer root sheath (ORS) and dermal papilla (DP). Sequencing results identified eight SNPs in KRT32, and association analyses revealed that the variations were significantly associated with multiple traits in wool (p < 0.05), including MFD, CF and MFC. The constructed haplotype combination H2H3 has higher CF and smaller MFD than other haplotype combination (p < 0.05). In conclusion, KRT32 can be used as a candidate gene for molecular genetic improvement of wool in Gansu Alpine Fine-wool sheep.


Introduction
Wool is one of the most moisture-sensitive natural fibers [1] and is unique in nature in its ability to both absorb moisture and wick away perspiration.According to X-ray diffraction analysis, the structure of keratin consists of an α-helical structure and a β-folded lamellar structure [2].There are two main types of wool keratin: one is the alpha-helixbased keratin intermediate filament proteins (IFPs), which act as the wool backbone and account for 58% of the total wool proteins; the other is keratin-associated proteins (KAPs), which serve as matrix components.These two proteins are cross-linked [3].Keratin is particularly chemically stable and possesses high mechanical strength [4] because it contains a high level of cystine and a particularly high level of disulfide bonds, which act as cross-links in the peptide chain of the protein.The physical properties of wool fibers are a very important factor in the grading and classification of wool for different textile applications [5].Wool diameter is a key factor in determining the quality and value of wool, and almost all properties of wool are related to or directly determined by fineness, such as strength, elongation, elasticity, bending stiffness and moisture absorption among other properties.Wool is a naturally curly animal fiber and its curvature is an important criterion for assessing the quality and craft value of wool.The natural curvature of wool's fibers defines it from other fibers.It is precisely because of its natural curl that wool can form many intervals of immobile space as a barrier to warmth.Therefore, in order to increase the economic value of wool, it is important to improve wool traits.
Keratins (KRTs) and Keratin-associated proteins (KRTAPs) are the most strongly expressed gene families in the hair follicle bulb [6].During hair follicle development, IFPs are expressed first, followed by the KRTAPs, and the total protein expressed by the two gene families undergoes a complex process of keratinization to form wool.The amount of IFPs in different qualities of wool is relatively stable and it is difficult to distinguish differences in IFPs from differences in the protein fractions of different wools [7]; therefore, there are fewer studies based on the relationship between IFP protein differences and gene expression.Although there are many keratin genes that regulate the same wool trait [8], they are different at the DNA level.Changes in keratin gene expression and spatio-temporal regulation determine the compositional content and cellular localization of wool keratin, thereby influencing traits such as curvature, fiber diameter, strength and length of wool.Related studies have suggested that variation in wool traits is based on changes in genes and the proteins they encode, and that changes in genes lead to phenotypic differences in wool [9].KRTs play an essential role in follicle development and in the formation of quality traits in wool [10], and most their variations correlate with wool structure and fiber properties, making them ideal candidates for improving wool traits.In this study, a penta-primer amplification refractory system (PARMS) was used for sheep genotyping to investigate the effect of keratin gene variation on wool traits in Gansu Alpine Fine-wool sheep.Genetic polymorphisms of the KRT32 (exons 1, 6 and 7) were analyzed and genotypes and haplotype combinations were analyzed for association with wool traits.The results may provide a theoretical basis for improving wool quality.

Collection of Wool, Blood and Tissue Samples from Sheep Samples
Two hundred and forty-two Gansu Alpine Fine-wool sheep ewe lambs from the same growing environment were selected for genetic variation for the KRT32.Blood was collected from the neck of the sheep and a portion of it was loaded onto an FTA card for DNA extraction using a two-step process [11].At the first shearing (12 months of age), wool samples were collected from the sheep at the position of the posterior margin of the left scapula and wool traits were measured by New Zealand Pastoral Measurements LTD (Ahuriri, Napier, New Zealand), including mean fiber diameter (MFD), fiber diameter standard deviation (FDSD), coefficient of variation of fiber diameter (CVFD); mean staple length (MSL), mean fiber curvature (MFC), mean staple strength (MSS) and comfort factor (CF).
Three newborn Gansu Alpine Fine-wool sheep ewes were selected and skin samples were collected from the posterior margin area of the shoulder joint of the left front leg on days 1, 30, 60, 90, 180, and 270 using a skin sampler with a radius of 0.44 cm; one part was placed in a freezing tube containing 4% paraformaldehyde and stored into a 4 • C ice box for immunofluorescence analysis, and the other part was placed in a freezing tube containing RNA protection solution for subsequent RT-qPCR analysis.In addition, three adult ewes in good physical condition were selected, tissue samples were taken and placed in cryopreservation tubes containing RNA-protecting solution for subsequent RT-qPCR analysis.

RT-qPCR Analysis
Total RNA was extracted from the collected tissues using Trizol reagent (Shanghai Yuanye Bio-Technology Co., Ltd., Shanghai, China) and cDNA was obtained by reverse transcription using an EvoM-MLV Reverse Transcription Kit (Accurate Biology, Wuhan, China), −20 • C ambient storage.The Beta-actin(ACTB) gene was used as a reference to calibrate the level of gene expression, and primers for the target gene were designed using Primer software (Primer Premier 5.0, Prime Corporation, Vancouver, BC, Canada) based on the gene data published in NCBI, and the primer sequences and PCR conditions are shown in Table 1.We analyzed the concentration, purity, and integrity of the RNA before doing RT-qPCR to ensure that the quality of the RNA supports the subsequent results.

Immunofluorescence Analysis
The primary antibody used in the experiment was K32 (TD9003, rabbit antibody) and the secondary antibody was CY3 (GB21303, goat anti-rabbit IgG); primary antibody was provided by Abmart PharmaTech Ltd. (Shanghai, China) and secondary antibody was provided by Servicebio Technology Co., Ltd.(Wuhan, China).Skin samples in good basic structural condition were selected, trimmed, dehydrated, paraffin-embedded, and sectioned.Paraffin sections were dewaxed using gradient biodewax and clear solution, anhydrous ethanol and distilled water.Antigen repair was performed using ethylene diamine tetra-acetic (EDTA) acid antigen retrieval buffer antigen repair solution (pH = 8).The slides were placed in PBS (pH 7.4) on a decolorizing shaker and shaken and washed 3 times for 5 min each.Sections were slightly shaken dry and circled around the tissue with a histochemical pen and sealed with a drop of 3% Bovine Serum Albumin (BSA) for 30 min.The primary antibody was added dropwise and the sections were incubated flat in a humid box at 4 • C overnight.The slides were washed in PBS (pH = 7.4) on a decolorizing shaker with shaking for 3 times, each time for 5 min, then the corresponding secondary antibody was added and incubated for 50 min at room temperature, protected from light.The slide wash was repeated under the same conditions and then the DAPI stain was dropped in and incubated for 10 min away from light.Then, we washed it again and added the tissue autofluorescence quench B solution for 5 min and rinsed with running water for 10 min.Sections were dried and sealed with an anti-fade sealant, and images were collected and analyzed.

PCR Amplification and Genotyping
According to the Ensembl database (https://www.ensembl.org(accessed on 2 March 2024)), polymorphism prediction showed that the polymorphisms of KRT32 gene were mainly concentrated in exons 1, 6 and 7, so these three exon fragments were selected as the analysis fragments.Primers for three exons were designed using the Primer software (Table 2).A total of 20 genomic DNAs of Gansu Alpine Fine-wool sheep were used to amplify three exons and all amplicons were then sequenced.The Primer synthesis, amplification, sequencing, and genotyping were performed by Gentides Biotech Co., Ltd.(Wuhan, China).Genotyping was performed using a penta-primer amplification refractory system (PARMS), a SNP PCR analysis technique that combines a pair of universal fluorescent primers, a pair of SNP allele-specific primers, and a reverse shared primer.Different primers (Table 3) bind specifically to FAM and HEX fluorescent markers in the PARMs master mix, respectively, and then fluorescence scanning is performed at the end of the reaction and the presence of alleles was determined based on the corresponding fluorescence signals.

Statistics and Analysis
After successful genotyping of KRT32, allele frequency, genotype frequency and polymorphism information content (PIC) were calculated using the formula described by Botstein et al. [12].Haploview 4.2 was used for linkage disequilibrium analysis and haplotype construction, after eliminating the results from sheep that could not be typed.One-way analysis of variance (ANOVA) was used to reveal the association between the frequency of Gansu Alpine Fine-wool sheep genotypes and wool traits.The results were output as mean ± standard error (S.E.), with p < 0.05 as the criteria for determining significant differences.The RT-qPCR results were calculated in excel using the 2 −∆∆CT method.

Expression of KRT32 mRNA in the Skin of Different Tissues and Stages of Gansu Alpine Fine-Wool Sheep
The OD260/OD280 of RNA were all between 1.8 and 2.0 and the concentrations were all greater than 100 ng/µL, which indicating that the RNA could be used for subsequent experiments.RT-qPCR results showed that KRT32 was expressed in heart, liver, spleen, lungs, kidneys and skin, and mRNA levels were significantly higher in skin tissues than in other tissues (p < 0.01).The KRT32 is expressed in the skin at all six periods of time, with expression increasing at birth, decreasing and stabilizing on days 30-90, increasing again on days 90-180, and decreasing again after day 180.In general, the expression was alternately increasing or decreasing.(p < 0.05; Figure 1).all greater than 100 ng/μL, which indicating that the RNA could be used for subsequent experiments.RT-qPCR results showed that KRT32 was expressed in heart, liver, spleen, lungs, kidneys and skin, and mRNA levels were significantly higher in skin tissues than in other tissues (p < 0.01).The KRT32 is expressed in the skin at all six periods of time, with expression increasing at birth, decreasing and stabilizing on days 30-90, increasing again on days 90-180, and decreasing again after day 180.In general, the expression was alternately increasing or decreasing.(p < 0.05; Figure 1).

Expression Localization and Distribution Density of KRT32 Encoding Protein in SKIN tissues of Gansu Alpine Fine-Wool Sheep of Various Ages
The expression and distribution of the KRT32-encoding protein was analyzed using immunofluorescence combined with the positive scoring of different parts of the hair follicle.Distribution of K32 in skin tissues of Gansu Alpine Fine-wool sheep of different ages was analyzed using immunofluorescence.The image on the right shows the H&E staining of the skin histology of Gansu alpine fine wool sheep to facilitate visualization of the entire hair follicle structure (Figure 2).The results showed that the cuticle, IRS, ORS, HM and DP were strongly positively expressed in the skin tissues of Gansu Alpine Fine-wool sheep at 1 day old.At 30, 60, 90, 180 and 270 days of age, skin tissue were high positive expression in the IRS, ORS, HM, sebaceous glands and DP (Table 4).

Expression Localization and Distribution Density of KRT32 Encoding Protein in SKIN Tissues of Gansu Alpine Fine-Wool Sheep of Various Ages
The expression and distribution of the KRT32-encoding protein was analyzed using immunofluorescence combined with the positive scoring of different parts of the hair follicle.Distribution of K32 in skin tissues of Gansu Alpine Fine-wool sheep of different ages was analyzed using immunofluorescence.The image on the right shows the H&E staining of the skin histology of Gansu alpine fine wool sheep to facilitate visualization of the entire hair follicle structure (Figure 2).The results showed that the cuticle, IRS, ORS, HM and DP were strongly positively expressed in the skin tissues of Gansu Alpine Fine-wool sheep at 1 day old.At 30, 60, 90, 180 and 270 days of age, skin tissue were high positive expression in the IRS, ORS, HM, sebaceous glands and DP (Table 4).

Genotyping and Polymorphism Analysis
The KRT32 exons 1, 6 and 7 were sequenced in 20 randomly selected DNA samples of Gansu Alpine Fine-wool sheep.Eight SNPs were detected: SNP1 (g.21455859), SNP2 (g.21455953), SNP3 (g.21455976), SNP4 (g.21456106), SNP5 (g.21460798), SNP6 (g.21460884), SNP7 (g.21463485) and SNP8 (g.21463503) (Figure 3).Of the 8 SNPs, three of them located in non-coding regions (SNP1, SNP7, SNP8) and five in coding regions (SNP2, SNP3, SNP4, SNP5 SNP6), of which three were synonymous (SNP2, SNP4, SNP6) and two were missense mutations (SNP3, SNP5).The GTG codon to GCG codon change in SNP3 results in the change of valine to alanine; hence, p.Val17Ala, and the GAC codon to AAC codon change in SNP5 results in the change of aspartic acid to asparagine, defined as p.Asp17Asn.The length of the DNAmarker for electrophoresis was 2000 bp.The PCR analysis produced a fragment of 670 bp, 612 bp and 468 bp, respectively, which were of the expected size (Figure 3).Genotyping using the PARMS technique revealed all the eight SNPs showed three genotypes; heterozygotes are located in the middle of the scatterplot and are indicated by red dots, while pure heterozygotes are located at the edges of the scatterplot, close to the x-axis and y-axis, and are indicated by blue and green dots, respectively.Grey dots are generally wells with no DNA samples or wells with no signal for PCR amplification (Figure 4).The horizontal coordinate is the FAM fluorescence value, and the vertical coordinate is the HEX fluorescence value.The dominant alleles of the eight SNPs in the KRT32 were C, A, C, A, A, G, T and C, corresponding to the genotypes with the highest frequencies as CC (0.402), GA (0.480), TC (0.452), GA (0.480), GA (0.452), AG (0.402), TT (0.712), CC (0.887).Population genetic analysis of eight positions in Gansu Alpine Fine-wool sheep revealed that all positions were moderately polymorphic (0.25 < PIC < 0.5), except SNP8 which was lowly polymorphic (PIC < 0.25) (Table 5).Genotyping using the PARMS technique revealed all the eight SNPs showed three genotypes; heterozygotes are located in the middle of the scatterplot and are indicated by red dots, while pure heterozygotes are located at the edges of the scatterplot, close to the x-axis and y-axis, and are indicated by blue and green dots, respectively.Grey dots are generally wells with no DNA samples or wells with no signal for PCR amplification (Figure 4).The horizontal coordinate is the FAM fluorescence value, and the vertical coordinate is the HEX fluorescence value.The dominant alleles of the eight SNPs in the KRT32 were C, A, C, A, A, G, T and C, corresponding to the genotypes with the highest frequencies as CC (0.402), GA (0.480), TC (0.452), GA (0.480), GA (0.452), AG (0.402), TT (0.712), CC (0.887).Population genetic analysis of eight positions in Gansu Alpine Finewool sheep revealed that all positions were moderately polymorphic (0.25 < PIC < 0.5), except SNP8 which was lowly polymorphic (PIC < 0.25) (Table 5).

Analysis of the Association of SNPs with Wool Traits
Eight SNPs for the KRT32 was analyzed for association with wool traits (Table 6).SNP1, SNP 3 and SNP 5 were significantly correlated with the coefficient of CVFD and MFC (p < 0.05); SNP 2, SNP 4 and SNP 6 were significantly correlated with MFC (p < 0.05); SNP7 was significantly correlated with the CVFD (p < 0.05); SNP8 was significantly associated with MFD and CF (p < 0.05).MFC was significantly higher in individuals with TT genotype than CC genotype on SNP1 (p < 0.05); MFC was significantly higher in individuals with GG genotype than AA genotype on SNP2 (p < 0.05); MFC was significantly higher in individuals with TT and TC genotypes than CC genotypes on SNP3 (p < 0.05);

Analysis of the Association of SNPs with Wool Traits
Eight SNPs for the KRT32 was analyzed for association with wool traits (Table 6).SNP1, SNP 3 and SNP 5 were significantly correlated with the coefficient of CVFD and MFC (p < 0.05); SNP 2, SNP 4 and SNP 6 were significantly correlated with MFC (p < 0.05); SNP7 was significantly correlated with the CVFD (p < 0.05); SNP8 was significantly associated with MFD and CF (p < 0.05).MFC was significantly higher in individuals with TT genotype than CC genotype on SNP1 (p < 0.05); MFC was significantly higher in individuals with GG genotype than AA genotype on SNP2 (p < 0.05); MFC was significantly higher in individuals with TT and TC genotypes than CC genotypes on SNP3 (p < 0.05); CVFD was significantly higher in individuals with a TC genotype than CC genotype, and other MFC was significantly higher in GG and GA genotypes than in AA genotypes on SNP4 (p < 0.05); MFC was significantly higher in individuals with GG and GA genotypes than in AA genotypes on SNP5, and CVFD was significantly higher in individuals with GA genotypes than in AA genotypes on SNP6 (p < 0.05); MFC was significantly higher in AG and AA genotypes than in GG genotypes on SNP6.In SNP7, the CVFD of TT genotypes was significantly higher than that of TC genotypes; in SNP8, the MFD of CG genotypes was significantly higher than that of GG genotypes, and the CF of CC genotypes was significantly higher than that of CG genotypes (p < 0.05).

Haplotype Reconstruction of Gene SNPs and Association Analysis of Haplotype Combinations with Wool Traits
Eight SNPs of the KRT32 were screened and analyzed for linkage disequilibrium and haplotype using the Haploview 4.2 software.The results showed that among the eight SNPs in KRT32, SNP1 and SNP8 were not involved in haplotype formation due to their weak association with other mutation sites, and the remaining six SNPs were screened to form two haplotype blocks, that is, Block1 (g.21455953 G>A, g.214555976 T>C, g.21456106 G>A and g.21460798 G>A) and Block2 (g.21460884 A>G and g.21463485 T>C), respectively, and both of them were in the strong interlocking state (Figure 5).In the Gansu Alpine Finewool sheep population, four haplotypes were constructed in Block1, and these haplotypes were combined to form six haplotype combinations with frequencies greater than 0.01 (Table 7).Three haplotypes were constructed in Block2, and combining these haplotypes together similarly resulted in six haplotype combinations with frequencies greater than 0.01 (Table 8).
Curr.Issues Mol.Biol.2024, 46, FOR PEER REVIEW 10 (SNP2, 3, 4 and 5) of the KRT32 were significantly correlated with MFC; different haplotype combinations of 2 SNPs (SNP6 and SNP7) of the KRT32 were significantly correlated with MFD, MFC and CF.The MFC of H2H2 and H2H3 individuals in Block1 of the KRT32 was significantly higher than that of H3H3 individuals (Table 9); In Block2, MFD was significantly higher in H1H3 individuals than in H2H3 individuals, while the CF was significantly lower in H1H3 individuals than in H2H3 individuals, MFC was lower in H1H3 individuals than in H1H2, H2H2, H2H3 and H3H3 individuals (Table 10).Associations with wool traits were analyzed for each of the six haplotype combinations of Block 1.The results showed that different haplotype combinations of four SNPs (SNP2, 3, 4 and 5) of the KRT32 were significantly correlated with MFC; different haplotype combinations of 2 SNPs (SNP6 and SNP7) of the KRT32 were significantly correlated with MFD, MFC and CF.The MFC of H 2 H 2 and H 2 H 3 individuals in Block1 of the KRT32 was significantly higher than that of H 3 H 3 individuals (Table 9

Discussion
Studies on the use of polymerase chain reaction-single-stranded conformational polymorphism (PCR-SSCP) to identify the effect of SNPs in KRT83 [13] and KRT31 [14] on wool traits have demonstrated nucleotide sequence variation in the sheep keratin genes.Analysis of KRT32 in Ensembl (version 110, 2023) showed that the gene is located on chromosome 11, and there is a sequence in the Rambouillet v1.0 genome construct CM008482.1 (a fine wool sheep breed), identified as ENSOARG00020008724 (GenBank GENE ID: 100526788), which is close to another keratin gene, KRT35 (ENSOART00020008687).The ATG start site of the KRT32 sequence is located in exon 1, with a predicted length of 401 amino acids.
Wool fibers consist mainly of keratin, encoded by the KRT and KRTAP genes, which are expressed in a highly regulated manner during hair follicle growth, as reflected by their ordered spatiotemporal distribution of mRNA expression.This study examined the expression levels of KRT32 mRNA in different other organs and in skin at different periods using RT-qPCR.The results showed that the expression level of KRT32 mRNA in skin was significantly higher than that in parenchyma organs, and it was expressed in skin at different developmental periods, with significant differences in the expression levels at various ages.It's related to the hair cycle.The expression of KRT32 increased to a peak during the anagen phase and then decreased to a steady state.During the catagen phase, the expression of KRT32 increased to a peak and then decreased, therefore, so it was inferred that KRT32 may regulate the maturation of hair follicles and wool growth mainly during the anagen and catagen phases.
Hair follicle growth is a complex cyclical physiological process influenced by a variety of factors.The feasibility of using fluorescent antibodies techniques to localise keratin in hair follicles was established as α-keratin has been shown to be antigenic [15][16][17].In this study, immunofluorescence was used to determine the location of protein expression, and the results showed that KRT32 mainly expressed in the DP, IRS, ORS and HM, with relatively weak expression in the cuticle, which is basically in line with the results of previous studies on the localization of mRNA expression [18].Positive scores showed that the protein encoded by the KRT32 is distributed at a relatively uniform density in all structures of the hair follicle except the stratum corneum, and since DP cells and ORS cells are the major cell types that regulate hair growth and regeneration [19], thus the KRT32 may have an important effect on follicle maturation and wool growth.
Chen et al. [20] elucidated that KRT75 plays a critical role in maintaining the hair shaft by constructing a mouse model.Chai et al. [21] revealed that variation in the promoter fragment of KRT34 affects MFD, FDSD, and MSL, and Bao et al. [10], in an experiment using yak as a test subject, discovered that KRTs showed a strong correlation throughout the hair follicle development cycle, suggesting a generalised effect of keratin on the hair growth process in different species.In this study it was found that KRT32 was significantly correlated with MFD, MFC, both of which are major traits determining the value of wool.We confirmed the genetic poly-morphism of the KRT32 in sheep, and a total of eight SNPs were screened in the three exons.Three of the variants occurred in the non-coding region, including one in the 5 ′ -UTR and two in the 3 ′ -UTR; and other five variants were identified in the coding region, including three synonymous and two missense mutations.SNPs in the coding region may alter the amino acid, affecting the structure and function of the corresponding proteins [22], ultimately leading to changes in the wool phenotype.In this study, we found that SNP3 and SNP5 in the coding region were missense mutations, the conversion of GTG codon to GCG codon in SNP3 resulted in the change of valine to alanine, thus p.Val17Ala, and the conversion of GAC codon to AAC codon in SNP5 resulted in the change of aspartic acid to asparagine, that is, p.Asp317Asn, and both SNPs were significantly associated with CVFD and MFC were significantly correlated, which may be due to the change in amino acid leading to the change in protein secondary structure, which in turn leads to the change in phenotype.SNP2, SNP4 and SNP6 were all located in the coding region and were synonymous mutations, but their variants are also significantly associated with wool traits, a result that is similar to the results of some previous studies [23][24][25].Theoretically, synonymous mutations in coding regions usually do not cause alterations in the function of coded proteins, but the results of studies have shown that synonymous mutations are not completely silent, and some of them can affect the formation or loss of variable spliceosomes in genes, mRNA stability, and translation efficiency [26,27].Some other synonymous mutations can affect protein folding and conformation, leading to alterations in gene function and phenotype [28].It is important to note that the analysis of the association results between genotypes and wool traits showed that SNPs in non-coding regions (SNP1, SNP7 and SNP8) were also significantly associated with wool traits, which may be due to the fact that SNP in non-coding regions can influence the binding efficiency of transcription factors [29], affect mRNA shearing or alter DNA methylation modifications, which are involved in the process of gene expression regulation and phenotype generation, and that the SNP in the 3 ′ UTR region of genes affects gene expression in multiple ways [30][31][32].
Haplotype analysis takes into account interactions between non-alleles and linkage disequilibrium between SNP, and is more comprehensive than analysis of individual SNPs [33].In this study, two blocks were constructed and the results of the association analysis of the haplotype combinations haplotype combinations of the two blocks with wool traits showed that the haplotype combinations H 2 H 3 and H 3 H 3 of Block1 were significantly associated with wool MFC, and that MFC was significantly higher in the haplotype combinations H 2 H 2 and H 2 H 3 populations than in the H 3 H 3 population, which suggests that the presence of the haplotype H2 in Block1 led to an increase in MFC; the haplotype combination H 2 H 3 of Block 2 was significantly associated with the MFD, CF and MFC, MFD was significantly higher in haplotype combination H 1 H 3 than in H 2 H 3 , while CF was significantly lower in diplotype H 1 H 3 than in H 2 H 3 .It should be noted that there is a correlation between fiber diameter and comfort factor, which is in agreement with the results of a previous study [34].Since the H 2 H 3 in Block2 increases MFC while decreasing MFD, which can substantially improve wool quality, and its frequency of occurrence is high in the population, haplotype combination H 2 H 3 can be used as a selection target.The results of the association analyses of haplotype combinations with wool traits were generally consistent with the results of the single locus association analyses, both of which significantly affected wool traits.The effects of KRT32 on wool traits suggest that this gene has a significant impact on the development of the hair follicle, and is an important gene in influencing wool phenotypes.It is important to note that we need to go further to understand the mechanism of the influence of other genes on wool traits as well as the interactions between the genes, to gain a deeper understanding of the formation of different types of wool, and better select and breed to meet human needs.

Conclusions
This study demonstrated for the first time the potential of the KRT32 gene as a genetic marker for wool traits.KRT32 was highly expressed in the skin and showed significant differences in expression at different stages of the skin, encoding proteins localized to the inner root sheath (IRS), outer root sheath (ORS) and dermal papilla (DP), and all eight SNPs in KRT32 were significantly correlated with wool traits, including MFD, CF and MFC.Constructed haplotype after haplotype combination H 2 H 3 has higher CF and smaller MFD.The combination of the above indicates that the sequence variation of KRT32 is significantly associated with economically important traits of wool (e.g., MFD and MFC), and it can be used as a candidate gene for wool improvement in Gansu Alpine Fine-wool sheep, which can provide a reference value for the improvement of wool quality.

Figure 1 .
Figure 1.(A,B) indicate the expression results of KRT32 mRNA in different tissues and different skin stages in Gansu Alpine Fine-wool sheep, respectively.The data are expressed as the mean ± S.E.Different letters indicate significant differences between different ages (p < 0.05).

Figure 1 .
Figure 1.(A,B) indicate the expression results of KRT32 mRNA in different tissues and different skin stages in Gansu Alpine Fine-wool sheep, respectively.The data are expressed as the mean ± S.E.Different letters indicate significant differences between different ages (p < 0.05).

Figure 2 .
Figure 2. Immunofluorescence staining for K32 in skin tissues of Gansu Alpine Fine-wool sheep at different times (5×).Blue (B,E,H,K,N,Q) and red (C,F,I,L,O,R) tissue in the figure show DAPI-labeled nuclear fluorescence staining and fluorescence staining for K32, respectively.ORS: outer root sheath; IRS: inner root sheath; HM: hair matrix; C: corneum; SG: sebaceous glands; DP: dermal papilla.H&E staining image are shown on the right.The first column (A,D,G,J,M,P) is a combined graph of the second and third columns

Figure 2 .
Figure 2. Immunofluorescence staining for K32 in skin tissues of Gansu Alpine Fine-wool sheep at different times (5×).Blue (B,E,H,K,N,Q) and red (C,F,I,L,O,R) tissue in the figure show DAPI-labeled nuclear fluorescence staining and fluorescence staining for K32, respectively.ORS: outer root sheath; IRS: inner root sheath; HM: hair matrix; C: corneum; SG: sebaceous glands; DP: dermal papilla.H&E staining image are shown on the right.The first column (A,D,G,J,M,P) is a combined graph of the second and third columns.
Curr.Issues Mol.Biol.2024, 46, FOR PEERREVIEW  7    results in the change of valine to alanine; hence, p.Val17Ala, and the GAC codon to AAC codon change in SNP5 results in the change of aspartic acid to asparagine, defined as p.Asp17Asn.The length of the DNAmarker for electrophoresis was 2000 bp.The PCR analysis produced a fragment of 670 bp, 612 bp and 468 bp, respectively, which were of the expected size (Figure3).

Figure 3 .
Figure 3. PCR amplification of exons 1, 6 and 7 of KRT32 and sequencing results, overlapping peaks indicate SNPs.The overlapping peaks in the direction of the arrows refer to mutation sites, and the different colored lines represent different bases.

Figure 3 .
Figure 3. PCR amplification of exons 1, 6 and 7 of KRT32 and sequencing results, overlapping peaks indicate SNPs.The overlapping peaks in the direction of the arrows refer to mutation sites, and the different colored lines represent different bases.

8 Figure 4 .
Figure 4. Genotyping results for eight positions of the KRT32 (A-H) gene in Gansu Alpine Finewool sheep.Different colors represent different genotypes, with green and blue representing heterozygotes and red representing purities

Figure 4 .
Figure 4. Genotyping results for eight positions of the KRT32 (A-H) gene in Gansu Alpine Fine-wool sheep.Different colors represent different genotypes, with green and blue representing heterozygotes and red representing purities.

Figure 5 .
Figure 5. Linkage disequilibrium analysis of SNPs in KRT32.The color of the square indicates the degree of interlocking, and the darker the color, the higher the degree of interlocking; The value represents the strength (percentage) of the correlation value between SNPs.

Figure 5 .
Figure 5. Linkage disequilibrium analysis of SNPs in KRT32.The color of the square indicates the degree of interlocking, and the darker the color, the higher the degree of interlocking; The value represents the strength (percentage) of the correlation value between SNPs.

Table 1 .
Gene primers and internal reference sequence information.

Table 2 .
Primer information for the three exons.

Table 3 .
PARMs-PCR primer sequences and fluorescent signal types.

Table 4 .
Distribution density of K32 in various parts of skin tissues at various ages.

Table 4 .
Distribution density of K32 in various parts of skin tissues at various ages.

Table 5 .
Genotype frequency and gene frequency of 8 SNPs of KRT32

Table 5 .
Genotype frequency and gene frequency of 8 SNPs of KRT32.

Table 6 .
Association analysis of KRT32 SNPs with wool traits.

Table 7 .
Haplotypes and haplotype combinations of block1 and their respective frequencies

Table 8 .
Association analysis of wool traits with haplotype combinations of Block1 MFD, mean fiber diameter; FDSD, fiber diameter standard deviation; CVFD, coefficient of variation of fiber diameter; CF, comfort factor; MSL, mean staple length; MSS, mean staple strength, MFC, mean fiber curvature.a,b Values with different superscripts within the same column differ significantly at p < 0.05.

Table 9 .
Haplotypes and haplotype combinations of block2 and their respective frequencies

Table 7 .
Haplotypes and haplotype combinations of block1 and their respective frequencies.

Table 8 .
Association analysis of wool traits with haplotype combinations of Block1.
diameter; CF, comfort factor; MSL, mean staple length; MSS, mean staple strength, MFC, mean fiber curvature.a,b Values with different superscripts within the same column differ significantly at p < 0.05.
); In Block2, MFD was significantly higher in H 1 H 3 individuals than in H 2 H 3 individuals, while the CF was significantly lower in H 1 H 3 individuals than in H 2 H 3 individuals, MFC was lower in H 1 H 3 individuals than in H 1 H 2 , H 2 H 2 , H 2 H 3 and H 3 H 3 individuals (Table 10).

Table 9 .
Haplotypes and haplotype combinations of block2 and their respective frequencies.

Table 10 .
Association analysis of wool traits with haplotype combinations of Block2.FDSD, fiber diameter standard deviation; CVFD, coefficient of variation of fiber diameter; CF, comfort factor; MSL, mean staple length; MSS, mean staple strength, MFC, mean fiber curvature.a,b Values with different superscripts within the same column differ significantly at p < 0.05.