利用连锁与连锁不平衡联合作图解析毛白杨重要性状的等位遗传变异

【摘要】森林作为陆地生态系统最重要的组成部分,蕴含了约82%的生物量和50%的生物多样性,是人类可持续发展与可再生利用的主要资源。但我国森林覆盖率低,国内所需木材原料大多依赖进口,木材资源的战略安全及国家林木育种工程面临严峻挑战；而“林以种为本,种以质为先”,常规的杂交育种策略虽为国家的速生丰产林建设提供了大批良种,但林木遗传杂合度高、育种周期长且重要经济性状多为数量性状等因素严重阻碍了林木遗传改良进程,当前林业生产所用品种的生产力难以满足工业发展对木材产量与质量的需求。因此,深入研究木材形成的生物学基础,解析木材形成的遗传调控机制,是加速林木良种化进程的根本。为此,本论文以我国北方重要的用材及生态型树种毛白杨(PopulustomentosaCarr.)的自然群体与杂交群体为材料,在毛白杨木材形成相关候选基因内微卫星标记资源大规模开发的基础上,首次开展了以连锁与连锁不平衡联合作图为策略的研究工作,剖析了控制林木生长及木材品质性状的等位遗传变异,并建立了林木基因标记辅助育种技术。主要研究内容及结论如下：1.利用直接测序法从138个毛白杨木材形成相关候选基因内开发出544个SSR标记位点,其中75%的位点位于基因的启动子或内含子等调控区域。利用这些标记在毛白杨自然群体中筛选出包含861个等位位点的188对多态性SSR引物,每一标记平均等位位点数目为2-7个。随机选择30对基因内SSR引物对杨树属内五个杨派26个种间个体(馒头柳作为类外群)进行种间转移检测发现,72%的基因内SSR标记可以在胡杨派中有效扩增,而在白杨派种间扩增率为100%；Neighbor-Joining(NJ)聚类分析确定了符合杨属传统分类的6个谱系分组。另外,通过比较220对包含基因内SSR位点的毛白杨与毛果杨同源序列,确定了种间SSR变异的主要原因是SSR位点侧翼序列的突变及重复基元类型的改变。进一步的杂交群体遗传分离检测及单标记关联分析表明开发的基因内SSR标记资源在关联作图、遗传连锁作图及近缘种比较基因组学等方面均具有较大的应用潜力。2.群体遗传多样性及群体结构评价是开展连锁不平衡作图及构建核心种质资源的前提。本文利用20对毛白杨特异性SSR标记对选择于毛白杨种质资源基因库的460株个体进行了遗传多样性评价,共确定了99个等位位点,其平均观测等位位点数目为4.05,显示了该群体具有中等遗传多样性水平。群体遗传结构分析表明,毛白杨种质资源可以划分为11个亚群体(K=11)。对STRUCTURE模型确定的K=3情况下的亚群体组成进行分析,表明了毛白杨三个气候区可以分别作为相对独立的遗传群体进行相关研究。对不同气候区的遗传多样性检测表明南部气候区具有最丰富的遗传多样性水平(N=93,Np=11,HE=0.445,F=-0.102),东北及西北气候区次之,表明南部气候区是当前毛白杨生长分布的中心区域。利用建立的“UPGMA-聚类分析逐次随机取样”与“群体结构重检测”联合分析法对毛白杨种质资源进行核心种质群体构建,确定15%(69)的比例选择出来的毛白杨核心种质群体是最理想的。3.进一步的连锁不平衡研究显示,在毛白杨候选基因内1,200bp距离内连锁不平衡水平已衰退至不明显,显示了基于候选基因的连锁不平衡作图策略在林木中具有更高的效力。利用混合线性模型进行单标记关联分析,确定了来自于39个候选基因的48个SSR标记位点与10个生长与木材品质性状显著连锁,构成80个关联(Q≥0.10),每一标记位点可解释表型变异的1.20%-14.5%；基于单倍型的关联分析确定了来自于32个高连锁区域(代表20个基因)的64个单倍型与10个性状的80个关联(Q≤0.10)。比较单标记及单倍型的关联结果显示单倍型效应主要来自于其包含的显著单标记位点,基于单倍型的关联分析并没有显现出更高的作图分辨率。4.利用1200株银腺杨×毛白杨杂交F1连锁作图群体构建了首张基于候选基因内共显性标记的毛白杨高密度、高解析度遗传连锁图谱。其中,在LOD≥5.0,r≤0.25条件下,确定了包含1270个标记位点(309个SSR、929个AFLP和32个InDel)的19个较大的连锁群,连锁位点覆盖毛白杨基因组总长约2758.6cM,标记间平均遗传距离为2.3cM。341个基因内共显性标记在连锁图谱的标定代表了候选基因的定位,为后期开展重要性状的QTL作图提供了科学的遗传平台。5.基于遗传图谱信息,利用复合区间作图法对光合相关性状、生物量、叶片性状和木材品质等28个数量性状进行了QTL分析。当LOD值≥2.5时,共检测到与25个表型性状相关的179个QTLs,每一标记可解释表型变异的1.27%-30.02%,显示不同复杂数量性状具有不同的遗传控制机制。所有179个QTL按照光合相关性状、叶片性状、生物量及木材品质分组后确定的QTL数目是10-84个,不均匀地分布到18个连锁群上,平均9.4个。LGVI和LGXII是控制这些性状的重要区域,LG_XIII上没有发现与所有28个数量性状相关的QTL。进一步研究发现80%的QTLs区间内包含基因内共显性标记,其隶属的基因主要包括与细胞生长分化相关的重要调控基因或转录因子,如类核酸剪切调控蛋白基因、蛋白酶体调控亚基S3蛋白家族基因和锌指蛋白家族基因等,另外一部分是控制木材形成及生长的功能基因,如纤维素合酶基因家族成员(CesA3、CesA7和CesA9),葡萄糖-6-磷酸脱氢酶基因与EXP1基因等。6.利用连锁与连锁不平衡联合作图策略,确定了来自于27个候选基因的29个SSR与8个生长与木材品质性状组成的39个显著关联(Q≤0.10)。进一步比较了每个显著标记在两个作图群体内对应的基因型效应,最终选择基因型效应一致的31个关联组合。随后,筛选出控制每个性状极高及极低表型值的优势等位位点联合体,并根据育种需要组合利用多个优势位点联合体建立基因标记辅助选择育种技术。利用该技术分别在毛白杨种质资源库(关联群体)及杂交作图群体内选择了10个具有潜在应用价值的优良单株。
【作者】杜庆章；
【导师】张德强；
【作者基本信息】北京林业大学，林木遗传育种，2014，博士
【关键词】毛白杨；基因内SSR标记；遗传连锁图谱；连锁与连锁不平衡联合作图；基因标记辅助育种；

【参考文献】
[1]邱胜根.数字图像处理技术在投影显示三色会聚中的应用研究[D].浙江大学,光学工程,2004,硕士.
[2]王维.黑果枸杞内生真菌次生代谢产物及其抑菌活性的研究[D].陕西科技大学,有机化学,2013,硕士.
[3]邱宗波,刘晓,李方民,田向军,岳明.He-Ne激光处理与干旱胁迫对小麦幼苗保护酶活性及脂质过氧化作用的影响[J].植物生态学报,2008,05:1002-1006.
[4]王珊珊.碳量子点的制备及性能研究[D].北京化工大学,2013.
[5]华正韬.土壤和石油烃性质对溶剂萃取法修复石油污染土壤的影响[D].天津大学,化学工程,2013,硕士.
[6]汪永红.浙江森林资源资产抵押贷款融资创新研究[D].浙江工业大学,2007.
[7]赵珊珊.基于现代物流的铁路零担货物运营网络优化设计[D].北京交通大学,2015.
[8]宋扬.不同心理弹性初中生对情绪面孔的注意偏向及情绪认知特点[D].吉林大学,基础心理学,2013,硕士.
[9]王璇.中国传统礼仪教育及其现代价值[D].兰州商学院,思想政治教育,2014,硕士.
[10]郭庆龙.故意毁坏财物罪中的毁坏行为研究[D].西南政法大学,刑法学,2012,硕士.
[11]孟双乐.钢—混凝土组合肋壳结构行波效应分析[D].西安建筑科技大学,防灾减灾工程及防护工程,2013,硕士.
[12]张煜.玻璃在时装设计中应用的探索[D].中央美术学院,艺术设计,2013,硕士.
[13]许鑫.司法回应民意的正当性研究[D].天津商业大学,法学理论,2013,硕士.
[14]徐苏男.水分胁迫及复水对结缕草生长生理和光合荧光特性的影响[D].辽宁大学,生态学,2012,硕士.
[15]王传禛.农产品资本化对其期货价格的影响[D].南京农业大学,金融学,2012,硕士.
[16]李科伟.四位优秀女子网球选手发球与接发球阶段线路变化特征的研究与分析[D].北京体育大学,体育教育训练学,2013,硕士.
[17]吴仲其,马庆莲,袁文林.船体分段吊装工艺[J].江苏船舶.1990(02)
[18]赵蓓蓓.树脂基固态胺吸附剂对低浓度二氧化碳的室温吸附行为研究[D].华东理工大学,材料化学工程,2014,硕士.
[19]尹盼盼.社会福利供给的绩效评价研究[D].上海工程技术大学,社会保障,2014,硕士.
[20]王绍臣.21世纪亚太经济社会和文化科技的发展[J].未来与发展.2001(06)
[21]沈夏燕.服务业前后台组织间知识转移机制研究[D].浙江工商大学,技术经济及管理,2013,硕士.
[22]王科俊,马慧,李雪峰.指纹及指静脉双模态识别的二级分类融合法研究[J].控制与决策,2011,08:1131-1135.
[23]王昱潭.基于机器视觉的灵武长枣定位与成熟度判别方法研究[D].北京林业大学,森林工程,2014,博士.
[24]徐洁.受限态高分子结构，动力学行为及其应用[D].南京大学,高分子化学与物理,2014,博士.
[25]张春燕.外商在华R&D投资的前向关联路径研究[D].华中科技大学,行政管理,2013,硕士.
[26]陈秋霞,汪卫国,秦军.中国IPTV的发展现状与前景分析[J].现代电信科技.2007(08)
[27]焦宁.吉林省棚膜蔬菜质量安全可追溯体系研究[D].吉林大学,农业经济管理,2014,硕士.
[28]梅帅.开放式结构拓扑优化软件设计与开发[D].大连理工大学,计算力学,2013,硕士.
[29]顾菁.菊属野生种抗旱生理机理及抗旱蛋白质组学研究[D].南京农业大学,园林植物与观赏园艺,2013,硕士.
[30]陈恩红,张振亚,王煦法.基于神经网络的增量式数据索引机制研究[A].中国自动化学会智能自动化专业委员会、昆明船舶设备集团公司、云南自动化学会、中国人工智能学会计算机视觉及智能控制学会、IEEE控制系统学会北京分会、中国自动化学会机器人竞赛工作委员会.2001年中国智能自动化会议论文集（上册）[C].中国自动化学会智能自动化专业委员会、昆明船舶设备集团公司、云南自动化学会、中国人工智能学会计算机视觉及智能控制学会、IEEE控制系统学会北京分会、中国自动化学会机器人竞赛工作委员会:,2001:6.
[31]崔张坤.电动汽车锂电池组高精度电压电流检测系统研究与设计[D].沈阳理工大学,检测技术与自动化装置,2012,硕士.
[32]王婷婷.现阶段中国农村合作社治理体系研究[D].内蒙古大学,公共管理,2014,硕士.
[33]曹秉基.肾上腺素α1受体在人食管下括约肌的表达研究[D].河北医科大学,外科学,2014,硕士.
[34]李静.一体式中空纤维帘式膜生物反应器性能的研究[D].天津工业大学,材料学,2004,硕士.
[35]高伟.基于特征知识库的遥感信息提取技术研究[D].中国地质大学,2010.
[36]葛阳.欣胃颗粒对慢性萎缩性胃炎大鼠血清胃泌素的影响[D].黑龙江中医药大学,中医内科学,2013,硕士.
[37]李志寒.从ErbB2阳性乳腺癌细胞代谢变化角度初探拉帕替尼诱导细胞耐受的分子机制[D].厦门大学,药理学,2014,硕士.
[38]谢雪宇.基于空气负离子资源的森林公园保健养生项目规划研究[D].中南林业科技大学,园林植物与观赏园艺,2014,硕士.
[39]王阳.供电企业信息资源整合平台研究与应用[D].电子科技大学,软件工程（专业学位）,2012,硕士.
[40]徐红如.对广义自缩序列的攻击算法研究[D].西安电子科技大学,密码学,2004,硕士.
[41]张生昌,张志鸿,马艺,邓鸿英,马中强.吸入压力对往复式油气混输泵外特性的影响[J].浙江工业大学学报,2015,01:34-38+120.
[42]武敏.分蘖葱头总黄酮纳米乳载药系统经皮靶向吸收研究[D].吉林大学,药物化学,2013,硕士.
[43]郁萌.复杂三角网格模型数控加工特征提取及刀具轨迹生成研究[D].华侨大学,材料加工工程,2013,硕士.
[44]刘峰雷.订单式生产企业预算管理研究[D].山东大学,工商管理（专业学位）,2012,硕士.
[45]陈黎.序列挖掘算法研究[D].重庆大学,2001.
[46]崔希利.介孔材料的模板合成及其应用研究[D].天津大学,化学工艺,2013,硕士.
[47]邓旭升.搅拌床快速热解反应器设计及其实验研究[D].北京林业大学,木材科学与技术,2013,硕士.
[48]陈秋兰.应用PAR指数评价骨性安氏Ⅱ类错（牙合）畸形矫治效果的临床研究[D].广西医科大学,口腔正畸学（专业学位）,2013,硕士.
[49]周超俊,蒋慰孙,叶银忠,藤井省三.一种基于贝叶斯分类原理的动态系统故障检测与诊断方法[J].控制与决策,1995,03:270-274+288.
[50]林淑萍.漆器行业发展的研究[D].厦门大学,工商管理,2014,硕士.