CN116790798A - 与花生含油量相关主效数量性状基因位点qCOA08_1及其应用 - Google Patents
与花生含油量相关主效数量性状基因位点qCOA08_1及其应用 Download PDFInfo
- Publication number
- CN116790798A CN116790798A CN202310859370.XA CN202310859370A CN116790798A CN 116790798 A CN116790798 A CN 116790798A CN 202310859370 A CN202310859370 A CN 202310859370A CN 116790798 A CN116790798 A CN 116790798A
- Authority
- CN
- China
- Prior art keywords
- dcaps
- oil content
- qcoa08
- seq
- peanut
- 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.)
- Pending
Links
- 235000019483 Peanut oil Nutrition 0.000 title claims abstract description 30
- 239000000312 peanut oil Substances 0.000 title claims abstract description 30
- 235000020232 peanut Nutrition 0.000 claims abstract description 41
- 235000017060 Arachis glabrata Nutrition 0.000 claims abstract description 32
- 235000010777 Arachis hypogaea Nutrition 0.000 claims abstract description 32
- 235000018262 Arachis monticola Nutrition 0.000 claims abstract description 32
- 239000003147 molecular marker Substances 0.000 claims abstract description 25
- 238000009395 breeding Methods 0.000 claims abstract description 23
- 230000001488 breeding effect Effects 0.000 claims abstract description 23
- 241001553178 Arachis glabrata Species 0.000 claims abstract 5
- 239000003921 oil Substances 0.000 claims description 49
- 235000019198 oils Nutrition 0.000 claims description 47
- 230000000694 effects Effects 0.000 claims description 14
- 239000002773 nucleotide Substances 0.000 claims description 8
- 125000003729 nucleotide group Chemical group 0.000 claims description 8
- 239000003550 marker Substances 0.000 claims description 7
- 108091028043 Nucleic acid sequence Proteins 0.000 claims description 6
- 238000003776 cleavage reaction Methods 0.000 claims description 6
- 230000007017 scission Effects 0.000 claims description 6
- 238000011144 upstream manufacturing Methods 0.000 claims description 5
- 108090000623 proteins and genes Proteins 0.000 description 52
- 244000105624 Arachis hypogaea Species 0.000 description 37
- 230000002068 genetic effect Effects 0.000 description 16
- 230000014509 gene expression Effects 0.000 description 12
- 101000988395 Homo sapiens PDZ and LIM domain protein 4 Proteins 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 8
- 210000000349 chromosome Anatomy 0.000 description 8
- 238000000034 method Methods 0.000 description 8
- 238000012216 screening Methods 0.000 description 6
- 101150100400 66 gene Proteins 0.000 description 4
- 108700028369 Alleles Proteins 0.000 description 4
- 108700039691 Genetic Promoter Regions Proteins 0.000 description 4
- 108700005075 Regulator Genes Proteins 0.000 description 4
- 108091023040 Transcription factor Proteins 0.000 description 4
- 102000040945 Transcription factor Human genes 0.000 description 4
- 238000010367 cloning Methods 0.000 description 4
- 238000012217 deletion Methods 0.000 description 4
- 230000037430 deletion Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 235000014113 dietary fatty acids Nutrition 0.000 description 4
- 229930195729 fatty acid Natural products 0.000 description 4
- 239000000194 fatty acid Substances 0.000 description 4
- 150000004665 fatty acids Chemical class 0.000 description 4
- 230000035772 mutation Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000011222 transcriptome analysis Methods 0.000 description 4
- 101100049737 Arabidopsis thaliana WRI1 gene Proteins 0.000 description 3
- 108091026890 Coding region Proteins 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 3
- 108700026244 Open Reading Frames Proteins 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000037433 frameshift Effects 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 230000004807 localization Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000001127 nanoimprint lithography Methods 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- 108020003589 5' Untranslated Regions Proteins 0.000 description 2
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 description 2
- 239000005977 Ethylene Substances 0.000 description 2
- 244000068988 Glycine max Species 0.000 description 2
- 235000010469 Glycine max Nutrition 0.000 description 2
- 108091092195 Intron Proteins 0.000 description 2
- 108091036066 Three prime untranslated region Proteins 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000008157 edible vegetable oil Substances 0.000 description 2
- 239000012636 effector Substances 0.000 description 2
- 230000001747 exhibiting effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012252 genetic analysis Methods 0.000 description 2
- 230000007614 genetic variation Effects 0.000 description 2
- 238000012253 re-sequencing analysis Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- WRIDQFICGBMAFQ-UHFFFAOYSA-N (E)-8-Octadecenoic acid Natural products CCCCCCCCCC=CCCCCCCC(O)=O WRIDQFICGBMAFQ-UHFFFAOYSA-N 0.000 description 1
- LQJBNNIYVWPHFW-UHFFFAOYSA-N 20:1omega9c fatty acid Natural products CCCCCCCCCCC=CCCCCCCCC(O)=O LQJBNNIYVWPHFW-UHFFFAOYSA-N 0.000 description 1
- QSBYPNXLFMSGKH-UHFFFAOYSA-N 9-Heptadecensaeure Natural products CCCCCCCC=CCCCCCCCC(O)=O QSBYPNXLFMSGKH-UHFFFAOYSA-N 0.000 description 1
- 229920000742 Cotton Polymers 0.000 description 1
- 108091029795 Intergenic region Proteins 0.000 description 1
- 239000005642 Oleic acid Substances 0.000 description 1
- ZQPPMHVWECSIRJ-UHFFFAOYSA-N Oleic acid Natural products CCCCCCCCC=CCCCCCCCC(O)=O ZQPPMHVWECSIRJ-UHFFFAOYSA-N 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000000692 Student's t-test Methods 0.000 description 1
- 238000000540 analysis of variance Methods 0.000 description 1
- 238000003975 animal breeding Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- 230000007849 functional defect Effects 0.000 description 1
- 238000003205 genotyping method Methods 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- QXJSBBXBKPUZAA-UHFFFAOYSA-N isooleic acid Natural products CCCCCCCC=CCCCCCCCCC(O)=O QXJSBBXBKPUZAA-UHFFFAOYSA-N 0.000 description 1
- 230000006372 lipid accumulation Effects 0.000 description 1
- 230000037356 lipid metabolism Effects 0.000 description 1
- 230000008774 maternal effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- ZQPPMHVWECSIRJ-KTKRTIGZSA-N oleic acid Chemical compound CCCCCCCC\C=C/CCCCCCCC(O)=O ZQPPMHVWECSIRJ-KTKRTIGZSA-N 0.000 description 1
- 230000008775 paternal effect Effects 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 230000022532 regulation of transcription, DNA-dependent Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 230000002103 transcriptional effect Effects 0.000 description 1
- 235000015112 vegetable and seed oil Nutrition 0.000 description 1
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
- C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/30—Detection of binding sites or motifs
-
- 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/13—Plant traits
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Analytical Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Organic Chemistry (AREA)
- Genetics & Genomics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Mycology (AREA)
- General Engineering & Computer Science (AREA)
- Botany (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一组与花生含油量相关主效数量性状基因位点及其应用,主效QTL命名为qCOA08_1,qCOA08_1位于花生基因组的第8个连锁群中,其标记区间为Chr.08_38056541‑Chr.08_40978133。本发明提供的技术能够用于花生含油量相关的分子标记辅助育种MAS。
Description
本申请为《与花生含油量相关主效数量性状基因位点及其应用》202211096080.6的分案申请。
技术领域
本发明涉及分子生物学技术领域,尤其是一种与花生含油量相关主效数量性状基因位点及应用分析。
背景技术
花生作为人类优质食用油的重要来源,已在100多个国家广泛种植。根据FAOSTAT(http://faostat3.fao.org/)公布的数据,全球花生荚果的年产量已经超过4400万吨,而中国占全球花生产量的38%,是世界上最大的花生生产国。即便如此,国内的花生产量仍然难以满足日益增长的食用花生油需求,而且由于中国现有耕地有限,通过增加耕地来增加花生产量是不现实的。据报道,花生含油量每增加1%,油脂加工的经济效益就会增加7%(Shasidhar et al.2017)。因此,提高花生单产和含油量以缓解不断增长的需求一直是全球花生育种计划的主要目标。
然而,几乎所有作物的产量和品质性状都是复杂的,数量性状受到遗传和环境因素的很大影响(Baring et al.2013;Wilson等,2013;Cui等,2020;Chen等人2019)。例如,Chen等人(2019)估计了多种环境下的几种产量相关性状。供试产量相关性状的遗传力在0.34~0.93之间,方差分析结果显示RIL群体间、环境间及RIL×环境互作均存在显著差异。在某些情况下,这些连续的数量性状和交替株系之间的边际变异会使田间试验中很难区分表现最好的株系。因此,弄清遗传力和遗传方式,有利于育种人员制定育种策略。目前,油脂含量的测定需要比较困难的过程,只有在收获后才能进行。在这种情况下,传统的提高含油量的育种策略被认为是一项成本高、效率低、耗时长的工作。在现代育种中,通过标记辅助选择(marker-assisted selection,MAS)可以在早期快速从种质资源中选择感兴趣的性状,且不受复杂环境的影响。MAS在培育具有一些复杂数量性状的优良作物方面取得了许多成功,如棉花、水稻和普通大豆(Tian et al.2021;Hulsekemp et al.2015;Shi等。2021)以及动物育种中(Chen等。2022)。然而,MAS的建立必须在特定的、有用的基因或QTL定位成功的条件下进行。
为了实现MAS育种目标,育种家们一直在努力寻找有价值的QTL。在世界各地的许多情况下,确定花生油含量的QTL并开展有针对性的育种工作,将是旨在确保人类食用油供应的项目的重要组成部分。遗憾的是,尽管育种家已经认识到MAS育种的重要性,但挖掘遗传资源提高花生油含量的研究明显滞后于大豆等其他油料作物,尤其是功能基因的克隆。由于花生品种缺乏相关的有利基因资源,利用MAS选育高油花生品种仍难以实现。因此,探索遗传资源,进一步通过图位克隆方式获得增加籽粒含油量有利基因,可以促进花生MAS育种,加快育种进程。
发明内容
本发明要解决的技术问题是提供一种一种与花生含油量相关主效数量性状基因位点及其应用。
为解决上述技术问题,本发明所采取的技术方案如下。
一种筛选与花生含油量相关的基因的方法,利用全基因组序列分析筛选候选基因,首先对两个亲本花生品种冀花5号JH5和开选01-6KX01-6进行基因组变异的比较,经过筛选在JH5和KX01-6两个候选区域之间识别出高质量的可变SNPs和Indels,去除分布在内含子、基因间区域和编码区下游及其他不能导致功能缺失的变异,保留分布在5'-UTR、3'-UTR、上游、剪接和外显子编码区的变异,所得变异对应的基因确定为候选基因。
一种筛选与花生含油量相关的基因的方法,通过转录组分析筛选候选基因,检测qCOA08_1和qCOA08_2两个候选区间中所有基因在油脂代谢关键期的表达模式,候选区间的全部基因作为候选基因,首先将在两个亲本花生品种冀花5号JH5和开选01-6KX01-6中均未表达的候选基因去除,这些基因不作为qCOA08_1和qCOA08_2的效应基因,剩余候选基因根据倍数变化和P值检测差异表达,检测标准为:倍数变化为log2fc≥1或≤-1,P指为P<0.001;即可得到qCOA08_1和qCOA08_2候选区间为数不多的差异表达基因DEGs,然后从所得差异表达基因DEGs中筛选强候选基因,经过试验验证或文献验证得到最终的候选基因或候选基因小组。
作为本发明的一种优选技术方案,所述强候选基因为花生Ah.UVRI66基因,其编码乙烯转录因子Wrinkled 1WRI1。
作为本发明的一种优选技术方案,所述Ah.UVRI66基因的FPKM值在低含油量品种KX01-6中为4.65,而Ah.UVRI66在高含油量品种JH5中表达量为22.03,为低油品种的5倍。
作为本发明的一种优选技术方案,所述Ah.UVRI66基因表达差异的分子来源在于启动子区主导的基因表达模式变异。
作为本发明的一种优选技术方案,所述Ah.UVRI66的启动子区在748和749核苷酸之间存在一个6bp ATATAT的插入,对应花生亲本之间含油量表达差异的分子来源或候选分子来源。
一种与花生含油量相关的主效QTL,所述主效QTL命名为qCOA08_1和/或qCOA08_2,所述qCOA08_1位于花生基因组的第8个连锁群中,其标记区间为Chr.08_38056541-Chr.08_40978133;所述qCOA08_2位于花生基因组的第8个连锁群中,其标记区间为Chr.08_44887537-Chr.08_46776535。
所述主效QTL的用途,用于开发与花生含油量相关的分子标记及其引物。
所述主效QTL的用途,用于花生含油量相关的分子标记辅助育种。
主效QTL的用途,基于所述主效QTL开发紧密连锁分子标记,通过专用的dCAPs标记识别所述紧密连锁分子标记,由此识别待测花生植株是否为对应的基因型,实现高含油量花生的分子标记辅助育种。
用于花生含油量性状识别的分子标记,为关联于主效QTL qCOA08_1下游侧翼序列的分子标记Chr.8_46776535,其核苷酸序列依次如SEQ ID NO.1所示。
用于鉴定或辅助鉴定花生含油量性状的dCAPs引物,包括:对应于分子标记Chr.8_46776535的dCAPs-F和dCAPs-R引物对,其核苷酸序列依次如SEQ ID NO.2和SEQ ID NO.3所示;对应的酶切位点为EcoRI,对应的识别序列为GAATTC。
用于花生含油量性状识别的分子标记,为关联于主效QTL qCOA08_1上游侧翼序列的分子标记Chr.8_44887537,其核苷酸序列依次如SEQ ID NO.4所示。
用于鉴定或辅助鉴定花生含油量性状的dCAPs引物,包括:对应于分子标记Chr.8_44887537的dCAPs-F和dCAPs-R引物对,其核苷酸序列依次如SEQ ID NO.5和SEQ ID NO.6所示;对应的酶切位点为BsrBI,对应的识别序列为GAGCGG。
用于花生含油量性状识别的分子标记,为关联于主效QTL qCOA08_2上游侧翼序列的分子标记Chr.8_38056541,其核苷酸序列依次如SEQ ID NO.7所示。
用于鉴定或辅助鉴定花生含油量性状的dCAPs引物,包括:对应于分子标记Chr.8_38056541的dCAPs-F和dCAPs-R引物对,其核苷酸序列依次如SEQ ID NO.8和SEQ ID NO.9所示;对应的酶切位点为FokI,对应的识别序列为GGATG。
用于花生含油量性状识别的分子标记,为关联于主效QTL qCOA08_2下游侧翼序列的分子标记Chr.8_40978133,其核苷酸序列依次如SEQ ID NO.10所示。
用于鉴定或辅助鉴定花生含油量性状的dCAPs引物,包括:对应于分子标记Chr.8_44887537的dCAPs-F和dCAPs-R引物对,其核苷酸序列依次如SEQ ID NO.11和SEQ ID NO.12所示;对应的酶切位点为HindIII,对应的识别序列为AAGCTT。
用于鉴定或辅助鉴定花生含油量性状的试剂盒,通过检测花生主效QTL qCOA08_1或qCOA08_2的多态性,用于花生含油量相关的分子标记辅助育种。
作为本发明的一种优选技术方案,所述试剂盒包含序列表中SEQ ID NO.2和SEQID NO.3组成的dCAPs引物对,和/或序列表中SEQ ID NO.5和SEQ ID NO.6组成的dCAPs引物对,和/或序列表中SEQ ID NO.8和SEQ ID NO.9组成的dCAPs引物对,和/或序列表中SEQ IDNO.11和SEQ ID NO.12组成的dCAPs引物对
所述的试剂盒、引物或分子标记,在花生分子育种中用于鉴定或辅助鉴定花生的含油量性状。
采用上述技术方案所产生的有益效果在于:本发明基于RIL群体构建了包含2840个多态性SNP标记的高分辨率遗传图谱,并在6个环境中稳定检测到2个效应较大的有价值QTL位点qCOA08_1和qCOA08_2。根据侧翼标记物的物理位置,qCOA08_1和qCOA08_2的物理间隔被定位在~2.9Mb和~1.7Mb之间。同时进行全基因组重测序和转录组分析,将候选基因的数量限定在40个以内。而一个编码WRI1转录因子并在亲本间表现出不同表达模式的强候选基因。本发明为花生优异等位基因的克隆提供了实用基础,与qCOA08_1和qCOA08_2紧密连锁的多态SNP标记或开发的dCAPs标记有助于加快MAS的育种进程。
附图说明
图1为6种环境下RIL种群油分的表型分布及差异分析相关图谱。
图2为20条染色体上多态性SNPs标记的物理位置(A)及对应的遗传图谱(B)。
图3为两个候选区变异的基因组分布图(a)及编码区突变类型图(b)。
图4为野外条件下通过NIL分析确认花生高油基因1-AhyHOF1的定位和功能;其中,测试性状的值以平均SE(n=6株植物)给出,星号表示Student’s t检验中p=0.05(*)、p=0.01(**)和p=0.001(**)水平的NILs之间存在显著差异,“ns”表示无显著差异。
具体实施方式
以下实施例详细说明了本发明。本发明所使用的各种原料及各项设备均为常规市售产品,均能够通过市场购买直接获得。在以下实施例的描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本申请实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本申请的描述。
如在本申请说明书和所附权利要求书中所使用的那样,术语“如果”可以依据上下文被解释为“当...时”或“一旦”或“响应于确定”或“响应于检测到”。类似地,短语“如果确定”或“如果检测到[所描述条件或事件]”可以依据上下文被解释为意指“一旦确定”或“响应于确定”或“一旦检测到[所描述条件或事件]”或“响应于检测到[所描述条件或事件]”。另外,在本申请说明书和所附权利要求书的描述中,术语“第一”、“第二”、“第三”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
实施例1、自然环境下RILs群体的表型评价
为挖掘能够提高花生含油量的基因资源,选用含油量表型差异较大的花生品种冀花5号(JH5)和开选01-6(KX01-6)构建重组自交系(RIL)(Table.1)。通过广义遗传力(h2 b)和含油量的分布,确定RIL群体是否适合进行油分含量的遗传分析。由表1可知,6种环境下的含油量均值大多数位于亲本均值之间,而RILs的最小值和最大值均超过亲本均值的极值,这表明JH5和KX01-6之间存在广泛的遗传变异。此外,6个环境中含油量的h2 b值较高,为0.81(表1),RILs直方图结果显示,含油量存在广泛的超亲遗传,这种超亲遗传力主要来源于遗传重组,这在QTL分析中是有利的。此外,尽管试验站环境对油分含量有显著影响,3502农场(NC)试验站油分含量平均值显著高于堤上(DS)试验站(图1),但峰度和偏度绝对值在0.02~0.66之间,表明RILs油分含量呈连续近似正态分布,符合典型的数量性状特征。综上所述,田间试验的含油量观察结果表明,由JH5和KX01-6杂交产生的RIL群体适合进行含油量的QTL分析。
表1.自然环境下含油性的遗传学分析
Note:DS,Dishang experiment station;NC,Nongchang experiment station;
实施例2、高分辨率遗传图谱构建
为了挖掘高含油量自交系中有利的等位基因,采用Affymetrix基因分型芯片对192个RILs家系进行了基因分型(Clevenger et.al,2017)。经过严格筛选,共鉴定出2840个多态性SNP标记,并将这些保留的标记定位到花生20条染色体上,构建了高分辨率的遗传图谱(图2)。结果显示,构建的遗传连锁图谱覆盖3706.382cM遗传距离,每个连锁群平均长度为185.32cM。在20个连锁群中,A07的LG Chr07连锁群长度最长为263.75cM,A02的LG Chr02连锁群最短,仅为77.34cM。每组平均有142个多态性标记,其中A03号染色体有309个多态性标记最多,A01号染色体有26个多态性标记最少。图谱相邻多态性标记之间的平均距离为1.31cM,在A02染色体上的平均距离为0.28cM,在B09染色体上的平均距离为4.44cM(表2)。为了检验构建的遗传图谱的准确性,将AhFAD2A和AhFAD2B基因控制的油酸含量性状(Varshney et al.2014;Kamdar等人,2021;Chi et al.2011)也进行了定位。结果表明,在09号染色体和19号染色体上分别定位到qFAD09和qFAD19两个主要位点,可解释54.4%的表型变异,LOD分别为14.54和11.80。AhFAD2A和AhFAD2B分别落在qFAD09和qFAD19的95%置信区间内,这强烈说明AhFAD2A和AhFAD2B是qFAD09和qFAD19的调控基因。这一结果验证了本研究构建的连锁图谱的准确性和实用性。
表2.用于鉴定与含油量相关的QTL的重组近交系(RIL)群体的SNP标记特征总结。
实施例3、自然环境下含油量QTL的鉴定
为进一步提高高油品种的标记辅助选择(MAS)的效率,利用2个试验站3年的RILs数据进行QTL分析。结果表明,在6个试验环境下,共检测到14个显著QTL,对192个F9:11花生RILs的表型变异贡献率为5.9%~20.7%,其LOD值为2.50~9.61。但除qCOA18_1和qCOA18_2外,大部分qtl不能在两个以上环境中同时检测到,这表明油分含量是由许多有效的次要QTL控制的,容易受到不同环境的影响。另外,由于qCOA08_1和qCOA08_2在5个或5个以上的环境中都能同时被检测到,其遗传变异最高值分别为16.1%和20.7%,增效等位基因均来自高油分品种冀花5号。qCOA08_1和qCOA08_2的最高LOD值分别为7.28和9.61,其峰值LOD值分别位于Chr.08_38056541至Chr.08_40978133和Chr.08_44936784至Chr.08_46776535的物理间隔内。总之,我们认为这两个QTL是两个有价值的主效应遗传位点,可用于提高花生油含量,但需要进一步确定更精确的物理区间或基因。
实施例4、利用全基因组序列确定候选基因
为了进一步筛选qCOA08_1和qCOA08_2的候选基因,对两个亲本进行了基因组变异的比较。经过严格筛选,在JH5和KX01-6两个候选区域之间共识别出912个高质量的可变SNPs和Indels(Table.S1)。其中超过87%(796)的变异分布在内含子、基因间区域和编码区下游,不太可能导致功能缺失,而分别有0.88%(8)、0.99%(9)、7.13%(65)、0.11%(1)和3.62%(33)分布在5'-UTR、3'-UTR、上游、剪接和外显子编码区(图3a)。在外显子编码区域的变异中,分别有18个和8个SNPs变异引起了非同义或同义突变,而其他7个SNPs或Indels变异引起了显著的结构突变,包括移码删除、移码插入、非移码删除、提前终止和终止子缺失,并可能导致功能缺陷(图3b)。进一步分析表明,所有这些变异共发生在24个预测基因中(表4),包括2个未知蛋白和22个根据Tifrunner v1.0注释基因组注释的基因。目前没有发现任何公开的信息可以支持这些注释基因之一可以增加油合成。因此假设qCOA08_1和qCOA08_2位点的功能变异也可能是由差异启动子引起的基因表达模式引起的。
表3.使用生长在六种环境中的192个F9花生RIL检测油含量的推定QTL
a标记或间隔表示链接图上LOD最大的区域。
b加值>0和<0分别代表来自JH5和KX01-6的QTL的递增效应。
表4.候选基因中出现由SNPs和Indels引起的突变类型.
实施例5、通过转录组分析确定候选基因
为了进一步确定qCOA08_1和qCOA08_2背后的调控基因,我们检测了qCOA08_1和qCOA08_2两个候选区间中所有基因在油脂代谢关键期的表达模式。结果qCOA08_1和qCOA08_2候选区间内分别注释了153和146个候选基因。其中,55个预测基因在两个亲本中均未表达,说明这些基因都不是qCOA08_1和qCOA08_2的候选基因,而其余144个预测基因根据倍数变化(log2fc≥1或≤-1)和P值(P<0.001)检测差异表达。结果显示,qCOA08_1和qCOA08_2候选区间内分别只有13个和6个差异表达基因(DEGs),但只有10个差异表达基因的启动子区存在差异,表明目标基因很可能隐藏在这些基因中。令人兴奋的是,在这10个预测基因中,有一个很强的候选基因,Ah.UVRI66编码一个乙烯转录因子Wrinkled 1(WRI1)引起了我们的高度关注。根据转录组分析(表5),Ah.UVRI66的FPKM值在低含油量品种KX01-6中约为4.65,而Ah.UVRI66在高含油量品种JH5中表达量约为22.03,为低油品种的5倍。此外,在Ah.UVRI66启动子区的在748和749核苷酸之间发现了一个6bp(ATATAT)的插入,暗示亲本之间表达的差异可能是由插入的6bp片段引起的。同时,根据前人的文献,我们发现WRI1家族基因已被广泛研究,并被证明是一个在油料植物中转录调控种子储油合成的主基因。因此,综合所有的结果,我们推测Ah.UVRI66是qCOA08_2背后的调控基因。无论如何,我们将qCOA08_2背后的调控基因命名为花生花生高油基因1(HighOilFavorablegene1inArachishypogaea AhyHOF1)。
表5.分别在qCOA08_1和qcoa 08_2的候选区间中差异表达基因的注释
实施例6、分子标记及其引物的开发。
课题组已经成功开发了dCAPs标记,具体引物和序列见下表。
实施例7、AhyHOF1对花生油品质影响的评价
为方便AhyHOF1功能分析,基于所开发的分子标记及其引物,从F5群体中筛选出一个高世代剩余杂合单株(RHP),并从其后代中构建了一组近等基因型(NILs),并对其脂肪酸组成进行了田间鉴定。此后,我们将与JH5母本表型一致的系记为#AhyHOF1,与KX01-6父本表型一致的等位系记为#AhyHOF1。通过观察油分含量和10种脂肪酸组成(FAC)对这些NILs进行了评估(图4)。正如预期的那样,#AhyHOF1的油分比#Ahyhof1高约1.6个百分点左右(图4a)。同时,在脂肪酸成分(FAC)评价中,10个FAC性状中有5个(图4b、c、e、f和j)在#AhyHOF1和#Ahyhof1之间差异不显著,2个FAC性状(包括C18:0和C20:0,图4d和h)在#AhyHOF1中显著高于#Ahyhof1,其他3个FAC在#AhyHOF1中均显著高于#Ahyhof1(图4j、i和k)。以上结果表明,AhyHOF1主要通过增加部分脂肪酸组成含量来提高花生油含量,在田间环境中AhyHOF1对花生油品质有一定影响。
综上实施例可见,本发明基于RIL群体构建了包含2840个多态性SNP标记的高分辨率遗传图谱,并在6个环境中稳定检测到2个效应较大的有价值QTL位点qCOA08_1和qCOA08_2。根据侧翼标记物的物理位置,qCOA08_1和qCOA08_2的物理间隔被定位在~2.9Mb和~1.7Mb之间。同时进行全基因组重测序和转录组分析,将候选基因的数量限定在40个以内。而一个编码WRI1转录因子并在亲本间表现出不同表达模式的强候选基因,很可能通过影响部分FAC合成来调控油脂积累。本研究为花生优异等位基因的克隆提供了有价值的信息。与qCOA08_1和qCOA08_2紧密连锁的多态SNP标记(或开发的dCAPs标记)有助于加快MAS的育种进程。本发明克服了以往研究中MAS育种的这些局限性,利用48K基因芯片构建了一个双亲本RIL群体的高分辨率遗传图谱,探索了可用于MAS育种的稳定有效的基因位点,开发了可用于育种的分子标记。这些信息对今后开展花生油含量育种工作具有重要科研和实用意义。
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (7)
1.一种与花生含油量相关的主效QTL,所述主效QTL命名为qCOA08_1,所述qCOA08_1位于花生基因组的第8个连锁群中,其标记区间为Chr.08_38056541-Chr.08_40978133。
2.用于花生含油量性状识别的分子标记,为关联于主效QTL qCOA08_1下游侧翼序列的分子标记Chr.8_46776535,其核苷酸序列依次如SEQ ID NO.1所示。
3.用于鉴定或辅助鉴定花生含油量性状的dCAPs引物,包括:对应于分子标记Chr.8_46776535的dCAPs-F和dCAPs-R引物对,其核苷酸序列依次如SEQ ID NO.2和SEQ ID NO.3所示;对应的酶切位点为EcoRI,对应的识别序列为GAATTC。
4.用于花生含油量性状识别的分子标记,为关联于主效QTL qCOA08_1上游侧翼序列的分子标记Chr.8_44887537,其核苷酸序列依次如SEQ ID NO.4所示。
5.用于鉴定或辅助鉴定花生含油量性状的dCAPs引物,包括:对应于分子标记Chr.8_44887537的dCAPs-F和dCAPs-R引物对,其核苷酸序列依次如SEQ ID NO.5和SEQ ID NO.6所示;对应的酶切位点为BsrBI,对应的识别序列为GAGCGG。
6.用于鉴定或辅助鉴定花生含油量性状的试剂盒,通过检测花生主效QTL qCOA08_1的多态性,用于花生含油量相关的分子标记辅助育种。
7. 根据权利要求6所述的试剂盒,其特征在于:所述试剂盒包含序列表中SEQ ID NO.2和SEQ ID NO.3组成的dCAPs引物对,和/或序列表中SEQ ID NO.5和SEQ ID NO.6组成的dCAPs引物对。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310859370.XA CN116790798A (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点qCOA08_1及其应用 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310859370.XA CN116790798A (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点qCOA08_1及其应用 |
CN202211096080.6A CN115862731B (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点及其应用 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211096080.6A Division CN115862731B (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点及其应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116790798A true CN116790798A (zh) | 2023-09-22 |
Family
ID=85660836
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310859370.XA Pending CN116790798A (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点qCOA08_1及其应用 |
CN202211096080.6A Active CN115862731B (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点及其应用 |
CN202310859371.4A Pending CN116694810A (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点qCOA08_2及其应用 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211096080.6A Active CN115862731B (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点及其应用 |
CN202310859371.4A Pending CN116694810A (zh) | 2022-09-08 | 2022-09-08 | 与花生含油量相关主效数量性状基因位点qCOA08_2及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (3) | CN116790798A (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108893551B (zh) * | 2018-07-16 | 2021-09-14 | 河北省农林科学院粮油作物研究所 | 一种检测花生高油酸含量的分子标记方法及应用 |
CN112080578B (zh) * | 2020-09-24 | 2022-11-18 | 中国农业科学院油料作物研究所 | 与花生含油量主效qtl位点连锁的分子标记及其应用 |
CN112063746B (zh) * | 2020-09-27 | 2022-05-17 | 中国农业科学院油料作物研究所 | 位于花生fad2a基因启动子区与花生油酸和亚油酸含量相关的snp分子标记及其应用 |
CN114736986B (zh) * | 2022-05-10 | 2024-04-30 | 贵州省油菜研究所 | 甘蓝型油菜种子含油量性状的主效qtl位点的snp分子标记相关检测引物或探针及应用 |
-
2022
- 2022-09-08 CN CN202310859370.XA patent/CN116790798A/zh active Pending
- 2022-09-08 CN CN202211096080.6A patent/CN115862731B/zh active Active
- 2022-09-08 CN CN202310859371.4A patent/CN116694810A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
CN115862731B (zh) | 2023-08-22 |
CN116694810A (zh) | 2023-09-05 |
CN115862731A (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10455783B2 (en) | Compositions and methods of plant breeding using high density marker information | |
Li et al. | Genome-wide association mapping of QTL underlying seed oil and protein contents of a diverse panel of soybean accessions | |
Ho et al. | Improvement of hybrid yield by advanced backcross QTL analysis in elite maize | |
Chai et al. | Validation of DGAT1-2 polymorphisms associated with oil content and development of functional markers for molecular breeding of high-oil maize | |
WO2021226806A1 (zh) | 甘蓝型油菜高密度全基因组snp芯片及其应用 | |
US20100037342A1 (en) | Methods and compositions for breeding plants with enhanced yield | |
Campbell et al. | Fast neutron-induced structural rearrangements at a soybean NAP1 locus result in gnarled trichomes | |
CN110157829B (zh) | 一种与油菜千粒重关联的分子标记snpa9-5及应用 | |
CN110656198B (zh) | 与小麦株高相关的kasp引物组及其应用 | |
Zhang et al. | Grain size selection using novel functional markers targeting 14 genes in rice | |
CN111187852B (zh) | 一个与小麦产量相关的snp位点及其应用 | |
Han et al. | QTL mapping pod dehiscence resistance in soybean (Glycine max L. Merr.) using specific-locus amplified fragment sequencing | |
Dang et al. | OsSYL2 AA, an allele identified by gene‐based association, increases style length in rice (Oryza sativa L.) | |
Li et al. | Molecular characterization of barley 3H semi-dwarf genes | |
US20170022574A1 (en) | Molecular markers associated with haploid induction in zea mays | |
Salem et al. | Assessment of genetic diversity and population structure on Kenyan sunflower (Helianthus annus L.) breeding lines by SSR markers | |
Qu et al. | Identification of QTL for seed coat colour and oil content in Brassica napus by association mapping using SSR markers | |
Zhou et al. | Mining favorable alleles for five agronomic traits from the elite rapeseed cultivar zhongshuang 11 by QTL mapping and integration | |
Labate et al. | Multilocus sequence data reveal extensive departures from equilibrium in domesticated tomato (Solanum lycopersicum L.) | |
Zhang et al. | Identification of a stable major QTL for fresh-seed germination on chromosome Arahy. 04 in cultivated peanut (Arachis hypogaea L.) | |
Mori et al. | Mapping of QTLs controlling epicotyl length in adzuki bean (Vigna angularis) | |
CN115862731B (zh) | 与花生含油量相关主效数量性状基因位点及其应用 | |
Okada et al. | Genetic factors associated with favourable pollinator traits in the wheat cultivar Piko | |
EA023910B1 (ru) | Способ идентификации растения подсолнечника или зародышевой плазмы подсолнечника, имеющего устойчивость к orobanche | |
Majid et al. | Marker-assisted identification of potential fertility restorers for the development of three-line hybrids in rice (Oryza sativa L.) |
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 |