CN113039288B - 一种蛋鸡全基因组snp芯片及其应用 - Google Patents
一种蛋鸡全基因组snp芯片及其应用 Download PDFInfo
- Publication number
- CN113039288B CN113039288B CN201880004561.5A CN201880004561A CN113039288B CN 113039288 B CN113039288 B CN 113039288B CN 201880004561 A CN201880004561 A CN 201880004561A CN 113039288 B CN113039288 B CN 113039288B
- Authority
- CN
- China
- Prior art keywords
- snp
- whole genome
- chip
- laying
- layer chicken
- 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
- 230000002068 genetic effect Effects 0.000 claims abstract description 25
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 208000035240 Disease Resistance Diseases 0.000 claims abstract description 10
- 108091028043 Nucleic acid sequence Proteins 0.000 claims abstract description 6
- 241000287828 Gallus gallus Species 0.000 claims description 67
- 238000009395 breeding Methods 0.000 claims description 16
- 230000001488 breeding effect Effects 0.000 claims description 16
- 239000003147 molecular marker Substances 0.000 claims description 12
- 230000036750 egg laying performance Effects 0.000 claims description 7
- 239000002773 nucleotide Substances 0.000 claims description 6
- 125000003729 nucleotide group Chemical group 0.000 claims description 6
- 230000017448 oviposition Effects 0.000 claims description 4
- 230000036039 immunity Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 abstract description 29
- 238000012098 association analyses Methods 0.000 abstract description 6
- 238000009394 selective breeding Methods 0.000 abstract description 3
- 235000013330 chicken meat Nutrition 0.000 description 48
- 238000001514 detection method Methods 0.000 description 10
- 210000000349 chromosome Anatomy 0.000 description 9
- 238000011160 research Methods 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 8
- 108020004414 DNA Proteins 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000000034 method Methods 0.000 description 6
- 108090000623 proteins and genes Proteins 0.000 description 6
- 241001465754 Metazoa Species 0.000 description 5
- 238000012163 sequencing technique Methods 0.000 description 5
- 239000013598 vector Substances 0.000 description 5
- 230000007614 genetic variation Effects 0.000 description 4
- 239000003550 marker Substances 0.000 description 4
- 238000003908 quality control method Methods 0.000 description 4
- 230000037396 body weight Effects 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 238000003205 genotyping method Methods 0.000 description 3
- 244000144977 poultry Species 0.000 description 3
- 235000013594 poultry meat Nutrition 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000010187 selection method Methods 0.000 description 3
- 108700028369 Alleles Proteins 0.000 description 2
- 208000006758 Marek Disease Diseases 0.000 description 2
- 241000220010 Rhode Species 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000007621 cluster analysis Methods 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 230000012447 hatching Effects 0.000 description 2
- 244000144972 livestock Species 0.000 description 2
- 150000007523 nucleic acids Chemical group 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000283690 Bos taurus Species 0.000 description 1
- 238000000018 DNA microarray Methods 0.000 description 1
- 239000003298 DNA probe Substances 0.000 description 1
- 241000287826 Gallus Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 208000027954 Poultry disease Diseases 0.000 description 1
- 101100495925 Schizosaccharomyces pombe (strain 972 / ATCC 24843) chr3 gene Proteins 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 241000282887 Suidae Species 0.000 description 1
- 238000010171 animal model Methods 0.000 description 1
- 238000005842 biochemical reaction Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 235000013365 dairy product Nutrition 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011143 downstream manufacturing Methods 0.000 description 1
- 239000000975 dye Substances 0.000 description 1
- 210000003746 feather Anatomy 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003100 immobilizing effect Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 108020004999 messenger RNA Proteins 0.000 description 1
- 238000002493 microarray Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000001259 photo etching Methods 0.000 description 1
- 102000054765 polymorphisms of proteins Human genes 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 230000001850 reproductive effect Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
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/6813—Hybridisation assays
- C12Q1/6834—Enzymatic or biochemical coupling of nucleic acids to a solid phase
- C12Q1/6837—Enzymatic or biochemical coupling of nucleic acids to a solid phase using probe arrays or probe chips
-
- 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
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
- C12N15/1034—Isolating an individual clone by screening libraries
- C12N15/1065—Preparation or screening of tagged libraries, e.g. tagged microorganisms by STM-mutagenesis, tagged polynucleotides, gene tags
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/11—DNA or RNA fragments; Modified forms thereof; Non-coding nucleic acids having a biological activity
-
- 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/6813—Hybridisation assays
- C12Q1/6827—Hybridisation assays for detection of mutation or polymorphism
-
- 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
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Genetics & Genomics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biotechnology (AREA)
- Molecular Biology (AREA)
- Analytical Chemistry (AREA)
- Biophysics (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Immunology (AREA)
- Plant Pathology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
提供了一种蛋鸡全基因组SNP芯片及其应用,该芯片上的SNP位点分别源于中国主要蛋鸡各品系共有的14624个SNP位点;与蛋鸡抗病性状相关的3677个SNP位点;与蛋鸡经济性状相关联的1.6万个SNP位点;弥补前述探针未覆盖的基因组区域的9358个SNP位点。该蛋鸡全基因组SNP芯片上的43681个SNP具有SEQ ID NO.1‑43681所示的DNA序列,该芯片可以特异性地鉴定商业化蛋鸡与地方品种蛋鸡的亲缘关系,还可以进行全基因组关联分析、全基因组选择育种,以及目标性状QTL定位分析和群体遗传学分析等应用。
Description
技术领域
本发明涉及分子生物学、功能基因组学、生物信息学和基因组育种领域,更具体涉及一种鸡全基因SNP芯片及其用途。
背景技术
SNP是指基因组水平上由单个核苷酸碱基的改变而导致核酸序列的多态性,其变异形式包括转换、颠换、缺失和插入等,具有数量多、分布广泛和易检测等特点。SNP作为一种遗传标记分子,对质量性状、经济性状以及一些复杂的疾病性状的遗传变异具有贡献,因此被广泛应用于遗传研究中。SNP芯片是将大量带有荧光标记的DNA探针固定在玻璃基片上通过光蚀刻而成,其原理是通过探针DNA链与检测目标基因组DNA进行杂交、结合,再经染色剂染色后利用激光扫描从而对SNP进行基因分型。SNP生物芯片的优势在于检测通量很大,一次可以检测几十万甚至几百万个SNP位点,其检测的准确性也很高。
全基因组SNP芯片即单核苷酸多态性(SNP)微阵列(Array),将数百万DNA标记序列排列在玻片或特殊硅片上,固定形成SNP探针阵列。它的工作原理是通过固定在芯片上的DNA标记序列与目标基因组发生碱基配对反应,从而精准鉴定基因信息。基因芯片运用先进的微生化反应技术、显微标记技术、微米级扫描分辨成像技术和生物信息计算机处理技术,准确鉴定检测生物个体中特定的mRNA或DNA序列。鸡的商业化SNP芯片主要有两款,Illumina公司的iSelect chicken 60K SNP芯片(已经下架,归美国科宝公司所有)和Affymetrix公司的600K高密度芯片。目前市场上仅售的Affymetrix 600K SNP芯片富含约56万个SNP,其位点主要来源于国外的商业化肉鸡品系、蛋鸡品系和一些近交系,相比之前Illumina 60K芯片,该芯片密度有了很大的提高,但是这些位点在中国地方品种中的应用却没有达到理想的效果。申请人前期的检测结果(约2,600个蛋鸡个体)发现Affymetrix600K高密度SNP芯片上约一半的位点多态性较差,再加上昂贵的成本也进一步阻碍了它在我国大规模群体中的实际应用。
畜禽全基因组SNP芯片由于标记密度高,均匀覆盖全基因组、具有测定准确率高、易实现标准化和自动化检测等诸多优势,可以应用于分子遗传研究和分子辅助育种等各个方面,可为育种上下游过程提供支持,包括种质资源遗传多样性分析、亲缘关系分析、全基因组关联分析、数量性状QTL定位分析和选择进化研究等。其中非常重要的应用是在全基因组选择育种(Genomic Selection)中的应用。全基因组选择是当前畜禽育种中主要选育技术,运用全基因组标记,估计所有可能的遗传效应,解释全部的遗传变异,通过统计标记的遗传效应,进行全基因组育种值的预测,具有育种值估计准确性高、遗传进展快的优点,已应用于奶牛、生猪、国外高产蛋鸡和快大型白羽肉鸡的商业化品系选育实践。
发明内容
本发明的目的在于提供一种蛋鸡全基因组SNP芯片及其应用。
本发明首先提供一种用于蛋鸡全基因组SNP分子标记组合,由43681个SNP分子标记组成,其核苷酸序列分别如SEQ ID NO.1~43681所示。
本发明提供一种蛋鸡全基因组SNP芯片,包含43681个SNP分子标记,其具有SEQ IDNO.1~43681所示的核苷酸序列。
本发明提供蛋鸡全基因组SNP芯片中,SNP位点包括五类探针:第一类是从中国国内7个品系(洛岛红Line1,洛岛白Line2、Line3,白来航Line4、Line5,贵妃鸡,东乡绿壳蛋鸡)共479个个体重测序数据筛选出的多态性良好的SNP位点,包含14,624个位点;第二类探针涉及申请人前期对蛋鸡抗病性状研究的发现的功能位点,包含3,677个SNP;第三类探针包含申请人前期通过全基因组关联分析(GWAS)技术挖掘出与蛋鸡重要经济性状(产蛋性能、蛋品质、饲料效率等)相关的显著位点,共2,677个SNP;第四类探针是根据已有的相关候选基因的文献和QTL数据库报道的跟蛋鸡经济性状相关的位点,共13,345个SNP;第五类是从SNP database数据库(ftp://ftp.ncbi.nlm.nih.gov/snp/organisms/archive/chicken_9031/)补充前几类未能覆盖的区域,包含9,358个SNP位点。以上五类探针共计包含43,681个SNP位点。
本发明的上述蛋鸡全基因组SNP分子标记组合及上述的蛋鸡全基因组SNP芯片中,所述分子标记的SNP位点位于SEQ ID NO.1~43681所示的核苷酸序列的第61位。
所述分子标记与蛋鸡抗病性状相关联、或与蛋鸡经济性状相关联;所述抗病为抗马立克氏病;所述蛋鸡经济性状为产蛋性能、产蛋品质、饲料报酬、免疫力。
本发明还提供了上述蛋鸡全基因组SNP芯片在检测蛋鸡DNA样品中的应用,该芯片可以对中外蛋鸡品种资源进行种质资源评价、对地方蛋鸡和商业化蛋鸡进行全基因组选择育种、目标性状QTL、关联位点及候选基因鉴定、亲缘关系鉴定等。
具体地,本发明提供了上述蛋鸡全基因组SNP分子标记组合和/或上述的蛋鸡全基因组SNP芯片在鸡种质资源改良中的应用。
本发明提供了上述蛋鸡全基因组SNP分子标记组合和/或上述的蛋鸡全基因组SNP芯片在鸡品种鉴定中的应用。
本发明提供了上述蛋鸡全基因组SNP分子标记组合和/或上述的蛋鸡全基因组SNP芯片在蛋鸡育种中的应用。
本发明提供了上述蛋鸡全基因组SNP分子标记组合和/或上述的蛋鸡全基因组SNP芯片在商品化蛋鸡与地方品种蛋鸡亲缘关系鉴定中的应用。
本发明提供了上述蛋鸡全基因组SNP分子标记组合和/或上述的蛋鸡全基因组SNP芯片在蛋鸡遗传多样性分析中的应用。
本发明提供了上述鸡全基因组SNP分子标记组合和/或上述的蛋鸡全基因组SNP芯片在蛋鸡目标性状QTL定位分析中的应用。
本发明蛋鸡全基因组SNP芯片上所有SNP标记根据SNP数据库名称命名,其核苷酸序列来源于鸡的5版本参考基因组(Gallus_gallus-5.0,2016)。
本发明的优点在于:1).优化设计:全基因组各染色体均匀分布,覆盖度高,保证了蛋鸡基因组遗传评估的准确性;2).功能相关性:该芯片收集了与蛋鸡产蛋性能、抗病性状、饲料利用效率和蛋品质等性状相关的显著位点,增加了芯片的实用性;3).适用性:芯片位点来源于自主选育的商业化高产蛋鸡品种和我国地方品种中多态性较高的SNP,既满足了科研工作者的科研需求也满足了蛋鸡企业进行全基因组分子育种的工作。
附图说明
图1为本申请的蛋鸡全基因组SNP芯片在鸡基因组分布情况,横坐标分别为20、40、60、80、100、120、140、160、180、200Mb,纵坐标由下到上分别表示chr1、chr2、chr3、chr4、chr5、chr6、chr7、chr8、chr9、chr10、chr11、chr12、chr13、chr14、chr15、chr16、chr17、chr18、chr19、chr20、chr21、chr22、chr23、chr24、chr25、chr26、chr27、chr28、chr33、chrLGE64和chrZ。
图2为本申请所有SNP位点在全基因组各染色体占有比例。
图3为各个品种鸡的主成分分析图。图上数字分别代表10个品种,1表示北京油鸡,2表示白来航鸡,3表示贵妃鸡,4表示东乡绿壳鸡,5表示丝羽乌骨鸡,6表示黄羽肉鸡,7表示洛岛白鸡,8表示藏鸡,9表示矮脚鸡,10表示洛岛红鸡。
图4为各个品种鸡的聚类分析图。枝条的大分叉代表每个品种都能独立区分开,小分叉代表每个个体。
图5为洛岛红鸡36周龄蛋重的曼哈顿图。
具体实施方式
以下实施例进一步说明本发明的内容,但不应理解为对本发明的限制。在不背离本发明精神和实质的情况下,对本发明方法、步骤或条件所作的修改或替换,均属于本发明的范围。
若未特别指明,实施例中所用的技术手段为本领域技术人员所熟知的常规手段。
实施例1蛋鸡全基因组SNP芯片的制备方法
本发明第一类探针获得:对中国国内自主选育的商业化品种以及地方品种蛋鸡进行全基因组重测序,其中包括洛岛红鸡、洛岛白鸡、白来航鸡、贵妃鸡和东乡绿壳蛋鸡等7个品系(洛岛白和白来航有两个品系)的479个个体(洛岛红Line1:92个,洛岛白Line2:74个,洛岛白Line3:91个,白来航Line4:60个,白来航Line5:69个,贵妃鸡:68个,东乡绿壳蛋鸡:25个),剔除每个品种MAF<0.05的SNP位点,同时筛选出各个品系共有的SNP以及以洛岛红鸡为主导的至少在其他一个品系中也共有的SNP位点,这样很大程度地保证了SNP位点的多态性和通用性。最终,各系共有的64,396个SNP以及以洛岛红鸡为主导的3,782,699个SNP位点作为本芯片的核心备选SNP。
本发明第二类探针获得:申请人前期对鸡的抗病性状研究,利用美国农业部下属的禽病与肿瘤实验室组建的抗马立克氏病的近交系白来航群体(其中抗病系和易感系各三个个体),通过全基因组重测序技术挖掘出各品系内存在的遗传变异,挑选出与抗病性状相关的系特异性SNP位点,共计1,554,188个,再通过与上一步所述的7个品系479个个体重测序得到的多态性较好的SNP进行overlapping,这样极大程度地保证了这些SNP位点在其他蛋鸡品种中的适用性,最终备选了480,341个抗病特异性的SNP。
本发明第三类探针获得:这类探针是根据申请人前期对蛋鸡重要经济性状的全基因组关联分析研究筛选得到,共两部分:第一部分是利用分别来自40个家系的385个白来航和361个矮脚蛋鸡对产蛋性能(产蛋数和蛋重)和蛋品质进行关联分析;第二部分是利用白来航和东乡绿壳蛋鸡正反交得到的F2代资源群体,对1,512个母鸡的产蛋性能,饲料利用效率,蛋品质和繁殖性能进行全基因组关联分析。以上获得的与蛋鸡重要经济性状相关联的SNP位点,再通过生物统计学方法分析去重后,得到2,677个功能性位点加入到芯片中。
本发明第四类探针获得:通过查找已发表文献的相关候选基因和在线QTLdatabase数据库(https://www.animalgenome.org/cgi-bin/QTLdb/GG/index)获得与鸡的蛋品质、产蛋性能、饲料效率和免疫相关的QTL,剔除置信区间较大和P值不显著(P>0.05)的QTL,最终得到了1,238个候选QTL区间,再利用重测序数据得到的SNP与候选基因和QTL进行overlapping,筛选出589,164个SNP作为性状相关的备选位点。
本发明第五类探针获得:通过NCBI(https://www.ncbi.nlm.nih.gov/snp)的鸡SNP database数据库下载得来。
SNP位点鉴定流程:考虑到鸡不同染色体片段大小差异较大,并且不同染色体上突变位点的重组率不一致,第一步先对每条染色体上的连锁不平衡水平进行检测,从而设置每条染色体的最优窗口/区间的大小,利用这些窗口/区间使得SNP在各条染色体上均匀分布,每个区间优先选择1个第一类各个品系共有的SNP探针,如果小于1个,则选择第二类探针中的SNP,第三类探针是强制放入芯片中的SNP,如果第一类探针和第二类探针选择完后,仍然不能覆盖的窗口,则选取第四类探针的SNP,第五类探针是补充上述几类探针不能覆盖的区间。
将鉴定好的SNP位点提交至Illumina公司,利用Infinium XT打分系统(http://www.illumina.com/)进行打分,将所有分值<0.8的位点去掉。对于删除的不合格SNP位点,选取距离其最近的SNP位点进行补充,并且再次进行打分。按照以上步骤鉴定和筛选,最后共获得50,000个标签SNP位点。由于Illumina Infinium XT平台利用这5万个SNP位点制作芯片时,存在成功转化率的问题,因此最终完成的蛋鸡全基因组50K SNP芯片上的有效SNP位点数为43,681个(见图1和图2)。
表1本发明蛋鸡全基因组SNP芯片位点在各染色体分布数量
实施例2基于本发明蛋鸡全基因组SNP芯片进行主成分和聚类分析
使用本发明提供的蛋鸡全基因组50K SNP芯片对20只洛岛红鸡(Rhode IslandRed)、20只洛岛白鸡(Rhode Island white)、19只北京油鸡(Beijing You)、19只白来航鸡(White Leghorn)、19只贵妃鸡(Houdan)、20只东乡绿壳蛋鸡(Dongxiang)、20只丝羽乌骨鸡(Silkie)、8只黄羽肉鸡(Yellow broiler)、20只藏鸡(Tibetan)和20只矮脚鸡(Dwarf)的DNA基因型进行检测,先剔除MAF小于0.01的SNP位点(641个),再分别利用PLINK和SNPhylo软件进行主成分和聚类分析。主成分分析和聚类进化树分析的结果见图3和图4所示,所有的品种都能单独分开,这也说明了本发明的芯片在不同品种中的适用性,其中我国地方品种的北京油鸡、东乡和丝羽乌骨鸡聚成一类,国外的商业品种洛岛白鸡、白来航和贵妃鸡聚成一类,特别注意的是洛岛红鸡单独聚成一类,且遗传距离相离其他品种较远,这也由于该品种已经经历了多个世代的高强度人工选择,基因组遗传结构以及核酸多态性发生了变化。上述的主成分和聚类分析结果都能够很好地证明芯片的适用性,其鉴定结果准确可靠。
实施例3本发明蛋鸡全基因组SNP芯片在全基因组关联分析中的应用
使用本发明提供的蛋鸡全基因组SNP芯片对北京市华都峪口禽业有限公司提供的680只洛岛红母鸡进行基因型检测,同时收集了这些母鸡36周龄的蛋重表型数据(EW36),SNP基因分型的质控标准:样本的检出率大于95%(sample Call rate>95%)、最小等位基因频率大于0.01(MAF>0.01)、哈迪温伯格平衡检验P值小于1e-06,质控后剩余673个个体和38,627个SNP用于接下来的全基因组关联(GWAS)分析。分析的模型采用线性混合模型,以SNP独立次数检验的P值8.1e-06和1.62e-04作为全基因组显著和潜在显著的阈值,GWAS结果的曼哈顿图见图5所示。结果显示,在一号染色体上存在一段与蛋重相关联的显著区域,这一结果也与之前研究报道的结果相一致(Liu et al.,2018,Frontiers in Genetics)。
实施例4本发明蛋鸡全基因组SNP芯片在基因组选择育种中的应用——遗传评估
使用本发明提供的蛋鸡全基因组SNP芯片进行基因组遗传评估。对北京市华都峪口禽业有限公司提供的某一品系三个世代的2,950只蛋鸡(285只公鸡,2,665只母鸡)的DNA进行SNP基因分型,同时使用基于系谱的选择方法和基于基因组的选择方法对这些个体达到28周龄的体重(BW28),母鸡28周龄的产蛋重(EW28)和38周龄的产蛋数(EN38)进行遗传评估。具体步骤如下:
(1)对2,950只鸡的分型检测结果进行质控,样本的检出率大于95%(sample Callrate>95%)、最小等位基因频率大于0.01(MAF>0.01)、哈迪温伯格平衡检验P值小于1e-06,质控后剩余2,950个个体和38,413个合格SNP用于全基因组选择分析;
(2)模型:BW28=Date+Hatch+Sex+Animal+Error
EW28/EN38=Date+Hatch+Animal+Error
√Date,Hatch,Sex分别为个体出生日期,孵化批次,以及性别的固定效应;
√Animal为动物个体的随机加性遗传效应,即个体育种值
√Error为随机残差
(3)分别使用基于系谱的选择方法(PBLUP模型)和基于基因组选择的方法(ssGBLUP:一步法基因组模型)进行遗传评估
PBLUP模型:常规育种值通过单性状动物模型估计:
y=Xb+Zu+e
其中y是表型值向量,X是固定效应设计矩阵,b是固定效应向量,包括世代(出生年)、孵化批次和性别,Z是随机效应关联矩阵,u是育种值向量并假设其服从其中A是分子血缘关系矩阵,是加性遗传方差,e是残差向量并假设其服从其中I是单位矩阵,是残差方差。
ssGBLUP模型:ssGBLUP模型在形式上与PBLUP模型相似,区别在于ssGBLUP模型中的随机效应向量服从其中H是将系谱信息和基因组信息整合的矩阵(Christensen and Lund,2010;Legarra,et al.,2009)。H矩阵构建如下:
下标1和2分别代表没有基因型的个体和有基因型的个体,H矩阵的逆矩阵在形式上相当简单(Aguilar,et al.,2010):
(4)结果:
表2不同遗传评估方法的准确性比较
(5)能够看出使用本发明芯片进行基因组选择相比传统选择方法,能够有效提高遗传评估的准确性。
工业实用性
本发明提供的蛋鸡全基因组SNP芯片上的SNP位点分别源于中国主要蛋鸡各品系共有的14624个SNP位点;与蛋鸡抗病性状相关的3677个SNP位点;与蛋鸡经济性状相关联的1.6万个SNP位点;弥补前述探针未覆盖的基因组区域的9358个SNP位点。本发明蛋鸡全基因组SNP芯片上的43681个SNP具有SEQ ID NO.1-43681所示的DNA序列,该芯片可以特异性地鉴定商业化蛋鸡与地方品种蛋鸡的亲缘关系,还可以进行全基因组关联分析、全基因组选择育种,以及目标性状QTL定位分析和群体遗传学分析等应用,具有中外鸡品种通用性,能够有助于加快蛋鸡行业的快速发展,具有很大的经济实用价值和科研价值。
Claims (8)
1.一种蛋鸡全基因组SNP芯片,其特征在于,由43681个SNP分子标记组成,其核苷酸序列分别如SEQ ID NO.1~43681所示;
所述分子标记的SNP位点位于SEQ ID NO.1~ 43681所示的核苷酸序列的第61位。
2.如权利要求1所述的蛋鸡全基因组SNP芯片,其特征在于,所述分子标记与蛋鸡抗病性状相关联、或与蛋鸡经济性状相关联;所述蛋鸡经济性状为产蛋性能、产蛋品质、饲料报酬、免疫力。
3.权利要求1所述的蛋鸡全基因组SNP芯片在商品化蛋鸡与地方品种蛋鸡亲缘关系鉴定中的应用。
4.权利要求1所述的蛋鸡全基因组SNP芯片在蛋鸡品种鉴定中的应用。
5.权利要求1所述的蛋鸡全基因组SNP芯片在蛋鸡育种中的应用。
6.权利要求1所述的蛋鸡全基因组SNP芯片在蛋鸡种质资源改良中的应用。
7.权利要求1所述的蛋鸡全基因组SNP芯片在蛋鸡群体遗传多样性分析中的应用。
8.权利要求1所述的蛋鸡全基因组SNP芯片在蛋鸡目标性状QTL定位分析中的应用。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2018/108741 WO2020062160A1 (zh) | 2018-09-29 | 2018-09-29 | 一种蛋鸡全基因组snp芯片及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113039288A CN113039288A (zh) | 2021-06-25 |
CN113039288B true CN113039288B (zh) | 2022-04-26 |
Family
ID=69952639
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880004561.5A Active CN113039288B (zh) | 2018-09-29 | 2018-09-29 | 一种蛋鸡全基因组snp芯片及其应用 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20210395803A1 (zh) |
CN (1) | CN113039288B (zh) |
WO (1) | WO2020062160A1 (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111690754B (zh) * | 2020-06-22 | 2022-09-09 | 浙江省农业科学院 | 一种鸭产蛋性能的全基因组关联分析方法 |
CN112725462B (zh) * | 2021-01-15 | 2022-07-05 | 广西大学 | 基于全基因组测序筛选的与龙胜凤鸡体重、体尺相关的snp分子标记组合及应用 |
CN115341035B (zh) * | 2021-05-12 | 2024-08-23 | 华中农业大学 | 用于选择母鸡产蛋蛋重的snp分子标记 |
CN113549699B (zh) * | 2021-07-01 | 2023-12-19 | 中国农业科学院北京畜牧兽医研究所 | 一种白羽肉鸡产蛋数的基因组选择方法 |
CN113667759B (zh) * | 2021-07-05 | 2022-06-14 | 江苏省家禽科学研究所 | 与鸡产蛋持续能力相关的snp遗传标记及其应用 |
CN113388687B (zh) * | 2021-08-06 | 2022-10-21 | 南昌师范学院 | 一种全血法鉴定东乡黑羽绿壳蛋鸡与黑羽乌鸡的方法 |
CN113699250B (zh) * | 2021-08-13 | 2023-11-03 | 中国农业科学院北京畜牧兽医研究所 | 与肉鸡饲料转化效率性状相关的分子标记及其应用 |
CN114480673B (zh) * | 2022-03-04 | 2022-12-23 | 江苏省家禽科学研究所 | 一种基于靶向捕获测序的鸡低密度snp液相芯片及其应用 |
CN115927647B (zh) * | 2022-07-13 | 2023-08-01 | 江苏省家禽科学研究所 | 与鸡蛋重曲线斜率相关的snp遗传标记及其应用 |
CN115261486B (zh) * | 2022-07-27 | 2023-06-23 | 中国农业科学院北京畜牧兽医研究所 | 一种华西牛全基因组选择育种芯片及其应用 |
CN116083597B (zh) * | 2022-12-01 | 2023-08-04 | 江苏省家禽科学研究所 | 影响鸡后期产蛋数的snp遗传标记及其应用 |
CN116072226A (zh) * | 2023-01-17 | 2023-05-05 | 中国农业大学 | 一种用于蛋鸡产蛋性状基因组选择的机器学习方法及系统 |
CN116377086B (zh) * | 2023-03-30 | 2024-03-15 | 山东省农业科学院家禽研究所(山东省无特定病原鸡研究中心) | 一种鸡全基因组低密度芯片及其制作方法和应用 |
CN116144796A (zh) * | 2023-04-06 | 2023-05-23 | 中国农业大学 | 一种与鸡蛋蛋壳重相关的snp分子模块及其应用 |
CN116515956A (zh) * | 2023-04-25 | 2023-08-01 | 重庆市畜牧科学院 | 基于gwas筛选鹅产蛋分子标记及选择方法 |
CN116452155A (zh) * | 2023-06-14 | 2023-07-18 | 内蒙古盛健农牧业工程技术研究有限公司 | 一种高产奶山羊快速育种的智能化构建系统 |
CN116790767B (zh) * | 2023-07-21 | 2024-02-20 | 华南农业大学 | 与鸡皮肤乌度相关的tyr基因分子标记及其应用 |
CN117070639B (zh) * | 2023-08-08 | 2024-02-02 | 江苏省家禽科学研究所 | 与鸡产蛋总重关联的snp遗传标记在鸡遗传育种中的应用 |
CN117187403B (zh) * | 2023-08-11 | 2024-05-17 | 江苏省家禽科学研究所 | 与鸡产蛋期剩余饲料消耗相关的snp遗传标记在鸡遗传育种中的应用 |
CN117646075A (zh) * | 2023-12-04 | 2024-03-05 | 江苏省家禽科学研究所 | 一种鸡繁殖性能分子微模块芯片及其应用 |
CN118222719B (zh) * | 2024-03-04 | 2024-09-13 | 江苏省家禽科学研究所 | Dis3l2基因中与鸡蛋长径关联的遗传标记在鸡遗传育种中的应用 |
CN118222720A (zh) * | 2024-03-04 | 2024-06-21 | 江苏省家禽科学研究所 | 影响鸡蛋黄比率的snp遗传标记在鸡遗传育种中的应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105349637A (zh) * | 2015-10-26 | 2016-02-24 | 中国农业大学 | 一种与鸡蛋蛋重相关的snp分子标记及其应用 |
CN107338317A (zh) * | 2017-08-24 | 2017-11-10 | 中国农业大学 | 一种与鸡蛋蛋壳强度性状相关的分子标记及其应用 |
CN107354220A (zh) * | 2017-08-24 | 2017-11-17 | 中国农业大学 | 一种与鸡产蛋数性状相关的分子标记及其应用 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105238859B (zh) * | 2015-10-13 | 2018-03-16 | 中国农业大学 | 一种获取鸡全基因组高密度snp标记位点的方法 |
-
2018
- 2018-09-29 US US17/281,223 patent/US20210395803A1/en active Pending
- 2018-09-29 WO PCT/CN2018/108741 patent/WO2020062160A1/zh active Application Filing
- 2018-09-29 CN CN201880004561.5A patent/CN113039288B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105349637A (zh) * | 2015-10-26 | 2016-02-24 | 中国农业大学 | 一种与鸡蛋蛋重相关的snp分子标记及其应用 |
CN107338317A (zh) * | 2017-08-24 | 2017-11-10 | 中国农业大学 | 一种与鸡蛋蛋壳强度性状相关的分子标记及其应用 |
CN107354220A (zh) * | 2017-08-24 | 2017-11-17 | 中国农业大学 | 一种与鸡产蛋数性状相关的分子标记及其应用 |
Also Published As
Publication number | Publication date |
---|---|
CN113039288A (zh) | 2021-06-25 |
WO2020062160A1 (zh) | 2020-04-02 |
US20210395803A1 (en) | 2021-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113039288B (zh) | 一种蛋鸡全基因组snp芯片及其应用 | |
CN110191965B (zh) | 猪全基因组50k snp芯片及应用 | |
US11578365B2 (en) | Chicken whole-genome SNP chip and use thereof | |
Fan et al. | Development and application of high-density SNP arrays in genomic studies of domestic animals | |
Wang et al. | High-density genetic linkage mapping in turbot (Scophthalmus maximus L.) based on SNP markers and major sex-and growth-related regions detection | |
CN107967409B (zh) | 一种猪全基因组低密度snp芯片及其制作方法和应用 | |
CA2673174C (en) | Artificial selection method and reagents | |
CN111278994B (zh) | 一种牙鲆抗病育种基因芯片及其应用 | |
CN110578008B (zh) | 一种预示和鉴定鸡腹脂重的分子标记方法 | |
CN109371143B (zh) | 与猪生长性状相关联的snp分子标记 | |
WO2023201950A1 (zh) | 一种用于北京黑猪基因分型的snp分子标记组合、芯片及其制备方法与应用 | |
US12054769B2 (en) | Sexual maturation in rainbow trout | |
CN109371144B (zh) | 一种与猪生长性状相关联的snp分子标记 | |
Lillie et al. | Genomic signatures of 60 years of bidirectional selection for 8-week body weight in chickens | |
CN114736974B (zh) | 与母猪产程性状相关的snp分子标记及其应用 | |
Liu et al. | Design and evaluation of a custom 50K Infinium SNP array for egg-type chickens | |
CN114736975B (zh) | 与母猪产仔数性状相关的snp分子标记及其应用 | |
Peripolli et al. | Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data | |
CN111926086A (zh) | 一种影响鸡体斜长的分子标记及其应用 | |
Bernard et al. | Development of a high-density 665 K SNP array for rainbow trout genome-wide genotyping | |
CN118480613A (zh) | 一种猪全基因组低密度5k snp芯片及应用 | |
CN110295236A (zh) | 猪饲料转化率的snp分子遗传标记 | |
CN115261486A (zh) | 一种华西牛全基因组选择育种芯片及其应用 | |
CN114317772A (zh) | 一种用于鉴定鸡高产性能的分子标记及其应用 | |
Goraga | Quantitative trait loci (QTL) and genetic parameters for economically important traits in chicken—a review |
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 |