Screening and identification of rice salt-tolerant mutants at germination stage induced by space mutagenesis and heavy ion radiation
-
摘要:目的
筛选萌发期耐盐突变体，为耐盐水稻新品种的创制与培育提供优异的种质资源。
方法对优质籼稻品种‘航聚香丝苗’干种子进行空间诱变和重离子诱变处理，并对诱变后代进行耐盐突变体的定向筛选。将野生型(WT)‘航聚香丝苗’种子用不同浓度NaCl溶液处理，探索萌发期耐盐突变体筛选的适宜浓度，以发芽率为指标对5 205份M2代诱变材料进行筛选，并在M3代进行验证。对WT和突变体进行耐盐相关生理指标的测定与分析，考查农艺性状；采用48个SSR标记对诱变后代与WT的遗传相似性进行评估。
结果萌发期耐盐突变体筛选的适宜NaCl溶液浓度为340 mmol/L，筛选得到12份萌发期耐盐突变体。生理指标测定结果显示，2个代表性突变体9-8和22-7的SOD、POD、CAT、APX活性均显著高于WT的，MDA、H2O2和O2−含量均显著低于WT的，耐盐突变体呈现出更强的抗氧化能力。突变体和野生型的农艺性状基本一致。SSR标记检测结果表明，WT ‘航聚香丝苗’与突变体19-49、20-15、20-49、20-62、15-6、16-65、22-7、20-19、20-38均无SSR差异位点，与突变体9-8、19-18、42-113的差异SSR位点均为1个。
结论筛选到12份萌发期耐盐突变体，且与WT ‘航聚香丝苗’的遗传背景高度一致，对水稻抗盐新品种的培育及抗盐性遗传机理的研究具有一定的理论和实际意义。
Abstract:ObjectiveTo screen salt-tolerant mutants during germination, provide excellent germplasm resources for creating and breeding new salt-tolerant rice varieties.
MethodThe dried seeds of the high quality indica rice variety ‘Hangjuxiangsimiao’, were treated with spatial mutagenesis and heavy ion mutagenesis, and the mutagenised progeny were subjected to targeted screening for salt-tolerant mutants. Wild type (WT) ‘Hangjuxiangsimiao’ seeds were treated with different concentrations of NaCl solution to explore the appropriate concentration for screening salt-tolerant mutants at the germination stage, and then 5 205 M2 mutant materials were screened using germination rate as an index and verified in the M3 generation. WT and mutants were measured and analysed for salt-tolerance-related physiological indices and agronomic traits. 48 SSR markers were used to assess the genetic similarity between the mutant progeny and WT.
ResultThe appropriate NaCl concentration for screening salt-tolerant mutants at the germination stage was 340 mmol/L, and 12 salt-tolerant mutants were screened at the germination stage. The results of physiological indices showed that the SOD, POD, CAT and APX activities of two representative mutants 9-8 and 22-7 were significantly higher than those of WT, and the MDA, H2O2 and O2− contents were significantly lower, the salt-tolerant mutants exhibited stronger antioxidant capacity. The results of SSR marker detection showed that there were no SSR differentiation locus between WT seedlings and mutants of 19-49, 20-15, 20-49, 20-62, 15-6, 16-65, 22-7, 20-19 and 20-38, and one differentiation SSR locus between them and mutants of 9-8, 19-18, 42-113.
ConclusionTwelve salt-tolerant mutants were screened at the germination stage, which were highly consistent with the genetic background of wild-type ‘Hangjuxiangsimiao’. It has certain theoretical and practical significance for breeding new salt-tolerant rice varieties and studing the genetic mechanism of salt tolerance.
-
Keywords:
- Rice /
- Space mutation /
- Heavy ion radiation /
- SSR marker /
- Mutant screening /
- Germination stage /
- Salt tolerance
-
杉木Cunninghamia lanceolata在我国已有上千年栽培历史，因其生长迅速、产量高、材质优良等特性[1]，长期以来深受人们的青睐而得以广泛种植。随着人们对杉木需求的日益高端化、高值化和生态化，杉木大径材培育对国家木材安全、生态安全、绿色发展越来越具有战略性意义[2-3]。近年来，许多学者致力于杉木大径材定向培育技术研究，提出了不同立地指数下的适宜保留密度范围[4-5]，针对不同遗传品质苗木及土壤条件提出了相应的抚育、施肥等经营措施[6-7]，运用林分经营模型对不同地区提出了间伐时间、间伐强度等优化管理模式[8]。但涉及杉木大径材材种形成的主要影响因子及其影响规律等方面的研究鲜有报道[9]。粤北、粤西北、粤东北等南岭山区丘陵地带是我国杉木四大中心产区之一，其地理位置优越，受季风环流影响，降水丰富，气候温暖湿润[10]，拥有培育优质高效杉木大径材林分得天独厚的自然条件。本研究基于广东省南岭山区特殊的自然生态条件，将现有杉木大径材林分特征因子数量化或标准化，形成综合反映林分控制因子和生态环境因子的指数评价体系，通过挖掘大径材形成主导影响因子并探究其影响规律，构建南岭山区杉木大径材出材量的预测模型，以期为广东省杉木大径材培育和经营提供理论指导，推动杉木由高密度、中小径材为经营目标向以生态为先导的大径材经营模式转变。
1. 研究区域自然概况
本研究在对粤北、粤西北、粤东北等南岭山区范围内现存具有杉木大径材的林分全面踏查的基础上，选择广东省肇庆市国有新岗林场，韶关市国有韶关林场、国有小坑林场、国有河口林场、国有大瑶山林场，广东省连山林场、广东连山笔架山省级自然保护区等7个地点，共计36块保存完好、长势优良的典型杉木大径材林分。研究区域位于111°57′~113°81′E，24°13′~25°55′N，地处南岭山脉西南至南麓，地貌以山地丘陵为主，土层深厚，土壤为赤红壤和黄壤；属中亚热带季风气候，年平均气温20 ℃左右，年均降雨量1600 mm，相对湿度80%。
2. 研究方法
2.1 标准样地设置与调查
2018—2019年，在上述7个地点的杉木大径材林分中选择不同坡向和坡位，设置400 m2(20 m×20 m)标准样地，4个角点埋设pvc管作为标志，共计设置标准样地36个。对标准样地内的杉木进行每木检尺，测量胸径、树高、冠幅、枝下高、干形等，同时按照梅花形选择5株优势木，测其树高，求优势木平均高。记录标准样地的郁闭度、经纬度、海拔、坡形、坡度、坡向、坡位等，并收集造林时间、造林密度、造林种苗情况以及抚育、施肥等林分经营历史和背景资料。标准样地基本情况见表1。
表 1 标准样地基本情况Table 1. Basic situation of standard plot地点
Site样地数
No. of
plots林龄
Stand
age保留密度/
(株·hm−2)
Reserved density平均胸径/cm
Average
DBH平均树高/m
Average
tree height优势木平均高/m
Average height of
dominant wood立地指数
Site
index国有新岗林场 Xingang State Forest Farm 9 23~28 1075~1626 21.7~23.4 15.0~17.2 17.9~23.8 16~22 国有韶关林场 Shaoguan State Forest Farm 3 20 1375~1405 20.5~21.6 12.5~14.8 16.1~17.2 16 国有小坑林场 Xiaokeng State Forest Farm 9 26~32 1150~1650 21.2~25.2 15.8~18.0 19.4~21.2 18~20 国有河口林场 Hekou State Forest Farm 3 26 1300 23.4~28.4 18.2~21.1 24.0~25.4 22 国有大瑶山林场
Dayaoshan State Forest Farm6 25 1176~1450 21.2~23.5 14.9~16.2 16.8~18.1 16 广东省连山林场
Guangdong Lianshan Forest Farm3 27 1150~1225 24.7~26.7 17.9~18.2 20.4~22.6 18~20 连山笔架山省级自然保护区
Bijiashan Provincial Nature Reserve
in Lianshan3 29 904~975 22.4~24.0 17.3~18.4 19.8~20.5 18 2.2 土样采集与理化性质测定
在标准样地中沿对角线布设取样点3个，每个样点分别在0~10、10~20和20~40 cm土层取样，采用多点取样，同层次混合，取土壤各层次混合样品，带回室内风干，拣去石砾、植物根系和碎屑，过2.000和0.149 mm土壤筛后，储藏于密封的塑料自封袋中，用于土壤化学性质(pH、有机质、全N、全P、全K、水解N、速效P和速效K等)测定[11]。同时，采用“环刀法”取各层次的原状土壤，带回实验室，测定其水分物理性质(土壤容重、自然含水量、最大持水量、最小持水量、田间持水量、总孔隙度、非毛管孔隙度、毛管孔隙度等)[11]。
2.3 林分数据计算
2.3.1 径阶平均胸径(
$ {\bar D}_{i} $ )和树高($ {\bar H}_{i} $ )以2 cm为1个径阶，将所测量的杉木8~42 cm胸径划分为18个径阶，计算各径阶平均胸径：
$$ {\bar D}_{i}=\sqrt{\frac{1}{N_i}\sum {{d}_{i}}^{2}}, $$ (1) 式中，di为第i径阶中实测胸径值，Ni为第i径阶总株数。
运用Excel软件对胸径和树高进行拟合，径阶平均树高最佳回归方程为：
$${\bar H_i} = 2.884\;8\bar D_i^{0.561\;0},\;\left( {R = 0.56} \right)\text{。}$$ (2) 2.3.2 单株材积(Vs)与径阶材积(Vi)
广东省杉木单株材积公式为：
$${V_{\rm{s}}} = 0.000\;058\;770\;42 \times {d^{1.969\;983\;1}} \times {h^{0.896\;461\;57}},$$ (3) 式中，d为单株胸径，h为单株树高。
径阶材积公式为：
$$\begin{split}{V_i} =& 0.000\;058\;770\;42 \times {\bar D_i}^{1.969\;983\;1} \times\\& {\bar H_i}^{0.896\;461\;57} \times {N_i}\text{。}\end{split}$$ (4) 2.3.3 林分蓄积量(V)
林分蓄积量为标准样地中各径阶活立木蓄积量之和，公式为：
$$V = \sum {V_i}\text{。}$$ (5) 2.3.4 大径材蓄积量(Vl)
参照《木材标准实用手册》[12]，根据过渡径阶大径材材种株数比例(Ai/%)，统计各标准样地大径材蓄积量之和，公式为：
$$\begin{split}{V_{\rm{l}}} = &\sum 0.000\;058\;770\;42 \times {\bar D_i}^{1.969\;983\;1} \times\\& {\bar H_i}^{0.896\;461\;57} \times {N_i} \times {A_i}\text{。}\end{split}$$ (6) 2.3.5 大径材出材量(Ol)
大径材出材量(规格材原条)公式[13]采用：
$$\begin{split}{O_{\rm{l}}} =& \sum 3.602\;437\;58 \times {10^{ - 5}} \times\\& {\bar D_i}^{1.947\;520\;76} \times {\bar H_i}^{1.007\;937\;69} \times {N_i} \times {A_i}\text{。}\end{split}$$ (7) 2.3.6 大径材株数比例(ND≥26)
统计各标准样地大径材(胸径D≥26 cm)株数占样地总株数的百分比。
2.3.7 大径材出材率(Yl)
统计大径材出材量占林分蓄积量的百分比，公式为：
$${Y_{\rm{l}}} = {O_{\rm{l}}}/V \times 100{\text{%}} \text{。}$$ (8) 2.4 评价指数
2.4.1 立地指数(SI)
根据标准样地内优势木平均高和林分年龄，查广东省杉木立地指数表[14]，确定各标准样地的立地指数。
2.4.2 土壤综合肥力指数(I)
采用隶属函数法对土壤8个物理性质和8个化学性质测定指标值进行标准化处理，其中，容重采用降型函数公式(9)，其他指标采用升型函数公式(10)：
$$ Q\left({X}_{i}\right)=\frac{{X}_{i\mathrm{m}\mathrm{a}\mathrm{x}}-{X}_{ij}}{{X}_{i\mathrm{m}\mathrm{a}\mathrm{x}}-{X}_{i\mathrm{m}\mathrm{i}\mathrm{n}}},$$ (9) $$ Q\left({X}_{i}\right)=\frac{{X}_{ij}-{X}_{i\mathrm{m}\mathrm{i}\mathrm{n}}}{{X}_{i\mathrm{m}\mathrm{a}\mathrm{x}}-{X}_{i\mathrm{m}\mathrm{i}\mathrm{n}}}, $$ (10) 式中，Q
$ \left({X}_{i}\right) $ 为各因子的隶属度值，$ {X}_{ij} $ 为第i项第j个指标的测定值，$ {X}_{i\mathrm{m}\mathrm{a}\mathrm{x}} $ 和$ {X}_{i\;\mathrm{m}\mathrm{i}\mathrm{n}} $ 分别为第$ i $ 项指标中的最大值和最小值。对隶属度值进行主成分分析，得到各指标的公因子方差，再计算各指标的权重，公式如下：
$${\omega _i} = {C_i}\Bigg/\sum\limits_{i = 1}^n {{C_i}} ,$$ (11) 式中，
$ {\omega }_{i} $ 为第$ i $ 项指标的权重，$ {C}_{i} $ 为第$ i $ 项指标的公因子方差，n为土壤理化性质测定指标的个数。应用模糊集加权综合法，计算土壤综合肥力指数(I)[15]，公式如下：
$$ I=\sum\limits_{i=1}^{n}\left[{\omega _i} \cdot Q\left({X}_{i}\right)\right], $$ (12) 式中，
${\omega _i}$ 为各单项指标权重。本研究先分别计算各标准样地3个层次土壤混合样的土壤综合肥力指数，再以3个层次平均土壤综合肥力指数作为该样地的土壤综合肥力指数(I)。2.4.3 微地形指数(M)
为了使一些相关性较大、性质相近的定性因素能进行数量统计分析，对坡向、坡位、坡度、坡形等指征林地微地形要素的指标进行数量化分级，分级标准参考《杉木大径材培育技术规程：LY/T 2809—2017》[16]，并按其对杉木生长的影响程度提出“林地微地形指标数量化表”(表2)，进行数据标准转化处理及计算各指标权重，根据数量化理论Ⅲ[17]计算数量化加权值并作为微地形指数(M)。
表 2 林地微地形指标数量化表Table 2. Quantization table of forest land micro-topography index指标分级
Grade of index坡向
Slope aspect坡位
Slope position坡度/(°)
Slope degree坡形
Slope shape1 南 South 山顶 Hilltop ＞35 凸 Convex 2 东南、西南 Southeast, Southwest 上坡 Upslope 26~35 直 Straight 3 东、西 East, West 中坡 Middle slope 16~35 凹 Concave 4 北、东北、西北 North, Northeast, Northwest 下坡 Downslope 10~16 2.4.4 地理指数(G)
根据GPS定位，将海拔和纬度数据进行几何平均数转化[14]，得到地理指数(G)。
2.5 数据分析
采用Microsoft excel 对所有数据进行统计及分析，用SAS V8.1进行显著性检验，用SPSS 19.0进行数据标准化转换及主成分分析等。
3. 结果与分析
3.1 立地指数(SI)对杉木大径材生长的影响
36个标准样地的SI范围为16~22，SI为16、18、20和22的标准样地分别占总样地数的38.9%、36.1%、13.9%和11.1%。不同SI与大径材蓄积量(Vl)、大径材出材量(Ol)、大径材株数比例(ND≥26)和大径材出材率(Yl)等指标进行相关性分析，结果(图1)显示，随着SI的增加，大径材成材各相关指标均呈增加的趋势，相关性均达到极显著水平(P<0.01)，相关系数为0.37~0.44。根据SI与Vl的二次项回归方程(y = 1.51x2 − 17.20 x + 111.16)，当SI为16、18、20和22时，Vl的理论值分别为223.12、291.56、372.10和464.74 m3·hm−2，相邻立地指数的Vl增幅为30%左右；根据立地指数与Ol的指数方程(y = 18.97e0.13x)，当SI为16、18、20和22时，Ol的理论值分别为163.67、214.27、280.52和367.25 m3·hm−2，相邻立地指数的Ol增幅均为31.0%，与Vl增幅基本一致；根据SI与ND≥26的二次项回归方程(y = −0.20 x2 + 12.53x − 120.74)，当立地指数为16、18、20和22时， ND≥26的理论值分别为20.1%、39.4%、49.1%和57.3%，SI为18的林地生产大径材株数比SI为16的提高近100.0%，SI为20、22的比18、20的增加幅度逐渐减缓；根据立地指数与Yl的二次项回归方程(y = −0.38 x2 + 17.57x−145.64)，当SI为16、18、20和22时，Yl的理论值分别为38.6%、48.0%、54.4%和57.7%。
![]()
图 1 不同立地指数(SI)下的杉木大径材生长相关指标分析Figure 1. Analysis of growth indexes for large-diameter wood of Cunninghamia lanceolataunder different site index (SI)3.2 土壤综合肥力指数(I)对杉木大径材生长的影响
36个标准样地的土壤综合肥力指数(I)为0.11~0.57，数值越大，土壤综合肥力越高。不同土壤综合肥力指数与大径材蓄积量(Vl)、大径材出材量(Ol)、大径材株数比例(ND≥26)和大径材出材率(Yl)等指标进行相关性分析，结果(图2)显示，随着土壤综合肥力指数的增加，大径材成材各相关指标均呈增加的趋势，其中，I与Vl、Ol、ND≥26等指标的相关性达显著水平(P<0.05)，相关系数为0.24~0.39；I与Yl相关性未达显著水平(P>0.05)。
![]()
图 2 不同土壤综合肥力指数(I)下的杉木大径材生长相关指标分析Figure 2. Analysis of growth indexes for large-diameter wood of Cunninghamia lanceolataunder different soil integrated fertility indexes (I)对标准样地土壤8个物理性质和8个化学性质测定指标值进行主成分分析，结果(图3)显示，土壤孔隙度和土壤持水量等土壤物理性质占较大权重，是土壤综合肥力的关键指标，其次是土壤N和有机质含量。可见，土壤物理结构、土壤N和有机质含量是杉木大径材形成的重要土壤指标。
![]()
图 3 土壤理化性质指标权重雷达图Figure 3. Weighted radar chart of indexes for soil physical and chemical properties3.3 微地形指数(M)对杉木大径材生长的影响
将坡向、坡位、坡度、坡形等指征林地微地形要素的指标数量化、综合化，以微地形指数(M)作为综合指标考察林地微地形对杉木大径材生长的影响。36个标准样地的M为−1.14~1.25，不同微地形指数与大径材蓄积量(Vl)、大径材出材量(Ol)、大径材株数比例(ND≥26)和大径材出材率(Yl)等指标进行相关性分析，结果(图4)显示，随着微地形指数的增加，大径材成材各相关指标均呈增加的趋势，相关性均达到极显著水平(P<0.01)，相关系数为0.22~0.33。根据林地微地形指标数量化表(表2)，微地形指数最大的东北坡向、下坡位、坡度低于35°的凹坡有利于南岭山区杉木大径材的生长与形成。
![]()
图 4 不同微地形指数(M)下的杉木大径材生长相关指标分析Figure 4. Analysis of growth indexes for large-diameter wood of Cunninghamia lanceolataunder different micro-topography indexes (M)3.4 地理指数(G)对杉木大径材生长的影响
树木的生长也受海拔和纬度的影响[18]，因此本研究将纬度和海拔转换为地理指数，以研究其对杉木大径材生长的影响。从图5可知，地理指数(G)与大径材成材各指标相关性均不显著(P>0.05)，随地理指数的增加，整体呈现先下降再上升的U型曲线趋势。南岭山区地貌特征以海拔200~1000 m的低山丘陵为主[10]，纬度范围为N 24°~26°，本研究选择的36个杉木大径材标准样地的海拔均在200~1000 m范围内，基本囊括了南岭山区所有海拔梯度级，具有代表性，可见，海拔和纬度对南岭山区杉木大径材形成的影响不明显。
![]()
图 5 不同地理指数(G)下的杉木大径材生长相关指标分析Figure 5. Analysis of growth indexes for large-diameter wood of Cunninghamia lanceolata under different geographical index (G)3.5 林龄(A)对杉木大径材生长的影响
36个标准样地林龄为20~32年，参照南方杉木用材林龄组划分标准[19]，杉木近熟林(20~25年)占41.7%，杉木成熟林(26~32年)占58.3%。从图6可知，林龄(A)与大径材成材各指标相关性均不显著(P>0.05)，随林龄的增大整体呈现先上升再下降的倒U型曲线趋势。由于本研究对象为南岭山区具有杉木大径材的林分，所有标准样地的杉木林龄均在近熟至成熟期10年左右的区间内，因此林龄效应及对大径材生长的影响趋势体现不充分，而图6显示，杉木大径材等指标后期出现下降的趋势，主要是标准样地中30~32年生的林分立地指数多为16或18，且土壤综合肥力指数、微地形指数值相对较低，其大径材成材各指标均低于立地指数为18以上且林龄较小的林分。可见，南岭山区杉木大径材普遍在近熟林至成熟林时期(23~30年)就已形成，不具杉木大径材生产潜力的林分即使生长期更长，也难以形成大径材。
![]()
图 6 不同林龄(A)下的杉木大径材生长相关指标分析Figure 6. Analysis of growth indexes for large-diameter wood of Cunninghamia lanceolataunder different stand age (A)3.6 保留密度(D)对杉木大径材生长的影响
南岭山区杉木造林密度普遍为2 m×2 m(2500株·hm−2)，本研究调查的36个标准样地林分保留密度(D)为900~1700株·hm−2，保留密度与大径材蓄积量(Vl)、大径材出材量(Ol)、大径材株数比例(ND≥26)和大径材出材率(Yl)等指标进行相关性分析，结果(图7)显示，Vl和Ol均随保留密度的增加呈先缓慢上升后缓慢下降的平缓二次项曲线趋势，相关性达显著水平(P<0.05)；ND≥26和Yl随保留密度的增加先平缓下降然后下降幅度逐渐增大，但相关性不显著。保留密度为1100~1300株·hm−2时，大径材成材各指标均处于较大值区间。
![]()
图 7 不同保留密度(D)下的杉木大径材生长相关指标分析Figure 7. Analysis of growth indexes for large-diameter wood of Cunninghamia lanceolataunder different reserved density (D)3.7 杉木大径材出材量预测
根据36个标准样地调查统计数据，选取林分特征因子：立地指数(SI)，林龄(A)、保留密度(D)，以Schumacher收获模型[3, 20]为基础，构建杉木大径材出材量(Ol)预测模型。该模型基本形式为：
$$\ln \left( {{O_{\rm{l}}}} \right) = {\beta _0} + {\beta _1}{\rm{SI}} + {\beta _2}{A^{ - 1}} + {\beta _3}\ln D + {\beta _4}{\rm{S}}{{\rm{I}}^{ - 1}},$$ 通过SPSS软件求解参数β0~β4，获得模型如下：
$$\begin{split}\;&\\ {O_{\rm{l}}}{\rm{ = }}&{{\rm{e}}^{{\rm{5}}{\rm{.781 + 0}}{\rm{.004SI + 5}}{\rm{.004}}{A^{{\rm{ - 1}}}}{\rm{ + 0}}{\rm{.261}}\ln D{\rm{ - 45}}{\rm{.151S}}{{\rm{I}}^{{\rm{ - 1}}}}}},\\& \qquad\qquad \left( {R = 0.62,P < 0.01} \right)\text{。}\end{split}$$ 选取不同立地指数、不同林龄、不同保留密度样地的实测数据，将模型的大径材预测产量与实际产量进行比较(表3)，平均预测精度90%，因此，可利用该模型初步预测评估南岭山区20~30年生杉木林分的大径材出材量。
表 3 杉木大径材出材量实际值与预测值Table 3. The actual value and predicted value of large-diameter wood output of Cunninghamia lanceolata立地指数
Site index林龄/年
Stand age保留密度/(株·hm−2)
Reserved density大径材出材量/(m3·hm−2) Output of large-diameter wood 实际值 Actual value 预测值 Predicted value 16 25 1176.0 187.7254 158.2417 16 28 1075.5 139.8997 151.3147 18 25 1176.0 206.9185 216.5177 18 28 1626.0 258.4599 230.6247 20 23 1176.0 294.6004 283.1329 20 23 1300.5 330.7747 290.6677 22 26 1350.0 316.8053 351.4470 4. 讨论与结论
杉木是广东省南岭山区传统优势造林树种，但现有林分以中幼林为主，且经营模式单一、集约化程度低，与周边省份相比，杉木林分生产力与自然资源优势潜力尚未充分发挥[21]。本研究对广东省南岭山区范围内仅存为数不多、保存完好的杉木大径材林分进行调查，以空间代替时间，基于现有样本数据挖掘大径材形成的主导生态环境因子和林分控制因子并探究其影响规律，以期为广东省杉木大径材培育和经营提供理论指导。
立地指数是指林分中在基准年龄时自由生长的优势木平均高度值，可考察林地生产某一类型林木的生产潜力，常作为指导人工林培育最直观的指标[16]。惠刚盈等[9]认为在材种形成过程中，无论有何种保留密度，只要立地指数过低，在目的材种的培育期限内，就无法产生需要的材种。立地指数越高，林木生长速度越快，林分内活立木间的竞争也越大，成熟时期林分内保留木的株数也相对越低，大径材出材率也相对越高[22-23]。现有研究普遍以立地指数16、18作为培育不同杉木目标材种的划分标准[4, 6, 9, 24-25]。在南岭山区，杉木大径材成材各相关指标也是随立地指数的增加均呈显著增加的趋势(P<0.01)，在立地指数16~22范围内，相邻立地指数的大径材蓄积量和大径材出材量增幅为30%左右，立地指数18的林地生产大径材株数比立地指数16的提高近100%，大径材出材率超过40%。可见，立地指数是南岭山区杉木大径材材种形成的重要限制因子，立地指数18及以上适宜培育杉木大径材。
土壤肥力是土壤的本质属性，对林木的长势和品质有直接影响[26-27]。土壤肥力评价能够揭示土壤利用与功能有关的土壤内在属性之间的异同，科学地确定森林土壤肥力指标并进行评价，对立地生产力和多目标森林经营的研究也具有重要价值[28]。本研究采用土壤综合肥力指数作为土壤整体响应特征考察不同土壤对杉木大径材生长的影响。研究发现，大径材成材各相关指标随土壤综合肥力指数的增加均呈增加的趋势，与大径材蓄积量、大径材出材量、大径材株数比例等指标的相关性达显著水平(P<0.05)。据主成分分析结果，土壤综合肥力指数的权重贡献率较大的是土壤物理性质、土壤N含量和有机质含量。土壤物理性质密切影响根系的生长、最大扎根深度和水分的有效吸收[28]，是与南岭山区杉木大径材的形成紧密相关的土壤指标，在进行杉木大径材经营时，应首先选择土壤孔隙度和持水量较高的立地。土壤有机质在提供林木生长所需要养分的同时，能有效改良土壤结构，使土壤变得疏松，促进土壤形成水稳性团聚体，增加土壤持水和保肥能力，有利于根系的生长、土壤微生物和土壤动物活动[29-30]；土壤全氮和水解氮绝大部分是以有机态的形式积累和贮存在土壤有机质中，与土壤有机质相关性更密切，相关系数分别为9.07和0.87[31]，因此，土壤有机质是与杉木大径材成材密切相关的又一重要土壤指标。
坡向、坡位、坡度、坡形等是形成林地微地形的主要因素，南岭山区杉木大径材成材各相关指标随着微地形指数的增加，均呈显著增加的趋势(P<0.01)。不同坡向、坡位、坡度、坡形等微地形要素会形成不同的山区小气候，东北坡向等阴坡接受太阳直接辐射较少，主要以地面散射为主，日温差较小、林分相对湿度高[32]；坡位和坡度对水土流失和集聚都有影响，坡形则与风速和蒸腾强度有关[32]，都直接影响到土层的厚度和含水量，故微地形指数最大的东北坡向、下坡位、坡度低于35°凹坡的土层较深厚，排水良好，所形成的小气候最有利于净光合产物长期积累。杉木是亚热带树种，温暖湿润、多雾静风的气候环境以及深厚、排水良好的土壤条件有利于杉木生长，而水湿条件对杉木径向生长的影响非常大，因此杉木大径材的形成对小气候的要求十分严格。南岭山区属于杉木产区分区中产力水平最高的中带东区[16]，气候温暖而不炎热，太阳辐射和水湿条件等微地形因素对杉木大径材形成是比温度更重要的限制因子。
林分密度是通过影响林分直径结构进而影响林分材种结构的[33-34]，保持合理的林分密度也是杉木大径材培育的关键。林分密度与林分直径呈负相关关系已达成共识[35]，关于如何对林分密度进行合理的调控，从而获得较高大径材出材量和出材率，并平衡单位面积整体木材产量，已有很多学者对此展开了广泛的研究。佟金权[24]、刘春燕等[36]认为保留密度为1200~1300 株·hm−2时，林分整体胸径最大，较为适合杉木大径材的培育；惠刚盈等[9]、涂育合等[4]根据不同立地指数提出了不同的保留密度范围，认为立地指数16~18、保留密度1050~1350株·hm−2时适宜培育中大径材，立地指数越大保留密度应越低。本研究发现，大径材蓄积量和大径材出材量均随保留密度的增加呈先缓慢上升后缓慢下降的平缓二次项曲线趋势，相关性达显著水平(P<0.01)；大径材株数比例和大径材出材率随保留密度的增加，先平缓下降，然后下降幅度逐渐增大，但相关性不显著。南岭山区18以上立地指数，保留密度在1100~1300株·hm−2范围内时，大径材成材各指标均达到较大值区间。
本研究以Schumacher收获模型为基础，构建杉木大径材出材量(Ol)预测模型：Ol=
${{\rm{e}}^{{\rm{5}}{\rm{.781 + 0}}{\rm{.004SI + 5}}{\rm{.004}}{A^{{\rm{ - 1}}}}{\rm{ + 0}}{\rm{.261}}\ln D{\rm{ - 45}}{\rm{.151S}}{{\rm{I}}^{{\rm{ - 1}}}}}}$ (R=0.62，P<0.01)，可利用该模型初步预测评估南岭山区20~30年生杉木林分的大径材收获量，据此指导南岭山区杉木大径材培育和经营。 -
![]()
图 1 不同浓度NaCl溶液处理下‘航聚香丝苗’干种子的萌发表型
Figure 1. The germination phenotypes of ‘Hangjuxiangsimiao’ dried seeds under different NaCl concentrations
![]()
图 2 不同浓度NaCl溶液处理下‘航聚香丝苗’干种子的萌发情况
Figure 2. The germination state of ‘Hangjuxiangsimiao’ dried seeds under different NaCl concentrations
![]()
图 4 M3代耐盐突变体盐处理7 d后的发芽率、发芽势和芽长对比
“**”表示突变体与野生型在0.01水平差异显著(t检验)；其中，红色柱子表示极耐盐突变体，蓝色柱子表示中度耐盐突变体。
Figure 4. Comparison of germination rate, germination potential and bud length of M3 generation salt-tolerant mutants after 7 days of salt treatment
“**” indicates that the mutants differed from the wild-type respectively at 0.01 levels of significant difference (t test); Red bars represent highly salt-tolerant mutants, while blue bars indicate mutants with moderate salt tolerance.
![]()
图 5 盐胁迫对2个代表性耐盐突变体抗逆相关生理指标的影响
各小图中不同小写字母表示差异显著(P<0.05，Duncan’s 法)。
Figure 5. Effect of salt stress on stress-related physiological indexes of two representative salt-tolerant mutants
Different lowercase letters in each figure indicate significant differences (P<0.05, Duncan’s method).
![]()
图 6 部分SSR引物对野生型材料及突变体的电泳图谱
Figure 6. Electrophoretic profile of partial SSR primers in the WT and the mutants
1~13：WT、9-8、20-62、42-113、20-38、15-6、16-65、19-49、20-15、20-19、20-49、22-7、19-18；红色箭头所指为差异位点。The red arrow indicates the differential loci.
表 1 空间诱变和重离子诱变的水稻突变体M3代株系的主要农艺性状1)
Table 1 The main agronomic characters of the rice mutants derived from space mutation and heavy ion mutation
株系
Line株高/cm
Plant
height剑叶长/cm
Flag-leaf
length剑叶宽/cm
Flag-leaf
width单株穗质量/g
Panicle weight
per plant穗长/cm
Panicle
length十粒长/cm
10-grain
length十粒宽/cm
10-grain
width千粒质量/g
1000-grain
weightWT 112.67±2.52b 34.60±2.04a 1.83±0.21a 22.30±2.82a 24.16±1.57a 11.13±0.16a 1.87±0.26a 21.16±1.50a 9-8 110.14±1.11b 33.75±2.02a 1.85±0.05a 22.24±0.15a 24.36±0.05a 10.59±0.02b 1.87±0.02a 18.52±0.43b 15-6 117.89±1.45a 34.05±0.32a 1.78±0.06a 21.81±1.26a 24.22±0.23a 11.08±0.07a 1.88±0.05a 21.27±0.32a 16-65 113.26±1.56b 33.71±2.18a 1.79±0.14a 22.72±1.31a 23.98±1.32a 11.12±0.09a 1.90±0.04a 20.89±1.03a 22-7 112.03±1.25b 34.19±1.65a 1.82±0.02a 22.99±1.06a 24.10±1.04a 11.13±0.06a 1.89±0.12a 21.25±0.14a 19-18 109.79±1.80b 34.04±1.93a 1.79±0.11a 22.57±1.47a 23.92±1.45a 10.71±0.12b 1.88±0.13a 21.03±0.36a 19-49 112.15±0.32b 33.92±1.21a 1.19±0.12a 22.40±2.03a 24.05±0.96a 11.09±0.11a 1.86±0.18a 20.84±1.21a 20-19 112.98±1.23b 33.92±2.08a 1.78±0.15a 22.16±1.65a 24.31±0.32a 11.14±0.03a 1.87±0.02a 20.96±1.07a 20-15 117.65±1.75a 35.24±1.15a 1.80±0.12a 22.90±1.54a 24.14±1.17a 11.12±0.12a 1.87±0.06a 21.08±0.35a 20-38 110.42±2.11b 34.13±1.33a 1.81±0.04a 22.62±1.63a 24.27±1.26a 11.10±0.11a 1.87±0.11a 21.23±0.22a 20-49 112.67±2.30b 33.68±1.82a 1.83±0.19a 23.14±1.24a 23.89±1.03a 11.03±0.21a 1.88±0.07a 21.35±0.53a 20-62 110.67±2.52b 34.15±1.63a 1.82±0.15a 22.09±2.12a 24.03±1.12a 11.19±0.03a 1.87±0.16a 21.57±0.28a 42-113 110.14±2.11b 30.50±1.29b 1.82±0.04a 22.36±1.36a 23.79±1.39a 11.01±0.22a 1.89±0.08a 21.24±0.55a 1) WT：航聚香丝苗；9-8为重离子诱变后代；其他株系为空间诱变后代；表中数据为平均值±标准误，n=5，同列数据后不同小写字母表示不同株系之间差异显著(P<0.05，Duncan’s法)。
1) WT: Hangjuxiangsimiao; 9-8 represents heavy ion-induced mutant; Other strains represent space-induced mutants; The data in the table are mean value ± standard error, n=5, different lowercase letters after the data in the same column indicate significant differences among different lines (P<0.05，Duncan’s method).
下载: 导出CSV
-
[1] 杨建昌, 李超卿, 江贻. 稻米氨基酸含量和组分及其调控[J]. 作物学报, 2022, 48(5): 1037-1050. [2] PONCE K, GUO L, LENG Y, et al. Advances in sensing, response and regulation mechanism of salt tolerance in rice[J]. International Journal of Molecular Sciences, 2021, 22(5): 2254. doi: 10.3390/ijms22052254
[3] CHENG X, HE Q, TANG S, et al. The miR172/IDS1 signaling module confers salt tolerance through maintaining ROS homeostasis in cereal crops[J]. New Phytologist, 2021, 230(3): 1017-1033. doi: 10.1111/nph.17211
[4] YOU J, ZONG W, HU H, et al. A STRESS-RESPONSIVE NAC1-regulated protein phosphatase gene rice protein Phosphatase18 modulates drought and oxidative stress tolerance through abscisic acid-independent reactive oxygen species scavenging in rice[J]. Plant Physiology, 2014, 166(4): 2100-2114. doi: 10.1104/pp.114.251116
[5] ZHANG J, FLOWERS T J, WANG S. Mechanisms of sodium uptake by roots of higher plants[J]. Plant and Soil, 2010, 326(1): 45-60.
[6] 黄洁, 白志刚, 钟楚, 等. 水稻耐盐生理及分子调节机制[J]. 核农学报, 2020, 34(6): 1359-1367. doi: 10.11869/j.issn.100-8551.2020.06.1359 [7] 曲梦宇, 许惠滨, 陈静, 等. 水稻耐盐分子机制与分子改良[J]. 植物遗传资源学报, 2022, 23(3): 644-653. [8] 姚栋萍, 吴俊, 胡忠孝, 等. 水稻耐盐碱的生理机制及育种策略[J]. 杂交水稻, 2019, 34(4): 1-7. [9] CYRANOSKI D. Satellite will probe mutating seeds in space[J]. Nature, 2001, 410(6831): 857-858.
[10] 陈淳, 严贤诚, 罗文龙, 等. 水稻空间诱变与重离子诱变生物学效应及突变体定向筛选[J]. 华南农业大学学报, 2021, 42(1): 49-60. doi: 10.7671/j.issn.1001-411X.202004012 [11] 俞法明, 严文潮, 毛雪琴, 等. 利用空间诱变技术进行早籼稻新品种的改良[J]. 核农学报, 2014, 28(6): 949-954. doi: 10.11869/j.issn.100-8551.2014.06.0949 [12] 王伟荣, 孙超才, 蒋美艳, 等. 甘蓝型油菜品种航天诱变后代性状分析[J]. 核农学报, 2015, 29(2): 215-220. doi: 10.11869/j.issn.100-8551.2015.02.0215 [13] 刘录祥, 赵林姝, 郭会君. 作物航天育种研究现状与展望[J]. 中国农业科技导报, 2007, 9(2): 26-29. doi: 10.3969/j.issn.1008-0864.2007.02.005 [14] 陈志强, 郭涛, 刘永柱, 等. 水稻航天育种研究进展与展望[J]. 华南农业大学学报, 2009, 30(1): 1-5. doi: 10.3969/j.issn.1001-411X.2009.01.001 [15] 周利斌, 李文建, 曲颖, 等. 重离子束辐照育种研究进展及发展趋势[J]. 原子核物理评论, 2008, 25(2): 165-170. doi: 10.11804/NuclPhysRev.25.02.165 [16] FU H, LI Y, SHU Q. A revisit of mutation induction by gamma rays in rice (Oryza sativa L. ): Implications of microsatellite markers for quality control[J]. Molecular Breeding, 2008, 22(2): 281-288. doi: 10.1007/s11032-008-9173-7
[17] 李瑞清, 武立权, 舒庆尧, 等. 一个新的水稻白化转绿突变体G9的特性研究[J]. 核农学报, 2010, 24(5): 881-886. doi: 10.11869/hnxb.2010.05.0881 [18] 曾晓珊, 彭丹, 石媛媛, 等. 利用SSR标记构建水稻核心亲本指纹图谱[J]. 作物研究, 2016, 30(5): 481-486. [19] 许绍斌, 陶玉芬, 杨昭庆, 等. 简单快速的DNA银染和胶保存方法[J]. 遗传, 2002, 24(3): 335-336. doi: 10.3321/j.issn:0253-9772.2002.03.027 [20] MANSOUR F M M, HASSAN F A. How salt stress-responsive proteins regulate plant adaptation to saline conditions?[J]. Plant Molecular Biology, 2022, 108(3): 175-224. doi: 10.1007/s11103-021-01232-x
[21] ZHAO S, ZHANG Q, LIU M, et al. Regulation of plant responses to salt stress[J]. International Journal of Molecular Sciences, 2021, 22(9): 4609. doi: 10.3390/ijms22094609
[22] CHOUDHURY F K, RIVERO R M, BLUMWALD E, et al. Reactive oxygen species, abiotic stress and stress combination[J]. The Plant Journal, 2017, 90(5): 856-867. doi: 10.1111/tpj.13299
[23] 中国水稻研究所, 农业部科技发展中心. 水稻品种鉴定技术规程 SSR标记法: NY/T 1433—2014[S]. 北京: 中国农业出版社, 2014. [24] 符德保, 李燕, 肖景华, 等. 中国水稻基因组学研究历史及现状[J]. 生命科学, 2016, 28(10): 1113-1121. [25] 郑静. 一份水稻白化突变体的遗传分析及其基因定位[D]. 雅安: 四川农业大学, 2008. [26] 魏溪涓. 水稻闭花授粉突变体fod的突变基因克隆与功能分析[D]. 金华: 浙江师范大学, 2015. [27] 钟琪. 组织培养技术筛选植物耐盐突变体的研究进展[J]. 安徽农业科学, 2016, 44(5): 145-148. doi: 10.3969/j.issn.0517-6611.2016.05.053 [28] 颜佳倩, 顾逸彪, 薛张逸, 等. 耐盐性不同水稻品种对盐胁迫的响应差异及其机制[J]. 作物学报, 2022, 48(6): 1463-1475. [29] LIU J, FU C, LI G, et al. ROS homeostasis and plant salt tolerance: Plant nanobiotechnology updates[J]. Sustainability, 2021, 13(6): 3552. doi: 10.3390/su13063552
[30] LIU C, MAO B, YUAN D, et al. Salt tolerance in rice: Physiological responses and molecular mechanisms[J]. The Crop Journal, 2022, 10(1): 13-25. doi: 10.1016/j.cj.2021.02.010
[31] VAN ZELM E, ZHANG Y, TESTERINK C. Salt tolerance mechanisms of plants[J]. Annual Review of Plant Biology, 2020, 71: 403-433. doi: 10.1146/annurev-arplant-050718-100005
[32] FOYER C, NOCTOR G. Redox homeostasis and antioxidant signaling: A metabolic interface between stress perception and physiological responses[J]. The Plant Cell, 2005, 17(7): 1866-1875. doi: 10.1105/tpc.105.033589
[33] KHARE T, SRIVASTAVA A K, SUPRASANNA P, et al. Individual and additive stress impacts of Na+ and Cl‾ on proline metabolism and nitrosative responses in rice[J]. Plant Physiology and Biochemistry, 2020, 152: 44-52. doi: 10.1016/j.plaphy.2020.04.028
[34] 温贤芳, 张龙, 戴维序, 等. 天地结合开展我国空间诱变育种研究[J]. 核农学报, 2004, 18(4): 241-246. doi: 10.3969/j.issn.1000-8551.2004.04.001 [35] 陈立凯, 黄明, 刘永柱, 等. 水稻开颖半不育突变体的观察、遗传分析和基因定位[J]. 中国农业科学, 2016, 49(1): 1-13. doi: 10.3864/j.issn.0578-1752.2016.01.001 [36] CALDECOTT R S, STEVENS H, ROBERTS B J. Stem rust resistant variants in irradiated populations mutations or field hybrids? [J]. Agronomy Journal, 1959, 51(7): 401-403. doi: 10.2134/agronj1959.00021962005100070011x
[37] 张福彦, 张建伟, 程仲杰, 等. 航天诱变技术在小麦育种上的应用[J]. 核农学报, 2019, 33(2): 262-269. doi: 10.11869/j.issn.100-8551.2019.02.0262 -
期刊类型引用(3)
1. 罗锡文，谷秀艳，胡炼，赵润茂，岳孟东，何杰，黄培奎，汪沛. 大田无人化智慧农场农田边界识别技术研究现状与展望. 农业机械学报. 2025(02): 1-18 . 
百度学术
2. 冯国练. 无人驾驶农机路径跟踪控制方法探索研究. 现代农机. 2025(02): 68-70 . 百度学术
3. 寇浩，曲双为. 基于AHP-FCE的花椒采摘机器人设计. 工业设计. 2024(06): 146-150 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 138
- HTML全文浏览量: 21
- PDF下载量: 56
- 被引次数: 4
目录
1. 研究区域自然概况
2. 研究方法
2.1 标准样地设置与调查
2.2 土样采集与理化性质测定
2.3 林分数据计算
2.3.1 径阶平均胸径(
$ {\bar D}_{i} $ )和树高($ {\bar H}_{i} $ )2.3.2 单株材积(Vs)与径阶材积(Vi)
2.3.3 林分蓄积量(V)
2.3.4 大径材蓄积量(Vl)
2.3.5 大径材出材量(Ol)
2.3.6 大径材株数比例(ND≥26)
2.3.7 大径材出材率(Yl)
2.4 评价指数
2.4.1 立地指数(SI)
2.4.2 土壤综合肥力指数(I)
2.4.3 微地形指数(M)
2.4.4 地理指数(G)
2.5 数据分析
3. 结果与分析
3.1 立地指数(SI)对杉木大径材生长的影响
3.2 土壤综合肥力指数(I)对杉木大径材生长的影响
3.3 微地形指数(M)对杉木大径材生长的影响
3.4 地理指数(G)对杉木大径材生长的影响
3.5 林龄(A)对杉木大径材生长的影响
3.6 保留密度(D)对杉木大径材生长的影响
3.7 杉木大径材出材量预测
4. 讨论与结论
