种群内的遗传多样性是种质资源保存的依据。遗传多样性是指种群内基因和基因型的多样性，它体现了物种基因的丰富程度，是研究生态系统多样性和物种多样性的基础^[1-5]。遗传多样性研究通常包括表型、生化、染色体、DNA和碱基序列等多层次遗传变异研究^[6]。尽管表型性状不能完全等同于基因位点，但利用表型性状研究种群的遗传多样性具有简便、高效、经济、直观和实用等优点^[7-12]。目前，国内外学者利用表型变异已在许多珍稀濒危植物上开展了研究^[13-15]，为保护或遗传改良珍稀濒危植物资源奠定了理论基础。

华重楼Paris polyphylla Smith var. chinensis (Franch.) Hara属百合科重楼属的多年生草本植物，作为重楼药材的基源植物收载于2020年版《中华人民共和国药典》中，是著名中成药云南白药、宫血宁胶囊、季德胜蛇药等的主要成分，具有非常高的药用和经济价值^[16-19]。华重楼主要分布在我国江苏、浙江、江西、福建、广东、广西、湖北、湖南和四川等南方省份，自然生长于海拔600~1350 m的林下荫处或沟谷边的草丛中^[20]。华重楼作为一种珍贵的中药材，主要以根茎入药，药用活性成分为重楼皂苷，具有抗癌、止血、消炎、镇痛、止咳平喘和免疫调节等作用，常用于治疗蛇毒、乳腺炎、跌打损伤等病症^[21-22]。虽然华重楼具有很高的药用和经济价值，在长江以南省份也偶见人工栽培，但研究主要集中在组培技术^[23-24]、人工栽培^[25]和药用成分分析^[26-27]等方面，而遗传育种的基础性工作地理种群遗传变异少见有相关报道。本文以华重楼7个地理种群为研究对象，系统测定叶片、茎、花萼、果球和种子等的17个表型性状，运用统计学分析方法，深入研究种群间和种群内的变异规律，以期为华重楼种质资源保护、药材资源开发利用和繁育改良提供理论支撑。

通过查阅中国植物数字标本馆( http://www.cvh.ac.cn/)和电子版《中国植物志》( http://www.iplant.cn/frps)初步确定华重楼的分布区域后，进一步通过文献检索、咨询当地相关专家和野外实地调查，最终确定华重楼的自然地理分布地，并于2018和2019年的9—11月间进行野生华重楼植株、叶片、花萼、果球和种子的样品采集和数据调查、测量工作。野生华重楼的样品数据采集点分别位于江西省新余市(JXXY)、浙江省丽水市(ZJLS)、四川省雅安市(SCYA)、广东省清远市(GDQY)、福建省三明市(FJSM)、湖南省娄底市(HNLD)和湖北省宜昌市(HBYC)，共计7个地理种群。数据采集时，选择生长情况正常、无明显的病虫害发生、无明显的生长畸形、达到生殖成熟、有正常饱满的果球的单株样品进行调查测量。鉴于华重楼天然种质资源分布稀少，调查时各地理种群单株保有量只需达到30株以上、样品单株间隔10 m即可。因野生华重楼多分布于深山沟壑，数据收集工作比较困难，故在选择样品单株时，各地理种群只调查15株，7个地理种群共计105株。调查测量时，现场记录样点的地理经纬度数据、海拔高度数据，其他气候因子数据通过查询资源学科领域基础科学数据整合与集成应用人地系统主题数据库( http://www.data.ac.cn)获得，具体数据见表1。

表 1 Table 1

表 1 7个野生华重楼地理种群的地理及主要气候因子情况 Table 1 Sample collection information of Paris polyphylla var. chinensis from different distribution regions and climatic data at the sampling sites 样品采集点 Sample site 种群 Population 经度 Longitude 纬度 Latitude 海拔/m Altitude 年均温度/℃ Mean annual temperature 无霜期/d Frost-free period 年日照/h Annual sunshine hour 年降水量/mm Annual precipitation 江西新余 Xinyu in Jiangxi JXXY 114°37′22″E 27°31′45″N 376 17.7 270 1 655.4 1 594.8 浙江丽水 Lishui in Zhejiang ZJLS 119°04′18″E 28°38′28″N 1 103 11.3 240 1 635.1 1 598.9 四川雅安 Yaan in Sichuan SCYA 103°01′57″E 30°30′24″N 1 450 17.0 295 1 519.2 1 350.0 广东清远 Qingyuan in Guangdong GDQY 112°20′11″E 25°10′10″N 648 18.0 210 2 290.0 1 329.0 福建三明 Sanming in Fujian FJSM 117°29′18″E 26°33′24″N 548 18.2 260 1 727.1 1 688.0 湖南娄底 Loudi in Hunan HNLD 112°02′62″E 27°48′05″N 290 17.1 270 1 726.0 1 455.6 湖北宜昌 Yichang in Hubei HBYC 111°21′19″E 30°27′44″N 400 16.9 275 1 710.5 1 215.5

表 1 7个野生华重楼地理种群的地理及主要气候因子情况 Table 1 Sample collection information of Paris polyphylla var. chinensis from different distribution regions and climatic data at the sampling sites

统计待调查表型性状的野生华重楼单株的叶片数和花萼片数，摘取其所有叶片和花萼片，平铺在硬纸板上，并利用数显游标卡尺测量其叶长、花萼长、叶宽、花萼宽，同时测量地径、果球横径和果球纵径等指标，精确到0.01 mm，测量结束后，计算叶片长宽比、花萼长宽比和果球横纵比；利用便携式叶面积测定仪(LI-3000C，LI-COR，Lincoln，USA)测定所有叶面积和花萼面积，精确到0.01 cm²；利用不锈钢卷尺测量每株野生华重楼的株高和花梗长，精确到0.1 cm；采集成熟种子，人工搓掉红色种皮，用清水冲洗干净，置于通风阴凉处，晾干后及时利用电子分析天平测定每个种群种子的千粒质量，精确到0.01 g。以上数据测量，每个种群选取5个样品单株，重复3次，每个种群测量15株，共计105株。

利用社会统计学分析软件SPSS25.0对数据进行处理分析。应用表型分化系数、单因素方差分析以及多因素方差分析方法分析种群间和种群内的变异情况；利用相关分析来探讨表型因子与地理因子的关系。

野生华重楼种群间和种群内的表型性状变异的方差分析结果见表2。所选取的17个表型性状除叶片数和花萼片数外，其余15个表型性状指标在种群间和种群内都存在显著(P<0.05)或极显著(P<0.01)的差异，这说明野生华重楼绝大多数表型性状在种群间和种群内存在广泛的变异。

表 2 Table 2

表 2 野生华重楼种群间和种群内全株、果、叶及种子表型性状的方差分析结果1) Table 2 Variance analysis of phenptypic traits of wholeplant, leaf, fruit and seeds among and within Paris polyphylla var. chinensis 性状 Trait 均方 Mean square 随机误差 Random error F 种群间 Among populations 种群内 Within populations 种群间 Among populations 种群内 Within populations 株高 Height 1 060.94 226.26 266.13 4.69** 5.97** 地径 Diameter 6.39 1.32 1.62 4.83** 7.13** 叶片数 No. of leaf 0.53 0.78 0.08 0.68 1.07 叶长 Leaf length 5 728.49 601.72 264.04 9.52** 7.18** 叶宽 Leaf width 1 043.13 43.62 51.48 23.91** 20.99** 叶面积 Leaf area 2 831.39 142.52 138.48 19.87** 3.56* 叶片长宽比 Blade aspect ratio 3.73 0.84 0.15 4.45** 5.39** 花萼片数 No. of calyx 4.52 0.55 1.22 0.82 0.79 花萼长 Calyx length 762.13 49.75 46.97 15.32** 10.32** 花萼宽 Calyx width 119.22 9.90 6.86 12.04** 11.31** 花萼面积 Calyx area 46.00 3.32 2.81 13.86** 7.15** 花萼长宽比 Calyx aspect ratio 0.55 0.13 0.03 4.12** 5.01** 花梗长 Pedicel length 152.91 70.67 25.22 3.16* 6.25** 果球横径 Fruit transverse diameter 89.74 18.79 5.95 4.78** 6.17* 果球纵径 Fruit longitudinal diameter 256.18 25.35 19.35 10.10** 8.33** 果球横纵比 Fruit transverse diameter to longitudinal diameter ratio 0.24 0.11 0.02 21.60** 17.05** 种子千粒质量 Seed 1 000 -grain weight 843.51 69.95 229.20 12.06** 10.87** 　1)“*”和“**”分别表示在0.05和0.01水平差异显著(单因素方差分析) 　1)“*” and “**” mean significant differences at 0.05 and 0.01 levels respectively(One-way ANOVA)

表 2 野生华重楼种群间和种群内全株、果、叶及种子表型性状的方差分析结果1) Table 2 Variance analysis of phenptypic traits of wholeplant, leaf, fruit and seeds among and within Paris polyphylla var. chinensis

野生华重楼7个地理种群17个表型性状的平均值、标准差及多重比较分析结果见表3。结果显示，除叶片数和花萼片数在各种群间差异不显著外，其余15个表型性状指标在各种群间均存在显著差异。其中，四川雅安种群的叶长、花萼长、叶宽、花萼宽、叶面积、花萼面积、株高和地径8个表型性状指标数值最大，与其他6个地理种群大多为差异显著；浙江丽水种群的果球纵径最大，叶片长宽比和果球横纵比最小；广东清远种群的花梗长、叶片长宽比和花萼长宽比最大，果球纵径最小；湖北宜昌种群的花梗长、叶长、花萼长、花萼宽、花萼长宽比、叶面积、种子千粒质量和地径8个表型性状指标为所有种群中最小，果球横纵比最大；湖南娄底种群果球横径和种子千粒质量最大；福建三明种群的果球横径最小，显著小于湖南娄底和浙江丽水种群。

表 3 Table 3

表 3 野生华重楼7个地理种群20个表型性状的平均值、标准差及多重比较1) Table 3 Mean, standard deviation and multiple comparison of 20 phenotypic traits in seven geographical populations of Paris polyphylla var. chinensis 种群2) Population 株高/cm Height 地径/mm Diameter 叶片数 No. of leaf 叶长/mm Leaf length 叶宽/mm Leaf width 叶面积/cm2 Leaf area JXXY 53.67±17.64ab 8.22±1.17ab 6.33±0.58a 134.91±21.64ab 44.89±10.01b 48.84±17.70b ZJLS 66.43±8.46a 7.93±0.70abc 6.67±1.15a 121.90±13.43bc 48.31±4.10ab 46.48±7.92b SCYA 69.20±12.51a 8.63±1.33a 7.00±1.00a 148.38±12.29a 51.85±4.90a 60.28±5.99a GDQY 49.63±24.01abc 6.47±0.84bcd 6.67±1.53a 114.58±55.20c 32.73±7.84d 30.50±21.37cd FJSM 35.33±5.81bc 5.29±0.78d 6.33±0.58a 123.02±10.82bc 35.80±3.96cd 34.73±5.96c HNLD 22.13±2.25c 6.05±0.84cd 6.67±0.58a 112.69±10.17c 38.91±5.91c 34.57±6.84c HBYC 25.83±2.25c 5.02±0.59d 5.67±0.58a 93.57±6.93d 35.70±7.53cd 26.35±6.28d 种群2) Population 叶片长宽比 Blade aspect ratio 花萼片数 No. of calyx 花萼长/mm Calyx length 花萼宽/mm Calyx width 花萼面积/cm2 Calyx area 花萼长宽比 Calyx aspect ratio JXXY 3.06±0.37bc 5.67±0.58a 49.91±7.87a 19.13±4.13a 7.68±2.56a 2.66±0.34ab ZJLS 2.53±0.24c 6.67±1.15a 49.34±5.34a 18.46±1.74ab 7.20±1.29a 2.68±0.26a SCYA 2.89±0.41b 5.80±0.84a 52.06±8.90a 19.23±3.43a 8.03±2.56a 2.73±0.35a

表 3 野生华重楼7个地理种群20个表型性状的平均值、标准差及多重比较1) Table 3 Mean, standard deviation and multiple comparison of 20 phenotypic traits in seven geographical populations of Paris polyphylla var. chinensis

表 3 Table 3

续表 3　Continued table 3 种群2) Population 叶片长宽比 Blade aspect ratio 花萼片数 No. of calyx 花萼长/mm Calyx length 花萼宽/mm Calyx width 花萼面积/cm2 Calyx area 花萼长宽比 Calyx aspect ratio GDQY 3.79±2.18a 5.33±0.58a 38.80±8.31b 13.79±1.90c 4.29±1.37b 2.81±0.42a FJSM 3.46±0.38ab 4.67±0.58a 36.71±8.41b 13.42±2.52c 3.97±1.59b 2.76±0.52a HNLD 2.95±0.45bc 5.67±0.58a 39.23±5.24b 16.56±2.48b 5.18±1.42b 2.39±0.25bc HBYC 2.72±0.51c 4.33±0.58a 34.19±2.86b 13.04±4.75c 4.26±1.16b 2.26±0.48c 种群2) Population 花梗长/cm Pedicel length 果球横径/mm Fruit transverse diameter 果球纵径/mm Fruit longitudinal diameter 果球横纵比 Fruit transverse diameter to longitudinal diameter ratio 种子千粒质量/g Seed 1 000-grain weight JXXY 18.17±5.03ab 30.45±5.25bc 35.93±7.28ab 0.86±0.10d 100.16±3.00a ZJLS 19.60±2.15ab 34.09±2.84ab 39.94±3.37a 0.85±0.04d 72.25±2.55b SCYA 18.60±3.47ab 29.50±4.41c 34.40±4.34bc 0.86±0.08d 73.25±2.55b GDQY 23.00±15.86a 29.79±3.46bc 26.33±4.40d 1.14±0.12b 77.49±10.10b FJSM 7.83±0.91ab 28.92±2.65c 30.47±3.37cd 0.96±0.88cd 95.76±4.17a HNLD 10.43±4.29ab 36.83±6.65a 36.06±5.95ab 1.02±0.11c 109.99±17.00a HBYC 5.80±1.06b 32.68±2.79abc 26.64±3.94d 1.25±0.19a 68.14±3.37b 　1)同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’s法)；2)JXXY：江西新余；ZJLS：浙江丽水；SCYA：四川雅安；GDQY：广东清远；FJSM：福建三明；HNLD：湖南娄底；HBYC：湖北宜昌 　1) Different lowercase letters in the same column indicated significant differences among different geographical populations (P<0.05, Duncan’s method)；2) JXXY: Xinyu in Jiangxi; ZJLS: Lishui in Zhejiang; SCYA: Ya’an in Sichuan; GDQY: Qingyuan in Guangdong; FJSM: Sanming in Fujian; HNLD: Loudi in Hunan; HBYC: Yichang in Hubei

野生华重楼7个地理种群的17个表型性状的方差分量及其种群间表型分化系数见表4。17个表型性状的平均方差分量在种群间占比26.39%、种群内占比52.19%、随机误差占比21.42%，就平均方差分量所占百分比而言，种群间只是种群内的一半，这说明野生华重楼生长表型在种群间的变异远小于种群内的变异；17个表型性状的平均分化系数种群间的贡献值为34.27% (即种群内的贡献值为65.73%)，这就进一步表明，野生华重楼的表型变异主要源于种群内部。

表 4 Table 4

表 4 野生华重楼表型性状的方差分量及种群间表型分化系数 Table 4 Variance component and differentiation coefficients of phenotypic traits among and within populations in Paris polyphylla var. chinensis 性状 Trait 方差分量 Variance component 方差分量占比/% Percentage of variance component 表型分化 系数/% Phenotypic differentiation coefficient 种群间 Among populations 种群内 Within populations 随机误差 Random error 种群间 Among populations 种群内 Within populations 随机误差 Random error 株高 Height 352.96 464.74 266.13 32.57 42.88 24.55 43.16 地径 Diameter 2.13 2.77 1.62 32.67 42.48 24.85 43.47 叶片数 No. of leaf 0.18 0.72 0.08 18.37 73.47 8.16 20.00 叶长 Leaf length 301.64 829.57 264.04 21.62 59.46 18.92 26.67 叶宽 Leaf width 52.36 88.04 51.48 27.29 45.88 26.83 37.29 叶面积 Leaf area 144.13 262.02 138.48 26.46 48.11 25.43 35.49 叶片长宽比 Blade aspect ratio 0.19 0.97 0.15 14.50 74.05 11.45 16.38 花萼片数 No. of calyx 1.47 1.63 1.22 34.03 37.73 28.24 47.42 花萼长 Calyx length 53.20 90.08 46.96 27.96 47.35 24.68 37.13 花萼宽 Calyx width 7.75 15.81 6.86 25.48 51.97 22.55 32.89 花萼面积 Calyx area 3.14 5.74 2.81 26.86 49.10 24.04 35.36 花萼长宽比 Calyx aspect ratio 0.04 0.16 0.03 17.39 69.57 13.04 20.00 花梗长 Pedicel length 44.11 93.10 25.22 27.16 57.32 15.53 32.15 果球横径 Fruit transverse diameter 8.46 23.80 5.95 22.14 62.29 15.57 26.22 果球纵径 Fruit longitudinal diameter 26.53 41.64 19.35 30.31 47.58 22.11 38.92 果球横纵比 Fruit transverse diameter to longitudinal diameter ratio 0.02 0.03 0.02 28.57 42.86 28.57 40.00 种子千粒质量 Seed 1 000-grain weight 272.61 271.75 229.20 35.24 35.13 29.63 50.08 平均值 Mean 26.39 52.19 21.42 34.27

表 4 野生华重楼表型性状的方差分量及种群间表型分化系数 Table 4 Variance component and differentiation coefficients of phenotypic traits among and within populations in Paris polyphylla var. chinensis

野生华重楼各地理种群表型性状的变异系数列于表5。野生华重楼17个表型性状在种群内的平均变异系数为17.72%，变异幅度为6.83%~29.95%。花梗长、叶面积、花萼面积和株高4个表型性状指标的平均变异系数均大于20%，超过所有表型性状指标的平均变异系数(17.72%)；叶片数、花萼片数、果球横径、果球纵径、果球横纵比和种子千粒质量6个表型性状指标的平均变异系数均小于14.47%，小于所有表型性状指标的平均变异系数(17.72%)，表明华重楼叶片数量、花萼数量、果球和种子的性状稳定性大于其他性状指标的稳定性。进一步比较7个地理种群野生华重楼17个表型性状指标的平均变异系数可以得到，广东清远、湖南娄底和江西新余3个地理种群的变异系数均超过平均变异系数，表明广东清远、湖南娄底和江西新余野生华重楼种群的表型变异和多样性较为丰富；湖北宜昌、浙江丽水、四川雅安和福建三明4个地理种群的变异系数均小于平均变异系数，说明湖北宜昌、浙江丽水、四川雅安和福建三明地理种群的表型变异和多样性较广东清远、湖南娄底和江西新余地理种源弱。

表 5 Table 5

表 5 野生华重楼种群17个表型性状变异系数 Table 5 Variation coefficients of 17 phenotypic traits in Paris polyphylla var. chinensis populations % 性状 Trait 种群1) Population 平均值 Mean JXXY ZJLS SCYA GDQY FJSM HNLD HBYC 株高 Height 32.86 12.73 18.08 48.40 16.44 72.02 8.73 29.89 地径 Diameter 14.22 8.86 15.42 27.17 14.79 13.89 11.80 15.16 叶片数 No. of leaf 9.12 17.32 14.29 17.28 9.12 9.12 10.19 12.35 叶长 Leaf length 16.04 11.02 8.28 48.18 8.79 9.03 7.40 15.53 叶宽 Leaf width 22.31 8.48 9.46 23.96 11.05 15.18 21.10 15.93 叶面积 Leaf area 36.24 17.04 9.93 70.08 17.15 19.78 23.81 27.72 叶片长宽比 Blade aspect ratio 12.06 9.35 14.31 57.54 10.96 15.21 18.87 19.76 花萼片数 No. of calyx 10.19 17.32 10.83 14.43 15.75 10.19 17.32 13.72 花萼长 Calyx length 15.76 10.81 17.10 21.42 22.90 13.35 8.38 15.67 花萼宽 Calyx width 21.60 9.42 17.83 13.77 18.75 15.00 36.47 18.98 花萼面积 Calyx area 33.36 17.87 31.85 31.98 39.97 27.35 27.26 29.95 花萼长宽比 Calyx aspect ratio 12.86 9.70 12.68 15.02 19.14 10.61 21.25 14.47 花梗长 Pedicel length 27.71 10.98 18.63 68.72 11.58 41.07 18.25 28.13

表 5 野生华重楼种群17个表型性状变异系数 Table 5 Variation coefficients of 17 phenotypic traits in Paris polyphylla var. chinensis populations

表 5 Table 5

续表 5　Continued table 5 性状 Trait 种群1) Population 平均值 Mean JXXY ZJLS SCYA GDQY FJSM HNLD HBYC 果球横径 Fruit diameter 17.26 8.34 14.96 11.63 9.16 18.07 8.55 12.57 果球纵径 Fruit longitudinal diameter 20.25 8.45 12.62 16.72 11.05 16.49 14.79 14.34 果球横纵比 Fruit aspect ratio 11.71 4.73 9.25 10.52 9.26 10.85 15.01 10.19 种子千粒质量 Seed 1 000-grain weight 3.00 3.52 3.52 13.03 4.35 15.45 4.94 6.83 平均值 Mean 18.62 10.94 14.06 29.99 14.72 19.57 16.12 17.72 　1) JXXY：江西新余；ZJLS：浙江丽水；SCYA：四川雅安；GDQY：广东清远；FJSM：福建三明；HNLD：湖南娄底；HBYC：湖北宜昌 　1) JXXY: Xinyu in Jiangxi；ZJLS: Lishui in Zhejiang；SCYA: Ya’an in Sichuan；GDQY: Qingyuan in Guangdong；FJSM: Sanming in Fujian；HNLD: Loudi in Hunan；HBYC: Yichang in Hubei

对野生华重楼全株、茎、叶、花、果球和种子17个表型性状与样品采集点地理气象因子之间的相关性进行了分析，结果表明，叶长、叶宽和叶面积与样品采集点的经度呈现显著或极显著的负相关关系，说明随着经度的不断增加，野生华重楼叶长、叶宽和叶面积将不断减小；叶宽、叶面积、花萼长、花萼宽、花萼面积和果球纵径与纬度呈现显著或极显著的正相关关系，叶片长宽比和果球横纵比与纬度呈现极显著的负相关关系，这就说明随着纬度的增加，叶宽、叶面积、花萼长、花萼宽、花萼面积、果球纵径将呈现上升趋势，而叶片长宽比、果球横纵比将呈现下降趋势；叶长、叶宽、叶面积、花萼长、花萼宽、花萼面积、花萼长宽比、株高和地径与海拔呈现显著或极显著的正相关关系，果球横纵比和种子千粒质量与海拔呈现极显著的负相关关系，说明海拔越高，叶片、花萼、株高、地径也随之越大，果球横纵比和种子越小(表6)；叶宽、叶面积、花萼长、花萼面积、果球横径和果球纵径与年平均温度呈显著或极显著的负相关关系，叶片长宽比、花萼宽与年平均温度呈极显著的正相关关系，说明随着年平均温度的升高，叶片将变得小而长，种球也随之变小；叶长、叶宽、叶面积、花萼长、花萼宽、花萼面积和果球纵径与无霜期呈显著或极显著的正相关关系，果球横纵比则与无霜期呈现极显著的负相关关系，说明无霜期越长，叶片和花萼越大，果球越大；叶长、叶面积、花萼长、花萼宽、花萼面积与年日照时数呈显著或极显著负相关关系，说明日照时间越长，叶片和萼片越小；叶长、花萼长、花萼宽、果球纵径和种子千粒质量与年降水量呈显著或极显著正相关，果球横纵比则与年降水量呈极显著负相关，这说明年降水量越大，叶片越长，种子越大(表6)。对表6的综合分析发现，野生华重楼的叶片数、花萼片数和花梗长与7个生态地理环境因子无显著相关关系，说明这些生态地理环境因子对野生华重楼这一物种的种群形态影响不大。

表 6 Table 6

表 6 野生华重楼17个表型性状与生态地理因子的相关系数1) Table 6 Correlation coefficients between 17 phenotypic traits of Paris polyphylla var. chinensis and ecological geography factors 性状 Trait 经度 Longitude 纬度 Latitude 海拔 Altitude 年均温度 Mean annual temperature 无霜期 Frost-free period 年日照时数 Annual sunshine hour 年降水量 Annual precipitation 株高 Height −0.123 0.085 0.671** −0.353 −0.072 −0.109 0.135 地径 Diameter −0.213 0.172 0.530* −0.293 0.135 −0.267 0.157 叶片数 No. of leaf −0.134 −0.074 0.286 −0.077 −0.041 0.007 0.088 叶长 Leaf length −0.206* 0.072 0.354** 0.010 0.200* −0.288** 0.188* 叶宽 Leaf width −0.222** 0.444** 0.530** −0.365** 0.373** 0.535** 0.133 叶面积 Leaf area −0.277** 0.301** 0.512** −0.179* 0.333** −0.425** 0.154 叶片长宽比 Blade aspect ratio 0.010 −0.340** −0.114 0.284** −0.219* 0.324** 0.008 花萼片数 No. of calyx 0.078 −0.230 0.299 −0.317 −0.344 0.152 0.151 花萼长 Calyx length −0.108 0.267** 0.449** −0.291** 0.199* −0.372** 0.210* 花萼宽 Calyx width −0.075 0.199* 0.314** 0.256** 0.226* −0.385** 0.255** 花萼面积 Calyx area −0.136 0.343** 0.386** −0.268** 0.287** −0.435** 0.170 花萼长宽比 Calyx aspect ratio 0.022 −0.186 0.229* −0.011 −0.176 0.118 0.146 花梗长 Pedicel length −0.071 −0.274 0.312 −0.114 −0.377 0.311 −0.036 果球横径 Fruit transverse diameter 0.059 0.081 −0.176 −0.230* 0.037 −0.063 −0.053 果球纵径 Fruit longitudinal diameter −0.027 0.328** 0.115 −0.362** 0.390** −0.495** 0.266* 果球横纵比　Fruit transverse diameter to longitudinal diameter ratio 0.098 −0.337** −0.317** 0.210 −0.490** 0.594** −0.423** 种子千粒质量 Seed 1 000-grain weight 0.232 −0.272 −0.535** 0.306 0.178 −0.135 0.498* 　1)“*”和“**”分别表示在0.05和0.01水平(双侧)上显著相关 　1)“*”and“**”indicated significant correlations at 0.05 and 0.01 levels (double-tail), respectively

表 6 野生华重楼17个表型性状与生态地理因子的相关系数1) Table 6 Correlation coefficients between 17 phenotypic traits of Paris polyphylla var. chinensis and ecological geography factors

通过对野生华重楼7个地理种群的17个性状和生长表型指标的深入研究发现，无论是在种群间还是在种群内，除叶片数和花萼片数外的其他指标均存在显著或极显著的变异，这与多花黄精Polygonatum cyrtonema ^[2,28]、人参Panax ginseng ^[29]、毕节小檗Berberis guizhouensis^[30]、八角莲Dysosma versipellis^[31]等林下经济植物的研究结果类似。通过种群间的表型分化系数分析可知，种群间的变异远小于种群内的变异，这就说明种群内的变异是野生华重楼表型变异的主要来源。华重楼17个表型性状的变异幅度为16.38%~50.08%，平均表型分化系数为34.27%，种群间的变异仅贡献34.27%，远小于种群内的变异贡献65.73%，这一现象出现的可能原因：一是各种群的自然分布有着极强的区域性，这就增加了不同种群间的基因交流难度；二是野生华重楼的生长发育与林下小生境关系密切。综上所述，野生华重楼种群间的基因交流机会要小于种群内部的基因交流，说明华重楼的表型变异主要源自于各地理种群内部。

变异系数在统计学中又称离散系数，是衡量离散程度的一个归一化量度指标，离散系数越大，则说明观测对象的变异越大，则对应的种群表型变异和多样性越丰富。通过对野生华重楼各地理种源表型指标的变异系数分析发现，所有指标的变异系数幅度在6.83%~29.95%之间，平均变异系数为17.72%，变异系数的变化幅度较大，种群内部的表型变异程度高，说明华重楼能够通过种群内部的遗传变异，更好地适应环境。

野生华重楼叶面积、花萼面积、株高、地径、果球横径、果球纵径、果球横纵比和种子千粒质量8个表型指标的变异系数分别为27.72%、29.95%、29.89%、15.16%、12.57%、14.34%、10.19%和6.83%，前3个指标远大于平均变异系数，后5个指标远小于平均变异系数，这些指标的稳定性排名依次为种子千粒质量>果球横纵比>果球横径>果球纵径>地径>叶面积>株高>花萼面积，说明野生华重楼种子的稳定性最高，果球的稳定性其次，茎的稳定性排第3，叶片和花萼的稳定性最差，这与多花黄精Polygonatum cyrtonema这方面性状的研究结果^[2]是类似的。

本研究发现，所选取的野生华重楼17个表型性状指标中，除叶片数、花萼片数和花梗长3个表型性状指标与所有地理气象因子无显著相关性外，其余性状指标均与部分地理气象因子存在显著或极显著的相关性，这说明地理气象因子对野生华重楼的种群表型影响较大。经度的增加会抑制叶长、叶宽和叶面积的生长，纬度的增加会促进叶片和果球的生长；海拔与野生华重楼的叶、茎营养生长呈现正相关关系，与果球和种子的生殖生长呈现负相关关系，说明海拔增加会促进营养生长，抑制生殖生长；年均温度、无霜期、年日照时长和年降雨量等环境因子均与绝大多数表型指标存在显著或极显著相关性，说明其表型变异受多个环境因子的共同控制。

华重楼是重楼药材的法定基原，是许多知名中成药不可或缺的主要组分材料，随着制药工业的不断发展，人们对华重楼的需求越来越大。因华重楼人工栽培种植周期长、生境要求苛刻、种植技术难，人们为了满足需求对其进行掠夺性、毁灭性的采挖，导致野生资源枯竭，处于濒危状态。

因此，开展野生华重楼资源保护的研究是非常有必要的，笔者结合研究结果，提出以下保护策略：首先，建议将遗传变异多样性丰富的华重楼地理种群原生地——广东清远设定为重点育种保护地；其次，严禁破坏野生华重楼的原生环境，为华重楼种群的恢复提供生境支撑；第三，及时开展华重楼人工繁育研究，争取在短时间内解决种苗繁育问题，为市场提供充足的人工种苗；最后，收集全分布区华重楼种质资源，建立保存圃，做好基因资源的抢救性保存工作。