-
丰富的种质资源为植物育种、遗传研究和新品种选育提供了条件和物质基础[1-2]。为更好地应对当前气候条件变化,急需培育出适宜在不同条件下生长的新品种,同时合理保护和利用植物种质资源是一项复杂的任务。在原始种质资源中挑选少部分有代表性的种质构建核心种质库可以最大限度地保存其表型或遗传多样性,从而有效地解决种质资源保存中成本高等问题,缩短育种进程,促进种质资源的合理开发和利用[3-4]。
与作物种质资源核心种质库构建相关的研究已有众多报道[5-8],其构建方法、策略和评价标准均相对成熟。参考作物核心种质库的构建策略,多种木本植物核心种质库已构建成功,包括葡萄(Vitis vinifera L.)[9]、油橄榄(Olea europaea L.)[10]和角豆(Ceratonia siliqua L.)[11]等经济树种,欧洲黑杨(Populus nigra L.)[12]、木荷(Schima superba Gardn. et Champ.)[13]、杉木(Cunninghamia lanceolata (Lamb.) Hook.)[14]和马尾松(Pinus massoniana L.)[15]等用材林树种。当前,构建核心种质库所选用的数据因研究目的不同,主要有表型数据或基因型数据[16]。植株的表型数据测量方便,在核心种质库的构建研究中应用广泛,同时随着分子标记技术和基因组测序技术在核心种质库构建研究中的应用,相关学者越发认识到植物表型研究的重要性和必要性。育种研究最终目的是获得具有特定性状和表型的新品种[7, 17],因此,基于表型数据构建核心种质库仍是未来种质资源研究的重要内容。由于不同物种具有不同的表型特征和群体结构,在构建核心种质库的研究中采用的取样策略不同,均具有一定的优缺点,如随机取样策略可以更好的保存原始种质的群体结构,但构建的核心种质库的代表性和多样性相对较差;偏离度取样可以保留较大的变异和多样性;位点优先取样策略不仅可以保留具有性状极值的种质,同时可以较好的保留群体结构[18-19]。因此,研究者应根据所研究种质资源的特点选择合适的取样策略。
美洲黑杨(Populus deltoides Marsh.)种质资源在杨树育种研究工作中占有重要地位[20-21],合理保护其表型和遗传多样性是高效开发和利用美洲黑杨种质资源的前提,本课题组通过前期对其表型和遗传多样性研究,发现原始种质资源库具有丰富的多样性[22]。为了更好地保存美洲黑杨种质资源,有研究者利用分子标记技术构建其核心种质库[23-24],但基于表型数据构建美洲黑杨核心种质库的研究仍未见报道。本研究基于美洲黑杨种质资源表型和生理性状,分析其群体结构,对种质资源进行分组,利用分组逐步聚类取样的方法构建表型核心种质库;通过与3种常规取样策略进行比较,首次提出性状频率取样策略并用于构建美洲黑杨表型核心种质库;在表型核心种质库的基础上通过补充部分种质构建表型优化核心种质库,使其表型保留比例达到100.00%,最大化保存种质资源的表型多样性信息。本研究有助于研究者更好的了解、保护和利用美洲黑杨种质资源,挖掘资源中的核心种质,为杨树育种研究及其他物种表型核心种质库构建奠定基础和提供参考。
HTML
-
本研究中的试验材料包括美洲黑杨6个种源群体、27个采样点的258个无性系(表1)。种质来源地包括哥伦比亚河流域(美国华盛顿州,Was)、北美洲圣劳伦斯河流域(加拿大魁北克省,Que)和密西西比河流域(美国密苏里州,Mis;艾奥瓦州,Iow;田纳西州,Ten;路易斯安那州,Lou)。试验材料与前期研究中的材料来源信息一致[22]。
种源群体
Provenances纬度
Latitude (N)经度
Longitude (W)采样点代码
Sample code无性系数量
Number of clones密西西比河流域
Mississippi River basin美国艾奥瓦州
Iowa, America (Iow)41°52′48″ 93°06′00″ I 19 美国密苏里州
Missouri, America (Mis)37°18′11″~38°49′08″ 89°32′40″~90°07′53″ M1 3 M2 8 小计 Subtotal 11 美国路易斯安那州
Louisiana, America (Lou)30°30′00″~32°03′00″ 91°27′00″~92°34′48″ L1 3 L3 2 L4 4 L5 8 L6 4 L7 7 L8 5 L9 4 小计 Subtotal 37 美国田纳西州
Tennessee, America (Ten)35°51′36″~36°24′36″ 89°12′00″~89°36′00″ T1 10 T2 10 T3 12 T4 6 T5 13 小计 Subtotal 51 哥伦比亚河流域
Columbia River basin美国华盛顿州
Washington, America (Was)46°13′07″ 119°04′49″ W 15 圣劳伦斯河流域
Saint Lawrence River basin加拿大魁北克省
Quebec, Canada (Que)46°02′07″~46°34′31″ 71°51′47″~73°02′03″ Q1 30 Q2 4 Q3 2 Q4 3 Q5 48 Q6 8 Q7 23 Q8 2 Q9 3 Q10 2 小计 Subtotal 125 总计 Total 258 Table 1. Source information of P. deltoides germplasm resources
-
将美洲黑杨植株1年生枝条剪成插穗,扦插于营养钵中(口径18 cm,高度25 cm),每个无性系扦插10株,在中国林业科学研究院通州试验基地(39°73′35″ N,116°75′18″ E)温室大棚进行培养,培养基质为黄土:粗砂石:草炭土 = 6:1:1。
在植株生长季,同时利用3台Li-6400便携式光合测定仪(Li-6400,LI-COR,美国)测定植株叶片的净光合速率(Pn)、胞间CO2浓度(Ci)、气孔导度(Gs)和蒸腾速率(Tr),测定叶片为上部第5~7片成熟叶片,每个无性系测定3株。测量过程中充分考虑仪器、测量时间和环境等因素对测量指标的影响,将258个无性系分为9组,每组28~29个无性系,测量时间分为3个时间段(9:00—9:40、9:40—10:20和10:20—11:00),每台仪器在1个时间段内可以完成1组无性系1株植株的1个叶片光合参数的测量,最终使每个无性系测定的3株的9片叶片中每台仪器在每个时间段各测定1片,1个生物学重复的测量周期为3 d,共测量9 d。使用叶绿素相对含量测定仪(SPAD-502,日本)测定植株上部第5~7片成熟叶片的SPAD值,每个叶片选取6个测定点,每个无性系测定3株。
每个无性系选取3株生长性状相似的植株,各采集3~5片成熟功能叶,称鲜质量,使用叶面积仪(Yaxin-1241,北京)测量叶片形态性状(长度、宽度和叶面积),然后使其风干称风干质量,最后烘干并称烘干质量。利用元素分析仪(Thermo Fisher Scientific,Inc.,美国),测定烘干叶片中的碳、氮元素的含量。
在植株生长后期,测量植株的株高、地径和叶片数(叶痕数),同时对根系和茎段的生长量(干质量和鲜质量)进行测定。
-
利用叶片的鲜质量、风干质量和烘干质量计算叶片自由水和结合水的含量,即自由水含量 = 鲜质量−风干质量;结合水含量 = 风干质量−烘干质量。
表型多样性指数(Shannon-Wiener指数,H')的计算公式为:
式中:H'为Shannon-Wiener指数;ln表示自然对数;n为性状指标总数;Pi表示某性状指标第i种级别的无性系数占总无性系数的比率,各性状指标被划分为10个级别,数值小于
$ {\overline{x}}-2 $ σ为第1级别,数值大于等于$ {\overline{x}}+2$ σ为第10级别,其中,每0.5σ划分为一个级别。无性系性状频率的计算公式为:
式中:Fj为第j个无性系的性状频率;n为性状指标总数;Pij表示种质库中第j个无性系的第i个性状指标所在的性状级别的频率。
-
利用表型和生理指标数据对6个种源、258个美洲黑杨无性系进行分组,Was种源的无性系单独为一组;Que和Iow种源无性系为一组;Mis、Lou和Ten种源无性系为一组。基于无性系标准化的表型和生理指标数据计算每组内无性系间的欧式距离,采用可变类平均法对无性系进行聚类分析,在最低聚类水平上分别采用随机取样、偏离度取样和位点优先取样的策略构建表型核心种质库。同时,本研究首次提出性状频率取样策略,即在聚类的最低级别的两个或多个无性系中优先选择具有性状极值的无性系,如果均具有性状极值均进入下一轮聚类分析,如果均没有性状极值则选择性状频率低的无性系进入下一轮聚类分析,以此类推,直到构建的种质库不满足核心种质库构建标准。
-
参照前期研究中表型核心种质库的检测方法,通过计算构建的种质库相对于原始种质库的极差符合率(CR)、变异系数变化率(VR)、均值差异百分率(MD)、方差差异百分率(VD)和表型保留比例(RPR)对种质库的代表性进行评价,各参数的计算公式[5, 18]如下:
式中:RC(i)是核心种质第i个性状的极差,RI(i)是原始种质第i个性状的极差,n是性状总数。
式中:CVC(i)是核心种质第i个性状的变异系数,CVI(i)是原始种质第i个性状的变异系数,n是性状总数。
式中:St是核心种质与原始种质进行t测验得到的均值差异显著(α = 0.05)的性状数,n是性状总数。
式中:SF是核心种质与原始种质进行F测验得到的方差差异显著(α = 0.05)的性状数,n是性状总数。
式中:MC(i)是核心种质第i个性状的表现型个数,MI(i)是原始种质第i个性状的表现型个数,n是性状总数。
通过计算各个种质库各个性状指标的表型多样性指数(H')对种质库的多样性进行评价。基于主成分分析构建的前2个主成分因子的得分,绘制原始种质库和表型核心种质库的分布图,对表型核心种质库进行确认。
表型核心种质库构建的最低标准为:MD < 20.00%且CR > 80.00%,RPR经常作为评价表型核心种质库有效性和取样比例必不可少的参数,本研究要求核心种质库的RPR > 90.00%,同时认为MD越小,CR、VR和VD越大表型核心种质库越能代表原始种质库的多样性[18]。
-
在原始种质库中挑选出表型核心种质库后剩余的种质称为保留种质。基于表型保留比例最大化的思想(RPR = 100.00%),分析并统计表型核心种质库相对于原始种质库中各性状指标丢失的表型类型,然后在保留种质中挑选最少数量的具备相应表型的无性系补充到表型核心种质库中形成表型优化核心种质库,当多个无性系均具有相应表型类型时,选择性状频率低的无性系,该部分无性系称为表型补充种质,最后对表型优化核心种质库进行评价。
1.1. 试验材料
1.2. 植株培养及表型生理指标的测定
1.3. 数据处理
1.4. 表型核心种质库构建策略
1.5. 表型核心种质库代表性检验方法
1.6. 优化核心种质库构建策略
-
对美洲黑杨3个类群的表型和生理指标进行对比分析(表2),发现第1类群(Was种源)无性系的株高(H)、叶片数(LN)、叶片碳元素含量(C)、叶片氮元素含量(N)、Pn和Tr高于另外2个类群的无性系;第2类群(Iow和Que种源)无性系的根系鲜质量(RFW)、根系干质量(RDW)和叶绿素相对含量(SPAD)最大;第3类群(Mis、Lou和Ten种源)无性系的地径(GD)、茎段和叶片的生长量、叶片形态性状、叶片自由水和结合水含量、Gs和Ci均高于另外2个类群。3个类群的表型和生理特征均比较明显,在该分组的基础上分别在组内筛选代表性的无性系构建核心种质库可以更好的代表美洲黑杨群体的表型和生理特征。
性状
Traits第1类群
Group 1 (Was)第2类群
Group 2 (Iow、Que)第3类群
Group 3 (Mis、Lou、Ten)H/cm 76.93 ± 16.17 a 62.93 ± 13.94 b 65.78 ± 17.92 b LN 42.32 ± 7.40 a 31.10 ± 6.46 c 34.27 ± 5.97 b C/(mg·g−1) 453.40 ± 12.68 a 430.84 ± 16.28 b 446.95 ± 16.74 a N/(mg·g−1) 33.55 ± 6.66 a 28.61 ± 4.58 b 29.03 ± 4.66 b Pn/(μmol·m−2·s−1) 10.94 ± 1.98 a 10.03 ± 2.37 ab 10.83 ± 2.53 a Tr/(mmol·m−2·s−1) 5.33 ± 1.37 a 5.01 ± 1.46 a 5.26 ± 1.45 a RFW/g 11.49 ± 2.31 ab 12.70 ± 3.88 a 9.50 ± 4.57 b RDW/g 4.92 ± 1.04 a 5.56 ± 1.78 a 3.60 ± 1.80 b SPAD 35.58 ± 3.82 ab 38.56 ± 5.41 a 35.37 ± 4.67 b GD/mm 7.38 ± 0.68 ab 7.02 ± 1.07 b 7.84 ± 1.25 a SFW/g 10.84 ± 4.09 ab 8.64 ± 4.44 b 12.36 ± 6.21 a SDW/g 5.44 ± 1.95 ab 4.31 ± 2.19 b 5.68 ± 2.93 a LFW/g 0.79 ± 0.14 c 1.00 ± 0.19 b 1.27 ± 0.29 a LDW/g 0.23 ± 0.04 c 0.30 ± 0.06 b 0.36 ± 0.08 a LA/cm2 35.32 ± 6.03 c 45.19 ± 9.00 b 55.17 ± 12.18 a LL/cm 7.33 ± 0.64 c 8.47 ± 0.99 b 9.49 ± 1.11 a LW/cm 7.53 ± 0.61 b 8.05 ± 0.79 ab 8.30 ± 0.97 a FWC/g 0.55 ± 0.11 c 0.67 ± 0.13 b 0.87 ± 0.21 a BWC/mg 23.43 ± 6.58 b 29.46 ± 9.95 b 39.68 ± 12.64 a Gs/(mol·m−2·s−1) 0.25 ± 0.08 ab 0.24 ± 0.09 b 0.27 ± 0.10 a Ci/(μmol·mol−1) 287.06 ± 21.47 a 288.01 ± 31.26 a 292.48 ± 25.42 a 注:H:株高;LN:叶片数;C:叶片碳元素含量;N:叶片氮元素含量;Pn:净光合速率;Tr:蒸腾速率;RFW:根系鲜质量;RDW:根系干质量;SPAD:叶绿素相对含量;GD:地径;SFW:茎段鲜质量;SDW:茎段干质量;LFW:单叶鲜质量;LDW:单叶干质量;LA:单叶面积;LL:叶片长度;LW:叶片宽度;FWC:单叶自由水含量;BWC:单叶结合水含量;Gs:气孔导度;Ci:胞间CO2浓度;DW:干重。不同字母表示同一性状不同类群在0.05水平上差异显著,下同。
Notes: H: height; LN: number of leaves; C: Carbon content of leaves; N: Nitrogen content of leaves; Pn: net photosynthetic rate; Tr: transpiration rate; RFW: fresh weight of root; RDW: dry weight of root; SPAD: relative content of chlorophyll; GD: ground diameter; SFW: fresh weight of stem; SDW: dry weight of stem; LFW: fresh weight of single leaf; LDW: dry weight of single leaf; LA: single leaf area; LL: leaf length; LW: leaf width; FWC: free water content of single leaf; BWC: bound water content of single leaf; Gs: stomatal conductance; Ci: intercellular CO2 concentration. Different letters represent significant differences at 0.05 level, the same as below.Table 2. Phenotypic and physiological characteristics of different groups
-
在构建美洲黑杨表型核心种质库时,基于不同的取样策略得到的结果存在较大差异(表3)。采用随机取样策略时,构建的第6个种质库(R6)的RPR为87.57%,低于90.00%,初步认定第5个种质库(R5)为核心种质库,包括32个无性系,取样比例为12.40%,MD、VD、CR、VR和RPR分别为0.00%、9.52%、85.69%、116.31%和91.90%;采用偏离度取样策略时,共构建了5个种质库,第4个种质库(D4)被初步认定为核心种质库,取样比例为19.38%,MD、VD、CR、VR和RPR分别为9.52%、28.57%、93.17%、113.44%和94.71%;采用位点优先取样策略时,共构建了8个种质库,第8个种质库(P8)的RPR低于90.00%,初步确定第7个种质库(P7)为核心种质库,包括32个无性系,MD、VD、CR、VR和RPR分别为0.00%、76.19%、100.00%、141.66%和91.69%;采用性状频率取样策略时,共构建了10个种质库,第10个种质库(F10)的RPR为89.89%,第9个种质库(F9)被初步认定为核心种质库,取样比例为10.47%,MD、VD、CR、VR和RPR分别为0.00%、80.95%、100.00%、148.94%和90.85%。在随机取样、偏离度取样和位点优先取样策略下形成的种质库之间H'没有显著的变化,而在性状频率取样策略下,H'在不同种质库间存在明显差异,随着取样比例的减少,H'先增大后减小,但在4种不同的取样策略下构建的核心种质库的H'与原始种质库均没有显著差异,均能代表原始种质库表型和生理性状的多样性。
取样策略
Sampling strategies种质库
Banks样本数(比例)
Number of
samples(ratio)均值差异
百分率
MD/%方差差异
百分率
VD/%极差
符合率
CR/%变异系数
变化率
VR/%表型保留
比例
RPR/%表型多样
性指数
H'随机取样
Random samplingR1 165(63.95) 0.00 0.00 94.20 99.70 100.00 2.018 a R2 108(41.86) 0.00 0.00 91.58 103.58 99.58 2.020 a R3 71(27.52) 0.00 0.00 89.85 108.17 99.05 2.045 a R4 49(18.99) 0.00 0.00 88.31 112.18 98.10 2.040 a R5 32(12.40) 0.00 9.52 85.69 116.31 91.90 2.001 a R6 22(8.53) 0.00 9.52 85.06 126.84 87.57 1.980 a 偏离度取样
Deviation samplingD1 165(63.95) 0.00 0.00 97.19 100.62 99.52 2.027 a D2 110(42.64) 0.00 0.00 95.87 104.36 99.05 2.043 a D3 74(28.68) 4.76 4.76 94.92 108.55 97.57 2.059 a D4 50(19.38) 9.52 28.57 93.17 113.44 94.71 2.060 a D5 35(13.57) 19.05 52.38 89.40 115.25 89.37 2.033 a 位点优先取样
Allele preferred samplingP1 167(64.73) 0.00 0.00 100.00 103.16 100.00 2.031 a P2 114(44.19) 0.00 0.00 100.00 108.31 100.00 2.037 a P3 79(30.62) 0.00 9.52 100.00 114.44 99.52 2.058 a P4 62(24.03) 0.00 33.33 100.00 120.40 97.57 2.056 a P5 48(18.60) 0.00 61.90 100.00 126.22 97.09 2.069 a P6 36(13.95) 0.00 66.67 100.00 136.74 95.66 2.045 a P7 32(12.40) 0.00 76.19 100.00 141.66 91.69 2.020 a P8 28(10.85) 0.00 85.71 100.00 148.65 87.83 1.992 a 性状频率取样
Trait frequency samplingF1 167(64.73) 0.00 0.00 100.00 105.02 100.00 2.059 abcd F2 106(41.09) 0.00 19.05 100.00 111.93 100.00 2.100 abc F3 75(29.07) 0.00 52.38 100.00 118.74 99.52 2.125 ab F4 57(22.09) 0.00 66.67 100.00 124.97 99.52 2.136 a F5 44(17.05) 0.00 90.48 100.00 132.41 99.05 2.130 ab F6 34(13.18) 0.00 76.19 100.00 140.11 95.61 2.079 abcd F7 29(11.24) 0.00 71.43 100.00 146.03 92.75 2.044 bcd F8 28(10.85) 0.00 71.43 100.00 146.92 91.80 2.035 cd F9 27(10.47) 0.00 80.95 100.00 148.94 90.85 2.015 cd F10 26(10.08) 0.00 71.43 100.00 148.80 89.89 2.007 d 表型优化核心种质库
Phenotypic optimize core bank32(12.40) 0.00 71.43 100.00 142.19 100.00 2.097 原始种质库
Original bank258 2.028 Table 3. Comparative analysis of banks with different sampling strategies
通过对比不同取样策略下构建的核心种质库的评价参数可以看出:通过性状频率取样策略构建的核心种质库在取样比例最低的前提下,MD为0.00%,VD、CR和VR均高于另外3种取样策略,代表性最好,而且在构建核心种质库的过程中,性状频率可以更好的对无性系进行区分,表明性状频率取样策略在构建美洲黑杨表型核心种质库的研究中是可行的,而且优于随机取样、偏离度取样和位点优先取样的策略。
-
通过对比不同取样策略构建的核心种质库的有效性,最终采用性状频率取样策略构建美洲黑杨表型核心种质库。核心种质库共包括27个无性系,在第1、第2、第3类群中分别选取了2、15、10个无性系,在每个种源中均有无性系被选择(表4),说明取样策略具有一定的代表性和广泛性。
种质库
Banks第1类群 Group 1 第2类群 Group 2 第3类群 Group 3 合计
TotalWas Que Iow 合计
TotalLou Ten Mis 合计
Total原始种质库 Original bank 15 125 19 144 37 51 11 99 258 表型核心种质库 Phenotypic core bank (F9) 2 14 1 15 5 4 1 10 27 表型优化核心种质库 Phenotypic optimize core bank 2 16 1 17 5 7 1 13 32 Table 4. The composition of original, phenotypic core and phenotypic optimize core banks
表型核心种质库的RPR为90.85%(表3),通过与原始种质库进行对比发现:其在14个表型和生理性状上共缺失了19个表现型,为了使RPR达到100.00%,在保留种质中挑选最少数的无性系,使其表型和生理性状包括缺失的19个表现型,最终筛选出5个美洲黑杨无性系为补充种质,补充到表型核心种质中,形成表型优化核心种质库。表型优化核心种质库的取样比例为12.40%,共32个美洲黑杨个体,其中,Was、Que、Iow、Lou、Ten和Mis种源各选择了2、16、1、5、7、1个无性系(表4)。表型优化核心种质库的MD、VD和VR分别为0.00%、71.43%和142.19%,CR和RPR均达到了100.00%,平均表型多样性指数为2.097(表3),表明表型优化核心种质库保留了原始种质库21个表型和生理性状的极值和10个等级的表现型,极大限度地去除了冗余,可以有效地保存美洲黑杨群体的表型多样性。
-
为了进一步对核心种质库进行确认,基于21个表型和生理性状进行主成分分析提取的第1、第2和第3主成分因子得分绘制保留种质、核心种质和补充种质的分布图,从图1可看出:保留种质的分布存在大量的重叠现象,表明种质间表型相似性高,而表型核心种质均匀分布在原始种质范围内,没有重叠现象,说明核心种质库去除了冗余,具有良好的代表性。补充种质进一步的丰富了核心种质的分布范围,形成的表型优化核心种质库具有更好的代表性。