CN114223453B - 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法 - Google Patents

一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法 Download PDF

Info

Publication number
CN114223453B
CN114223453B CN202111560662.0A CN202111560662A CN114223453B CN 114223453 B CN114223453 B CN 114223453B CN 202111560662 A CN202111560662 A CN 202111560662A CN 114223453 B CN114223453 B CN 114223453B
Authority
CN
China
Prior art keywords
sequencing
apple
rootstock
mutagenesis
whole genome
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
Application number
CN202111560662.0A
Other languages
English (en)
Other versions
CN114223453A (zh
Inventor
张世忠
张步航
郑成超
李菡
齐琦
姜舒纳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Agricultural University
Original Assignee
Shandong Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong Agricultural University filed Critical Shandong Agricultural University
Priority to CN202111560662.0A priority Critical patent/CN114223453B/zh
Publication of CN114223453A publication Critical patent/CN114223453A/zh
Application granted granted Critical
Publication of CN114223453B publication Critical patent/CN114223453B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G17/00—Cultivation of hops, vines, fruit trees, or like trees
    • A01G17/005—Cultivation methods
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G2/00—Vegetative propagation
    • A01G2/10—Vegetative propagation by means of cuttings
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G31/00—Soilless cultivation, e.g. hydroponics
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01H—NEW PLANTS OR NON-TRANSGENIC PROCESSES FOR OBTAINING THEM; PLANT REPRODUCTION BY TISSUE CULTURE TECHNIQUES
    • A01H5/00—Angiosperms, i.e. flowering plants, characterised by their plant parts; Angiosperms characterised otherwise than by their botanic taxonomy
    • A01H5/06—Roots
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01H—NEW PLANTS OR NON-TRANSGENIC PROCESSES FOR OBTAINING THEM; PLANT REPRODUCTION BY TISSUE CULTURE TECHNIQUES
    • A01H6/00—Angiosperms, i.e. flowering plants, characterised by their botanic taxonomy
    • A01H6/74—Rosaceae, e.g. strawberry, apple, almonds, pear, rose, blackberries or raspberries
    • C—CHEMISTRY; METALLURGY
    • C07—ORGANIC CHEMISTRY
    • C07K—PEPTIDES
    • C07K14/00—Peptides having more than 20 amino acids; Gastrins; Somatostatins; Melanotropins; Derivatives thereof
    • C07K14/415—Peptides having more than 20 amino acids; Gastrins; Somatostatins; Melanotropins; Derivatives thereof from plants
    • 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/63—Introduction of foreign genetic material using vectors; Vectors; Use of hosts therefor; Regulation of expression
    • C12N15/79—Vectors or expression systems specially adapted for eukaryotic hosts
    • C12N15/82—Vectors or expression systems specially adapted for eukaryotic hosts for plant cells, e.g. plant artificial chromosomes (PACs)
    • C12N15/8201—Methods for introducing genetic material into plant cells, e.g. DNA, RNA, stable or transient incorporation, tissue culture methods adapted for transformation
    • C12N15/8202—Methods for introducing genetic material into plant cells, e.g. DNA, RNA, stable or transient incorporation, tissue culture methods adapted for transformation by biological means, e.g. cell mediated or natural vector
    • C12N15/8205—Agrobacterium mediated transformation
    • 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/63—Introduction of foreign genetic material using vectors; Vectors; Use of hosts therefor; Regulation of expression
    • C12N15/79—Vectors or expression systems specially adapted for eukaryotic hosts
    • C12N15/82—Vectors or expression systems specially adapted for eukaryotic hosts for plant cells, e.g. plant artificial chromosomes (PACs)
    • C12N15/8241—Phenotypically and genetically modified plants via recombinant DNA technology
    • C12N15/8261—Phenotypically and genetically modified plants via recombinant DNA technology with agronomic (input) traits, e.g. crop yield
    • 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
    • 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)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Biotechnology (AREA)
  • Wood Science & Technology (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Botany (AREA)
  • Zoology (AREA)
  • Environmental Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Microbiology (AREA)
  • Developmental Biology & Embryology (AREA)
  • Analytical Chemistry (AREA)
  • Cell Biology (AREA)
  • Plant Pathology (AREA)
  • Physiology (AREA)
  • Gastroenterology & Hepatology (AREA)
  • Medicinal Chemistry (AREA)
  • Natural Medicines & Medicinal Plants (AREA)
  • Mycology (AREA)
  • Immunology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法，包括步骤：苹果无融合生殖砧木全基因组测序；测序木本材料单株去掉柱头后采收种子；采收的种子低温冷藏层积；不同浓度EMS诱变，种子萌发并播种；收集种子进行层积萌发；水培并鉴定表型；差异材料植物组织培养，诱导生根鉴定根系表型；差异材料进行全基因组重测序和转录组测序；寻找诱变基因；发根农杆菌转化体系进行诱变基因功能验证。本发明利用苹果无融合生殖特性、根系表型分析及转录组高通量快速筛选鉴定新型无融合砧木，极大的节省了筛选的时间与空间，且方法流程简单、操作技术容易掌握，具有广阔的应用前景。

Description

一种基于全基因组诱变创制苹果无融合异源多倍体砧木的 方法
技术领域
本发明属于植物组织培养技术领域，尤其涉及一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法。
背景技术
无融合生殖是指不经过精子与卵细胞的受精结合而直接产生无性胚和无性种子的一种特殊的无性繁殖，其子代在遗传上是其母本的克隆，保留着母本的基因型，其后代不产生分离现象，从而表现出高度的整齐一致性及其他的优良性状，并在繁殖过程中能够保留下来。与现行苹果栽培应用的有性种子作基砧和无性系矮砧相比，无融合生殖砧木具有许多优点：通过无性种子繁殖，后代整齐一致；繁殖系数高，与扦插、压条等繁殖方式相比，繁殖效率高，根系壮，立地性好，适应性强；不带病毒，适应现代果树无毒栽培的发展趋势。
无融合生殖实生砧木首先要具备无融合生殖特征，并且一别以往的无性繁殖方法，而采用种子播种进行实生繁殖。大多数果树的遗传物质都是高度杂合的，实生繁殖将产生大幅度的分离，而使后代参差不齐。为防止这种现象的发生，实生砧木必须具有无融合生殖特征。
但目前对于具有无融合生殖特性的实生砧木存在创制困难的问题，制约了无融合生殖在育种程序和无性种子生产中的应用。尤其是苹果砧木选育上，存在童期长、植株占空间大、筛选周期长，高度的无融合生殖性状不能通过有性杂交传递给后代的问题。
发明内容
本发明提供了一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法。本发明克服了苹果无融合生殖砧木的增殖特征，实现其快速高目的性增殖，节省了极大的时间与空间，可以大规模快速获得新型异源多倍体砧木并进行诱变基因功能鉴定。
为了实现上述发明目的，本发明提供以下技术方案：
本发明提供了一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法，包括以下步骤：
（1）对苹果无融合生殖砧木进行全基因组测序；
（2）在苹果无融合生殖砧木开花时期，去除柱头采收种子；
（3）采收的种子低温冷藏层积后，进行诱变处理、消毒，然后再次低温冷藏层积，进行播种，待得到的萌发幼苗长出根系移栽水培；
（4）鉴定水培幼苗的根系表型；
（5）选取根系表型有差异的幼苗扦插生根，再次鉴定根系表型；
（6）选取扦插生根的根系表型有差异的幼苗进行组织培养，诱导生根后鉴定根系表型，得到差异诱变砧木；
（7）将差异诱变砧木进行全基因组测序和转录组测序，与步骤（1）全基因组测序结果进行比较获得诱变基因；
（8）对诱变基因进行功能验证。
进一步的，所述步骤（1）中采用PacBio Hifi测序、Oxford Nanopore Technology测序和Illumina 二代测序平台对苹果无融合生殖砧木进行全基因组测序；然后利用Canuhifi 软件对得到的原始测序数据进行组装，组装结果用NextPolish 软件和二代Illumina数据进行校正并组装高质量基因组，利用已知苹果基因组序列和Hi-C数据，采用ALLHiC软件进行Contig的分相和挂载，得到挂载到全部染色体的苹果无融合生殖砧木基因组序列，对挂载到染色体的基因组利用NECAT校正的ONT测序数据进行Gap Closing, 然后利用PacBio hifi 数据进行Gap Closing，得到最终的染色体水平的分相的苹果无融合生殖砧木基因组。
进一步的，所述步骤（3）低温冷藏层积的步骤为：将采收的种子与潮湿沙土混匀置于4℃冰箱中层积25-50天。
进一步的，所述诱变处理为化学诱变，所述化学诱变的试剂包括EMS、秋水仙素、5-溴尿嘧啶。
进一步的，所述步骤（5）中根系表型的差异包括侧根数量增多。
进一步的，所述步骤（7）中采用Solexa平台和SOLiD平台对差异诱变砧木进行全基因组测序和转录组测序，通过比较基因组学和转录组WGCNA分析，获得诱变基因的核苷酸序列。
进一步的，所述步骤（8）中利用诱变基因来构建目的基因超表达载体，转入发根农杆菌后，侵染野生型组培苗茎尖1.0 cm，经共培养和筛选培养后获得转基因根系，进行功能验证。
本发明还提供了利用所述的方法获得的诱变基因，所述诱变基因的核苷酸序列如SEQ ID No.1所示。
进一步的，所述诱变基因编码的蛋白质的氨基酸序列如SEQ ID No.2所示。
本发明还提供了所述的诱变基因在用于鉴定无融合诱变系苹果砧木中的应用。
本发明与现有技术相比，具有以下优点和有益效果：
本发明对平邑甜茶进行全基因组测序；后代无融合种子用化学试剂诱变；将诱变后的种子继续层积萌发并播种；收集后代种子进行层积萌发；利用根系水培扫描系统分析根系表型；根系构型显著差异材料进行组织培养培养繁殖；组培苗生根验证表型；差异材料进行全基因组重测序和转录组测序并寻找差异基因；利用快速的发根农杆菌转化体系对突变基因进行基因功能验证。本发明基于异源三倍体分相全基因组解析，利用苹果无融合生殖特性、根系表型分析及转录组高通量快速筛选鉴定新型无融合砧木，极大的节省了筛选的时间与空间，使优秀无融合生殖苹果砧木的繁殖成为现实。本发明所述方法所需流程简单、操作技术容易掌握，且能够适用于所有的苹果无融合生殖砧木品种，具有广适性，并实现了优秀无融合砧木的快速大量的增殖，因此具有广阔的应用前景。
附图说明
图1为实施例2中平邑甜茶去除柱头前的花；
图2为实施例2中平邑甜茶经EMS诱变后的种子；
图3为实施例2中水培获得的平邑甜茶差异幼苗根系；
图4为实施例2中继代组织培养得到的平邑甜茶组培苗；
图5为实施例2中生根诱导得到的诱变平邑甜茶幼苗（左）和对照组（右）平邑甜茶幼苗；
图6为实施例2中侧根数量差异显著性分析；
图7为实施例2中基因超表达（左）和对照组（右）平邑甜茶幼苗。
具体实施方式
结合以下具体实例对本发明的技术方案作进一步详细的说明。
下述实施例中，如无特殊说明，所使用的实验方法均为常规方法，所用材料、试剂等均可从生物或化学试剂公司购买。
实施例1
本发明的一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法，包括以下步骤：
（1）苹果无融合生殖砧木全基因组测序：对苹果无融合生殖目的材料采用但不限于PacBio Hifi测序、Oxford Nanopore Technology 测序和Illumina 二代测序平台等综合手段进行全基因组测序，利用Canu hifi 软件对得到的原始数据进行组装，组装结果用NextPolish 软件和二代Illumina 数据进行校正并组装高质量基因组，利用已发表的苹果基因组序列如寒富、金冠等和Hi-C数据，采用ALLHiC软件进行Contig的分相和挂载，最终得到挂载到全部染色体的无融合生殖目的材料基因组序列，对挂载到染色体的基因组利用NECAT校正的ONT测序数据进行Gap Closing, 然后再次利用PacBio hifi 数据进行GapClosing，得到最终的染色体水平的分相的苹果无融合生殖目的材料基因组；
（2）测序木本材料后代种子进行诱变：在无融合生殖目的材料开花时期，去除柱头，成熟后收集种子；将种子与潮湿沙土混匀（沙土潮湿程度为手握成团不滴水），置于4℃冰箱中层积，25-50天后获得即将萌发的种子（不同目的材料可根据自身特点适当调整）；清水冲洗出种子，0.5% EMS诱变溶液避光处理种子16小时（可根据不同的实验要求调整EMS浓度和诱变时间），0.1M硫代硫酸钠溶液消毒处理，继续4℃层积20天，进行大田播种；其中，诱变溶液可以选用秋水仙素、5-溴尿嘧啶等。
（3）水培苗观察根系表型：收集诱变株系后代种子，消毒处理后，置于潮湿沙土中，于22-24℃，培养至幼苗长出根系即可移栽水培，2周后观察确定差异根系；
（4）根系构型有差异的幼苗进行扦插生根：根系构型有差异的幼苗进行土壤移栽，2个月后进行扦插，待长出足够的根后统计生根率及表型；
（5）组培苗进行根系构型表型及分析：扦插生根有差异的株系取幼嫩的组织进行组织培养，然后转移到继代培养基（每升MS培养基中含0.5mg 6-BA + 0.05mg IBA，30g蔗糖，7.5g琼脂，pH为5.8-5.85）中，于22-24℃，进行扩增繁殖获得更多幼苗；2个月后选取继代培养基中长势良好的幼苗转移至生根培养基（每升1/2MS培养基中含0.5 mg IBA，25 g/L蔗糖，7.5g琼脂，pH为5.8-5.85）诱导生根，再次确定表型；
（6）诱变材料进行全基因组重测序及转录组测序寻找差异基因：采用Solexa平台和SOLiD平台对差异诱变材料进行全基因组测序和转录组测序，通过比较基因组学和转录组WGCNA分析，获得差异基因的核苷酸序列及其编码的氨基酸序列；
（7）进行转基因根系验证：构建目的基因超表达载体并转入发根农杆菌K599, 侵染野生型组培苗茎尖1.0 cm，经共培养和筛选培养后获得转基因根系。
本发明是基于苹果无融合生殖砧木的特性，并利用化学诱变与组织培养技术而进行的高通量高规模快速筛选和鉴定无融合砧木。上述所述的方法能够适用于所有的无融合生殖砧木品种，具有广适性，并实现了优秀无融合砧木的快速大量的增殖。
实施例2
现以平邑甜茶为例，对本发明所述的基于全基因组诱变创制苹果无融合异源多倍体砧木的方法，进行详细介绍，该方法包括以下步骤：
（1）对单株沂蒙山区野生百年‘平邑甜茶’老树幼嫩树叶采样后用液氮迅速冷冻，送往北京百迈客生物科技有限公司测序，基因组提取由百迈客公司完成；
（2）采用PacBio Hifi测序、Oxford Nanopore Technology 测序和Illumina 二代测序平台进行全基因组测序，并利用Canu hifi 软件对 PacBio Hifi数据进行组装，组装结果用NextPolish 软件和二代Illumina 数据进行校正，利用寒富的基因组序列和Hi-C数据，采用ALLHiC软件进行Contig的分相和挂载，最终得到挂载到51条染色体的平邑甜茶基因组序列，染色体挂载率为94.33%，对挂载到染色体的基因组利用NECAT校正的ONT测序数据进行Gap Closing, 然后再次利用PacBio hifi 数据进行Gap Closing，得到最终的染色体水平的分相的平邑甜茶基因组共1.88Gb；
（3）在平邑甜茶开花时期（图1），去除柱头，成熟后收集种子；
（4）将种子与潮湿沙土混匀（沙土潮湿程度为手握成团不滴水），置于4℃冰箱中层积，30天后获得即将萌发的种子；
（5）清水冲洗出种子，0.5% EMS诱变溶液（每100mL诱变溶液中含有EMS 0.5mL，10M磷酸二氢钠溶液9.2 mL，1M磷酸氢二钠溶液0.8mL，DMSO溶液5mL,去离子水84.5mL）避光处理种子16小时，0.1M硫代硫酸钠溶液消毒处理（图2），继续4℃层积20天，进行大田播种；
（6）收集诱变株系后代种子，采用步骤（4）方法层积种子40天，获得萌发幼苗；
（7）将萌发的幼苗转移到Hoagland营养液中，2周后使用扫描仪（EpsonPerfection V850 Pro）扫描幼苗根系，WinRHIZO软件对根系进行侧根数量统计，GraphPadPrism软件对数据进行差异显著性分析（图6 A）；
（8）选择侧根数量明显增多的幼苗（图3），取幼嫩的组织在超净工作台中用75%的酒精处理30 s进行表面消毒，再用无菌水冲洗3次；然后将外植体放进0.1%升汞溶液中浸泡9 min进行最后一步消毒，再用无菌水冲洗3次；最后，将消毒好的外植体接种到生芽培养基（每升MS培养基中含3mg TDZ，0.5mg NAA，30g蔗糖，7.5g琼脂，pH为5.8-5.85）中暗培养14天，光培养30天；
（9）将组培苗转移到继代培养基（每升MS培养基中含0.5mg 6-BA + 0.05mg IBA，30g蔗糖，7.5g琼脂，pH为5.8-5.85）中，于22-24℃条件下进行扩增繁殖获得更多组培苗（图4）；
（10）选取继代培养基中长势良好的幼苗切取茎尖1.5cm转移至生根培养基（每升1/2MS培养基中含0.5 mg IBA，25 g/L蔗糖，7.5g琼脂，pH为5.8-5.85）中于22-24℃生根培养1个月，得到平邑甜茶幼苗（图5），使用步骤（7）中方法再次鉴定材料侧根明显增多（图6B）；
（11）将步骤（10）中变异材料送至北京百迈客生物科技有限公司进行全基因组重测序，基因组提取由公司完成；
（12）将步骤（10）中变异材料送至武汉贝纳科技服务有限公司，将步骤（2）获得的全基因组作为参考进行有参转录组测序，RNA提取由公司完成；利用双端测序技术对mRNA进行测序，原始 fastq数据经过fastp过滤得到剩余的有效数据如reads和碱基数等（表1）；通过Star将过滤后的转录组序列与参考基因比对，并进行统计，获得测序数据经过滤后单端的 clean reads 总数和参考序列上有唯一比对位置的 reads 所占比例等（表2）；利用表达定量得到的基因在各个样品中的表达量 reads count 数据，进行差异表达分析，差异表达分析的软件为 DESeq2，筛选阈值为 padj < 0.05 且 |log2FoldChange| > 1，如若显著差异基因数目过少时筛选阈值为 p value < 0.05 且 |log2FoldChange| > 1，本次分析使用的是FDR值进行显著差异过滤，获得经鉴定的差异表达基因数目（表3）；
表1：测序数据量信息统计
Sample Raw reads Raw bases Clean reads Clean bases Q20 rate Q30 rate GC content
CK1 47369523 7186254800 47074423 7069620678 98.152% 94.523% 46.325%
CK2 41852790 6527369100 41592057 6421428803 98.437% 94.635% 47.258%
Line1 45706349 6810493300 45421610 6699957850 97.583% 93.986% 47.947%
Line2 39479586 5870469600 39233638 5775190892 98.428% 94.471% 46.956%
注：Sample 为样本名称；Raw reads 为原始未过滤数据的 Reads 数；Raw bases为原始未过滤数据的碱基数； Clean reads 为过滤后剩余的 Reads 数；Clean bases 过滤后剩余的碱基数； Q20_rate 为过滤后，总序列中质量值大于 20（错误率小于1%）的碱基数的比例，该值越大说明测序质量越好； Q30_rate 为过滤后，总序列中质量值大于 30（错误率小于0.1%）的碱基数的比例，该值越大说明测序质量越好；GC_content 为过滤后数据GC含量。
表2：数据比对结果统计
Sample Total_reads(bp) Unique_mapped(%) Multiple_mapped(%) Unmapped(%)
CK1 21328786 35.62 54.19 10.19
CK2 22802512 35.59 53.47 10.63
Line1 20654880 33.37 52.03 14.60
Line2 19.801722 32.58 51.26 16.16
注：上表中 Sample 为样品名称，Total_reads 为测序数据经过滤后单端的clean reads 总数，Unique_mapped 为参考序列上有唯一比对位置的 reads 所占比例，Multiple_mapped 为在参考序列上有多个比对位置的 reads 所占比例，Unmapped 为未比对上的 reads 所占比例。
表3：差异表达基因数目
Comparison Significant_Diff_Number Up Down
CK1_vs_ Line1 10612 4788 5383
CK2_vs_ Line2 11533 5351 5840
注：上表为每个比较组显著差异基因数目，以及上调和下调基因数目。
（13）通过KEGG和GO等手段进行转录本功能注释，共得到涉及34个主要通路的2764个突变功能基因，预测其包含参与调控主根、侧根和根毛的基因；选择其中一个调控侧根数量的主要通路中的主效基因进行功能验证，即选择核苷酸序列如SEQ ID No.1（MA01：Mh04G130085）所示的差异基因，其编码得到的氨基酸序列如SEQ ID No.2（MA01：Mh04G130085）所示；
（14）为了验证定向改造该砧木的可行性，对选择的诱变基因进行基因功能验证。首先构建目的基因过表达载体，载体选择购自BioVector NTCC典型培养物保藏中心的pROK2，酶切位点选择BamH1和Kpn1，使用苏州新赛美生物科技有限公司的EasyFusionAssembly Master Mix定向同源重组克隆试剂盒连接目的基因，连接条件为50℃、30min；大肠杆菌为实验室保存的感受态TOP10，卡那霉素抗性；将连接好的质粒转化进感受态TOP10（冰水浴30 min，37℃热激90s，冰水浴5 min）中，涂板后37℃过夜培养，挑取单菌落，经菌液PCR和载体酶切鉴定后送生工生物工程有限公司测序，测序结果与步骤（2）获得的全基因组中序列比对结果相同，表示载体构建成功；
表4：目的基因扩增所用引物序列
引物名称 引物序列
MA01-F acgggggactctagaggatccATGGAGAAAGAAAAAGAAAGAGCC
MA01-R gcccttgctcaccatggtaccTTATTTGCATGAGAGAGTAGCAATAAAT
（15）转入发根农杆菌K599。感受态K599购买自上海唯地生物技术有限公司，利福平抗性；将构建好的载体转化后涂布固体YEP抗性板（每升YEP含100μL利福平，50μL卡那霉素，7.5g琼脂），28℃培养48h后挑取单菌落于液体YEP抗性培养基（每升YEP含100μL利福平，50μL卡那霉素）中，振荡培养6h后进行菌液PCR，选取转化成功的保存菌种。
（16）取少量菌种于液体YEP抗性培养基中过夜活化，第二天早上用移液枪吸取100μL菌液加入到20mL液体YEP抗性培养基中28℃振荡培养4-6小时，用分光光度计测量农杆菌菌液浓度，选取OD值为0.4-0.6的菌液离心（5000 rpm，5 min）收集菌体，无菌水清洗后离心倒掉上清，取少量菌在20mL侵染液（1/2 MS+15 g/L蔗糖，pH为5.20）中重悬；切取野生型组培苗茎尖1.0 cm，蘸取侵染液后于共培养培养基（1/2 MS+15 g/L蔗糖，8.5 g/L琼脂，pH为5.20）中22℃暗培养3天，转移到筛选培养基（1/2 MS+15 g/L蔗糖，6.0 g/L琼脂，pH为5.90，特美汀500 mg/L）,每10天更换一次筛选培养基；诱变基因过表达的平邑甜茶（左）和对照组的平邑甜茶（右）如图7所示，相较于对照组来说，过表达的材料侧根数量多（图6C）且密度大，单体长度稍短，主根增长；主根长度增加可以更好的固定植株、分生侧枝，吸收深层土壤中的水分及微量元素，侧根密度大扩大了根系的分布范围、增加了根系吸收土壤养分的面积，增强了综合吸收水分与养分的能力，可更有效地实现早期丰产优质的目标。
通过上述实验，可以看出，利用本发明所述的方法能够快速筛选诱变获得的无融合砧木，相对于现有技术显著的缩短了筛选时间，缩小了筛选空间，使无融合苹果的增殖成为可能。
以上实施例仅用以说明本发明的技术方案，而非对其进行限制；尽管参照前述实施例对本发明进行了详细的说明，对于本领域的普通技术人员来说，依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明所要求保护的技术方案的精神和范围。
序列表
<110> 山东农业大学
<120> 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法
<160> 2
<170> SIPOSequenceListing 1.0
<210> 1
<211> 1371
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 1
atggagaaag aaaaagaaag agccaataga atccacatac ttgtggttcc tcttcctgtc 60
caaggccaca taaacccaat gctccaattc tcaaaacgct tagcatccaa aggtttaaga 120
gtgacacttg tcacaatctc atccaacact gaagccattg aaactcaatt gtgtaaggtc 180
aagattgagc caattcatgc aacatctgaa gaagaagaag aagccggcaa taatgtggag 240
aacaaggttc tgtggttcca aaccattttg tcacggagat tgccacagct gatactgaag 300
caagagaatg atggctaccc aatcatctgt cttgtctatg actcagccat gccttgggtt 360
ttagacacag ccaaagggct aggaattttt ggggcttcat ttttcactca ctcttgtgct 420
gttggagctg tgtactataa tgtgcatgaa gggttgctta gggttcctct gaaagagcct 480
acaatatcat tacctgggtt gcccttactt gagccacaag atctgccttc tttcatctat 540
gatcccgggt catatccgtc cttcttgaac ttggtggtta gtaggttttc aaatatcagt 600
ggcgctgatt ggatcttctt caacacattc catgggctgg aaaaagaggt ggtgaactgg 660
atgagaacgc gatggccgat caagacaata ggaccaaccc ttccatcaat gtacctagac 720
aagagattag aagatgacaa agattatggc ttcagcattt tcaactccaa cactgaagct 780
tgcaagaagt ggctagaatc aaaagaaaca ggcacagtcg tatatgtatc atttggaagc 840
atggcaaatt taggagaaaa acaaatggag gaaatagcat tgggattaaa agggagcaac 900
gctaacttct tgtgggtagt tagagaatcc gaaactcaga agcttccaag caattttgaa 960
gagcaaacat cagagaaggg tttggttgta aattggtgtc ctcaattaga agttttgggt 1020
cacagggccc tcggttgctt catgacacat tgcggttgga actcaacgct cgaggcattg 1080
agctctggag tccccatggt ggcaatgcca cagtggacag accaaacgac taatgccaag 1140
tttgtagaag atgaatggaa agtaggagta agggtcaagg ttgaccaaat gggaattgtc 1200
actaaagaag aaatagaaag atgtatagca gaagtgattg aaggagagag agggaaggag 1260
attaagagga atttaatgag atggagagaa ctggctaagg aggccatgga tgaaggagga 1320
agctctgata agaacattga ggaatttatt gctactctct catgcaaata a 1371
<210> 2
<211> 456
<212> PRT
<213> 人工序列(Artificial Sequence)
<400> 2
Met Glu Lys Glu Lys Glu Arg Ala Asn Arg Ile His Ile Leu Val Val
1 5 10 15
Pro Leu Pro Val Gln Gly His Ile Asn Pro Met Leu Gln Phe Ser Lys
20 25 30
Arg Leu Ala Ser Lys Gly Leu Arg Val Thr Leu Val Thr Ile Ser Ser
35 40 45
Asn Thr Glu Ala Ile Glu Thr Gln Leu Cys Lys Val Lys Ile Glu Pro
50 55 60
Ile His Ala Thr Ser Glu Glu Glu Glu Glu Ala Gly Asn Asn Val Glu
65 70 75 80
Asn Lys Val Leu Trp Phe Gln Thr Ile Leu Ser Arg Arg Leu Pro Gln
85 90 95
Leu Ile Leu Lys Gln Glu Asn Asp Gly Tyr Pro Ile Ile Cys Leu Val
100 105 110
Tyr Asp Ser Ala Met Pro Trp Val Leu Asp Thr Ala Lys Gly Leu Gly
115 120 125
Ile Phe Gly Ala Ser Phe Phe Thr His Ser Cys Ala Val Gly Ala Val
130 135 140
Tyr Tyr Asn Val His Glu Gly Leu Leu Arg Val Pro Leu Lys Glu Pro
145 150 155 160
Thr Ile Ser Leu Pro Gly Leu Pro Leu Leu Glu Pro Gln Asp Leu Pro
165 170 175
Ser Phe Ile Tyr Asp Pro Gly Ser Tyr Pro Ser Phe Leu Asn Leu Val
180 185 190
Val Ser Arg Phe Ser Asn Ile Ser Gly Ala Asp Trp Ile Phe Phe Asn
195 200 205
Thr Phe His Gly Leu Glu Lys Glu Val Val Asn Trp Met Arg Thr Arg
210 215 220
Trp Pro Ile Lys Thr Ile Gly Pro Thr Leu Pro Ser Met Tyr Leu Asp
225 230 235 240
Lys Arg Leu Glu Asp Asp Lys Asp Tyr Gly Phe Ser Ile Phe Asn Ser
245 250 255
Asn Thr Glu Ala Cys Lys Lys Trp Leu Glu Ser Lys Glu Thr Gly Thr
260 265 270
Val Val Tyr Val Ser Phe Gly Ser Met Ala Asn Leu Gly Glu Lys Gln
275 280 285
Met Glu Glu Ile Ala Leu Gly Leu Lys Gly Ser Asn Ala Asn Phe Leu
290 295 300
Trp Val Val Arg Glu Ser Glu Thr Gln Lys Leu Pro Ser Asn Phe Glu
305 310 315 320
Glu Gln Thr Ser Glu Lys Gly Leu Val Val Asn Trp Cys Pro Gln Leu
325 330 335
Glu Val Leu Gly His Arg Ala Leu Gly Cys Phe Met Thr His Cys Gly
340 345 350
Trp Asn Ser Thr Leu Glu Ala Leu Ser Ser Gly Val Pro Met Val Ala
355 360 365
Met Pro Gln Trp Thr Asp Gln Thr Thr Asn Ala Lys Phe Val Glu Asp
370 375 380
Glu Trp Lys Val Gly Val Arg Val Lys Val Asp Gln Met Gly Ile Val
385 390 395 400
Thr Lys Glu Glu Ile Glu Arg Cys Ile Ala Glu Val Ile Glu Gly Glu
405 410 415
Arg Gly Lys Glu Ile Lys Arg Asn Leu Met Arg Trp Arg Glu Leu Ala
420 425 430
Lys Glu Ala Met Asp Glu Gly Gly Ser Ser Asp Lys Asn Ile Glu Glu
435 440 445
Phe Ile Ala Thr Leu Ser Cys Lys
450 455

Claims (8)

1.一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法，其特征在于，包括以下步骤：
（1）对苹果无融合生殖砧木进行全基因组测序；
（2）在苹果无融合生殖砧木开花时期，去除柱头采收种子；
（3）采收的种子低温冷藏层积后，进行诱变处理、消毒，然后再次低温冷藏层积，进行播种，待得到的萌发幼苗长出根系移栽水培；所述低温冷藏层积的步骤为：将采收的种子与潮湿沙土混匀置于4℃冰箱中层积25-50天；
（4）鉴定水培幼苗的根系表型；
（5）选取根系表型有差异的幼苗扦插生根，再次鉴定根系表型；所述根系表型的差异为侧根数量增多；
（6）选取扦插生根的根系表型有差异的幼苗进行组织培养，诱导生根后鉴定根系表型，得到差异诱变砧木；
（7）将差异诱变砧木进行全基因组测序和转录组测序，与步骤（1）全基因组测序结果进行比较获得诱变基因；
（8）对诱变基因进行功能验证。
2.根据权利要求1所述的方法，其特征在于，所述步骤（1）中采用PacBio Hifi测序、Oxford Nanopore Technology 测序和Illumina 二代测序平台对苹果无融合生殖砧木进行全基因组测序；然后利用Canu hifi 软件对得到的原始测序数据进行组装，组装结果用NextPolish 软件和二代Illumina 数据进行校正并组装高质量基因组，利用已知苹果基因组序列和Hi-C数据，采用ALLHiC软件进行Contig的分相和挂载，得到挂载到全部染色体的苹果无融合生殖砧木基因组序列，对挂载到染色体的基因组利用NECAT校正的ONT测序数据进行Gap Closing, 然后利用PacBio hifi 数据进行Gap Closing，得到最终的染色体水平的分相的苹果无融合生殖砧木基因组。
3.根据权利要求1所述的方法，其特征在于，所述诱变处理为化学诱变，所述化学诱变的试剂为EMS或秋水仙素或5 -溴尿嘧啶。
4.根据权利要求1所述的方法，其特征在于，所述步骤（7）中采用Solexa平台和SOLiD平台对差异诱变砧木进行全基因组测序和转录组测序，通过比较基因组学和转录组WGCNA分析，获得诱变基因的核苷酸序列。
5.根据权利要求1所述的方法，其特征在于，所述步骤（8）中利用诱变基因来构建目的基因超表达载体，转入发根农杆菌后，侵染野生型组培苗茎尖1.0 cm，经共培养和筛选培养后获得转基因根系，进行功能验证。
6.利用权利要求1-5任一项所述的方法获得的诱变基因，其特征在于，所述诱变基因的核苷酸序列如SEQ ID No.1所示。
7.根据权利要求6所述的诱变基因，其特征在于，所述诱变基因编码的蛋白质的氨基酸序列如SEQ ID No.2所示。
8.权利要求6所述的诱变基因在用于鉴定无融合诱变系苹果砧木中的应用。
CN202111560662.0A 2021-12-20 2021-12-20 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法 Active CN114223453B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111560662.0A CN114223453B (zh) 2021-12-20 2021-12-20 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111560662.0A CN114223453B (zh) 2021-12-20 2021-12-20 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法

Publications (2)

Publication Number Publication Date
CN114223453A CN114223453A (zh) 2022-03-25
CN114223453B true CN114223453B (zh) 2022-12-27

Family

ID=80759199

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111560662.0A Active CN114223453B (zh) 2021-12-20 2021-12-20 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法

Country Status (1)

Country Link
CN (1) CN114223453B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115119744A (zh) * 2022-07-26 2022-09-30 山东农业大学 一种苹果无融合一系四群体育种方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101124887A (zh) * 2007-09-28 2008-02-20 南开大学 丹参三倍体的培育方法
CN101974629A (zh) * 2010-10-26 2011-02-16 西南大学 一种虚拟合成物种调查异源多倍体物种起源的方法
CN103026912A (zh) * 2013-01-08 2013-04-10 山东农业大学 一种提高苹果砧木无融合生殖率和苗木整齐度的方法
CN105925680A (zh) * 2016-05-06 2016-09-07 中国农业科学院蔬菜花卉研究所 一种四倍体马铃薯高通量测序开发标记的方法及其应用
CN111411107A (zh) * 2020-03-27 2020-07-14 武汉古奥基因科技有限公司 一种多倍体基因组survey的方法
CN111583995A (zh) * 2020-05-12 2020-08-25 西藏自治区农牧科学院水产科学研究所 一种多倍体生物基因组二倍化程度量化评估的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150181822A1 (en) * 2013-12-31 2015-07-02 Dow Agrosciences Llc Selection based on optimal haploid value to create elite lines

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101124887A (zh) * 2007-09-28 2008-02-20 南开大学 丹参三倍体的培育方法
CN101974629A (zh) * 2010-10-26 2011-02-16 西南大学 一种虚拟合成物种调查异源多倍体物种起源的方法
CN103026912A (zh) * 2013-01-08 2013-04-10 山东农业大学 一种提高苹果砧木无融合生殖率和苗木整齐度的方法
CN105925680A (zh) * 2016-05-06 2016-09-07 中国农业科学院蔬菜花卉研究所 一种四倍体马铃薯高通量测序开发标记的方法及其应用
CN111411107A (zh) * 2020-03-27 2020-07-14 武汉古奥基因科技有限公司 一种多倍体基因组survey的方法
CN111583995A (zh) * 2020-05-12 2020-08-25 西藏自治区农牧科学院水产科学研究所 一种多倍体生物基因组二倍化程度量化评估的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
桑树多倍体鉴定和育种研究现状;焦锋等;《西北农业学报》;20151231(第12期);第1-13页 *
金冠苹果成熟胚离体诱导四倍体的研究;李林光等;《园艺学报》;20071231(第5期);第1129-1134页 *

Also Published As

Publication number Publication date
CN114223453A (zh) 2022-03-25

Similar Documents

Publication Publication Date Title
CN107475210B (zh) 一种水稻白叶枯病抗性相关基因OsABA2及其应用
CN112725360B (zh) 棉花GhHDA6基因在调控植物开花期中的应用
CN109797157B (zh) 一种抗非生物逆境转录因子PbrbHLH92及其引物、编码的蛋白和应用
CN107177610B (zh) 一种调控种子大小的拟南芥mpk基因及增加种子大小的方法
CN109879944B (zh) 一种与植物抗旱性相关的ear1蛋白及其编码基因与应用
CN109423492B (zh) SlTOE1基因在调控番茄开花时间和产量中的应用
CN114223453B (zh) 一种基于全基因组诱变创制苹果无融合异源多倍体砧木的方法
CN108486149B (zh) 一种黄瓜CsWRKY50基因在增强黄瓜霜霉病抗性中的应用
CN108676081B (zh) 紫云英leafy基因及其应用
CN108342396B (zh) 玉米基因ZmEREB180在植物耐渍中的应用
CN110938122A (zh) 雄性不育基因OsNIN5及其应用和育性恢复的方法
CN114250233B (zh) 拟南芥钙离子通道基因AtCNGC3在菌核病防控中的应用
CN114164229B (zh) 利用FvePILS5基因的CRISPR/Cas9基因敲除载体获得再生效率高的草莓新种质的方法及应用
CN113215172B (zh) 雄性不育基因MsJMT及其应用
CN116064571A (zh) 一种芒果耐盐基因及其鉴定方法和应用
CN114703198A (zh) 一种番茄转运蛋白SlZIF1的克隆及其应用
CN114214333A (zh) 一种调控植物叶表皮毛发育和次生壁厚度的基因及其应用
CN113337534A (zh) 一种提高矮牵牛遗传转化效率的组培方法
CN100360672C (zh) 棉花金属硫蛋白基因GhMT1序列及其克隆与应用
CA2547514A1 (en) Method for screening genomic dna fragments
CN116445523B (zh) 一种小麦抗逆性基因TaBSL3与应用
CN114958866B (zh) 调控大豆分枝数的基因及其用途
CN116004672B (zh) 提高植物生物量和产量的磷酸甘油酸激酶基因及其应用
CN113249393B (zh) 大豆GmPCBER4基因、蛋白及应用
CN116254288B (zh) 一种春兰MIR156b基因在调控植物开花时间中的应用

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