论文部分内容阅读
木材产量和品质改良是林木育种研究的关键科学问题,木材形成是一个多基因控制的复杂过程,深入解析木材形成的关键因子和调控机理,有助于借助分子设计手段加快育种进程。日本落叶松(Larix kaempferi)是我国东北、华北和西南亚高山地区的重要针叶造林树种之一,也是我国主要针叶纸浆材用树种。本研究以17年生日本落叶松种源/家系试验林为研究对象,选取具有代表性的603个基因型个体建立关联群体,在分析表型多样性、遗传多样性、核苷酸多样性及连锁不平衡水平基础上,开展100个SNP标记和20对SSR标记与21个生长、形质、木材纤维及木材化学等表型性状的关联研究,旨在筛选出与木材生长和品质相关联的功能标记,为日本落叶松分子设计育种和标记辅助育种奠定基础。主要研究结论如下:1、日本落叶松关联群体表型多样性分析显示:日本落叶松具有丰富的表型变异,表型变异系数为3.08%-54.84%,不同性状的表型变异大小为:生长性状>纤维性状>形质性状>木材化学性状。关联群体各性状间及与地理环境因子间的相关性分析共检测到134对显著或极显著正负相关关系。主成分分析有效验证了性状间的相关性,并将21个表型性状降维成10个主成分。关联群体中6个自然群体的表型方差分析显示,除冠幅和早晚材微纤丝角的群体间差异不显著外,其他性状均存在显著或极显著差异;表型变异系数及表型分化系数分析表明日本落叶松自然群体具有丰富的表型变异;表型性状聚类分析将6个自然群体划分为两个类群。关联群体和人工群体表型性状的比较分析显示,大部分表型性状平均值均较自然群体显著提高,但部分木材纤维性状(如晚材纤维长和纤维长宽比等)平均值反而降低。2、生长及材性相关基因的SNP分析表明日本落叶松具有较高的核苷酸多样性水平。在扩增的27301 bp序列中共发现1698个SNP位点,SNP发生频率为1/16;总核苷酸多样性水平θw和πT变化范围分别为0.0035-0.2167和0.00224-0.03045。中性检验结果表明大部分候选基因在进化过程中主要受纯化选择作用。进一步的LD分析显示对日本落叶松开展基于候选基因的关联研究是可行的。20个候选基因中仅有Lac2和AE2基因在1500 bp时,R2值才降至0.1以下,其余18个基因的LD在R2<0.1,核苷酸序列长度约为250 bp-1000bp时就已消失。定量PCR分析表明20个生长及材性相关基因在日本落叶松根、成熟木质部、未成熟木质部、形成层、韧皮部和叶片等6个组织均有表达,其中生长相关基因在成熟叶片表达量最高,材性相关基因在成熟木质部和未成熟木质部高表达。3、利用SSR和SNP标记对日本落叶松自然群体进行遗传多样性分析,鉴于日本落叶松本身的生物学特性(物候差异、花粉流和种子流)、分布区域的复杂地形变化及生境片段化影响,日本落叶松在针叶树种中具有中等水平的遗传多样性(He=0.626);人工群体遗传多样性分析显示经过高强度的生长性状人工选择后,子代依然保存着较高的遗传多样性(He=0.613)。分子方差和遗传分化系数分析表明日本落叶松遗传变异主要以群体内变异为主(遗传距离:0.05-0.093)。基于SSR、SNP和SSR+SNP的自然群体遗传结构分析可将6个自然群体划分为2个类群。以自然群体结构划分结果为基础,采用SSR+SNP的标记组合可将关联群体划分为2个类群。类群Ⅰ主要由草津和浅间等自然群体构成,类群Ⅱ主要由富士和日光等自然群体及部分人工群体构成。4、以关联群体的遗传结构矩阵(Q矩阵)和亲缘关系矩阵(K矩阵)作为协变量,采用MLM模型进行关联研究。20个SSR位点中只有BC241与3年生胸径和材积显著关联。在20个生长及材性相关基因的100个候选SNP中,共检测到25个SNP位点与19个表型性状显著关联,对表型性状的解释率为1.09%-4.75%。生长相关基因中共有10个SNP位点与13个表型性状显著关联,其中AE基因家族中的AE2和AE3基因与部分形质性状和木材纤维性状显著关联。10个木材性状相关基因中,共检测到12个SNP位点与12个表型性状显著关联,对表型性状的解释率为1.09%-4.75%。转录及调控因子中,mi R397a有1个位点与早材微纤丝角和4个晚材性状(晚材壁腔比、晚材纤维长、晚材纤维宽和晚材率)显著关联;MYB转录因子中有两个位点(MYB1-337和MYB1-392)与密度显著关联。20个候选基因中,AE2、AE3和POD等3个基因与生长、形质、早晚材纤维性状和木材化学性状均显著关联,可作为木材性状改良的重点关注对象;MIR397a、CESA2和CSL基因主要与晚材发育显著相关,可作为晚材性状改良的重要标记;MYB转录因子主要与密度显著关联,是进行密度改良的重要标记。