桃种质资源群体遗传分析及果实数字基因表达谱构建

【摘要】桃（Prunuspersica[L.]Batsch）是蔷薇科李属重要的核果类果树,因果实口感好,营养物质丰富,深受消费者喜欢。过去20年,通过传统育种方法国内外已经选出2000多个优良品种,但大多来源于少数亲本材料,遗传背景狭窄,使桃遗传多样性没有随品种资源数量的增加而增加,不利于一些优良性状的改良。国内外应用分子标记已经对桃遗传多样性进行了大量研究,但不同研究组所选标记、植物材料及分析方法不同,难以进行数据整合,所得结果主要适用于所测样本,缺乏一个全景式品种群体分析图。作为桃的原产国,我国拥有4000多年的栽培历史及丰富多样的种质资源,如何选择优良性状亲本,对桃育种进程具有重要意义。为此,本研究采用ABI3130分析系统及48个高多态性SSR标记对我国的434份桃种质资源进行基因分型,并结合前期西班牙研究组用相同标记、相同方法下所得到的224份欧美桃种质资源基因型数据,对世界范围桃种质资源进行遗传对比分析,明确世界范围内桃品种资源的遗传多样性、品种群结构、连锁不平衡程度和在连锁群上分布。选择出中国核心桃种质资源材料采用新型高密度SNP芯片检测遗传变异,并进行简单性状的全基因组关联分析。桃果实品质研究方面,果实香气及硬度是评价果实品质的重要指标,两者均是由多基因控制的复杂性状,前人采用候选基因克隆、EST文库及微阵列芯片在香气、质地基因功能研究中取得了很大进展,但受基因数量有限,其分子调控机理仍不明确。本研究采用高通量IlluminaRNA-seq技术,对‘湖景蜜露’、‘玉露’、‘锦绣黄桃’和‘中华寿桃’4个优良桃栽培品种及‘湖景蜜露’的三个成熟期进行数字基因表达谱构建,从转录组水平全面系统地分析了果实香气物质合成及果实软化的分子调控机理。主要结果如下：1.653份桃种质资源的遗传多样性分析48个SSR位点在653份桃种质资源中共扩增得到588个等位基因,平均每位点扩增得到12.25个等位基因,观察杂合度（Ho）则从0.13（PMS02）到0.63（UDP96-005）,平均0.47。东亚桃亚组的遗传多样性显著高于欧美桃亚组,引进东亚品种的可以提高欧美品种的遗传多样性水平。聚类分析可将658份材料按地理起源、驯化进程和果实性状显著区分：野生近缘种与普通桃显著分开；古老地方品种与栽培种显著分开,东亚桃和欧美桃品种群可显著分开；欧美桃和油桃显著分开,选出281份中国桃种质资源用于高密度SNP芯片构建。不同群体LD分析发现东亚品种群LD遗传衰减距离（3.36cM）比欧美品种群（5.50cM）小,东亚品种群在应用高密度SNP全基因组关联分析中有更大的优势。‘玉露’亚群的LD水平最高,地方品种群LD水平最低。东亚品种群、欧美品种群及地方品种群在第1、第4连锁群的LD位点组合比例较高,重要的果肉色泽和果肉质地、香气分别位于这两的连锁群上。2.桃品种特异性荧光SSR分子标记数据库建立及应用从全基因组范围的48个荧光SSR分子标记中选出15个高多态性标记,经不同荧光标识,最终分为3组在ABI3130测序仪测试。目前数据库包含669份材料。445个中国桃品种中有426（96.6%）个具有特异基因型组合,研究初步建立了159份中国桃核心种质资源分子标记数据库,并提出8个中国品系用作今后资源评价的参照。对17个桃品种进行系谱验证发现,我国资源圃所育成的7个品种系谱信息正确,但日本桃品种系谱信息与文献出现偏差。3.基于9K高密度SNP芯片的全基因组关联分析对281份中国桃种质资源进行9KSNP芯片构建,共得到3781个多态性SNP,最小等位基因频率均大于0.03,并选出159个中国桃核心种质资源用于资源保存、管理及全基因组关联分析。邻接法聚类结果显示：3781个SNP与48个SSR聚类结果一致,且与群体结构分析结果一致。叶腺、花型、果肉颜色、果实质地、果皮茸毛、粘离核及果形7个质量性状进行关联分析,Scaffold5上的SNPIGA544512标记与叶腺性状呈显著相关性,Scaffold1上的SNPIGA103771和Scaffold6上的SNPIGA686125与花形呈显著相关,Scaffold1上的SNPIGA93479和SNPIGA93768与果肉颜色显著相关,Scaffold4上的SNPIGA477941、SNPIGA477945、SNPIGA477951与果肉溶质／不溶质呈显著相关,Scaffold5上的SNPIGA596063与果皮茸毛显著相关,Scaffold6和Scaffold1上的SNPIGA610889釉SNPIGA124466与果实粘离核显著相关,Scaffold6上的SNPIGA698099、SNPIGA697021和SNPIGA697026与果形显著相关。4.应用数字基因表达谱挖掘桃果实香气物质合成及果实质地相关基因12个数字基因表达谱经过滤共产生89881875测序序列,平均序列长度100bp,78.5%的测序片段比对到了桃参考基因组上。对15个果实香气、软化相关基因进行Q-PCR表达水平验证,结果与RNA-seq所测结果呈现出显著相关性。等级聚类分析将‘湖景蜜露’果实成熟过程中4431个差异表达基因分为24个组,将不同品种间5541个差异表达基因分为14组。‘湖景蜜露’果实成熟过程中,与内酯类物质合成相关的脂质转移酶基因（LTP）、酰基辅酶A氧化酶基因（ACXI）及环氧化物水解酶基因（EPH）表达水平变化显著,其中LTP表达水平为先升高后降低趋势,而ACX2对C1s脂肪酸具底物特异性,推断与γ-十内酯合成相关。不同品种和不同成熟期的差异表达基因在黄酮类物质合成途径中的富集程度在所有KEGG代谢途径中最高。
【作者】李雄伟；
【导师】高中山；
【作者基本信息】浙江大学，果树学，2013，博士
【关键词】桃；SSR；遗传多样性；连锁不平衡；单核苷酸多态性；全基因组关联分析；数字基因表达谱；香气；

【参考文献】
[1]吴悠黄艳.石墨烯产业化渐行渐近[N].中国有色金属报,2014-12-20007.
[2]韩家新,王家华.一种以相关性确定条件属性的决策树[J].微机发展,2003,05:38-39+42.
[3]崔恒进.β-Catenin调控间充质干细胞向HGF的定向迁移[D].苏州大学,细胞生物学,2013,硕士.
[4]余伟林.距舟楔关节融合治疗Muller-Weiss病的疗效分析[D].苏州大学,临床医学（专业学位）,2014,硕士.
[5]罗少宏.天然棕色棉色素组分、结构及性能研究[D].浙江理工大学,纺织品材料与纺织品设计,2014,硕士.
[6]王新辉.基于Web的法院综合管理系统的设计与实现[D].中南大学,2004.
[7]易巍.锁胸筋膜及其毗邻结构的临床解剖学结构[D].吉林大学,人体解剖与组织胚胎学,2012,硕士.
[8]李玥.纳米硅薄膜晶体管压力传感器制作及特性研究[D].黑龙江大学,微电子学与固体电子学,2013,硕士.
[9]吴保伟.网络体育新闻专题策划的共性特征研究[D].首都体育学院,体育新闻,2013,硕士.
[10]吴统顺.散装炸药在地下金属矿山的应用[J].长沙矿山研究院季刊.1983(03)
[11]林伦毅.葛根汤配合推拿手法治疗颈型颈椎病之临床研究[D].北京中医药大学,中医临床基础,2014,硕士.
[12]杨金标.系统仿真的基本概念和基本技术(下)——电子计算机原理与应用讲座之八[J].冶金自动化,1982,04:44-48.
[13]刘爽.多元智能理论和小学英语作业设计研究[D].杭州师范大学,学科教育（专业学位）,2012,硕士.
[14]施英吉.装饰工程施工中的项目管理研究[D].吉林大学,项目管理,2014,硕士.
[15]乔林波.面向高速网络的复杂网络威胁模式提取方法研究与实现[D].国防科学技术大学,计算机科学与技术,2012,硕士.
[16]邵燕.部分运动队90后运动员榜样教育教育的现状与对策研究[D].北京体育大学,思想政治教育,2013,硕士.
[17]谢少强.重构的天人：中国前卫艺术的现代性演进[D].东北师范大学,美术学,2012,硕士.
[18]李大伟.风景艺术的审美价值与真实性研究[D].河北师范大学,学科教学,2012,硕士.
[19]覃斌健.汉英小说翻译中语用充实的研究[D].广东外语外贸大学,翻译学,2014,硕士.
[20]钟婉茹.中国内地与香港食品安全法律制度比较研究[D].兰州大学,经济法学（专业学位）,2013,硕士.
[21]张哲宁.基于有限元理论的枢纽内部行人仿真研究[D].北京建筑大学,道路与铁道工程,2013,硕士.
[22]高磊.基于光流的动态场景中运动车辆检测与跟踪算法研究[D].中国科学技术大学,机械电子工程,2014,硕士.
[23]肖志远,郭凡良.社会主义核心价值观的培育路径——基于“场域-惯习”理论的视角分析[J].中共南昌市委党校学报,2014,06:42-45+58.
[24]冯科.辽圣宗朝统治阶层中的契汉矛盾[D].内蒙古大学,中国古代史,2014,硕士.
[25]马文博.安康市中心城区中小学教师媒介素养研究[D].西北大学,新闻学,2013,硕士.
[26]李扬.基于任务驱动法的导学案设计与实施[D].首都师范大学,现代教育技术,2014,硕士.
[27]王素青.《儿子与情人》中女性的觉醒[D].武汉理工大学,英语语言文学,2012,硕士.
[28]徐冰.基于下肢姿态传感器检测的人体防摔策略研究[D].吉林大学,生物医学工程,2013,硕士.
[29]焦理解.世界末日的救赎之火[D].中国矿业大学,英语语言文学,2014,硕士.
[30]张明媛.三种化学药物中有关物质的分析方法研究[D].第二军医大学,药物分析学,2013,硕士.
[31]孙洋.论浮动抵押中消极担保条款的效力[D].吉林大学,法律,2014,硕士.
[32]陈学荣.基于OpenGL的树木三维形态模拟[D].西安电子科技大学,机械制造及其自动化,2012,硕士.
[33]亓立敏.基于声强法的高速列车车内噪声源识别技术研究[D].西南交通大学,车辆工程,2014,硕士.
[34]李卓莉.城市基层管理体制改革背景下的新型城市社区区域化党建工作研究[D].天津大学,公共管理,2013,硕士.
[35]王丽.九节龙皂苷Ⅰ衍生物的合成及活性筛选[D].湖南大学,药物化学,2012,硕士.
[36]杨占宇.速度轮滑作为短道速度滑冰辅助训练的可行性研究[D].东北师范大学,体育,2012,硕士.
[37]左伟.基于人居环境科学视角下的徽州新聚落规划设计研究[D].安徽建筑大学,城市规划与设计,2013,硕士.
[38]江涛.刑事执行若干问题探讨[D].湘潭大学,法律,2003,硕士.
[39]刘聪.非债务税盾对上市公司融资行为影响的实证研究[D].陕西科技大学,会计学,2013,硕士.
[40]许益铵.舟山附近海域大黄鱼遗传多样性及家系的生长性状研究[D].浙江海洋学院,养殖（专业学位）,2014,硕士.
[41]李金祥.企业一体化管理体系研究[D].武汉理工大学,产业经济学,2003,硕士.
[42]赵倩.论民间借贷存在的问题及解决思路[D].烟台大学,法律,2014,硕士.
[43]殷惠君.膜结构风荷载的数值模拟研究[D].同济大学,2006.
[44]孙煜.小型水电站计算机监控系统设计[D].内蒙古大学,电子与通信工程,2013,硕士.
[45]唐哲,丁二玉,骆斌,陈世福.基于WebMining的推荐系统[J].计算机科学,2005,12:193-196.
[46]金丹.爱尔兰语与爱尔兰国家认同建构[D].北京外国语大学,英语语言文学,2014,硕士.
[47]张鹏飞.网络著作权合理使用制度中的一般条款研究[D].宁波大学,民商法学,2013,硕士.
[48]马飞.山西能源竞争力研究[D].天津财经大学,数量经济学,2012,硕士.
[49]彭可.中国农批公司发展战略研究[D].首都经济贸易大学,工商管理（专业学位）,2013,硕士.
[50]彭谊.跌创散联合透皮剂外敷促进闭合性骨折模型愈合的实验研究[D].辽宁医学院,中西医结合临床,2012,硕士.