1.本发明涉及技术领域，特别是一种用于烟草根黑腐病主效抗性位点筛选/检测的功能分子标记及其应用。

背景技术：

2.烟草根黑腐病(black root rot,brr)是由土壤真菌根串珠霉菌(thielaviopsis basicola)引起的一种常见土传病害。该真菌在世界主要烟叶产区广泛分布，在我国云南、贵州等主要烟叶主产区均发现有此病的发生。烟草根黑腐病的特征是根和幼茎基部等部位的黑色病变，进而导致营养缺乏，发育迟缓，晚熟，植株矮化且生长不均匀，叶片皱缩或萎蔫，严重时烟苗倒伏、死亡，从而严重影响烟草产量和质量。尽管通过轮作及化学药剂等手段，也能较好的防控根黑腐病对烟叶生产的危害，但上述手段不仅会造成环境污染且会增加烟叶的种植成本，而选育具有抗根黑腐病基因的品种是最经济有效且可持续的策略。3.迪勃纳氏烟草(nicotiana debneyi，2n＝4x＝48)是源于澳大利亚的一种野生烟草，它不仅是研究杂交特性的理想材料，而且携带有烟草根黑腐病、霜霉病和白粉病等多种病害抗性基因，已成为烟草性状改良的有效遗传资源。通过种间杂交和回交，其根黑腐病的单显性抗病基因rbrr1已成功被转育到白肋烟与烤烟等栽培烟草中。4.已知rbrr1对烟草根黑腐病菌株具有完全抗性，并在烟草各生育阶段均对根串珠霉菌具有较高的免疫能力。最新研究表明，烟草抗性位点rbrr1被定位在烟草17号染色体(k326参考基因组)上，ss219555和ss192650是与其共分离的显性caps标记。然而，caps标记必须使用限制性内切酶，需要较冗长的分析步骤，并受限于酶切效率不能适应高通量自动化的问题。另外，各种突变会引起酶切位点的增加或消失，极大限制了两个功能标记在分子标记辅助育种中的应用。5.因此，重新开发一种方便检测且能适应高通量要求的根黑腐病抗性功能标记，是开展烟草根黑腐病品种选育的迫切需求。

技术实现要素：

6.本发明提供了一种用于烟草根黑腐病主效抗性位点筛选/检测的功能分子标记及其应用，用于解决现有技术中存在的caps标记必须使用限制性内切酶，需要较冗长的分析步骤，并受限于酶切效率不能适应高通量自动化的技术问题。7.本发明提供了一种用于抗烟草根黑腐病烟草植株选育的功能分子标记，扩增所述功能分子标记所用引物对选自引物对组a或b，第一引物对由seq id no:1、seq id no:2组成；8.第二引物对由seq id no:3、seq id no:4组成；9.第三引物对由seq id no:5、seq id no:6组成；10.第四引物对由seq id no:7、seq id no:8组成；11.引物对组b为选自以下引物对中的任意两个：12.第五引物对由seq id no:9、seq id no:10组成；13.第六引物对由seq id no:11、seq id no:12组成；14.第七引物对由seq id no:13、seq id no:14组成。15.对应的引物对组a中：16.第一引物对中seq id no:为ind192650‑f、seq id no:2为ind192650‑r；17.第二引物对seq id no:3为ind219555‑3f、seq id no:4为ind219555‑3r；18.第三引物对seq id no:5为brr_lg2_ind‑2‑f、seq id no:6为brr_lg2_ind‑2‑r；19.第四引物对seq id no:7为brr‑lg2‑ind‑5‑2f、seq id no:8王brr‑lg2‑ind‑5‑4r；20.引物对组b中：21.第五引物对为seq id no:9为dm192650‑2f、seq id no:10为dm192650‑2r；22.第六引物对seq id no:11为dm219555‑2f、seq id no:12为dm219555‑2r；23.第七引物对为seq id no:13为dm5‑1f、seq id no:14为dm5‑2r。24.本发明提供功能分子标记，仅需以烟草dna为模块，采用高通量的扩增设备，以上述核苷酸序列为引物进行扩增后，即可得到所需的功能分子标记。采用该功能分子可有效指示含有对烟草根黑腐病明确抗性的rbrr1抗性位点，能有效辅助烟草植株育种后的筛选工作，同时检测结果不受突变后酶切准确率影响，能有效避免假阳性检测结果的出现，检测结果较准确，检测效率提高。25.从自然环境中筛选得到具有该抗性位点的烟草植株后，进行推广种植后，能有效降低烟叶主产区烟草根黑腐病的发病比例和程度。26.优选的，采用引物对组a扩增后指示烟草rbrr1抗性位点共显性标记。采用该组引物对进行扩增后，能准确有效指示rbrr1抗性位点的共显性标记。27.优选的，采用引物对组b扩增后指示烟草rbrr1抗性位点显性功能标记。采用该组引物对进行扩增后，能准确有效指示rbrr1抗性位点的显性功能标记。28.优选的，功能分子标记是由所述引物对以烟草植株基因组dna为模板经pcr扩增得到。根据已获得引物对，可采用高通量扩增设备进行扩增得到功能分子标记，操作效率高，准确性高。扩增后还需进行电泳检测，具体按现有制备分子标记常用方法进行，在此不累述。29.优选的，所述引物对组a选自所述第一引物对、第二引物对中的至少一对。所述第一引物对、第二引物对在扩增中均显示较优良的多态性。30.优选的，所述引物对组a中至少包括：第四引物对。通过采用第四引物对进行扩增后，扩增效率高，在检测流程不精细的情况下也能很好的指示检测结果，适合相对较简陋的检查环境中使用。31.本发明的另一面还提供了一种功能分子标记在烟草辅助育种或者抗烟草根黑腐病烟草rbrr1抗性位点基因检测中的用途，所述功能分子标记为如上述。32.采用上述功能分子标记点，能有效检测烟草植株中是否含有该抗病位点。所用检测方法中其他未详述步骤和操作按现有基因位点检测常规pcr扩增操作进行，在此不累述。33.本发明的另一方面还提供了一种烟草rbrr1抗性位点基因的检测方法，包括以下步骤：34.以待检测烟草植株基因组dna为模板，采用如上述的引物对对所述模块进行扩增，扩增结果通过电泳显示指示待检测烟草植株是否含有烟草rbrr1抗性位点基因。35.所用检测方法中其他未详述步骤和操作按现有基因位点检测常规pcr扩增操作进行，在此不累述。36.本发明的另一方面还提供了一种烟草辅助育种方法，方法包括：使用上述的功能分子标记对待筛选烟草植株进行rbrr1抗性位点基因检测的步骤。上述功能分子标记可用于烟草辅助育种中，对烟草植株中是否含有该抗病位点进行准确、高效、有效的检测，从而便于从大量植株中快速筛选出具有抗病基因的植株。37.本发明“功能分子标记”是指实现100％鉴定抗病基因，能有效指示共分离，并能有效消除假阳性检测结果。38.本发明能产生的有益效果包括：39.1)本发明所提供的用于烟草根黑腐病主效抗性位点筛选/检测的功能分子标记及其应用，利用实验室现有基因组和重测序数据，以云晒1号(g306)基因组为参考，通过对携带rbrr1抗病基因材料tn90和感病材料云烟87的重测序数据分析，开发出了方便检测的indel功能标记。40.2)本发明所提供的用于烟草根黑腐病主效抗性位点筛选/检测的功能分子标记及其应用，采用本发明提供的功能标记进行了分子标记辅助抗病选育，鉴定了rbrr1抗病位点导入的遗传效应，并分析了rbrr1位点在烟草核心种质资源中的分布，为下一步开展烤烟抗根黑腐病品种选育及rbrr1基因的克隆奠定基础。附图说明41.图1为抗感材料及g306的rbrr1抗性位点共分离标记ss192650和ss219555两翼序列对比示意图；其中(a)ss192650位点的序列对比；(b)ss219555位点的序列对比，图中红色框和箭头分别显示转化indel功能标记(ind192650‑f、ind219555‑3f、ind219555‑3r)所在位置。tkf2002和tkf7002分别是前人报道的抗性和感病材料；42.图2为rbrr1抗病位点不同共显性功能标记的多态性检测结果示意图；其中各功能标记点中泳道1和2为云烟87；泳道5和6为tn90；泳道3和4为f1。ind‑2，‑5，和‑6分别为brr_lg2_ind‑2，‑5，和‑6。ind192650和ind219555‑3f/3r具有良好的多态性；43.图3为rbrr1抗病位点brr_lg2_ind‑5标记优化检测结果示意图；泳道a、f1和b分别为云烟87、杂交基因型和tn90；泳道1～15为bc3f2单株；图中箭头指示云烟87扩增条带所在位置，检测结果显示44.brr_lg2_ind‑5‑2f/2r，‑2f/4r和‑2f/5r均具有良好的多态性，其中brr_lg2_ind‑5‑2f/4r的条带比较清晰；45.图4为云烟87和tn90以及bc3f2群体中rbrr1+和rbrr1‑基因型根黑腐抗性对比；图中(a)rbrr1位点在染色体上的位置(以g306为参考基因组)；(b)rbrr1位点的遗传定位，数据参考(qin et al.2018)；(c)rbrr1位点定位区间各分子标记的物理位置锚定；(d)bc3f2世代rbrr1位点的两种重组类型；(e)bc3f2世代rbrr1+和rbrr1‑基因型单株接菌根串珠霉菌后表型(左)和根黑腐抗性对比(右)。比例尺＝2cm；**p;0.01，t测验；46.图5为rbrr1抗病位点显性功能标记在不同材料中的扩增结果；其中，sota2、burley21和kentucky14为白肋烟，hicks broad leaf、virginia gold和yellow special为国外烤烟。具体实施方式47.本发明公开的功能标记分子可通过采用上述引物对，以烟草植株dna基因序列为模块采用pcr高通量扩增后，指示烟草植株中是否包含烟草rbrr1抗性位点。48.本发明中引物对为选自引物对组a或b，引物对组a为选自以下引物对中的任意两个：49.第一引物对由ind192650‑f、ind192650‑r组成；50.第二引物对由ind219555‑3f、ind219555‑3r组成；51.第三引物对由brr_lg2_ind‑2‑f、brr_lg2_ind‑2‑r组成；52.第四引物对由brr‑lg2‑ind‑5‑2f、brr‑lg2‑ind‑5‑4r组成；53.引物对组b为选自以下引物对中的任意两个：54.第五引物对由dm192650‑2f、dm192650‑2r组成；55.第六引物对由dm219555‑2f、dm219555‑2r组成；56.第七引物对由dm5‑1f、dm5‑2r组成。57.本领域技术人员熟知，在上述引物对所示核苷酸序列中,可在其5'端或3'端分别增加1～10个碱基,所增加的碱基类型可根据烟草基因组dna上与上述引物对相匹配区域的碱基类型，并依据碱基配对原则来确定。58.1.1材料59.所用材料为携带有rbrr1抗性位点的普通烟草tn90与烤烟主栽品种云烟87。60.前期研究中，以tn90为供体亲本、以云烟87为受体亲本，通过连续回交已经获得74个bc3f2家系。61.tn90/云烟87的bc3f2群体用于新开发/转化标记的验证，并结合分子标记辅助选育(molecular marker assisted selection，mas)方法筛选bc3f2群体中纯合抗性位点基因型和非抗性纯合基因型单株用于抗性表型鉴定。62.以337份不同烟草类型种质资源为研究材料，分析抗性位点的等位基因频率分布。63.本实施例中所用两份亲本材料和337份种质资源来源于云南省烟草农业科学研究院种质库。64.病原菌根串珠霉菌(thielaviopsis basicola)由西南大学植物保护学院窦彦霞老师提供。65.1.2引物设计及优化66.ss219555和ss192650标记是rbrr1抗性位点的两个共分离caps标记。根据前人的研究结果(qin et al.2018)，下载根黑腐病抗病和感病材料中ss219555和ss192650标记所在位点的扩增序列。67.利用软件dnaman(version 9.0.1.116)分别对每个标记所在位点抗、感材料中的等位序列进行序列对比分析，根据序列的插入缺失(indel)差异设计开发标记，将两个共分离caps标记转化成免酶切的indel类型标记。68.根据实验要求，可按现有技术中常用基因序列测定方法对所用云晒1号(g306)、普通烟草tn90和云烟87进行的基因测序。本发明采用已获取的云晒1号(g306)的基因组序列、普通烟草tn90和云烟87的重测序序列。69.以云晒1号为参考序列，通过blast方法锚定ss219555和ss192650标记之间的物理区段，分析tn90和云烟87重测序数据以获取两者之间的序列差异信息。随后根据两个caps标记所在物理区间内的indel差异位点设计开发标记。70.参照[1]吴迷,汪念,沈超,等.基于重测序的陆地棉indel标记开发与评价[j].作物学报,2019,045(002):196‑203.中公开的方法进行indel分子标记开发，并结合此部分内容进行操作：即：筛选出插入/缺失大于等于5bp的indel位点，提取indel位点上下游125bp序列并生成fasta文件，使用batchprimer3(v1.0,https://wheat.pw.usda.gov/demos/batchprimer3/)批量设计引物，引物长度范围为18～27bp、最适长度21bp，gc含量范围在40～60％之间，扩增片段150‑200bp、最优扩增片段180bp。之后，检查返回的设计结果，并对需要的地方进行手动调整。[0071]为了提高特异性，本发明通过对brr_lg2_ind‑5标记进行pcr扩增测序和引物设计优化。[0072]1.3dna提取、pcr扩增和电泳检测[0073]本试验选用天根公司植物基因组提取试剂盒(目录号：dpyc3‑50，tiangen biotech(beijing)co.，ltd)提取供试材料基因组dna，具体操作步骤按说明书进行。[0074]pcr反应体系：dreamtaq green pcr master mix(2×)(dntp浓度为0.4mmol/l,mg2+浓度为4mmol/l)5μl，正、反引物(2μmol/l)各1μl，dna模板2μl(浓度50～100ng/μl)，补充ddh2o到10μl。pcr程序设置：95℃变性30s，退火温度是55℃持续30s，72℃延伸30s，设置35个循环，最后72℃延伸5min。pcr产物用8％的聚丙烯酰胺凝胶或1％的琼脂糖电泳检测。[0075]1.4烟草根黑腐抗性鉴定[0076]本论文按照chen等描述的根部扩散法(root irradiation)进行烟草根黑腐病抗性鉴定(chen et al.2020)。在人工气候室内采用灭菌土培养烟草供试材料，设定30±1℃/28±1℃的日/夜温度、85‑90％的相对湿度以及14h的光照时间，直至第四片真叶出现。[0077]同时，将单孢分离得到的串珠霉(t.basicola)病原菌接种到马铃薯葡萄糖琼脂(potato dextrose agar,pda)平板上，在25℃下暗箱培养10～15天，收获分生孢子后，用无菌水配成孢子浓度为107孢子/ml的悬浮液。[0078]利用无菌剪刀对四叶期烟苗的茎基两侧根系进行轻微伤害，立即灌根接种20ml的孢子悬浮液于烟株茎基部，并随后将供试烟苗放置于人工气候室进行培养；在接种后第15天(days post‑inoculation,dpi)，此时第一片枯萎叶片出现，用清水洗净烟株根系，进行病害等级鉴定。[0079]所有对照处理用灭菌水进行。[0080]参照国家烟草行业标准gb/t 23222‑2008对烟草根黑腐分级标准中根部症状的描述，主要针对根系表型观察，将病害分为0、1、3、5、7、9共六级；其中，[0081]0级：无病，主、侧根系均无明显坏死；[0082]1级：少数主根和个别侧根坏死呈特异黑色；[0083]3级：半数主根及少数侧根坏死呈黑色；[0084]5级：大部分根系坏死呈黑色；[0085]7级：根系全部坏死呈黑色，茎基部受害明显；[0086]9级：病株基本枯死。[0087]根据评级结果计算病情指数(disease severity index,dsi)，代表根黑腐病的发病程度。计算公式如下：[0088]dsi＝[∑(各级代表值×各级病株数)]/(最高级值×调查总株数)×100[0089]试验设计采取随机区组设计，共设置三个生物学重复，每个重复每个基因型包含15个单株。每个重复独立计算病情指数。[0090]2结果与分析[0091]2.1rbrr1抗性位点共显性功能标记的转化、开发及优化[0092]rbrr1是一个单显性广谱性根黑腐抗病基因，前人已完成其两个共分离caps标记ss219555和ss192650的两翼序列测序(qin et al.2018)。[0093]在此基础上，本发明首先利用ss219555和ss192650标记所在位点在抗病和感病材料中的dna扩增序列以及云晒1号参考基因组对应序列进行序列对比分析。[0094]结果显示，在ss192650标记位点817bp的序列区间内共检测到50个snp和5个indel差异位点，其中有两个indel位点序列差异大于5bp；在ss219555标记位点629bp的序列区间内共检测到9个snp和4个indel差异位点，其中有两个indel位点序列差异大于5bp(图1)。[0095]根据ss192650位点和ss219555位点的indel差异信息分别设计1和3对引物(图1；表1)。利用含有rbrr1抗病基因的材料tn90和不含该抗病位点的材料云烟87对引物进行差异分析，其中ind192650和ind219555‑3f/3r具有良好的多态性(图2)。[0096]以云晒1号为参考序列，通过blast方法将ss219555和ss192650标记锚定到其chr2染色体上，并将两者之间的物理区段(191,732,059‑193,313,373)作为目标区间。[0097]利用tn90和云烟87重测序数据以及云晒1号参考基因组序列对目标区间开发indel分子标记，在目标区间内共获得6个indel标记，分别命名为bbr‑lg2‑indel‑1～6(表1)。[0098]对这些标记进行筛选，共有3个标记具有多态性，分别是brr_lg2_ind‑2，‑5，和‑6，其中brr_lg2_ind‑2和‑5呈现较好的多态性(图2)。[0099]在后续的标记使用过程中，brr_lg2_ind‑5在分离群体中表现出三种以上带型，不利于准确读取基因型。因此，进一步对brr_lg2_ind‑5标记进行优化，在多态性位点两侧分别增设1条正引物和4条反引物(表1)。[0100]表1转化与开发的rbrr1抗性位点共显性标记[0101][0102][0103]通过brr_lg2_ind‑5标记多条正反引物的组合及筛选，检测结果显示brr_lg2_ind‑5‑2f/2r，‑2f/4r和‑2f/5r均具有良好的多态性，其中brr_lg2_ind‑5‑2f/4r的条带比较清晰(参见图3)。然而，在tn90/云烟87的分离群体中进行检测时，虽然三个标记能够检测出rbrr1抗性位点，但是三个标记均表现为显性标记特征，即所有单株基因型均具有云烟87的带型(图3)。[0104]综上所述，将ind192650，brr_lg2_ind‑2，brr_lg2_ind‑5‑2f/4r和ind219555‑3f/3r作为rbrr1抗性位点共显性功能标记用于后续分析。[0105]2.2分子标记辅助选择选育[0106]本研究院前期研究已获得以tn90为供体亲本、云烟87位轮回亲本的74个bc3f2家系。选取每个家系27个单株混合取样用于基因型检测，并利用上述转化或开发的标记进行分子标记辅助选择，共发现14个bc3f1世代单株含有抗病位点(杂合基因型)。同样，通过前景选择，从这14个衍生的bc3f2家系中获得含有纯合rbrr1抗病位点的单株，这些高世代改良株系可用作后续的抗病育种。[0107]2.3抗性位点的遗传效应分析[0108]为了分析rbrr1抗病位点在导入云烟87后的遗传效应，选取云烟87和tn90衍生的一个含有抗病位点的bc3f2群体进行研究。选用上述4个功能标记进行基因型分析，在285个bc3f2单株中，目标区间基因型为纯合rbrr1抗病位点基因型(记作rbrr1+)、纯合感病基因型的材料(记作rbrr1‑)及杂合基因型的个体数目分别为60、76和149个，且没有检测到遗传交换单株。[0109]选取rbrr1+和rbrr1‑基因型单株在温室条件下进行根黑腐抗病性鉴定。同时，在同样条件下考察tn90和云烟87的根黑腐病抗性。[0110]结果表明，云烟87平均病情指数为59.26％，tn90所有测试单株均未发病，其病情指数为0。对应的rbrr1+和rbrr1‑在根黑腐病的抗性上表现出极显著差异；rbrr1+材料的根黑腐抗病平均病情指数为20.39％，rbrr1‑材料的根黑腐抗病平均病情指数为72.35％(图4)。[0111]2.4根黑腐抗性位点的单倍型分析[0112]为了进一步探究rbrr1抗病位点在自然群体中的分布，利用4个indel标记对包含不同调制类型的普通烟草(sstt)材料进行基因分型，以检测rbrr1抗病位点的等位变异频率。[0113]结果表明，在337份普通烟草核心种质资源中，共检测到1个烟草材料(sota2)的基因型同tn90基因型一致(图5)，基因型频率为0.59％(表2)。[0114]表2不同材料中rbrr1抗病位点的等位变异频率[0115][0116]表性鉴定结果显示，sota2与tn90具有相似的根黑腐抗性表型。由此推测，该抗病位点尚未在我国烟草生产中普及应用。[0117]由于上述4个indel共显性标记在自然群体中均表现为只有两种基因型，因此，将其中3个indel标记转化为检测rbrr1抗病位点的特异显性标记，分别为dm192650‑2f/2r、dm219555‑2f/2r和dm5‑1f/2r，为kasp转化及高通量应用奠定基础(表3；图5)检测rbrr1抗病位点的特异显性标记，分别为dm192650‑2f/2r、dm219555‑2f/2r和dm5‑1f/2r。[0118]表3 rbrr1抗性位点显性功能标记[0119][0120]3讨论[0121]3.1rbrr1位点功能标记转化、开发与利用[0122]分子标记辅助选择回交育种已经在农作物抗病育种、品质育种中得到广泛应用(singh et al.2015；vida et al.2009；zhao et al.2012；zhou et al.2003)。在分子标记辅助选择过程中，分子标记的可靠性、检测时所需dna的质和量、分子标记的检测流程、其多al.1981)，而紧密连锁标记将有助于识别重组事件，打破rbrr1和n.debneyi邻近基因之间潜在的不利连锁。理论上，位于n.debneyi易位/插入片段的标记在纯合抗病材料中的扩增带型应该一致，而在普通烟草中可能具有扩增多态性。[0133]本研究结果表明，ind192650，brr_lg2_ind‑2，brr_lg2_ind‑5‑2f/4r和ind219555‑3f/3r标记在自然群体中扩增具有二元性，只有两种pcr扩增表现。[0134]因此，又将其中三个标记转化为检测rbrr1抗病位点的显性标记，为后续的kasp标记转化及高通量基因型鉴定奠定基础。[0135]为使本发明实施方式的目的、技术方案和优点更加清楚，下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整地描述。[0136]实施例[0137]以下实施例中所用物料及仪器如无特殊说明，均通过商业渠道购得。[0138]1.1材料[0139]所用材料为携带有rbrr1抗性位点的普通烟草tn90与烤烟主栽品种云烟87。[0140]前期研究中，以tn90为供体亲本、以云烟87为受体亲本，通过连续回交已经获得74个bc3f2家系。[0141]tn90/云烟87的bc3f2群体用于新开发/转化标记的验证，并结合分子标记辅助选育(molecular marker assisted selection，mas)方法筛选bc3f2群体中纯合抗性位点基因型和非抗性纯合基因型单株用于抗性表型鉴定。[0142]以337份不同烟草类型种质资源为研究材料，分析抗性位点的等位基因频率分布。[0143]本实施例中所用两份亲本材料和337份种质资源来源于云南省烟草农业科学研究院种质库。[0144]病原菌根串珠霉菌(thielaviopsis basicola)由西南大学植物保护学院窦彦霞老师提供。[0145]1.2引物设计及优化[0146]ss219555和ss192650标记是rbrr1抗性位点的两个共分离caps标记。根据前人的研究结果(qin et al.2018)，下载根黑腐病抗病和感病材料中ss219555和ss192650标记所在位点的扩增序列。[0147]利用软件dnaman(version 9.0.1.116)分别对每个标记所在位点抗、感材料中的等位序列进行序列对比分析，根据序列的插入缺失(indel)差异设计开发标记，将两个共分离caps标记转化成免酶切的indel类型标记。[0148]根据实验要求，可按现有技术中常用基因序列测定方法对所用云晒1号(g306)、普通烟草tn90和云烟87进行的基因测序。本发明采用已获取的云晒1号(g306)的基因组序列、普通烟草tn90和云烟87的重测序序列。[0149]以云晒1号为参考序列，通过blast方法锚定ss219555和ss192650标记之间的物理区段，分析tn90和云烟87重测序数据以获取两者之间的序列差异信息。随后根据两个caps标记所在物理区间内的indel差异位点设计开发标记。[0150]参照[1]吴迷,汪念,沈超,等.基于重测序的陆地棉indel标记开发与评价[j].作物学报,2019,045(002):196‑203.中公开的方法进行indel分子标记开发，并结合此部分内容进行操作：即：筛选出插入/缺失大于等于5bp的indel位点，提取indel位点上下游125bp序列并生成fasta文件，使用batchprimer3(v1.0,https://wheat.pw.usda.gov/demos/batchprimer3/)批量设计引物，引物长度范围为18～27bp、最适长度21bp，gc含量范围在40～60％之间，扩增片段150‑200bp、最优扩增片段180bp。之后，检查返回的设计结果，并对需要的地方进行手动调整。[0151]为了提高特异性，本发明通过对brr_lg2_ind‑5标记进行pcr扩增测序和引物设计优化。[0152]1.3dna提取、pcr扩增和电泳检测[0153]本试验选用天根公司植物基因组提取试剂盒(目录号：dpyc3‑50，tiangen biotech(beijing)co.，ltd)提取供试材料基因组dna，具体操作步骤按说明书进行。[0154]pcr反应体系：dreamtaq green pcr master mix(2×)(dntp浓度为0.4mmol/l,mg2+浓度为4mmol/l)5μl，正、反引物(2μmol/l)各1μl，dna模板2μl(浓度50～100ng/μl)，补充ddh2o到10μl。pcr程序设置：95℃变性30s，退火温度是55℃持续30s，72℃延伸30s，设置35个循环，最后72℃延伸5min。pcr产物用8％的聚丙烯酰胺凝胶或1％的琼脂糖电泳检测。[0155]1.4烟草根黑腐抗性鉴定[0156]本论文按照chen等描述的根部扩散法(root irradiation)进行烟草根黑腐病抗性鉴定(chen et al.2020)。在人工气候室内采用灭菌土培养烟草供试材料，设定30±1℃/28±1℃的日/夜温度、85‑90％的相对湿度以及14h的光照时间，直至第四片真叶出现。[0157]同时，将单孢分离得到的串珠霉(t.basicola)病原菌接种到马铃薯葡萄糖琼脂(potato dextrose agar,pda)平板上，在25℃下暗箱培养10～15天，收获分生孢子后，用无菌水配成孢子浓度为107孢子/ml的悬浮液。[0158]利用无菌剪刀对四叶期烟苗的茎基两侧根系进行轻微伤害，立即灌根接种20ml的孢子悬浮液于烟株茎基部，并随后将供试烟苗放置于人工气候室进行培养；在接种后第15天(days post‑inoculation,dpi)，此时第一片枯萎叶片出现，用清水洗净烟株根系，进行病害等级鉴定。[0159]所有对照处理用灭菌水进行。[0160]参照国家烟草行业标准gb/t 23222‑2008对烟草根黑腐分级标准中根部症状的描述，主要针对根系表型观察，将病害分为0、1、3、5、7、9共六级；其中，[0161]0级：无病，主、侧根系均无明显坏死；[0162]1级：少数主根和个别侧根坏死呈特异黑色；[0163]3级：半数主根及少数侧根坏死呈黑色；[0164]5级：大部分根系坏死呈黑色；[0165]7级：根系全部坏死呈黑色，茎基部受害明显；[0166]9级：病株基本枯死。[0167]根据评级结果计算病情指数(disease severity index,dsi)，代表根黑腐病的发病程度。计算公式如下：[0168]dsi＝[∑(各级代表值×各级病株数)]/(最高级值×调查总株数)×100[0169]试验设计采取随机区组设计，共设置三个生物学重复，每个重复每个基因型包含15个单株。每个重复独立计算病情指数。[0170]2结果与分析[0171]2.1rbrr1抗性位点共显性功能标记的转化、开发及优化[0172]rbrr1是一个单显性广谱性根黑腐抗病基因，前人已完成其两个共分离caps标记ss219555和ss192650的两翼序列测序(qin et al.2018)。[0173]在此基础上，本发明首先利用ss219555和ss192650标记所在位点在抗病和感病材料中的dna扩增序列以及云晒1号参考基因组对应序列进行序列对比分析。[0174]结果显示，在ss192650标记位点817bp的序列区间内共检测到50个snp和5个indel差异位点，其中有两个indel位点序列差异大于5bp；在ss219555标记位点629bp的序列区间内共检测到9个snp和4个indel差异位点，其中有两个indel位点序列差异大于5bp(图1)。[0175]根据ss192650位点和ss219555位点的indel差异信息分别设计1和3对引物(图1；表1)。利用含有rbrr1抗病基因的材料tn90和不含该抗病位点的材料云烟87对引物进行差异分析，其中ind192650和ind219555‑3f/3r具有良好的多态性(图2)。[0176]以云晒1号为参考序列，通过blast方法将ss219555和ss192650标记锚定到其chr2染色体上，并将两者之间的物理区段(191,732,059‑193,313,373)作为目标区间。[0177]利用tn90和云烟87重测序数据以及云晒1号参考基因组序列对目标区间开发indel分子标记，在目标区间内共获得6个indel标记，分别命名为bbr‑lg2‑indel‑1～6(表1)。[0178]对这些标记进行筛选，共有3个标记具有多态性，分别是brr_lg2_ind‑2，‑5，和‑6，其中brr_lg2_ind‑2和‑5呈现较好的多态性(图2)。[0179]在后续的标记使用过程中，brr_lg2_ind‑5在分离群体中表现出三种以上带型，不利于准确读取基因型。因此，进一步对brr_lg2_ind‑5标记进行优化，在多态性位点两侧分别增设1条正引物和4条反引物(表1)。[0180]表1转化与开发的rbrr1抗性位点共显性标记[0181][0182]通过brr_lg2_ind‑5标记多条正反引物的组合及筛选，检测结果显示brr_lg2_ind‑5‑2f/2r，‑2f/4r和‑2f/5r均具有良好的多态性，其中brr_lg2_ind‑5‑2f/4r的条带比较清晰(参见图3)。然而，在tn90/云烟87的分离群体中进行检测时，虽然三个标记能够检测出rbrr1抗性位点，但是三个标记均表现为显性标记特征，即所有单株基因型均具有云烟87的带型(图3)。[0183]综上所述，将ind192650，brr_lg2_ind‑2，brr_lg2_ind‑5‑2f/4r和ind219555‑3f/3r作为rbrr1抗性位点共显性功能标记用于后续分析。[0184]2.2分子标记辅助选择选育[0185]本研究院前期研究已获得以tn90为供体亲本、云烟87位轮回亲本的74个bc3f2家系。选取每个家系27个单株混合取样用于基因型检测，并利用上述转化或开发的标记进行分子标记辅助选择，共发现14个bc3f1世代单株含有抗病位点(杂合基因型)。同样，通过前景选择，从这14个衍生的bc3f2家系中获得含有纯合rbrr1抗病位点的单株，这些高世代改良株系可用作后续的抗病育种。[0186]2.3抗性位点的遗传效应分析[0187]为了分析rbrr1抗病位点在导入云烟87后的遗传效应，选取云烟87和tn90衍生的一个含有抗病位点的bc3f2群体进行研究。选用上述4个功能标记进行基因型分析，在285个bc3f2单株中，目标区间基因型为纯合rbrr1抗病位点基因型(记作rbrr1+)、纯合感病基因型的材料(记作rbrr1‑)及杂合基因型的个体数目分别为60、76和149个，且没有检测到遗传交换单株。[0188]选取rbrr1+和rbrr1‑基因型单株在温室条件下进行根黑腐抗病性鉴定。同时，在同样条件下考察tn90和云烟87的根黑腐病抗性。[0189]结果表明，云烟87平均病情指数为59.26％，tn90所有测试单株均未发病，其病情指数为0。对应的rbrr1+和rbrr1‑在根黑腐病的抗性上表现出极显著差异；rbrr1+材料的根黑腐抗病平均病情指数为20.39％，rbrr1‑材料的根黑腐抗病平均病情指数为72.35％(图4)。[0190]2.4根黑腐抗性位点的单倍型分析[0191]为了进一步探究rbrr1抗病位点在自然群体中的分布，利用4个indel标记对包含不同调制类型的普通烟草(sstt)材料进行基因分型，以检测rbrr1抗病位点的等位变异频率。[0192]结果表明，在337份普通烟草核心种质资源中，共检测到1个烟草材料(sota2)的基因型同tn90基因型一致(图5)，基因型频率为0.59％(表2)。[0193]表2不同材料中rbrr1抗病位点的等位变异频率[0194][0195]表性鉴定结果显示，sota2与tn90具有相似的根黑腐抗性表型。由此推测，该抗病位点尚未在我国烟草生产中普及应用。[0196]由于上述4个indel共显性标记在自然群体中均表现为只有两种基因型，因此，将其中3个indel标记转化为检测rbrr1抗病位点的特异显性标记，分别为dm192650‑2f/2r、dm219555‑2f/2r和dm5‑1f/2r，为kasp转化及高通量应用奠定基础(表3；图5)检测rbrr1抗病位点的特异显性标记，分别为dm192650‑2f/2r、dm219555‑2f/2r和dm5‑1f/2r。[0197]表3 rbrr1抗性位点显性功能标记[0198][0199]3讨论[0200]3.1rbrr1位点功能标记转化、开发与利用[0201]分子标记辅助选择回交育种已经在农作物抗病育种、品质育种中得到广泛应用(singh et al.2015；vida et al.2009；zhao et al.2012；zhou et al.2003)。在分子标记辅助选择过程中，分子标记的可靠性、检测时所需dna的质和量、分子标记的检测流程、其多态性水平和成本被认为是限制分子标记广泛应用的五个因素(collard and mackill 2008)。[0202]然而，酶切步骤使得caps标记使用流程更加冗长，增加了caps标记的使用成本和试验难度，并限制其在高通量基因型检测领域里面的应用(an et al.2021；shavrukov 2016)。caps标记的这些特点使其特别适合于相对较小的实验室来定位克隆基因组目标区段或目标基因或多态性分析研究(shavrukov 2016)。[0203]本论文在前人研究基础上，锚定rbrr1抗性位点的两个共分离caps标记ss219555和ss192650的基因组所在物理区间，并随后根据ss219555和ss192650标记的两翼序列以及为云烟87和tn90的重测序数据，基于抗、感材料序列的indel多态性转化、开发获得多个共显性、免酶切的分子标记。[0204]进而，利用ind192650，brr_lg2_ind‑2，brr_lg2_ind‑5‑2f/4r和ind219555‑3f/3r进行分子标记辅助选择，筛选出根黑腐抗性显著提高的bc3f2世代rbrr1位点纯合单株，这些材料可用于后续的抗根黑腐品种选育。[0205]这些结果说明，本发明通过转化或开发的4个indel标记可用于根黑腐抗性育种的分子标记辅助选择和后续的抗病基因定位。[0206]易位系，尤其是小片段易位系，是目前将野生资源优异基因导入生产作物最为理想的方式。一般来说，具有较小的外来片段的易位系在遗传上更稳定，产生有害影响的可能性更小。[0207]例如，在主粮作物小麦的研究中发现，冰草6p染色体小片段的插入易位系pubing3035具有显著提高的粒重和穗长，并根据连锁图谱将易位断点定位在小麦1a染色体短臂的着丝粒附近(zhang et al.2015)。与此类似，前人同样通过遗传连锁图谱将rbrr1基因定位在17号染色体(k326基因组)的ss127301和ss166247标记之间(qin et al.2018)；同时，本研究将rbrr1基因共分离标记ss219555和ss192650锚定到chr2染色体(g306基因组，chr2长＝201,712,017)长臂末端，表明rbrr1基因是以片段易位而非整臂易位的形式被导入到普通烟草中。通常情况下，位于外源染色体片段的标记表现为显性标记。[0208]然而，在本研究的标记开发鉴定过程中，所转化和开发的标记几乎均为共显性标记。前人测序还发现携带rbrr1基因的抗性材料中ss219555和ss192650标记两翼序列与n.debneyi一致(qin et al.2018)。同时，在对brr_lg2_ind‑5标记的改良过程中，brr_lg2_ind‑5‑2f/2r，‑2f/4r和‑2f/5r在分离群体的所有单株中均能够扩增云烟87的带型。综上，推测rbrr1基因是以一个外源片段插入易位的形式导入到普通烟草，才能使得这些标记既能够在扩增出外源片段条带的同时也能够扩增出普通烟草的条带。然而，关于携带rbrr1基因的抗性材料的基因组组成鉴定尚需要进一步研究。[0209]3.2rbrr1位点在烟草育种中的利用价值[0210]由于早期栽培所用普通烟草品种对根黑腐病没有免疫能力，烟草根黑腐的影响给全世界烟草种植者造成了巨大的损失。rbrr1抗病位点在白肋烟上的首先应用使得根黑腐病情得到了很大的控制(朱贤朝2002)。本研究发现，在我国主栽烤烟品种、国内烤烟、国内晒烟和香料烟中均未携带rbrr1抗病位点，仅在引进的白肋烟中发现两个材料携带有该抗病位点，表明rbrr1抗病位点在我国尚未进行广泛利用。[0211]有报道称rbrr1抗病位点和烤烟的不良农艺性状及化学特质相关(legg et al.1981)，而紧密连锁标记将有助于识别重组事件，打破rbrr1和n.debneyi邻近基因之间潜在的不利连锁。理论上，位于n.debneyi易位/插入片段的标记在纯合抗病材料中的扩增带型应该一致，而在普通烟草中可能具有扩增多态性。[0212]本研究结果表明，ind192650，brr_lg2_ind‑2，brr_lg2_ind‑5‑2f/4r和ind219555‑3f/3r标记在自然群体中扩增具有二元性，只有两种pcr扩增表现。[0213]因此，又将其中三个标记转化为检测rbrr1抗病位点的显性标记，为后续的kasp标记转化及高通量基因型鉴定奠定基础。[0214]尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。