极小种群大别山五针松种群生态学及遗传多样性研究

大别山五针松(Pinus dabeshanensis C.Y.Cheng&Y.W.Law)隶属于松科(Pinaceae)松属(Pinus),是我国特有的珍稀濒危植物,仅局限分布在大别山地区.本文利用种群生态学和遗传学的方法,从群落物种多样性、区系特点、种群数量动态、种内和种间竞争、种群生殖力调查、转录组测序及微卫星引物开发、种群遗传多样性等方面对其进行了研究,以期为该濒危植物的保护和利用提供依据.主要研究结果如下:1 大别山五针松群落物种多性及区系特点大别山五针松分布区的群落内,共有维管植物72科141属218种.其中,蕨类植物4科4属6种,裸子植物3科3属4种,被子植物65科134属208种.蔷薇科、樟科、菊科、杜鹃花科、豆科和禾本科等科为优势科.区系分析显示,该群落区系起源古老,地理成分复杂,68个科中,世界分布、热带分布和温带分布的科分别有21、32和15个,显示由热带向温带过渡的特性.137个属中,热带成分共40属,占总属数(不包括15个世界广布属)的32.79%;温带成分共有82属,占总属数的67.21%,并以北温带为主.群落垂直结构分为乔木层、灌木层和草本层3层,其中,Margalef指数、Simpson指数和Shannon-Wiener指数均以灌木层最大,其次是乔木层,草本层最低.灌木层与草本层的Pielou指数相近,且均大于乔木层.2 大别山五针松种群数量动态用空间序列代替时间的方法分析种群结构,编制种群特定时间生命表,绘制死亡率曲线和消失率曲线,并用4个生存函数进行种群的生存分析;同时结合谱分析方法,分析了大别山五针松种群数量的动态变化.结果表明,大别山五针松种群数量少,结构存在波动性.幼苗阶段个体较丰富,幼树阶段个体较少,种群趋于衰退.生命表和存活曲线表明,大别山五针松存活曲线介于Deevey Ⅱ和Deevey Ⅲ型之间,种群死亡率和消失率曲线变化趋势基本一致,在第5龄级出现峰值.4个生存函数曲线表明,大别山五针松具有前期快速减少、中期稳定和后期衰退的特点.谱分析显示,大别山五针松种群动态除受基波影响外,还具有明显的小周期波动.3 大别山五针松种内、种间竞争关系用Hegyi单木竞争模型分析了大别山五针松的种内和种间竞争强度.结果表明,竞争主要来自种间.使大别山五针松受压的植物主要有短柄袍、黄山松、满山红和茅栗等.大别山五针松在胸径小于25 cm之前,受到的竞争压力较大;当胸径大于25 cm后,竞争强度变化不明显,竞争强度与胸径符合幂函数关系(CI=AD-B).4 大别山五针松种群生殖力大小生殖力调查表明,自然种群种子的败育率为61.08%,饱满种子的活力为100%;栽培种群种子的败育率高达89.88%,饱满种子的活力为89.76%.低温冷藏更有利于种子活力的保存.种子可在较宽的温度范围内萌发,且室温储藏的种子萌发率比低温冷藏的要高;破壳处理有利于种子吸胀但对种子萌发没有明显的促进作用;自然种群种子萌发率显著高于栽培种群.5 大别山五针松的转录组测序及微卫星引物开发通过Illumina HiSeqTM 2000测序平台,对大别山五针松针叶进行转录组测序.从所获得的unigene中共检测到SSR 1966个,平均每23.08 kb含一个SSR位点.SSR序列中最丰富的重复基元为三核苷酸,其次是二核苷酸、单核苷酸、六核苷酸、五核苷酸和四核苷酸,所占的比例分别为41.60%、23.25%、21.41%、7.83%、4.43%和1.53%.三核苷酸和二核苷酸中最丰富的基元类型分别是AGC/CTG和AT/AT.基于2-6核苷酸重复序列,共成功设计出引物431对.随机挑选82对扩增,共获得19对多态性引物,这些引物在24个个体中共检测到57个等位基因,平均每个位点的等位基因数NA=3.0个.引物的预期杂合度和观测杂合度范围分别为0-0.958和0.082-0.759,PIC值变化范围为0.077-0.707,平均值为0.389.6 大别山五针松的遗传多样性和遗传结构13对EST-SSR引物在233个个体中共检测到50个等位基因数,平均每个位点3.85个,自然种群和栽培种群的预期杂合度HE分别为0.458和0.423,观察杂合度HO分别为0.481和0.510,表明大别山五针松种群的遗传多样性属于中等水平.当代基因流和历史基因流均显示种群间存在较大基因流,种群间遗传分化极小(FST=0.028).MSVAR分析表明,有效种群从古代的大约104-105株下降至现在的102-103株,种群大小发生变化的时间大约在10000年前.基于BOTTLENECK软件和TPM及SMM两种模型检测结果显示,在TPM模型下,2个自然种群可能经历了近期的瓶颈效应.STRUCTURE聚类分析显示,自然种群和栽培种群被分为两组.其中,栽培种群主要含1个基因池,而自然种群包含2个.3个自然种群均含有稀有等位基因,而栽培种群没有,表示栽培种群不足以代表自然种群所包含的遗传信息.更新世气候的回升可能是导致种群大小发生变化的主要原因.7大别山五针松的濒危机制与保护策略鉴于上述研究,大别山五针松的濒危主要由以下几个方面导致:1)种群结实率低、种子败育率高,导致种群生殖力较低;2)幼苗死亡率高,难以向幼树转化;3)种群竞争力弱,在竞争中处于劣势;4)小种群效应导致近交衰退现象严重,难以抵御各类风险.针对以上研究结果,本文提出如下保护策略:1)就地保护,适当间伐杂木,为幼苗、幼树的生长提供光、营养等资源,促使幼苗形成幼树;同时减少生存竞争压力;2)人工控制授粉,提高种子品质;3)保护大王沟的3个种群,尤其注意保护含稀有等位基因的个体;4)从自然种群采集种子,并多采集大王沟1和大王沟3种群及含有稀有等位基因个体的种子,人工育苗扩大种群.