CN114990226B - 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 - Google Patents
一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 Download PDFInfo
- Publication number
- CN114990226B CN114990226B CN202210602891.2A CN202210602891A CN114990226B CN 114990226 B CN114990226 B CN 114990226B CN 202210602891 A CN202210602891 A CN 202210602891A CN 114990226 B CN114990226 B CN 114990226B
- Authority
- CN
- China
- Prior art keywords
- deoxynucleotide
- chr
- wool
- sheep
- diameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6858—Allele-specific amplification
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/124—Animal traits, i.e. production traits, including athletic performance or the like
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/16—Primer sets for multiplex assays
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P60/00—Technologies relating to agriculture, livestock or agroalimentary industries
- Y02P60/80—Food processing, e.g. use of renewable energies or variable speed drives in handling, conveying or stacking
- Y02P60/87—Re-use of by-products of food processing for fodder production
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明属于生物技术领域,具体涉及一种与细毛羊羊毛纤维直径相关的SNP位点组合及应用。本发明提供了与细毛羊羊毛纤维直径相关的33个SNP位点组合,所述SNP位点基于绵羊v4.0基因组序列比对确定;本发明提供了检测SNP位点组合的基因型的试剂在检测细毛羊羊毛纤维直径或细毛羊分子标记辅助育种中的应用;以及基于该位点组合形成的分子探针组合、基因芯片、试剂盒,利用本发明提供的位点组合及分子探针组合、基因芯片、试剂盒能对细毛羊个体的羊毛纤维直径进行分析,对早期难以度量的细毛羊羊毛纤维直径性状进行个体选择,缩短世代间隔,加速育种进程,节约育种成本,为今后细毛羊的鉴定、保种和分子遗传育种提供支撑。
Description
技术领域
本发明属于基因育种技术领域,具体涉及一种与细毛羊羊毛纤维直径相关的SNP位点组合及其应用。
背景技术
绵羊(Ovis aries)是一种具有重要农业和生物学意义的家畜,它作为最早被驯化的动物之一,为人类提供肉、奶、羊毛和羔皮等,在全球农业经济中发挥着至关重要的作用。羊毛是优质纺织原料的来源,对国民经济具有重大作用。
我国拥有丰富的毛用绵羊品种资源。羊毛是一种天然的高性能材料,具有防污、柔软、防晒、保暖和透气等多种特性。因此,羊毛在纺织加工中占有重要地位。羊毛质量由纤维直径、纤维长度、卷曲、颜色和髓质百分比决定。纤维直径通常与羊毛加工性能相关,并决定其最终用途。上述羊毛性状受遗传和非遗传因素的共同影响。纤维直径作为细毛羊重要经济性状之一,通常占羊毛单位价值的75%–80%。
纤维直径是羊毛质量和价值的主要决定因素。平均纤维直径是可测量的最重要的原毛性能之一。因此,平均纤维直径是含脂羊毛价格的重要决定因素。纤维直径也是少数在加工过程中基本保持不变的原毛参数之一。直径限制了可以从给定原材料纺出的纱线厚度(支数),对于给定的纱线支数,纱线的各种物理特性,例如弯曲刚度和伸长率,均取决于其组成纤维直径。为了探究细毛羊纤维直径相关的遗传特征,本发现选择采用应用广泛的分子遗传标记方法。
分子遗传标记方法(Molecular Genetic Markers)是以个体遗传物质内核苷酸序列变异为基础的遗传标记,是DNA水平遗传多态性的直接反应。它具有显著的优越性,如生物发育的不用阶段,不同组织的DNA均可用于遗传标记分析;基因组变异丰富;检测手段简单快捷,易实现自动化处理等。分子遗传标记目前应用广泛的有RFLP(RestrictionFragment Length Polymorphism,限制性片段长度多态分析技术)、RAPD(RandomAmplified Polymorphism DNA,随机引物扩增多态性DNA技术)、AFLP(Amplified FragmentLength Polymorphism,扩增片段长度多态性分析技术)以及SNP(Single-nucleotidepolymorphisms,单核苷酸多态性)等。其中,SNP是研究人类家族和动植物品系遗传变异的重要依据,因此常被广泛应用于群体遗传学研究和疾病相关的基因探究。近年来,SNP常在动物遗传学分析和遗传育种中发挥着重要作用,因此,家畜遗传育种也常借助SNP来加速传统育种技术的革新,建立创新型的育种理论与体系。
目前,可以进行SNP基因型分型的技术有很多,本发明主要选用了博瑞迪生物技术公司的GenoBaits(基于液相探针杂交的靶向基因捕获技术方案)和GenoPlexs(基于多重PCR的靶向基因捕获技术方案)的靶向基因捕获测序技术体系。GenoBaits主要应用于全基因选择、全基因组关联分析和遗传图谱构建。GenoPlexs的主要应用方向是种质资源分析、分子标记辅助回交和分子标记辅助早筛。该技术体系主要通过降低文库丰度实现仅对目标为点进行深度重测序,具备稳定性较高、对序列保守型和GC含量容忍度高等优势。目前绵羊SNP研究中缺乏对细毛羊羊毛纤维直径相关的SNP位点结合研究。
本申请首先提供了与细毛羊羊毛纤维直径相关的33个SNP位点组合,所述SNP位点基于绵羊v4.0基因组序列比对确定;分别为:位于chr 1第203825947位,其脱氧核苷酸为C或A;位于chr 1第226733906位,其脱氧核苷酸为C或T;位于chr 3第45470146位,其脱氧核苷酸为A或C;位于chr 4第68142771位,其脱氧核苷酸为C或G;位于chr 5第93335425位,其脱氧核苷酸为C或T;位于chr 5第93344882位,其脱氧核苷酸为C或T;位于chr 5第93387255位,其脱氧核苷酸为C或T;位于chr 5第93391985位,其脱氧核苷酸为G或A;位于chr 5第93392877位,其脱氧核苷酸为C或T;位于chr 5第93393426位,其脱氧核苷酸为A或G;位于chr 5第93507537位,其脱氧核苷酸为A或C;位于chr 6第25952072位,其脱氧核苷酸为T或A;位于chr 6第37126564位,其脱氧核苷酸为T或C;位于chr 10第51739659位,其脱氧核苷酸为G或A;位于chr 11第8917643位,其脱氧核苷酸为C或A;位于chr 12第25119445位,其脱氧核苷酸为G或A;位于chr 12第25120732位,其脱氧核苷酸为C或G;位于chr 12第25135944位,其脱氧核苷酸为A或G;位于chr 12第25149517位,其脱氧核苷酸为C或A;位于chr 12第25152554位,其脱氧核苷酸为T或A;位于chr12第25154575位,其脱氧核苷酸为C或T;位于chr 12第25155325位,其脱氧核苷酸为T或C;位于chr 12第36292909位,其脱氧核苷酸为G或A;位于chr 12第78576808位,其脱氧核苷酸为T或C;位于chr 15第68745093位,其脱氧核苷酸为C或T;位于chr 20第25185724位,其脱氧核苷酸为C或T;位于chr 20第37863763位,其脱氧核苷酸为C或T;位于chr 20第44119346位,其脱氧核苷酸为G或A;位于chr 22第46988971位,其脱氧核苷酸为G或A;位于chr 22第47002481位,其脱氧核苷酸为A或C;位于chr 22第50203275位,其脱氧核苷酸为G或A;位于chr 25第13791395位,其脱氧核苷酸为T或C;位于chr26第40609491位,其脱氧核苷酸为C或T;其次,本发明通过GenoPlexs(基于多重PCR的靶向基因捕获技术方案)和GenoBaits(基于液相探针杂交的靶向基因捕获技术方案)技术,获得目标SNP的基因型,实现了细毛羊羊毛纤维直径的快速有效检测,对绵羊的分子育种以及种质资源的保护与改造均有着重要的意义。
发明内容
为满足我国当前育种生产上对生产性状方向的芯片位点功能检测。功能研究方面的需求,本发明提供了一种中国具有代表性的四个细毛羊品种(中国美利奴羊、高山美利奴羊、敖汉细毛羊和青海细毛羊)的高深度全基因组重测序数据,以绵羊v4.04基因组为参考,结合已有的与绵羊生产性状相关的研究,获得了一种检测准确,使用方便,市场前景广阔的一种细毛羊羊毛纤维直径性状的SNP位点组合,所述位点能够用于绵羊品种的选育、保护以及改良。具体包括以下内容:
第一方面,本发明提供了一种与细毛羊羊毛纤维直径相关的33个SNP位点组合,所述33个SNP位点组合基于绵羊v4.0基因组序列比对确定;分别为:位于chr 1第203825947位,其脱氧核苷酸为C或A;位于chr 1第226733906位,其脱氧核苷酸为C或T;位于chr 3第45470146位,其脱氧核苷酸为A或C;位于chr 4第68142771位,其脱氧核苷酸为C或G;位于chr 5第93335425位,其脱氧核苷酸为C或T;位于chr 5第93344882位,其脱氧核苷酸为C或T;位于chr 5第93387255位,其脱氧核苷酸为C或T;位于chr 5第93391985位,其脱氧核苷酸为G或A;位于chr 5第93392877位,其脱氧核苷酸为C或T;位于chr 5第93393426位,其脱氧核苷酸为A或G;位于chr 5第93507537位,其脱氧核苷酸为A或C;位于chr 6第25952072位,其脱氧核苷酸为T或A;位于chr 6第37126564位,其脱氧核苷酸为T或C;位于chr 10第51739659位,其脱氧核苷酸为G或A;位于chr 11第8917643位,其脱氧核苷酸为C或A;位于chr 12第25119445位,其脱氧核苷酸为G或A;位于chr12第25120732位,其脱氧核苷酸为C或G;位于chr 12第25135944位,其脱氧核苷酸为A或G;位于chr 12第25149517位,其脱氧核苷酸为C或A;位于chr 12第25152554位,其脱氧核苷酸为T或A;位于chr 12第25154575位,其脱氧核苷酸为C或T;位于chr 12第25155325位,其脱氧核苷酸为T或C;位于chr 12第36292909位,其脱氧核苷酸为G或A;位于chr 12第78576808位,其脱氧核苷酸为T或C;位于chr 15第68745093位,其脱氧核苷酸为C或T;位于chr 20第25185724位,其脱氧核苷酸为C或T;位于chr 20第37863763位,其脱氧核苷酸为C或T;位于chr 20第44119346位,其脱氧核苷酸为G或A;位于chr22第46988971位,其脱氧核苷酸为G或A;位于chr 22第47002481位,其脱氧核苷酸为A或C;位于chr 22第50203275位,其脱氧核苷酸为G或A;位于chr 25第13791395位,其脱氧核苷酸为T或C;位于chr26第40609491位,其脱氧核苷酸为C或T。
第二方面,本发明提供了检测上述第一方面所述与细毛羊羊毛纤维直径相关的33个SNP位点组合的试剂在检测与细毛羊羊毛纤维直径中的应用。
优选地,所述试剂包括用于检测所述SNP位点组合的引物,本领域技术人员根据本发明提供的细毛羊羊毛纤维直径相关的SNP位点组合中的每一个位点的序列信息设计引物,可以在同一反应条件下实现检测目的的引物。其中,引物的设计为常规方法,根据本申请提供的细毛羊羊毛纤维直径相关的SNP位点组合中的位点信息,不再需要付出创造性劳动的情况下就能获得,因此,根据本申请提供的细毛羊羊毛纤维直径相关生物SNP位点组合获得引物也属于本发明的保护范围。
优选地,所述试剂包括用于检测所述SNP位点组合的分子探针组合。分子探针的设计为常规方法,根据本申请提供的细毛羊羊毛纤维直径相关的SNP位点组合中的位点信息,不再需要付出创造性劳动的情况下就能获得,因此,根据本申请提供的细毛羊羊毛纤维直径相关生物SNP位点组合获得分子探针也属于本发明的保护范围。
优选地,所述分子探针组合如表1所示。
表1与细毛羊羊毛纤维直径相关的SNP位点组合的分子探针组合
优选地,所述试剂包括基因芯片,所述基因芯片采用常规方法将获得的引物或探针固定在聚合物基片上,例如尼龙膜、硝酸纤维膜、塑料、硅胶晶片、微型磁珠等,或将探针固定在玻璃板上,或在玻璃等硬质表面上直接合成获得的引物或探针,本申请的SNP基因芯片的使用方法与常规方法相同。
第三方面,本发明提供了一种分析细毛羊羊毛纤维直径性状的分子探针组合,所述分子探针组合检测上述第一方面中所述的与细毛羊羊毛纤维直径相关的33个SNP位点组合。
优选地,所述分子探针组合如上表1所示。
第四方面,本发明提供了一种分析细毛羊羊毛纤维直径性状的基因芯片,所述基因芯片负载有上述第三方面所述的分析细毛羊羊毛纤维直径性状的分子探针组合。
第五方面,本发明提供了一种分析细毛羊羊毛纤维直径性状的试剂盒,所述试剂盒包括上述第三方面所述的分析细毛羊羊毛纤维直径性状的分子探针组合或第四方面所述的分析细毛羊羊毛纤维直径性状的基因芯片。
第六方面,本发明提供了上述第三方面所述的分子探针组合,或上述第四方面所述的基因芯片,或上述第五方面所述的试剂盒在细毛羊羊毛纤维直径性状评价,或在细毛羊品种筛选,或在细毛羊品种鉴定,或在细毛羊分子标记辅助育种中的应用。
第七方面,本发明提供了一种分析细毛羊羊毛纤维直径性状的方法,所述方法为:检测待测细毛羊的基因组DNA中如上述第一方面所述的与细毛羊羊毛纤维直径相关的33个SNP位点基因型;对照细毛羊基因组DNA的所述33个SNP位点基因型进行比较,根据基因型检测结果判断细毛羊的羊毛纤维直径性状。
本发明的有益效果是:本发明提供了与细毛羊羊毛纤维直径相关的33个SNP位点组合,所述SNP位点基于绵羊v4.0基因组序列比对确定;其次,本发明发现通过分子探针或基因芯片等方式待测细毛羊基因组DNA中的与细毛羊羊毛纤维直径相关的33个SNP位点组合的基因型,能够用于细毛羊羊毛纤维直径性状分析,用于细毛羊早期育种选择,实现对早期难以度量的净毛率性状进行个体选择,缩短世代间隔,加速育种进程,节约大量的育种成本,为今后细毛羊的鉴定、保种、以及遗传育种提供支撑;而且基于本发明提供的细毛羊羊毛纤维直径相关的33个SNP位点组合形成的分子探针组合、基因芯片、试剂盒与现有的高密度芯片相比,通量小、成本低,分析更容易,普适性广,市场前景广阔。
附图说明
图1对与细毛羊羊毛纤维直径性状相关的SNP数据进行GWAS中GLM模型下实施例1中计算获得的p值取-log10以后绘制的曼哈顿图;
图2对与细毛羊羊毛纤维直径性状相关的SNP数据进行GWAS中GLM模型下实施例1中计算获得的p值取-log10以后绘制的Q-Q图。
具体实施方式
下面结合实施例对本发明的技术方案进行详细说明。需要指出的是以下实施例的给出仅是为了起到说明的目的,并不是用于对本发明的范围进行限制。本领域的技术人员在不背离本发明的宗旨和精神的情况下,可以对本发明进行各种修改和替换。
下述实施例中所有试验的实验方法如无特殊说明,均为常规方法。
下述实施例中所有试验的实验条件如无特殊说明,均为常规条件,如Sambrook等分子克隆实验手册,或按照制造厂商说明书建议的条件。
本发明所述的SNP是单核苷酸多态性的简称,指在基因组水平上由单个核苷酸的变异所引起的DNA序列多态性。
实施例1细毛羊羊毛纤维直径相关的SNP位点
1、总SNP集合的获取
对我国具有代表性的四个细毛羊品种中460个细毛羊个体进行了全基因组重测序,平均深度为5X,应用重测序分析流程,与2015年发布的绵羊v4.0参考基因组(从NCBI中获取)进行比对,两种方式比对获得的共同结果形成一个SNP集合。
具体而言,对多个细毛羊个体进行高深度重测序是由生物测序公司完成,生物测序公司完成的测序结果均能够实现本发明的技术目的,本发明不做限制。本申请将测序公司返回的Fastq文件通过BAM文件比对到参考基因组绵羊v4.0后的得到BAM文件,利用SAMtools和GATK软件对样本BAM文件分析得到的包含群体SNP分型信息的VCF文件,将两种方式获得的VCF文件结果合并,经过质量筛选后获得了包含33个SNP位点的SNP集合。
具体的,本发明所使用的细毛羊品种为中国四个具有代表性的细毛羊品种,分别是中国美利奴羊、高山美利奴羊、敖汉细毛羊和青海细毛羊。
2、候选基因与所在功能区域的筛选
根据中国具有代表性的四个细毛羊品种(中国美利奴羊、高山美利奴羊、敖汉细毛羊和青海细毛羊)在羊毛纤维直径性状上表现出来的显著性差异,首先利用自写的perl脚本对其进行标记质控,去除此等位基因频率小于0.05,缺失率大于20%,杂合比例大于80%以及非二等位的位点。其次借助五部分的群体分析,其中包括由MEGA-X软件完成的系统发育树的构建、由Admixture软件(v1.3)软件完成的群体结构分析、由gcta(v1.92.2)软件完成的PCA分析和亲缘关系分析以及由软件HaploviewLD完成的衰减分析,可以综合评判材料的遗传多样性以及遗传背景是否存在较大的差异,揭示非家系群体或系谱不明确的群体材料的遗传相似度以及各个亚群和整体材料的受选择程度,由此确定调整GWAS(全基因组关联分析)所采用的模型。最终通过GLM(简单线性模型),以曼哈顿图与QQ图展示定位到了与细毛羊羊毛纤维直径相关的SNP位点以及候选基因,以阈值为0.01筛选出了其显著结果,确定了30个与细毛羊羊毛纤维直径相关的,功能确定的候选基因或标记SOX2、DNAJC19、MFSD1、RARRES1、EHBP1、TMEM17、JAZF1、CAST、ERAP1、ERAP2、TSPAN5、FAM184B、LOC101103163、KCTD12、RNF43、CAPN2、PRRX1、TNNT2、LOC101112664、LOC101108158、ELOVL5、ID4、RNF144B、ELOVL2、MKI67、MGMT、LOC101110287、HNRNPF、BICC1和UBE2E1。
3、功能基因位点对应SNP位点的提取
利用GWAS模型的统一表达式:y=Xα+Qβ+Kμ+e,其中,y为表型向量,X为基因型矩阵,α为基因型效应向量,Q为固定效应矩阵(可以为群体结构/性别/地点/场次等信息),β为固定效应向量,K为随机效应矩阵,主要指亲缘关系矩阵,μ为随机效应向量,e为残差向量。针对每个SNP位点,都检验α是否为0,α为0的概率值p用于衡量标记基因型与表型的关联程度,p值越小,α为0的概率越小,该标记越可能与性状关联。从而有步骤2确定的候选基因所在的功能区域对应的SNP位点,得到SOX2、DNAJC19、MFSD1、RARRES1、EHBP1、TMEM17、JAZF1、CAST、ERAP1、ERAP2、TSPAN5、FAM184B、LOC101103163、KCTD12、RNF43、CAPN2、PRRX1、TNNT2、LOC101112664、LOC101108158、ELOVL5、ID4、RNF144B、ELOVL2、MKI67、MGMT、LOC101110287、HNRNPF、BICC1和UBE2E1共30个与纤维直径相关联的功能基因或标记,且仅包含33个SNP位点的位点组合。
所述33个SNP位点组合的物理信息具体如下表2所示。
表2细毛羊羊毛纤维直径性状SNP位点组合的物理信息
number | Mark | Chr | Position | ref | mut |
SNP00001 | 1_203825947 | 1 | 203825947 | G | A |
SNP00002 | 1_226733906 | 1 | 226733906 | C | T |
SNP00003 | 3_45470146 | 3 | 45470146 | A | C |
SNP00004 | 4_68142771 | 4 | 68142771 | C | G |
SNP00005 | 5_93335425 | 5 | 93335425 | C | T |
SNP00006 | 5_93344882 | 5 | 93344882 | C | T |
SNP00007 | 5_93387255 | 5 | 93387255 | C | T |
SNP00008 | 5_93391985 | 5 | 93391985 | G | A |
SNP00009 | 5_93392877 | 5 | 93392877 | C | T |
SNP00010 | 5_93393426 | 5 | 93393426 | A | G |
SNP00011 | 5_93507537 | 5 | 93507537 | A | C |
SNP00012 | 6_25952072 | 6 | 25952072 | T | A |
SNP00013 | 6_37126564 | 6 | 37126564 | T | C |
SNP00014 | 10_51739659 | 10 | 51739659 | G | A |
SNP00015 | 11_8917643 | 11 | 8917643 | C | A |
SNP00016 | 12_25119445 | 12 | 25119445 | G | A |
SNP00017 | 12_25120732 | 12 | 25120732 | C | G |
SNP00018 | 12_25135944 | 12 | 25135944 | A | G |
SNP00019 | 12_25149517 | 12 | 25149517 | C | A |
SNP00020 | 12_25152554 | 12 | 25152554 | T | A |
SNP00021 | 12_25154575 | 12 | 25154575 | C | T |
SNP00022 | 12_25155325 | 12 | 25155325 | T | C |
SNP00023 | 12_36292909 | 12 | 36292909 | G | A |
SNP00024 | 12_78576808 | 12 | 78576808 | T | C |
SNP00025 | 15_68745093 | 15 | 68745093 | C | T |
SNP00026 | 20_25185724 | 20 | 25185724 | C | T |
SNP00027 | 20_37863763 | 20 | 37863763 | C | T |
SNP00028 | 20_44119346 | 20 | 44119346 | G | A |
SNP00029 | 22_46988971 | 22 | 46988971 | G | A |
SNP00030 | 22_47002481 | 22 | 47002481 | A | C |
SNP00031 | 22_50203275 | 22 | 50203275 | G | A |
SNP00032 | 25_13791395 | 25 | 13791395 | T | C |
SNP00033 | 26_40609491 | 26 | 40609491 | C | T |
实施例2细毛羊羊毛纤维直径相关的SNP的panel制备
基于实施例1获得的SNP位点组合,本申请委托博瑞迪生物科技有限公司进行了纤维直径相关的SNP的panel制备。在定量质检合格的DNA中加入多重PCR Panel mix和多重PCR扩增酶体系,置于PCR仪上完成PCR反应。PCR产物利用羧基磁珠进行纯化后,再次计入带有Barcode的测序引物和高保真PCR反应体系进行PCR扩增,不同的Barcode用于区分不同的样品。经过羧基磁珠纯化后扩增产物,即完成多重PCR捕获及文库建库。本领域技术人员根据本发明提供的细毛羊净毛率相关的SNP位点组合中的每一个位点的序列信息设计引物,为常规方法,不再需要付出创造性劳动的情况下就能获得;而且,panel制备也是依据本发明提供的细毛羊净毛率相关的SNP位点组合能够常规制备的。
实施例3 437个细毛羊个体羊毛纤维直径的检测
基于实施例1获取到的SNP位点与实施例2的1K penal制备对细毛羊的个体进行检测,在本发明的一个实施例中,选择了采用博瑞迪自主研发的GenoBaits(基于液相探针杂交的靶向基因捕获技术方案)对其个体进行羊毛纤维直径的检测。该技术的工作原理是基于目标探针与靶向序列互补结合进行定点捕获,对捕获的靶点序列进行洗脱、靶点扩增、建库和测序,最终获得目标SNP的基因型,在经济有效的条件下,所能检测的靶向位点及其标记数目在检测密度和通量上等同于高密度固相芯片。通过该技术从而获得目标样本的结果值。细毛羊羊毛纤维直径相关位点多态性检测结果如表3所示。
表3细毛羊羊毛纤维直径相关位点多态性
细毛羊不同基因型与羊毛纤维直径之间的相关分析结果如表4所示。
表4细毛羊不同基因型与羊毛纤维直径之间的相关分析结果
上述结果表明,通过检测本发明所述的33个细毛羊羊毛纤维直径SNP位点组合的基因型,可以对细毛羊羊毛纤维直径进行分析,其中33个SNP位点组合分别为:分别为:位于chr 1第203825947位,其脱氧核苷酸为C或A;位于chr 1第226733906位,其脱氧核苷酸为C或T;位于chr 3第45470146位,其脱氧核苷酸为A或C;位于chr 4第68142771位,其脱氧核苷酸为C或G;位于chr 5第93335425位,其脱氧核苷酸为C或T;位于chr 5第93344882位,其脱氧核苷酸为C或T;位于chr 5第93387255位,其脱氧核苷酸为C或T;位于chr 5第93391985位,其脱氧核苷酸为G或A;位于chr 5第93392877位,其脱氧核苷酸为C或T;位于chr 5第93393426位,其脱氧核苷酸为A或G;位于chr 5第93507537位,其脱氧核苷酸为A或C;位于chr 6第25952072位,其脱氧核苷酸为T或A;位于chr 6第37126564位,其脱氧核苷酸为T或C;位于chr 10第51739659位,其脱氧核苷酸为G或A;位于chr 11第8917643位,其脱氧核苷酸为C或A;位于chr 12第25119445位,其脱氧核苷酸为G或A;位于chr 12第25120732位,其脱氧核苷酸为C或G;位于chr 12第25135944位,其脱氧核苷酸为A或G;位于chr 12第25149517位,其脱氧核苷酸为C或A;位于chr12第25152554位,其脱氧核苷酸为T或A;位于chr 12第25154575位,其脱氧核苷酸为C或T;位于chr 12第25155325位,其脱氧核苷酸为T或C;位于chr 12第36292909位,其脱氧核苷酸为G或A;位于chr 12第78576808位,其脱氧核苷酸为T或C;位于chr 15第68745093位,其脱氧核苷酸为C或T;位于chr 20第25185724位,其脱氧核苷酸为C或T;位于chr 20第37863763位,其脱氧核苷酸为C或T;位于chr 20第44119346位,其脱氧核苷酸为G或A;位于chr 22第46988971位,其脱氧核苷酸为G或A;位于chr 22第47002481位,其脱氧核苷酸为A或C;位于chr 22第50203275位,其脱氧核苷酸为G或A;位于chr25第13791395位,其脱氧核苷酸为T或C;位于chr26第40609491位,其脱氧核苷酸为C或T。
本领域技术人员基于本发明提供的仅由33个SNP位点组成的细毛羊羊毛纤维直径SNP位点组合可以制成的分析细毛羊羊毛纤维直径的SNP探针组合、基因芯片和试剂盒,能够在基因组水平上对细毛羊羊毛纤维直径进行分析,或遗传评估、品种筛选、品种鉴定,以获得更高的育种值估计准确性,控制育种进程,还能够应用于绵羊系谱重构、绵羊品种溯源、种质资源保护和种质资源改良。由于与细毛羊羊毛纤维直径相关的研究较为稀缺,因此本申请旨在加速细毛羊的分子育种进程,增强对该物种资源的保护与改良,节约大量的育种成本,由此提升细毛羊带来的毛用经济效益。
以上所述仅为帮助理解本发明的优选实例,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化,在不违背本发明的思想下,本领域技术人员在此基础上对本发明作出的各种改动或者修改,同样应属于本发明的范围。
Claims (3)
1.检测与细毛羊羊毛纤维直径相关的33个SNP位点组合的试剂在检测细毛羊羊毛纤维直径中的应用;所述33个SNP位点组合基于绵羊v4.0基因组序列比对确定;分别为:位于chr1第203825947位,其脱氧核苷酸为C或A;位于chr 1第226733906位,其脱氧核苷酸为C或T;位于chr 3第45470146位,其脱氧核苷酸为A或C;位于chr 4第68142771位,其脱氧核苷酸为C或G;位于chr 5第93335425位,其脱氧核苷酸为C或T;位于chr 5第93344882位,其脱氧核苷酸为C或T;位于chr 5第93387255位,其脱氧核苷酸为C或T;位于chr 5第93391985位,其脱氧核苷酸为G或A;位于chr 5第93392877位,其脱氧核苷酸为C或T;位于chr 5第93393426位,其脱氧核苷酸为A或G;位于chr 5第93507537位,其脱氧核苷酸为A或C;位于chr 6第25952072位,其脱氧核苷酸为T或A;位于chr 6第37126564位,其脱氧核苷酸为T或C;位于chr 10第51739659位,其脱氧核苷酸为G或A;位于chr 11第8917643位,其脱氧核苷酸为C或A;位于chr 12第25119445位,其脱氧核苷酸为G或A;位于chr 12第25120732位,其脱氧核苷酸为C或G;位于chr 12第25135944位,其脱氧核苷酸为A或G;位于chr 12第25149517位,其脱氧核苷酸为C或A;位于chr12第25152554位,其脱氧核苷酸为T或A;位于chr 12第25154575位,其脱氧核苷酸为C或T;位于chr 12第25155325位,其脱氧核苷酸为T或C;位于chr 12第36292909位,其脱氧核苷酸为G或A;位于chr 12第78576808位,其脱氧核苷酸为T或C;位于chr 15第68745093位,其脱氧核苷酸为C或T;位于chr 20第25185724位,其脱氧核苷酸为C或T;位于chr 20第37863763位,其脱氧核苷酸为C或T;位于chr 20第44119346位,其脱氧核苷酸为G或A;位于chr 22第46988971位,其脱氧核苷酸为G或A;位于chr 22第47002481位,其脱氧核苷酸为A或C;位于chr 22第50203275位,其脱氧核苷酸为G或A;位于chr25第13791395位,其脱氧核苷酸为T或C;位于chr26第40609491位,其脱氧核苷酸为C或T。
2.如权利要求1所述的应用,其特征在于,所述试剂包括用于检测所述SNP位点组合的探针组合。
3.检测与细毛羊羊毛纤维直径相关的33个SNP位点组合的试剂在制备检测细毛羊羊毛纤维直径性状的基因芯片中应用,所述33个SNP位点组合基于绵羊v4.0基因组序列比对确定;分别为:位于chr 1第203825947位,其脱氧核苷酸为C或A;位于chr 1第226733906位,其脱氧核苷酸为C或T;位于chr 3第45470146位,其脱氧核苷酸为A或C;位于chr 4第68142771位,其脱氧核苷酸为C或G;位于chr 5第93335425位,其脱氧核苷酸为C或T;位于chr 5第93344882位,其脱氧核苷酸为C或T;位于chr 5第93387255位,其脱氧核苷酸为C或T;位于chr 5第93391985位,其脱氧核苷酸为G或A;位于chr 5第93392877位,其脱氧核苷酸为C或T;位于chr 5第93393426位,其脱氧核苷酸为A或G;位于chr 5第93507537位,其脱氧核苷酸为A或C;位于chr 6第25952072位,其脱氧核苷酸为T或A;位于chr 6第37126564位,其脱氧核苷酸为T或C;位于chr 10第51739659位,其脱氧核苷酸为G或A;位于chr 11第8917643位,其脱氧核苷酸为C或A;位于chr 12第25119445位,其脱氧核苷酸为G或A;位于chr 12第25120732位,其脱氧核苷酸为C或G;位于chr12第25135944位,其脱氧核苷酸为A或G;位于chr 12第25149517位,其脱氧核苷酸为C或A;位于chr 12第25152554位,其脱氧核苷酸为T或A;位于chr 12第25154575位,其脱氧核苷酸为C或T;位于chr 12第25155325位,其脱氧核苷酸为T或C;位于chr 12第36292909位,其脱氧核苷酸为G或A;位于chr 12第78576808位,其脱氧核苷酸为T或C;位于chr 15第68745093位,其脱氧核苷酸为C或T;位于chr 20第25185724位,其脱氧核苷酸为C或T;位于chr 20第37863763位,其脱氧核苷酸为C或T;位于chr 20第44119346位,其脱氧核苷酸为G或A;位于chr 22第46988971位,其脱氧核苷酸为G或A;位于chr22第47002481位,其脱氧核苷酸为A或C;位于chr 22第50203275位,其脱氧核苷酸为G或A;位于chr 25第13791395位,其脱氧核苷酸为T或C;位于chr26第40609491位,其脱氧核苷酸为C或T。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210602891.2A CN114990226B (zh) | 2022-05-30 | 2022-05-30 | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 |
PCT/CN2023/096722 WO2023231940A1 (zh) | 2022-05-30 | 2023-05-29 | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 |
AU2023280943A AU2023280943A1 (en) | 2022-05-30 | 2023-05-29 | Snp locus combinations related to diameter of fine wool sheep wool fibers, and use |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210602891.2A CN114990226B (zh) | 2022-05-30 | 2022-05-30 | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114990226A CN114990226A (zh) | 2022-09-02 |
CN114990226B true CN114990226B (zh) | 2023-06-09 |
Family
ID=83031577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210602891.2A Active CN114990226B (zh) | 2022-05-30 | 2022-05-30 | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 |
Country Status (3)
Country | Link |
---|---|
CN (1) | CN114990226B (zh) |
AU (1) | AU2023280943A1 (zh) |
WO (1) | WO2023231940A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114990226B (zh) * | 2022-05-30 | 2023-06-09 | 中国农业科学院兰州畜牧与兽药研究所 | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103243167A (zh) * | 2013-05-15 | 2013-08-14 | 新疆农垦科学院 | 与绵羊羊毛纤维直径相关的分子标记及其应用 |
CN109385471A (zh) * | 2018-09-26 | 2019-02-26 | 塔里木大学 | 绵羊基因组26号染色体上绵羊脱毛症相关基因的分析方法 |
CN113278712A (zh) * | 2021-07-23 | 2021-08-20 | 中国农业大学 | 分析绵羊毛色的基因芯片、分子探针组合、试剂盒及应用 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2368428A1 (en) * | 2010-03-16 | 2011-09-28 | Institut National De La Recherche Agronomique | Obtention of a rex animal by molecular methods based on the alteration of the LIPH function |
CN103276098B (zh) * | 2013-06-14 | 2014-07-16 | 东北农业大学 | 一种预示和鉴定绵羊羊毛长度的分子标记方法 |
CN104894253B (zh) * | 2015-05-20 | 2018-04-03 | 东北农业大学 | 能预示和鉴定绵羊羊毛细度的分子标记方法及其引物对 |
CN107604078B (zh) * | 2017-10-23 | 2021-05-14 | 新疆畜牧科学院畜牧研究所 | 与绵羊羊毛纤维直径性状相关的分子标记及其特异性引物对和应用 |
CN110734984B (zh) * | 2019-10-15 | 2023-02-03 | 甘肃农业大学 | 与细毛羊羊毛纤维直径相关的遗传标记及其应用 |
CN111394479A (zh) * | 2020-04-28 | 2020-07-10 | 中国计量大学 | 一种用于细毛羊辅助选育的分子标记及应用方法 |
CN112048562B (zh) * | 2020-08-28 | 2021-05-25 | 中国农业科学院兰州畜牧与兽药研究所 | 一种影响高山美利奴羊羊毛纤维直径的snp分子标记及其应用 |
CN111996265B (zh) * | 2020-09-20 | 2021-07-06 | 中国农业科学院兰州畜牧与兽药研究所 | 一种影响细毛羊羊毛纤维直径的snp分子标记及其应用 |
CN112029872B (zh) * | 2020-09-22 | 2021-04-20 | 中国农业科学院北京畜牧兽医研究所 | 一种与细毛羊羊毛性状相关的snp标记及其检测引物组、试剂盒、检测方法和应用 |
CN114231642B (zh) * | 2022-01-07 | 2024-01-02 | 新疆畜牧科学院畜牧研究所 | 与鄂尔多斯细毛羊羊毛纤维直径性状相关的分子标记及特异性引物对和应用 |
CN114990226B (zh) * | 2022-05-30 | 2023-06-09 | 中国农业科学院兰州畜牧与兽药研究所 | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 |
-
2022
- 2022-05-30 CN CN202210602891.2A patent/CN114990226B/zh active Active
-
2023
- 2023-05-29 WO PCT/CN2023/096722 patent/WO2023231940A1/zh active Application Filing
- 2023-05-29 AU AU2023280943A patent/AU2023280943A1/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103243167A (zh) * | 2013-05-15 | 2013-08-14 | 新疆农垦科学院 | 与绵羊羊毛纤维直径相关的分子标记及其应用 |
CN109385471A (zh) * | 2018-09-26 | 2019-02-26 | 塔里木大学 | 绵羊基因组26号染色体上绵羊脱毛症相关基因的分析方法 |
CN113278712A (zh) * | 2021-07-23 | 2021-08-20 | 中国农业大学 | 分析绵羊毛色的基因芯片、分子探针组合、试剂盒及应用 |
Also Published As
Publication number | Publication date |
---|---|
AU2023280943A1 (en) | 2024-01-25 |
WO2023231940A1 (zh) | 2023-12-07 |
CN114990226A (zh) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114959059B (zh) | 一种与细毛羊羊毛纤维直径变异系数相关的snp位点组合及其应用 | |
WO2023001211A1 (zh) | 一种分析绵羊毛用性状的基因芯片、分子探针组合、试剂盒及应用 | |
CN114790483B (zh) | 一种与细毛羊净毛率相关的snp位点组合及其应用 | |
CN113278712B (zh) | 分析绵羊毛色的基因芯片、分子探针组合、试剂盒及应用 | |
CN113265476B (zh) | 分析绵羊产奶性能的基因芯片、分子探针组合、试剂盒及应用 | |
WO2023001209A1 (zh) | 分析绵羊脂尾的基因芯片、分子探针组合、试剂盒及应用 | |
CN114990226B (zh) | 一种与细毛羊羊毛纤维直径相关的snp位点组合及其应用 | |
US20220186289A1 (en) | Molecular marker related to wool yield of long-haired rabbit and use thereof | |
CN113718052B (zh) | 5000个snp位点组合的应用及小麦品种真实性身份鉴定的方法 | |
CN109706247B (zh) | 利用微卫星技术对近交系小鼠进行遗传质量监测的方法 | |
CN114752678B (zh) | 与猪达115公斤体重背膘厚相关的snp分子标记及其应用 | |
CN113293220B (zh) | 分析绵羊耳部大小的基因芯片、分子探针组合、试剂盒及应用 | |
CN107937557A (zh) | 一个与猪经济性状相关的snp位点及其应用 | |
CN117106936B (zh) | 分析山羊毛色性状的分子标记组合及应用 | |
CN117089633B (zh) | 分析山羊绒毛有无的分子标记组合及应用 | |
CN117089636B (zh) | 分析山羊肉用性能的分子标记组合及应用 | |
CN117089634B (zh) | 分析山羊奶用性能的分子标记组合及应用 | |
CN117089635B (zh) | 分析山羊繁殖性能的分子标记组合及应用 | |
CN114622020B (zh) | 一种与鸡生长性状相关的klhl31基因分子标记及应用 | |
CN117126948B (zh) | 分析山羊耳性状的分子标记组合及其应用 | |
CN117106935A (zh) | 分析山羊有无角性状的分子标记组合及应用 | |
CN112725464A (zh) | 与长白猪无效乳头数相关的snp分子标记及其获取方法 | |
CN117867132A (zh) | 一种bmp2基因下游snp标记在绵羊品种尾型选择中的应用 | |
CN112725463A (zh) | 与长白猪总乳头数相关的snp分子标记及其应用和获取方法 | |
CN116497125A (zh) | 渤海三疣梭子蟹原种的snp标记及其应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |