于2020年8—11月蒙古栎种子成熟期，在全国蒙古栎天然分布区调查并收集8省12个种源的蒙古栎种质资源，种源确定依据河北省洪崖山国有林场2018—2019年的全国蒙古栎种质资源踏查结果，结合厉月桥^[15]关于蒙古栎分布的研究，选择蒙古栎主要天然分布区且具有优良林分的蒙古栎天然林，所选种源基本涵盖了蒙古栎主要天然分布区，可代表全国蒙古栎分布情况。所选采种林分为当地分布的天然林分，林龄＞30 a，每个种源选取30株以上采种母树，母树间距50 m以上，收集的种子均已成熟饱满，无病虫害。蒙古栎不同种源的主要生态因子见表1。

8—9月，选取健壮植株，挂牌标记，拍摄植株照片，进行蒙古栎种质资源踏查，并准确记录采种地的经纬度、海拔及年均气温等；9—11月，收集挂牌标记且具代表性的蒙古栎植株种子。

分别对每个种源所采集的蒙古栎单株种子等量混合后，采用四分法随机选取300粒种子进行种长、种宽及种子单粒重、千粒重的测定，用电子游标卡尺进行蒙古栎种宽和种长的测定，精确到0.01 mm，用电子天平进行种子单粒重的测定，用百粒法测量种子千粒重，精确到0.01 g。

用Excel对测量的基础数据进行整理和统计；种形指数为种子的宽(直径)和长(纵径)之比^[16]；运用SPSS 11.5进行方差分析、相关分析(Pearson相关)、主成分分析；同时运用SPSS 11.5将原始数据进行标准化转换，采用欧式距离法进行聚类分析。

如表2所示：不同种源蒙古栎种长、种宽差异较大，种长为11.38~20.96 mm，种宽为11.07~17.80 mm。其中辽宁省抚顺市蒙古栎种子的种长和种宽最长，均值分别为20.96、17.80 mm，其种宽、单粒重显著高于其他地区(P＜0.05)；内蒙古赤峰市最低，种长和种宽均值分别为11.38、11.07 mm，显著低于其他地区(P＜0.05)。河北省承德市围场县与兴隆县蒙古栎种子的种长、种宽无显著差异。种形指数均存在一定差异，但差异较小，其中河北省秦皇岛市最高，山东省泰安市最低，均值分别为1.28、1.02。

单粒重和千粒重差异较大，其中各种源千粒重之间均存在显著差异(P＜0.05)，最高为辽宁省抚顺市(3 027.39 g)，最低为内蒙古赤峰市(1 244.40 g)，可以看出千粒重的变异幅度较大。

如表3所示：不同种源蒙古栎种子表型性状存在一定的相关性。其中蒙古栎种子的种长、种宽、单粒重、千粒重相互之间均存在极显著正相关(P＜0.01)，即种长越长，种宽越宽，单粒重和千粒重随着种长、种宽的增大而增加，表明这4个表型性状是相互影响的。种形指数与种长、种宽、单粒重和千粒重的相关性呈不显著正相关。

从不同种源蒙古栎种子表型性状与生态因子之间的相关性分析(表4)可知：海拔与种宽呈显著负相关(P＜0.05)，与单粒重、千粒重呈极显著负相关(P＜0.01)，相关系数达−0.67、−0.75、−0.75，即海拔越高，种宽、单粒重和千粒重越小。东经与千粒重、单粒重呈显著正相关(P＜0.05)，相关系数均达0.72，即随着经度的增加，单粒重、千粒重增大。无霜期与单粒重呈正相关，即无霜期越长，单粒重越大。种长、种形指数与生态因子之间也存在一定的相关性，但差异不显著，说明蒙古栎种子表型性状的差异与种源生态环境因子有一定的相关性，即存在地理变异特性。

以11个种源间蒙古栎种子相关特性为样本，将种子5个表型性状和6个生态因子数据进行主成分分析。由表5可知：特征值大于1有3个主成分，贡献率分别为43.78%、22.92%、12.57%，累计贡献率达79.27%。其中，第1主成分特征值为5.37，对应的特征向量千粒重和单粒重最大，均为0.41。可以把第1主成分称为种子质量因子。第2主成分特征值为2.52，对应的特征向量无霜期最大，为0.51，可以把第2主成分称为环境因子；第3主成分特征值为1.38，对应的特征向量北纬最小，可以把第3主成分称为地理因子。综上可知，影响蒙古栎种子特性的因素从大到小为单粒重、千粒重、无霜期、北纬。

由聚类分析结果(图1)可以看出：根据单粒重和千粒重，不同种源蒙古栎共分为4类，分别为Ⅰ、Ⅱ、Ⅲ和Ⅳ类种群，其中黑龙江省黑河市、辽宁省本溪市、辽宁省抚顺市为Ⅰ类，种子的单粒重和千粒重远高于所有种源的均值，为最优种源；北京市门头沟区、吉林省延边市为Ⅱ类，单粒重和千粒重略高于均值；河北省张家口市、河北省承德市兴隆县、山东省泰安市、河北省秦皇岛市为Ⅲ类，单粒重和千粒重略低于均值；河北省承德市围场县、内蒙古赤峰市为Ⅳ类，单粒重和千粒重最小。从聚类分析结果看，蒙古栎种子表型性状表现出了一定的区域性，天然种源基本上分布于东北部地区，Ⅰ类分布于东部偏冷地区，Ⅱ类分布于东部和北部地区，Ⅲ类分布于中部地区(华北和华东)，Ⅳ类分布于北部地区，表现出从东北向西南的走向，种源质量逐渐降低，这与相关性分析结果相吻合。

图  1  不同种源蒙古栎种子的聚类分析

Figure 1.  Cluster analysis for Q. mongolica seeds from different provenances

种质资源的调查和评价是研究品种起源、演化和驯化的基础，对优异资源筛选、品种选育具有重要意义^[17]。本研究通过对蒙古栎种子5种表型性状的测定和分析，得出不同种源蒙古栎种实表型性状之间存在差异，其中种长、千粒重均存在显著差异(P＜0.05)，这与厉月桥等^[18]对蒙古栎、李迎超^[19]对栓皮栎Q. variabilis、麻栎Q. acutissima的研究结果一致，且种长与种宽、单粒重、千粒重之间均存在极显著正相关(P＜0.01)，表明蒙古栎种质资源丰富，各生长性状在种源间存在较大的遗传改良潜力。

植物与其生存环境的关系一直是关注的热点^[20]，处于不同生态环境的同种物种，经过长期对当地生境的适应，发生了地理变异，形成了不同的地理种源^[21]。蒙古栎地理分布较为广泛，生存环境存在较大差异。气候条件，如气温、水分等对植物的有性繁殖有着重要的影响^[22]，海拔等对种子、果实的品质也具有影响^[23-27]，而气候条件的波动及其他环境条件的改变会影响物种更新和种群动态^[28-30]。本研究表明：蒙古栎不同种源种子表型特征与生态因子之间均存在相关性，其中海拔与种宽存在显著负相关，这与厉月桥等^[18]对蒙古栎研究结果相同，单粒重和千粒重受分布区海拔的影响存在垂直变异，说明蒙古栎具有广泛的适应性，随着环境条件的变化，出现遗传变异。这与李斌等^[31]对白皮松Pinus bungeana及张清等^[32]对铁橡栎Quercus cocciferoides的研究结果一致，即变异越大，适应环境的能力就越强。单粒重与东经存在显著正相关，随着经度的降低，单粒重和千粒重逐渐降低，种形指数无显著变化，种源质量相对降低，可能由于种子储藏营养物质的量直接影响幼苗建植，大种子有利于幼苗建植早期根系的发育和生物量的积累^[33]，从而具有更强的生长势、逆境耐受能力和对资源的竞争能力^[34]。种子的形状与种长和年均气温的相关性较大，说明蒙古栎种子形状受自身遗传物质控制，同时也受周边环境因子的影响，存在地理变异，由于长期生长在生态条件明显不同的地区，蒙古栎种源间可能产生遗传分化，这与相关的研究结果^{[15, 35]}一致。

本研究表明：影响种子特性的因素按照重要性从大到小依次为单粒重、千粒重、无霜期、北纬。主成分分析表明：贡献率最大的第1主成分表征的特征向量为单粒重和千粒重，因此可将单粒重和千粒重作为快速筛选种源优良与否的优先筛选指标，这与常恩福等^[36]对铁橡栎的研究结果一致。第2、3主成分主要为环境因子和地理因子，这与常恩福等^[36]、张清等^[32]的研究结果存在差异。常恩福等^[36]研究表明：第2、3、4主成分是种形指数、纬度和海拔，即形状因子和生态环境因子；张清等^[32]研究表明：第2主成分为种形指数，可能与树种特性及分布范围存在一定关系。

根据种子性状的综合聚类分析结果，将11个种源按单粒重和千粒重聚为4类，地理分布相近的种源并没有聚在一起，说明了种源间表型性状变异不具有连续性。Ⅰ类种群中，黑河同辽宁地理位置距离较大，但聚为同一类群，可能由于其经度相似，经度和千粒重、单粒重呈显著正相关有关，这与李文文^[37]对蒙古栎种源变异的研究结果一致。Ⅱ类和Ⅲ类种群中，河北和山东相隔距离较大，可能由于同一类群在不同种源之间地形起伏较大，形成了一定程度上的地理隔离，从而造成同一类群的种源，地理距离反而较大，这与李文文^[37]的研究结果一致。Ⅲ类种群除山东泰安外，均位于河北地区，地理距离相对集中，且与Ⅳ类种群相隔地理距离较小，单粒重和千粒重均低于均值，种源质量相对较差。

其中Ⅰ类种群在种源质量方面优于Ⅱ、Ⅲ和Ⅳ类，Ⅰ类分布区均为蒙古栎集中分布区^[15]，同时根据蒙古栎适生分布区的研究可知：不同种源间地理梯度变异明显，呈东北—西南走向，这与李斌等^[31]对白皮松的研究结果相同。

只通过5种表型性状特征确定各种源之间的表型差异存在一定的局限性，因此，还需进一步连续多年收集测定或继续开展分子水平方面的研究，深入分析蒙古栎不同种源之间的遗传变异。