Application of detection reagent of molecular marker related to wool length character of down producing goat

Abstract

The invention belongs to the technical field of animal breeding, and particularly relates to application of a detection reagent of a molecular marker related to the hair length character of a down producing goat. The molecular marker is located at 155, 344, 701-155, 346 and 600 positions on a No.1 chromosome of a goat genome, and the nucleotide sequence of the molecular marker is as shown in SEQ ID NO. 1. It is found that the wool length of an individual with a normal molecular marker area is larger than that of an individual with a deletion molecular marker area, and therefore the cashmere goat variety with the large wool length can be bred by detecting deletion mutation of the molecular marker on the chromosome 1 of the goat genome.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of animal breeding, and particularly relates to application of a detection reagent of a molecular marker related to a wool length trait of a cashmere goat. BACKGROUND

[0002] Inner Mongolia cashmere goat is an excellent local livestock breed formed through long-term natural selection and artificial systematic breeding. Copy number variation (CNV) is widely present in the genome, including deletion and duplication. With the deepening understanding of complex animal and plant genomes, more and more evidence shows that CNVs play an important role in determining breeding and agricultural important traits. Although CNV is an important source of genetic diversity in animals, the influence of CNV on the genome and wool traits of goats is still largely unknown. Genome-wide association study (GWAS) is a powerful tool for trait-related gene mapping and has been widely used to identify candidate genomic loci affecting phenotypic variation. However, the research of GWAS usually relies on the information of single nucleotide polymorphism (SNP) associated with phenotypes, while ignoring large fragment CNVs with larger phenotype effect size. Through CNV detection and CNV-GWAS analysis of Inner Mongolia cashmere goats, new insights can be provided for the genetic breeding of Inner Mongolia cashmere goats. It is found that with the increase of wool length, the wool yield and body weight also increase. Selecting long wool type individuals for breeding can accelerate the genetic progress of wool yield and body weight of cashmere goats, realize the indirect selection of wool on wool yield and body weight, and facilitate processing in subsequent application. Therefore, mining molecular markers affecting the wool length of Inner Mongolia cashmere goats at the genome level is conducive to breeding cashmere goat breeds with larger wool yield and body weight. SUMMARY

[0003] In view of the above technical problems, the application provides application of a detection reagent of a molecular marker related to a wool length trait of a cashmere goat. The molecular marker can be used to breed cashmere goats with long wool length and indirectly improve the wool yield and body weight of cashmere goats.

[0004] In a first aspect of the application, application of a detection reagent of a molecular marker related to a wool length trait of a cashmere goat in identifying the wool length trait of the cashmere goat is provided. The nucleotide sequence of the molecular marker is shown in SEQ ID NO. 1.

[0005] As a preferred embodiment of the application, the method for identifying the wool length trait of the cashmere goat comprises the following steps: Detecting deletion or insertion of the molecular marker on chromosome 1 of the genome of the goat to be tested. The wool length of an individual with insertion is greater than that of an individual with deletion.

[0006] In the second aspect of the present application, the application of a detection reagent of a molecular marker related to the hair length trait of cashmere goats in the genetic breeding of cashmere goats is provided, and the nucleotide sequence of the molecular marker is shown as SEQ ID NO. 1.

[0007] As a preferred embodiment of the present application, the genetic breeding is to select cashmere goats with long hair length.

[0008] Further preferably, the method for selecting cashmere goats with long hair length comprises the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome is detected, and an individual with insertion is selected as a parent to increase the hair length of the offspring cashmere goats.

[0009] As a preferred embodiment of the present application, the genetic breeding is to select cashmere goats with large wool yield or large body weight.

[0010] Further preferably, the method for selecting cashmere goats with large wool yield comprises the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome is detected, and an individual with insertion is selected as a parent to increase the wool yield of the offspring cashmere goats.

[0011] Further preferably, the method for selecting cashmere goats with large body weight comprises the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome is detected, and an individual with insertion is selected as a parent to increase the body weight of the offspring cashmere goats.

[0012] The goat genome of the present application is the reference genome ARS1 of goats.

[0013] The cashmere goat of the present application is Inner Mongolia cashmere goat.

[0014] Compared with the prior art, the present application has the following beneficial effects: The present application provides a molecular marker related to the hair length trait of cashmere goats, and it is found that by detecting the deletion mutation of the molecular marker on chromosome 1 of the goat genome, cashmere goat breeds with long hair length can be selected, specifically: by selecting an individual with insertion of the molecular marker on chromosome 1 of the goat genome as a father or a mother, the hair length of the offspring can be increased, and the wool yield and body weight of the offspring can be indirectly increased. BRIEF DESCRIPTION OF DRAWINGS

[0015] Figure 1 is the feature of CNV in IMCGs; A, the number distribution of CNV on 29 autosomes, B, the correlation between chromosome length and CNV number, C, the number distribution of CNV in different length intervals.

[0016] Figure 2are descriptive statistics of CNVRs in IMCGs; A, the proportion of CNVRs in three categories, B, the distribution of CNVRs on different chromosomes, C, the distribution of CNVRs in different length intervals, D, the distribution of CNVRs on autosomes based on 1 Mb windows.

[0017] Figure 3 are PCA visualizations; X and Y axes represent the first two explained variance percentages, respectively.

[0018] Figure 4 are TDRD3 qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0019] Figure 5 are KRT222 qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0020] Figure 6 are MDGA1 qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0021] Figure 7 are NOX4 qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0022] Figure 8 are DYRK1A qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0023] Figure 9 are SATB1 qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0024] Figure 10 are SH2D5 qPCR validation of important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents the relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0025] Figure 11 areMMADHC qPCR verification of top important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents relative copy number. Relative copy number < 2 represents that the CNV is a deletion type.

[0026] Figure 12 is TPK1 qPCR verification of top important candidate CNVs. The horizontal axis represents 20 individual Cashmere goats, and the vertical axis represents relative copy number. Relative copy number < 2 represents that the CNV is a deletion type. DETAILED DESCRIPTION

[0027] It should be noted that the following detailed description is illustrative only, and is intended to provide further description in order to provide a fuller enabling teaching of the exemplary embodiments according to the present application. Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application pertains.

[0028] It is also important to note that the terms "comprises", "comprising", "includes", "including", "contains", "containing" or variations thereof herein, do not specify an exhaustive or complete list of components or features as used in describing processes or devices. Unless otherwise expressed, the terms "comprises", "comprising", "includes", "including", "contains", "containing" or variations thereof are intended to be inclusive or open-ended and not exclude additional, unrecited components, features, steps, processes, devices, compositions or combinations thereof.

[0029] The application will be further described with reference to the following examples, which are intended to be purely exemplary of the application, and not limiting. The test samples and test procedures used in the following examples include the following (if the specific experimental conditions are not specified in the examples, the general conditions are usually used, or the conditions recommended by the reagent company; the reagents, consumables, etc. used in the following examples, if not specifically stated, can be obtained from commercial channels.

[0030] 1. Materials and methods 1.1. Experimental animals and samples The Inner Mongolian cashmere goats used in the present application are all from Erlangshan Ranch of Inner Mongolia Beiping Textile Co., Ltd. The phenotype data collection method used is consistent with the previous research published by the present application team, see "Rong Y, Wang X, Na Q, et al. Genome-wide association study for cashmere traits in Inner Mongolian cashmere goat population reveals new candidate genes and haplotypes. BMC genomics 2024; 25:658." Specifically, the cashmere trait phenotype data of 3842 individuals were collected between 2020 and 2023, including 1401 males and 2441 females, and the pedigree can be traced back to two generations. The phenotype records include 6966 cashmere yield, 6814 cashmere fineness, 5830 cashmere thickness and 5831 cashmere length. Among the 404 individuals used for resequencing, 84 were rams, of which 75 were one year old, 4 were two years old, and 5 were three years old. 320 were ewes, of which 92 were one year old, 54 were two years old, and 174 were three years old. For the correction of the phenotype value, the present application uses the method of "breeding value + residual", and the descriptive statistics of the corrected phenotype value are shown in Table 1.

[0031] Table 1 Descriptive statistics of four cashmere traits

[0032] 1.2, Quality control of whole genome sequencing data and sequence alignment Genomic DNA was extracted from ear tissue samples using the phenol-chloroform method, and the concentration and purity of the DNA sample were determined by ultraviolet spectrophotometry, A 260nm / A 280nmThe ratio was 1.8-2.0. Subsequently, high-quality DNA samples were submitted to Shijiazhuang Boruit Biotechnology Co., Ltd. for 25 X deep whole-genome resequencing using the DNBSEQ-T7 platform. The original reads were filtered using V0.20.0 fastp software to remove adapter sequences, remove paired reads with N content exceeding 10 bases, and remove paired reads with Q≤20 low-quality base number exceeding 40% of the length of the read. The reference genome was indexed, and the quality-controlled clean reads were aligned with the goat reference genome ARS1, GCF_001704415.1 using V0.7.17 Burrows-Wheeler Aligner software. The aligned sam file was converted into a bam file using V1.20 Samtools software, and the bam file was sorted. The sorted bam file was removed from the duplicate data to obtain the final bam file using the MarkDuplicates program in V3.8 Genome Analysis Toolkit software.

[0033] 1.3, Whole genome CNV detection CNV calling of Inner Mongolia Cashmere Goat population by V1.3.2 CNVpytor, which is a read-depth-based CNV detection method. Referring to the literature record "Abyzov A, Urban AE, Snyder M, Gerstein M. CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. Genome research 2011;21:974-84.", the bin size of all individuals is set to 100bp. In order to improve the accuracy of CNV detection, the following standards are used to filter false positive CNV calls: P-value <0.01, fraction of reads mapped with q0 quality in call segments >0.5, size of CNV >1 kb, fraction of reference genome gaps within call region >5%, calls close to gaps in reference genome <100 kb. This produces high-quality CNVs suitable for subsequent analysis. Among them, CNVs of the duplication type refer to regions that include more than two copy numbers, and CNVs of the deletion type represent regions that only contain 0 or 1 copy number. According to the following principles, a Python script is written: the calculation method of CNV length is termination position-starting position+1, and the length interval is set as: 1 Kb~2 Kb, 2 Kb~5 Kb, 5 Kb~10 Kb, 10 Kb~20 Kb, 20 Kb~50 Kb, 50 Kb~100 Kb, 100 Kb~300 Kb. The present application uses a self-programmed Python script to count the length and number of all CNVs on each chromosome. In addition, the present application respectively counts the length and number of deletion and duplication types. The Graphpad Prism 8.0.2 software is used to draw the distribution diagram of CNVs in different chromosomes and different length intervals, and the fitting diagram of the number of CNVs and the length of chromosomes.

[0034] 1.4, CNVR determination To obtain CNVR, the present application wrote a Python script according to the following principle: considering two CNVs, CNV1 starts at position A and ends at position C, CNV2 starts at position B and ends at position D, A < B < C < D, if the mutual overlap between the two CNVs is at least 1 bp, then they are merged into a CNVR on the genome from X to W. Using this script, two or more CNVs that overlap at least 1 bp are merged into one CNVR, and in order to reduce the false positive rate, only CNVRs found in three or more samples are considered, thereby generating a final CNVR list. This list contains the following columns: CNVR_ID, Chr, Start, End, n_Sample, Type, Length, n_CNV. CNVR_ID is the generated CNVR number; Chr is the chromosome number; Start is the start position of the generated CNVR; End is the end position of the generated CNVR; n_Sample is the number of samples involved in this CNVR; Type is the type of CNVR, CNVR is divided into two types of "Deletion" and "Duplication", overlapping "Deletion" and "Duplication" type CNVRs are merged into a single region, called "Mixed" CNVR; Length is the length of the generated CNVR; n_CNV is how many CNVs are included in this CNVR. In addition, the present application uses a self-programmed Python script to count the length and number of all CNVRs on each chromosome and the length and number of each of the three types. The length interval of CNVR is set as: 1Kb~2 Kb, 2Kb~5 Kb, 5Kb~10 Kb, 10Kb~20 Kb, 20Kb~50 Kb, 50Kb~2000 Kb. The total length of the goat ARS1 reference genome autosomes is about 2,466,191,353 bp, and the present application calculates the proportion of all CNVRs in the whole genome. Statistical analysis was performed on the three types of CNVRs, and Graphpad Prism 8.0.2 software was used to visualize the number of CNVRs in different chromosomes and different length intervals. The distribution of CNVRs in 29 pairs of autosomes was visualized using the R package CMplot.

[0035] 1.5, PCA analysis Based on the obtained whole genome CNVs data, the present application uses the " --pca 3" parameter of V1.90 Plink software to calculate the first three principal components PC, and uses V4.2.2 R to draw the principal component analysis PCA graph.

[0036] 1.6, whole genome association analysis CNV association analysis with CY, CD, CT, FL traits was performed using the fastGWA-mlm model in V1.94.1 GCTA software. The software will automatically filter CNVs with MAF < 0.0001 and missingness rate > 10%. As follows: ; where y is the n x 1 phenotype vector; X cnv is the genotype vector with effect β cnv ; g is the vector of total genetic effects captured by CNV-derived genetic relationship matrix (GRM), g ~ N(0, ); π is the CNV-derived GRM vector with all off-diagonal elements set to 0; e is the residual vector, e ~ N(0, ).

[0037] Since the Bonferroni correction method is too strict to determine the significance threshold of GWAS, the present application adjusts the threshold of genome-wide significant association to P = 1 x 10 -3 The genome inflation factor λ of the test statistic is calculated by the slope of the linear regression between the observed quantiles and the theoretical quantiles in V4.2.2 R. Significant SNPs are displayed with the threshold line in the Manhattan plot. The Manhattan plot and quantile-quantile (QQ) plot are generated by the CMplot package in R.

[0038] 1.7, Functional enrichment analysis of candidate genes overlapping with CNV To determine the candidate genes associated with the wool traits, the present application uses V2.30.0 Bedtools software to annotate the significantly associated CNVs based on the goat ARS1 version genome. The CNVs are located on the identified genes, only those genes located within 1 Mb ± 500 kb around each CNV are considered. In addition, the candidate genes are subjected to GO functional annotation and KEGG enrichment analysis using the clusterProfiler package in R. The KEGG pathways and GO biological processes of the candidate genes are manually searched to infer the potential gene functions.

[0039] 1.8, Real-time fluorescent quantitative PCR verification To verify the accuracy of CNVs significantly associated with wool traits, the present application used real-time quantitative PCR to verify 10 randomly selected deletion-type CNVs. These CNVs were verified by qPCR on 20 Inner Mongolia Cashmere Goat samples. After comparing the reference genome ARS1 and determining the sequences of the 10 CNVs, primers were designed using Primer premier 5.0 software, synthesized by Beijing Sinonow Gold Biotech Co., Ltd., and qPCR reactions were performed using TB Green ® Premix Ex Taq™ II kit. See Table 2 for detailed information on all CNVs and primers. The genomic DNA of the sample to be tested was used as the template for qPCR, and β-actin, which did not undergo copy number variation, was selected as the internal reference gene. The CNVs to be verified and the internal reference gene were each repeated three times. The formula for calculating the relative copy number is: ; where 2 indicates that a normal individual has 2 copies of a certain gene; 2 -ΔΔCt represents the fold change of a certain DNA fragment between a normal individual and a variant individual, and the base number 2 refers to the fold change in product quantity during PCR amplification, ΔCt = Ct value of the target fragment - Ct value of the corresponding internal reference gene, and ΔΔCt = ΔCt value of the sample with copy number variation - ΔCt value of the control sample. Finally, the calculated results are compared with the normal value 2, and values close to 2 are considered normal, while values of 3 or more indicate copy number duplication, and values of 1 or less indicate copy number deletion.

[0040] Table 2 List of qPCR verification primer information for CNVs

[0041] 2、Data analysis The association between the genotypes of 31 SNPs and wool traits was evaluated using independent sample t-test in V25.0 SPSS statistical software. When the variance analysis results of a certain group of genotypes showed significant differences, the statistical differences between two genotypes were evaluated using Bonferroni correction, and the data were expressed as "mean ± standard error".

[0042] 3、Results 3.1、Identification of copy number variation CNVs, and 400,578 high-quality CNVs were retained from the CNVs, with an average of 991.5 per animal. These CNVs included 193,741 deletions, accounting for 48.37%, and 206,837 duplications, accounting for 51.63%. The length of all CNVs ranged from 1000 bp to 255,600 bp, with an average length of 6,351.66 bp, of which the average length of deletions was 4,718.59 bp, and the average length of duplications was 7,881.34 bp, as shown in Table 3. Among the 29 pairs of autosomes, the most deletions were on chromosome 12, with 11,905, and the most duplications were on chromosome 3, with 12,678, as shown in Figure 1 Table 4. The number of CNVs was significantly positively correlated with the length of the chromosome, R 2 = 0.49, P = 0.0001, as shown in Figure 1 Table 3 and B and Table 4. Among all CNVs of deletion types, the CNVs with a length of 2Kb-5Kb accounted for the highest proportion, 47.16%. Among all CNVs of duplication types, the CNVs with a length of 5Kb-10Kb accounted for the highest proportion, 41.62%. When the length was 50Kb-300Kb, the number of CNVs decreased significantly, and the CNVs of deletion types accounted for only 0.23%, and the CNVs of duplication types accounted for only 0.63%, as shown in Figure 1 Table 5.

[0043] Table 3 Statistics of CNVs and CNVRs

[0044] Table 4 Distribution of CNVs on 29 pairs of autosomes of Inner Mongolia Cashmere Goat

[0045] Table 5 Distribution of CNVs in different length intervals

[0046] 3.2, Distribution of copy number variation regions All high-quality CNVs overlapping between samples were combined together, and a total of 12,609 CNVRs were obtained, with an average length of 12,258.59 bp, of which the average length of deletions was 6,211.79 bp, the average length of duplications was 17,534.39 bp, and the average length of mixtures was 26,248.65 bp. These CNVRs included 8,131 deletions, accounting for 64.49%, 1,547 duplications, accounting for 12.27%, and 2,931 mixtures, accounting for 23.24%, as shown in Table 3, andFigure 2 The number of CNVRs of deletion type on chromosome 5 and 6 was the most, which was 415; the number of CNVRs of duplication type on chromosome 1 was the most, which was 115; the number of CNVRs of mixed type on chromosome 5 was the most, which was 167, see Table 6 and Supplementary Table 6. Figure 2 Table 6 and Supplementary Table 6. The coverage of CNVRs on chromosome 18 was the highest, which was 14.26%. Compared with other chromosomes, the coverage of CNVRs of deletion, duplication and mixed type on chromosome 18 was the highest, which was 4.00%, 2.35% and 7.91%, respectively. The coverage of CNVRs on chromosome 20 was the lowest, which was 3.80%; compared with other chromosomes, the coverage of CNVRs of deletion type on chromosome 10 was the lowest, which was 1.26%; the coverage of CNVRs of duplication type on chromosome 12 was the lowest, which was 0.64%; the coverage of CNVRs of mixed type on chromosome 20 was the lowest, which was 1.38%, see Table 6 and Supplementary Table 6.

[0047] All CNVRs on 29 autosomes covered 154,568,600 bp, which was 6.27% of the reference genome. The length of 77.64% of CNVRs was 2 Kb~20 Kb, see Table 7. Figure 2 Table 7. Distribution of CNVRs on 29 autosomes of Inner Mongolia Cashmere Goat Figure 2 It can be seen from Table 7 that these CNVRs had good uniformity on 29 autosomes.

[0048] Table 6. Distribution of CNVRs length and number on 29 autosomes of Inner Mongolia Cashmere Goat

[0049] Supplementary Table 6. Distribution of CNVRs length and number on 29 autosomes of Inner Mongolia Cashmere Goat

[0050] Table 7. Distribution of CNVRs on different length intervals

[0051] 3.3. PCA analysis The first three principal components were calculated by Plink software, and the population structure clustered by the first and second principal components was as shown in Figure 2. Figure 3The individuals in the population were roughly grouped into 3 categories, indicating that there was a certain degree of population stratification in the population. For the correction of population stratification, the published "breeding value + residual" was used as the corrected phenotype value to ensure the accuracy of the association analysis results, see "Rong Y, Wang X, Na Q, et al. Genome-wide association study for cashmere traits in Inner Mongolia cashmere goat population reveals new candidate genes and haplotypes. BMC genomics 2024; 25:658."

[0052] 3.4 CNV affecting hair length traits Further study of CNVs reaching the genome-wide significance level found that deletion mutations at 155,344,701~155,346,600 on chromosome 1 of the goat genome can significantly affect the hair length traits of cashmere goats.

[0053] The association analysis between the deletion mutation at positions 155,344,701 to 155,346,600 on chromosome 1 of the goat genome and the hair length trait is shown in Table 8: Table 8 Polymorphisms at positions 155,344,701 to 155,346,600 on chromosome 1 of the goat genome

[0054] Note: Different letters indicate significant differences ( p <0.05), the same letter indicates no significant difference ( p >0.05) Table 8 shows that wool length is greatest for individuals with a normal insertion in this region, while wool length is smallest for individuals with a deletion in this region. It should be noted that normal individuals are those with an insertion at positions 155,344,701–155,346,600 on chromosome 1 of the goat genome.

[0055] When using a mixed linear model to perform genome-wide association analysis of CNV, the CNV marker at 155,344,701~155,346,600 on chromosome 1 of the goat genome reached a genome-wide significant level, indicating that the marker is significantly associated with the hair length trait of cashmere goats, and when the marker is normal, it is beneficial for cashmere goats to have longer wool length.

[0056] This suggests that the presence of a deletion mutation at positions 155,344,701 to 155,346,600 on chromosome 1 of the goat genome can be used to select cashmere goats with longer fleeces. By selecting normal individuals as either sires or dams, the fleece length of offspring can be increased. Furthermore, based on the prior art "Study on the Inheritance of Coat Types in Inner Mongolia Cashmere Goats and Their Indirect Selection for Important Economic Traits [D]. Li Xuewu, Inner Mongolia Agricultural University," which states that "selecting long-haired types can reduce cashmere fineness and increase cashmere length and body weight," the molecular markers associated with cashmere fleece length provided by this invention can also indirectly achieve the goal of increasing cashmere production and weight in offspring.

[0057] 3.6. Verification of Detected CNVs Using qPCR In order to verify the association analysis results in the present invention, 10 important candidate CNVs were randomly selected, which were located in 10 co-localized genes. DYRK1A, MMADHC, TPK1, TDRD3, KRT222, MDGA1, NOX4, SH2D5, MIR17 and SATB1 qPCR validation was performed.

[0058] The results are as follows Figure 4~Figure 12 As shown, the 10 CNVs were of deletion type in all individuals, proving that the CNVs significantly associated with villus traits obtained by the present invention were accurate.

[0059] The above merely describes the preferred embodiments of the present application, and it should be pointed out that those skilled in the art can make several improvements and refinements without departing from the principles of the present application, and these improvements and refinements should also be considered as falling within the protection scope of the present application.

Claims

1. Use of a detection reagent for molecular markers related to cashmere goat hair length traits in identifying cashmere goat hair length traits or cashmere goat genetic breeding, characterized in that: The nucleotide sequence of the molecular marker is shown in SEQ ID NO.

1.

2. Use of the molecular marker detection reagent related to cashmere goat hair length traits according to claim 1 in identifying cashmere goat hair length traits or cashmere goat genetic breeding, characterized in that: The method for identifying the hair length traits of cashmere goats includes the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome to be tested is detected, and the wool length of the insertion type individual is greater than the wool length of the deletion type individual.

3. Use of the molecular marker detection reagent related to cashmere goat hair length trait according to claim 1 in identifying cashmere goat hair length trait or cashmere goat genetic breeding, characterized in that: The genetic breeding is to select cashmere goats with long hair length.

4. Use of the molecular marker detection reagent related to cashmere goat hair length trait according to claim 3 in identifying cashmere goat hair length trait or cashmere goat genetic breeding, characterized in that: The method for breeding cashmere goats with large hair length comprises the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome is detected, and the insertion-type individuals are selected and used as parents to increase the wool length of the offspring cashmere goats.

5. Use of the molecular marker detection reagent related to cashmere goat hair length trait according to claim 1 in identifying cashmere goat hair length trait or cashmere goat genetic breeding, characterized in that: The genetic breeding is to select cashmere goats with large cashmere production or large body weight.

6. Use of the molecular marker detection reagent related to cashmere goat hair length traits according to claim 5 in identifying cashmere goat hair length traits or cashmere goat genetic breeding, characterized in that: The method for breeding cashmere goats with large cashmere production comprises the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome is detected, and the insertion-type individuals are selected and used as parents to increase the cashmere production of offspring cashmere goats.

7. Use of the molecular marker detection reagent related to cashmere goat hair length trait according to claim 5 in identifying cashmere goat hair length trait or cashmere goat genetic breeding, characterized in that: The method for breeding the heavy cashmere goat comprises the following steps: The deletion or insertion of the molecular marker on chromosome 1 of the goat genome is detected, and the insertion-type individuals are selected and used as parents to increase the weight of the offspring cashmere goats.

8. Use of the molecular marker detection reagent related to cashmere goat hair length traits according to claim 1 in identifying cashmere goat hair length traits or cashmere goat genetic breeding, characterized in that: The cashmere goats are Inner Mongolia cashmere goats.

Images