基于辣椒基因组重测序的InDel标记开发及应用.docx
《基于辣椒基因组重测序的InDel标记开发及应用.docx》由会员分享,可在线阅读,更多相关《基于辣椒基因组重测序的InDel标记开发及应用.docx(6页珍藏版)》请在冰豆网上搜索。
基于辣椒基因组重测序的InDel标记开发及应用
基于辣椒基因组重测序的InDel标记开发及应用
郭广君,孙茜,刘金兵,潘宝贵,刁卫平,戈伟,高长洲,王述彬
【摘要】摘要:
基于全基因组重测序进行分子标记开发是一种新的方法和趋势。
本研究基于辣椒全基因组重测序数据鉴定出的InDel位点,在2号染色体上随机选择并设计了40个InDel标记,以24份辣椒种质验证其有效性。
结果显示:
40个标记可获得扩增产物,其中具有多态性的标记35个,多态率为87.5%,同时35个标记可以有效地区分24份辣椒种质的多样性和特异性。
以上结果表明,基于重测序数据开发辣椒InDel标记切实可行。
【期刊名称】江苏农业学报
【年(卷),期】2015(031)006
【总页数】7
【关键词】关键词:
辣椒;重测序;InDel标记;种质鉴定
辣椒(Capsicumspp.)属于茄科(Solanaceae)辣椒属(Capsicum)作物,作为一种重要的蔬菜作物和调味品,广泛栽培于南北美、亚洲、欧洲和大洋洲等世界各地。
随着消费者需求的多样性,利用传统方法在辣椒的遗传改良上取得了很大进展,但是分子标记[1]和第二代测序技术[2]的出现为辣椒的遗传改良提供了更大的机遇[3]。
目前辣椒分子标记的开发和应用主要集中于利用公用数据或者基因组文库开发SSR标记[4-8]。
随着高通量测序技术的发展和辣椒基因组数据的发布,使得辣椒分子标记的开发和应用进入了一个新的阶段[9-10],第三代分子标记SNP标记和InDel标记的开发和应用在辣椒中逐渐广泛[11-14]。
插入/缺失多态性(Insertion/deletion,InDel)标记是由于等位基因位点处的DNA序列在不同个体间发生了核苷酸片段的插入/缺失而产生的长度多态性变异[15]。
根据目标位点两侧的序列设计特异引物进行PCR扩增,扩增片段的长度多态性即是InDel标记。
在整个基因组中,InDel标记的多态性频率仅次于SNP标记,远高于SSR标记。
InDel标记多态性可通过聚合酶链式反应(PCR)和琼脂糖凝胶电泳或非变性聚丙烯酰胺凝胶电泳等简单的步骤达到基因分型的目的[16-17],相对于SNP标记,其标记的设计和检测更为简易,应用也更为方便[18]。
此外,InDel标记的另一优点是不可能在同1个基因组位点上发生同样长度的2个InDel位点突变,所以具有共同InDel突变位点的2个物种间代表其具有血缘一致性(Identity-by-descent)[19]。
综上所述,In-Del标记为共显性标记,具有较好的稳定性和较为丰富的多态性,是高分辨率图谱构建,关联分析和基因定位等方面最佳的候选标记[20]。
但是,辣椒中已开发和应用的InDel标记数量非常少,本研究的目的在于基于辣椒基因组重测序数据开发InDel标记,利用不同遗传背景的辣椒种质对开发出的InDel标记进行验证,筛选获得有效的InDel标记,进一步探讨InDel标记在种质鉴定中的应用。
1材料与方法
1.1InDel位点的挖掘
本研究利用重测序数据与参考基因组进行比对查找InDel位点。
参考基因组版本为pepper.V.1.55,下载网址:
http:
//peppergenome.snu.ac.kr/download.php。
采用illuminaHiSeqTM2500测序平台对一年生辣椒(Capsicumannuum)“G29”和灌木状辣椒(Capsicum.frutescents)“G108”进行测序。
将测序得到的原始Reads进行质量评估并过滤得到CleanReads。
利用bwa软件比对定位CleanReads在参考基因组上的位置。
InDel的检测主要使用GATK软件工具包实现。
根据CleanReads在参考基因组的定位结果,使用Samtools进行去重复(Markduplicates),使用GATK软件进行局部比对(Localrealignment)和碱基质量值校正(Baserecalibration)等预处理,以保证检测结果的准确性,再使用GATK软件进行变异检测和校准,选取可靠的In-Del位点[21-22]。
利用SnpEff软件根据InDel位点在参考基因组上的位置信息,对比参考基因组的基因、CDS位置等信息对其进行注释和预测其影响。
1.2InDel引物设计和分析
根据重测序数据中预测到InDel位点,筛选出插入/缺失碱基数为5~10个碱基的位点。
基于辣椒基因组序列,将筛选出的位点定位在基因组上,取InDel位点两翼各200bp碱基长度,共401bp长度进行引物设计。
采用Primer3.0进行引物设计,上游引物设计范围在1~195bp,下游引物设计范围为205~401bp,产物大小为150~250bp,退火温度为52~60℃。
优化后的扩增反应采用10μl反应体系,包括20ng模板DNA、5μl2×Mix混合液(含Mg2+的10×PCRbuffer、2.5mmol/L的dNTPs、0.5UTaqDNA聚合酶)、每个引物终浓度为0.2μmol/L。
反应程序为94℃预变性4min;94℃变性30s,55℃复性30s,72℃延伸40s,循环30次;72℃延伸7min,置于4℃下保存。
扩增产物用10%非变性聚丙烯酰胺凝胶电泳,银染显色,统计分析条带类型。
1.3供试辣椒种质及其DNA提取
用于检测标记有效性的辣椒材料来源于国外和国内多个科研单位,其中包括一年生辣椒(Capsicumannuum)21份、灌木状辣椒(Capsicumfrutescents)2份,下垂辣椒(Capsicumbaccatum)1份。
除11份辣椒品种外,其余13份以G开头的材料为项目组引进的国外种质资源(表1)。
采用改良CTAB法提取叶片的基因组DNA[23],经琼脂糖凝胶电泳和微量分光光度计检测,DNA浓度调至50ng/μl备用。
2结果与分析
2.1辣椒InDel标记的鉴定
通过与参考基因组序列比对,一年生辣椒品种G29和灌木状辣椒品种G108的全基因组范围内共检测到的InDel数量分别为533523个和1664770个,其中编码区的InDel数量分别为1019个和2515个,具体结果见表2。
通过进一步比较分析,G29和G108之间的全基因组范围内InDel数量为1586427个。
2.2InDel标记有效性验证
根据预测的InDel位点,在2号染色体上随机选择了40个InDel位点设计引物,引物信息见表3。
以24份不同来源的辣椒材料验证标记的有效性。
经过PCR和产物检测发现,40对引物在24份辣椒材料中均能扩增出条带,扩增条带的大小与预测产物大小基本相同。
如图1为引物InDel-2-3对24份辣椒种质PCR扩增结果。
经统计,在24份材料中共扩增出70个清晰可辨的位点,其中35个标记表现出多态性,5个无多态性,多态率达到87.5%。
2.3InDel标记在品种和种质资源鉴定中的应用
利用Popgene32软件对40对引物在24份辣椒材料中的扩增结果进行分析发现,每对引物反映的基因多样性指数范围为0.08~0.50,其均值为0.20;Shannon′s多样性信息指数分布在0.17至0.80之间,其均值为0.34(表3)。
利用NTSYSpc2.10e中的UPGMA方法对24份材料进行聚类分析,由聚类结果(图2)可以看出,在遗传相似系数0.42处,24份辣椒种质分为2个大类,一类为21份一年生辣椒,另一大类为2份灌木状辣椒(C.frutescens)和1份下垂辣椒(C.baccatum)。
21份一年生辣椒中可以聚类为羊角形辣椒、牛角形(长灯笼形)辣椒和甜椒三大类别。
同时,聚类结果显示牛角形(长灯笼形)辣椒和甜椒相似度更高,同一单位品种间相似度高于不同单位间相似度。
3讨论
InDel标记为共显性标记,具有较好的稳定性和较为丰富的多态性。
虽然在整个基因组中,其多态性频率次于SNP标记[24],但是,相对于SNP标记,InDel标记在引物设计、多态性检测等方面更为简易,应用也更为方便[16-17,25]。
但是,目前辣椒中In-Del标记的开发数量较少,应用较为有限,这一现状急待改善[11]。
本研究基于全基因组重测序数据,通过与参考基因组序列比对,在一年生辣椒品种G29和灌木状辣椒品种G108的全基因组范围内检测到的InDel数量分别为533523个和1664770个。
以上结果说明利用全基因组重测序技术可以在辣椒中更好地挖掘InDel位点,为全基因组范围内的InDel设计提供了更优选择。
这一观点已经在包括辣椒在内的多个物种中得以验证[12,26-27]。
利用24份辣椒材料对2号染色体上的40个标记的有效性进行验证,发现24份材料在40个标记中共扩增出70个清晰可辨的位点,其中35个标记表现出多态性,多态率达到87.5%。
这一结果很好地验证了本研究开发出的InDel标记的有效性和较高的多态性。
目前辣椒分子标记的开发和应用主要集中于利用公用数据或者基因组文库开发SSR标记,但是大部分研究结果表明SSR标记的多态性及其应用率远低于InDel标记。
如Huang等人基于数据库中的辣椒序列开发出12个SSR标记,其中多态性标记为5个,多态率为41.7%[5];Lee等构建辣椒分子遗传连锁图谱时开发出76个SSR标记,其中具有多态性的标记为46个,多态率为60.5%[4];Yi等开发出1201个EST-SSR标记,其多态性标记只有150个,多态率仅为29.2%[6];Huang等依据EST序列开发出SSR标记755个,利用8份辣椒栽培种对其中的210个进行验证,其多态性标记为127个,多态率为60.5%[8];Shirasawa等基于辣椒EST序列设计了5751个SSR标记,利用辣椒材料对其中77个进行验证,其中60个表现出多态性,多态率为77.9%[7]。
但是,也有少数研究结果显示,InDel标记与SSR标记的多态率并无明显差异。
Li等利用InDel标记构建辣椒种内遗传图谱时显示,1000个InDel标记中251个标记可以成功应用于图谱构建,其多态率为25.1%[11]。
随后,该团队(2015年)基于SSR和InDel标记构建辣椒种间遗传图谱时,分别筛选了1038个InDel标记和674个SSR标记,其中多态性标记分别为140个和102个,多态率分别为13.5%和15.1%[12]。
利用NTSYSpc2.10e中的UPGMA方法对24份材料进行聚类分析,结果表明35个InDel标记不仅可以有效区分出21份一年生辣椒、2份灌木状辣椒和1份下垂辣椒,而且21份材料中又可以聚类为不同种类,其聚类结果与表型特征非常吻合。
这一结果表明InDel标记在种质资源的遗传分析、品种鉴定和品种保护方面具有极大的应用价值。
但同一单位育成品种的相似度很高,说明可利用的辣椒种质的遗传背景较窄,对于品种改良极为不利,需要开发和应用新的种质进行品种的选育和改良。
总的来说,基于全基因组重测序数据开发InDel标记可以有效地应用于辣椒种质鉴定,遗传图谱的构建,关联分析和基因定位等方面,是辣椒种质资源开发利用和分子标记辅助育种的高效工具。
参考文献:
[1]王述彬,刁卫平,刘金兵,等.辣椒胞质雄性不育相关基因的克隆及表达[J].江苏农业学报,2014,30(4):
833-838.
[2]徐婉莉,裴徐梨,荆赞革,等.辣椒actin基因电子克隆与生物信息学分析[J].江苏农业科学,2014,42(5):
46-48.
[3]RAMCHIARYN,KEHIEM,BRAHMAV,etal.ApplicationofgeneticsandgenomicstowardsCapsicumtranslationalresearch[J].PlantBiotechnologyReports