Factors affecting individual growth of Pinus koraiensis and Tilia amurensis in coniferous and broadleaved mixed forests in Jiaohe, Jilin Province of northeastern China
-
摘要:目的
以吉林蛟河针阔混交林优势树种红松、紫椴为研究对象，探究叶性状的种间差异，不同叶性状间的相关性，以及生物（邻域竞争、初始胸径、叶片性状）与非生物因素（地形因子）对树木生长速率的影响。
方法利用Welch-ANOVA方差分析量化红松、紫椴叶性状的种间差异；利用主成分分析及斯皮尔曼相关系数评估性状间的相关性；利用普通最小二乘模型计算不同因素对红松、紫椴生长速率的相对影响。
结果红松的光饱和净光合速率、暗呼吸速率显著高于紫椴；紫椴的光饱和气孔导度、比叶面积显著高于红松。红松的光饱和净光合速率与光饱和气孔导度、比叶面积呈正相关；紫椴的光饱和净光合速率与光饱和气孔导度呈正相关，与比叶面积呈负相关。初始胸径、邻域拥挤指数、坡度对红松年均胸径生长量有显著负效应。初始胸径对红松年均生物量生长量有显著正效应，邻域拥挤指数、坡度有显著负效应。初始胸径对紫椴年均胸径生长量有显著负效应，对紫椴年均生物量生长量有显著正效应。海拔对紫椴年均胸径生长量、年均生物量生长量均有显著负效应。
结论红松、紫椴叶性状表现出显著的种间差异，两物种叶片性状的相关性呈现不同规律。紫椴的生长受初始胸径、海拔的影响。初始胸径越大，年均胸径生长量越低，而年均生物量生长量越高；海拔越高，年均胸径生长量与年均生物量生长量越低。红松的生长受初始胸径、邻域拥挤指数和坡度的影响。初始胸径越大，年均胸径生长量越低，而年均生物量生长量越高；邻域拥挤指数、坡度越大，年均胸径生长量与年均生物量生长量越低。
Abstract:ObjectivePinus koraiensis and Tilia amurensis, the dominant tree species in coniferous and broadleaved mixed forests in Jiaohe, Jilin Province of northeastern China, were taken as the research object. We investigated the interspecific differences in leaf traits, correlations among leaf traits, and the effects of biotic factors (neighborhood competition, initial DBH, leaf traits) and abiotic factors (topographical factor) on growth rate of trees.
MethodWelch ANOVA was used to quantify the interspecific differences in leaf traits of Pinus koraiensis and Tilia amurensis. Principal component analysis (PCA) and Spearman correlation coefficient were used to evaluate the correlations among leaf traits. The relative effects of different factors on growth rate of Pinus koraiensis and Tilia amurensis were calculated by ordinary least squares model (OLS).
ResultThe light-saturated photosynthetic rate and leaf dark respiration of Pinus koraiensis were significantly higher than those of Tilia amurensis; light-saturated stomatal conductance and specific leaf area of Tilia amurensis were significantly higher than those of Pinus koraiensis. The light-saturated photosynthetic rate of Pinus koraiensis was positively correlated with light-saturated stomatal conductance and specific leaf area. The light-saturated photosynthetic rate of Tilia amurensis was positively correlated with light-saturated stomatal conductance and negatively correlated with specific leaf area. Initial DBH, neighborhood crowding index and slope had significant negative effects on average annual DBH growth of Pinus koraiensis. Initial DBH had significantly positive effect on average annual biomass growth of Pinus koraiensis, neighborhood crowding index and slope had significantly negative effect. Initial DBH had significantly negative effect on average annual DBH growth and positive effect on average annual biomass growth of Tilia amurensis. Altitude had significant negative effects on the average annual DBH growth and annual biomass growth of Tilia amurensis.
ConclusionThe leaf traits of Pinus koraiensis and Tilia amurensis have significantly interspecific differences, the correlation of leaf traits of two species has different patterns. The growth of Tilia amurensis is affected by the initial DBH and elevation. The larger the initial DBH is, the lower the average annual DBH growth is, and the higher the average annual biomass growth is. The higher the elevation is, the lower the average annual DBH growth and biomass growth is. The growth of Pinus koraiensis is affected by the initial DBH, neighborhood crowding index and slope. The larger the initial DBH is, the lower the average annual DBH growth is, while the higher the average annual biomass growth is. The larger the neighborhood crowding index and slope is, the lower the average annual DBH growth and average annual biomass growth are.
-
多倍化是一种普遍的自然现象，是植物进化和物种形成的重要推动力[1]。利用2n配子的有性多倍化是植物多倍化的主要途径之一。但自然界中天然2n配子发生频率一般较低，难以应用于植物多倍体育种实践。为了提高2n配子的发生频率，进而提高多倍体得率，育种学家们常常利用物理或化学方法进行2n配子人工诱导。
秋水仙碱作为一种重要的化学诱变剂，已成功应用于诱导林木大孢子、胚囊、合子以及体细胞染色体加倍、选育多倍体的育种实践，获得了数以百计的杨树（Populus spp.）[2-8]、刺槐（Robinia pseudoacacia）[9]、杜仲（Eucommia ulmoides）[10]、尾叶桉（Eucalyptus urophylla）[11]等林木三倍体与四倍体新种质。秋水仙碱也曾用于诱导杨树[12-13]、杜仲[14]以及桉树（Eucalyptus）[15]等树种的花粉染色体加倍研究，获得了20.0% ~ 88.0%的2n花粉。但是，秋水仙碱价格昂贵，且具有剧毒性，容易导致加倍处理操作人员中毒。
高温是一种有效的物理诱导方法，目前已成功应用于诱导杨树[16-19]、杜仲[20]、桉树[21]等树种的大孢子、胚囊以及合子染色体加倍、选育多倍体的育种实践，获得了大量的毛白杨（Populus tomentosa）[16-17]、‘折引3号’杨（P. pseudo-simonⅡ × P. nigra ‘Zheyin3#’）[18]、毛新杨（P. tomentosa × P. bolleana）[19]以及银腺杨（P. alba × P. glandulosa）[19]等杂种三倍体。高温也曾用于诱导毛新杨[22]、响叶杨（P. adenopoda） [23]、银白杨（P. alba）[24]及银灰杨（P. canescens）[25]的花粉染色体加倍研究，确定了高温诱导2n花粉的有效时期为终变期附近，且获得了60.0% ~ 80.0%的高频率2n花粉。但有研究表明：高温诱导获得的2n配子均为SDR型[26]，其遗传亲本的杂合性仅为FDR型2n配子的一半，严重限制了高温在多倍体育种中的应用。
赤霉素（gibberellins，GA3）是一类重要的内源植物激素，具有调控细胞分化、种子萌发、叶片扩展、毛状体发育、花粉成熟以及开花转换等方面的功能[27-30]。赤霉素还是一种非生物逆境应激信号分子，水稻（Oryza sativa）经历冷胁迫时，通过降低内源赤霉素水平以抑制根系生长，并阻滞花粉发育[31]。此外，赤霉素能够调节拟南芥（Arabidopsis thaliana）小孢子发生和花粉发育[32]，具有诱导拟南芥产生2n花粉的功能[33]。但赤霉素是否能够有效地诱导杨树产生2n花粉尚不清楚。
减数分裂是植物有性生殖过程中配子形成的必要阶段。在植物花粉母细胞减数分裂过程中，微管骨架参与纺锤体和成膜体的形成，以及核周和胞质微管网络的构建，从而决定同源染色体和姐妹染色单体分离以及细胞质的正常分裂[34-36]。有研究表明：植物天然2n花粉的产生与花粉母细胞减数分裂以及胞质分裂的异常密切相关[37]，而这些异常的现象均受细胞内微管网络系统的严格调控。研究赤霉素喷洒处理后新疆杨（Populus bolleana）花粉母细胞微管骨架的异常变化，可为赤霉素喷洒处理诱导新疆杨2n花粉的发生机制提供有力的细胞学依据。
新疆杨为杨柳科（Salicaceae）杨属（Populus）白杨派乔木，是银白杨的变种，广泛种植于新疆、内蒙、北京、陕西、辽宁等省市。新疆杨是分布区内生态防护林建设和城市绿化的主要树种之一，在防风固沙和涵养水源方面具有重要作用。但有关新疆杨多倍体育种研究的报道甚少。本研究以新疆杨花粉母细胞为研究对象，在掌握新疆杨花粉母细胞减数分裂规律的基础上，研究赤霉素喷洒处理诱导新疆杨产生2n花粉的可行性及其对花粉母细胞微管骨架的影响，旨在为开发诱导效果好且价格低廉的新型化学诱变剂、构建高效的林木多倍体育种技术体系奠定重要的理论和实践基础。
1. 材料和方法
1.1 试验材料
每年冬季落叶后，从内蒙古自治区包头市国家杨树良种基地采集新疆杨雄花枝，塑料布包裹后运输至北京林业大学温室待用。
1.2 研究方法
1.2.1 新疆杨花粉母细胞减数分裂进程观察
新疆杨雄花枝水培于10 ~ 20 ℃的温室，每隔2 ~ 3 h摘取2 ~ 3个花芽，去除鳞片后，将其放于5 mL的离心管内，利用新配的卡诺固定液进行固定。在4 ℃冰箱固定24 h后，摘取花药，利用醋酸洋红染色法观察花粉母细胞减数分裂进程。染色体制片具体操作参照鲁敏等[23]的方法。在Olympus BX-51光学显微镜下观察花粉母细胞减数分裂进程[38]，并运用Olympus DP70数码摄像系统进行照相。
1.2.2 新疆杨雄花芽的赤霉素喷洒处理
以育性良好的新疆杨雄花芽为材料，于花粉母细胞、第一次减数分裂和第二次减数分裂时期，使用浓度为50、100、200 μmol/L的赤霉素对花芽分别进行3、5、7次喷洒处理。当观察到大部分新疆杨花粉母细胞发育至第二次减数分裂时，使用水对花芽分别进行3、5、7次喷洒处理作为对照（表1）。处理后的雄花芽继续在温室中水培至花粉成熟，按照不同处理收集花粉，放置于1.5 mL的离心管中，添加适量的变色硅胶，密封保存于–20 ℃的冰箱。
表 1 新疆杨雄花芽的赤霉素喷洒处理Table 1. GA3 spraying treatment of flower buds in P. bolleana减数分裂时期
Meiotic stage赤霉素浓度
Concentration of
GA3/(μmol· L−1)喷洒次数
Spraying time花粉母细胞时期
Pollen mother cell (PMC) stage50 3 5 7 100 3 5 7 200 3 5 7 第一次减数分裂
Meiosis Ⅰ50 3 5 7 100 3 5 7 200 3 5 7 第二次减数分裂
Meiosis Ⅱ50 3 5 7 100 3 5 7 200 3 5 7 第二次减数分裂（对照）
Meiosis Ⅱ (control)0 3 5 7 1.2.3 新疆杨2n花粉发生频率的统计
用镊子分别夹取少量不同处理组合的花粉，均匀地洒在滴有醋酸洋红的载玻片上制作临时涂片，在Olympus BX-51光学显微镜下随机选取5个视野观察并进行拍照，进一步利用Image J软件测量花粉粒直径。以对照组花粉粒平均直径的1.28倍作为判别2n花粉的标准[39]，统计每个视野中总花粉粒数与2n花粉粒的数量，以2n花粉占花粉粒总数的百分比作为2n花粉发生频率，计算出不同处理组合的2n花粉发生频率。试验重复3次，每次重复测量200 ~ 300个花粉粒。
1.2.4 新疆杨成熟花粉精核和营养核的细胞学观察
用镊子夹取少量不同处理组合的花粉，均匀的洒在滴有蒸馏水的载玻片上水合2 h，用DAPI（4,6-diamidino-2-phenylindole）对细胞核进行染色，用Olympus BX-51荧光显微镜观察新疆杨单倍性与2n花粉粒精核与营养核的形态差异。
1.2.5 新疆杨花粉母细胞减数分裂微管骨架的免疫荧光标记
摘取赤霉素处理后的新疆杨雄花芽新鲜花药，迅速投入PEM缓冲液（50 mmol/L PIPES，5 mmol/L EGTA，1 mmol/L MgSO4，pH 6.8）配制的4%多聚甲醛固定液中固定45 min后，参考Brown等[40]的方法，开展赤霉素处理对新疆杨花粉母细胞减数分裂微管骨架的影响研究，利用激光共聚焦扫描显微镜观察并照相。
1.2.6 数据处理方法
2n花粉诱导率等百分比数据经反正弦变换[41]后，利用SPSS18.0统计软件进行差异显著性分析和多重比较。
2. 结果与分析
2.1 新疆杨雄花芽减数分裂进程的观察
温室水培条件下，新疆杨花粉母细胞由减数分裂细线期发育至四分体时期历时约4 d（图1，表2）。其花粉母细胞的减数分裂是一个连续而非同步的过程（表2）。温室水培96 h后，花粉母细胞开始进入减数分裂，大部分花粉母细胞发育至细线期（65.5%）（图1A），此时期染色体呈细丝状，分布杂乱无章。水培至124 h，减数分裂发育至粗线期（40.0%）（图1B），染色体收缩变粗，形态较细线期明显。水培136 h后，大部分花粉母细胞进入双线期（47.5%）（图1C），染色体进一步缩聚，可明显观察到同源染色体非姊妹染色单体间的交叉。水培144 h后，大部分花粉母细胞进入终变期（46.0%）（图1D），此时染色体长度继续缩短，核仁缩小乃至消失，染色体缩聚呈点状，形态清晰可见。水培至152 h时，花粉母细胞发育至中期Ⅰ（27.0%）（图1E），联会的同源染色体整齐地排列在赤道板中央。水培160 h后，大部分花粉母细胞进入后期Ⅰ（22.5%）和末期Ⅰ（37.5%），联会的同源染色体开始分离，移向细胞的两极（图1F）。当染色体抵达细胞的两极后，开始解聚（图1G），随后核仁重新出现，形成2个子核，标志着第一次减数分裂结束。新疆杨花粉母细胞完成第一次减数分裂大约需要64 h。
![]()
图 1 新疆杨花粉母细胞减数分裂进程A. 细线期；B. 粗线期；C. 双线期；D. 终变期；E. 中期Ⅰ；F. 后期Ⅰ；G. 末期Ⅰ；H. 前期Ⅱ；I. 中期Ⅱ；J. 后期Ⅱ；K ~ L. 末期Ⅱ；M. 末期Ⅱ多核仁现象；N. 四分体；O. 三分体。标尺 = 10 μm。A, leptotene; B, pachytene; C, diplotene; D, diakinesis; E, metaphase I; F, anaphase I; G, telophase I; H, prophase Ⅱ; I, metaphase Ⅱ; J, anaphase Ⅱ; K-L, telophase Ⅱ; M, multi-nucleoli phenomena at telophase Ⅱ; N, tetrad; O, triad. Bar = 10 μm.Figure 1. Meiotic process of PMC in P. bolleana表 2 新疆杨花粉母细胞的减数分裂进程Table 2. Meiotic process of PMC in P. bolleana% 减数分裂时期
Meiotic stage水培时间 hydroponics time/h 96 110 124 136 144 152 160 168 174 180 186 192 细线期 Leptotene 65.5 28.5 9.0 3.5 偶线期 Zygotene 25.5 51.5 36.5 14.5 5.0 粗线期 Pachytene 9.0 16.5 40.0 27.5 18.0 12.5 双线期 Diplotene 3.5 14.5 47.5 22.0 14.0 终变期 Diakinesis 7.0 46.0 20.5 中期Ⅰ Metaphase Ⅰ 9.0 27.0 8.5 6.0 后期Ⅰ Anaphase Ⅰ 10.0 22.5 8.5 末期Ⅰ Telophase Ⅰ 12.0 37.5 13.0 5.5 2.0 前期Ⅱ Prophase Ⅱ 4.0 21.0 34.0 28.0 11.0 5.5 0.5 中期Ⅱ Metaphase Ⅱ 10.5 30.5 37.0 25.5 6.5 2.0 后期Ⅱ Anaphase Ⅱ 8.0 21.0 39.0 10.5 2.5 末期Ⅱ Telophase Ⅱ 8.5 19.5 50.5 14.0 四分体 Tetrad 3.0 27.0 81.0 注：表中数据为所在时期花粉母细胞比例。 Note: data in the table are the proportion of PMCs in the period. 第一次减数分裂结束后，花粉母细胞未进行细胞质分裂，继续开始第二次分裂。当雄花枝水培至168 h时，花粉母细胞发育至前期Ⅱ（34.0%）（图1H），此时2个子核区的染色体重新缩聚，核仁逐渐消失。水培174 h后，大部分花粉母细胞进入中期Ⅱ（37.0%）（图1I），2组染色体整齐地排列在子细胞的赤道板上。水培180 h后，大部分花粉母细胞进入后期Ⅱ（39.0%）（图1J），姐妹染色单体开始分开，移向子细胞的两极。水培186 h后，花粉母细胞发育至末期Ⅱ（50.5%）（图1K），染色单体抵达子细胞的两极，染色体解螺旋，核仁重新出现（图1L），此时期部分花粉母细胞存在多核仁现象（图1M）。水培192 h后，绝大部分花粉母细胞进入四分体时期（81.0%）（图1N）。此时期还可以观察到少量的三分体存在（图1O）。第二次减数分裂经历的时间明显少于第一次减数分裂，仅需要32 h。
2.2 新疆杨2n花粉的赤霉素诱导
新疆杨雄花芽经不同浓度的赤霉素喷洒处理后，部分花芽因赤霉素的刺激而生长发育过快，导致干枯死亡或提早脱落（图2A），未收集到花粉。存活花芽继续水培于温室，直至花药成熟开裂，均可收集到一定量的花粉。花粉经2%的醋酸洋红染色后，在光学显微镜下进行花粉粒镜检，可观察大部分花粉粒被染成红色，且着色均匀，这说明该花粉粒具有较好的生活力。通过花粉粒形态比较可以发现：新疆杨对照组中存在少量的天然2n花粉（图2B），其发生频率为0.11% ~ 0.42%（表3）。
![]()
图 2 赤霉素喷洒处理诱导新疆杨2n花粉及其对花芽发育的影响A. 赤霉素处理后干瘪的花芽；B. 天然2n花粉（箭头指向）；C.赤霉素诱导的2n花粉（箭头指向）；D. DAPI染色后的单倍性花粉；E. DAPI染色后的赤霉素诱导的2n花粉。A, the dry flower buds after GA3 treatment; B, natural 2n pollen (arrow); C, induced 2n pollen by GA3 treatment (arrow); D, haploid pollen after DAPI staining; E. induced 2n pollen by GA3 treatment after DAPI staining.Figure 2. 2n pollen induced by GA3 and its effects on the development of male flower buds in P. bolleana表 3 不同赤霉素喷洒处理组合的新疆杨2n花粉诱导率Table 3. Induction rates of 2n pollen with different combinations of GA3 spraying treatments in P. bolleana减数分裂时期 Meiotic stage 赤霉素浓度
Concentration of GA3/（μmol· L−1）喷洒次数
Spraying time2n花粉诱导率
Induction rate of 2n pollen/%花粉母细胞时期 PMC stage 50 3 0.40 ± 0.26 5 1.03 ± 0.92 7 5.01 ± 1.74 100 3 3.01 ± 0.36 5 4.69 ± 0.44 7 5.04 ± 2.16 200 3 3.22 ± 0.40 5 4.89 ± 0.53 7 5.00 ± 3.63 第一次减数分裂 Meiosis Ⅰ 50 3 2.15 ± 1.18 5 2.12 ± 0.21 7 5.98 ± 1.35 100 3 2.67 ± 1.76 5 2.55 ± 0.80 7 4.98 ± 3.78 200 3 5.53 ± 1.64 5 4.03 ± 1.29 7 4.72 ± 0.44 第二次减数分裂 Meiosis Ⅱ 50 3 4.72 ± 1.83 5 5.50 ± 1.89 7 8.83 ± 3.10 100 3 4.08 ± 1.25 5 3.09 ± 0.96 7 3.04 ± 2.75 200 3 4.25 ± 2.13 5 5.66 ± 1.29 7 5.22 ± 1.35 第二次减数分裂（对照） Meiosis Ⅱ （control） 0 3 0.11 ± 0.04 5 0.42 ± 0.27 7 0.26 ± 0.11 不同赤霉素处理组合中，均可观察到一定比例的2n花粉（图2C）。花粉粒经DAPI染色后，无论是单倍性花粉粒（图2D）还是2n花粉粒（图2E），均可观察到1个较大的营养核和2个高度凝缩的精核。但2n花粉粒的营养核和精核明显大于正常的单倍性花粉。这说明赤霉素喷洒处理诱导产生的大花粉的确是2n花粉。
从表3可以看出：不同处理组合的2n花粉诱导率存在明显的差异，其中以50 μmol/L的赤霉素于第二次减数分裂时期喷洒处理7次的效果最好，平均2n花粉诱导率达到（8.83 ± 3.10）%；以50 μmol/L的赤霉素于花粉母细胞时期喷洒处理3次的效果最差，平均2n花粉诱导率仅为（0.40 ± 0.26）%。不同处理组合2n花粉诱导率的方差分析结果表明（表4）：处理次数（F = 9.600，P < 0.001）对新疆杨2n花粉诱导率具有极显著影响；处理时期（F = 4.426，P = 0.017）对2n花粉诱导率具有显著影响；赤霉素浓度（F = 2.538，P = 0.088）对2n花粉诱导率的影响不显著。处理时期和处理次数的多重比较结果表明（表5）：第二次减数分裂时期的2n花粉诱导率显著高于花粉母细胞时期和第一次减数分裂时期；赤霉素喷施7次的2n花粉诱导率显著高于3次和5次。因此，当花粉母细胞发育至第二次减数分裂时，以50 μmol/L的赤霉素进行7次喷洒处理是诱导新疆杨花粉染色体加倍的最佳处理组合。
表 4 不同处理组合新疆杨2n花粉诱导率的方差分析Table 4. Variance analysis on induction rates of GA3-induced 2n pollen in P. bolleana under different treatment combinations变异来源 Source of variation df F P 减数分裂时期 Meiotic stage 2 4.426* 0.017 喷洒次数 Spraying time 2 9.600** 0.000 赤霉素浓度 Concentration of GA3 2 2.538 0.088 减数分裂时期 × 喷洒次数 Meiotic stage × spraying time 4 0.940 0.448 减数分裂时期 × 赤霉素浓度 Meiotic stage × concentration of GA3 4 5.103** 0.001 喷洒次数 × 赤霉素浓度 Spraying time × concentration of GA3 4 3.565* 0.012 减数分裂时期 × 喷洒次数 × 赤霉素浓度
Meiotic stage × spraying time × concentration of GA38 0.562 0.804 注：**代表差异极显著（P < 0.01），*代表差异显著（P < 0.05）。Notes: ** represents extremely significant differences (P < 0.01), * represents significant differences (P < 0.05). 表 5 不同减数分裂时期和不同喷洒次数新疆杨2n花粉诱导率的多重比较Table 5. Multiple comparisons in induction rates of 2n pollen under different meiotic stages anddifferent spraying times in P. bolleana减数分裂时期
Meiotic stage平均2n花粉诱导率
Average induction rate of 2n pollen/%喷洒次数
Spraying time平均2n花粉诱导率
Average induction rate of 2n pollen/%花粉母细胞时期 PMC stage 3.59 ± 2.18b 3 3.34 ± 1.85b 第一次减数分裂 Meiosis Ⅰ 3.86 ± 2.03b 5 3.73 ± 1.78b 第二次减数分裂 Meiosis Ⅱ 4.93 ± 2.33a 7 5.31 ± 2.54a 注：不同小写字母代表差异显著（P < 0.05）。 Note: different lowercase letters represent significant differences (P < 0.05). 2.3 赤霉素喷洒处理对新疆杨花粉母细胞减数分裂微管骨架的影响
在新疆杨花粉母细胞减数分裂过程中，其微管骨架的分布总体表现出规律性变化（图3A ~ L）。花粉母细胞进入减数分裂前期Ⅰ时，网状核周微管分布于细胞质中，包裹着细胞核（图3A ~ B）。花粉母细胞发育至终变期时，微管发生方向改变，原来分布于胞质中的微管减少，逐步转向细胞核内，包围着细胞核内的染色体（图3C）。当花粉母细胞进入中期Ⅰ时，从两极发出的细而长的极丝微管和与着丝粒相连的短而粗的着丝粒微管聚集形成狭窄的纺锤体维管束，连接着排列在赤道板上的染色体（图3D）。花粉母细胞发育至后期Ⅰ时，纺锤体微管牵引着同源染色体移向细胞两极（图3E）。同源染色体抵达细胞两极后，花粉母细胞进入末期I，此时微管又一次转向，分布至整个细胞质区域（图3F）。
![]()
图 3 新疆杨花粉母细胞减数分裂过程中微管骨架的动态变化A. 细线期；B. 偶线期；C. 终变期；D. 中期Ⅰ；E. 后期Ⅰ；F. 末期Ⅰ；G. 前期Ⅱ；H. 中期Ⅱ；I. 后期Ⅱ；J. 末期Ⅱ；K. 四分体；L. 三分体；M. 赤霉素处理后的前期Ⅱ；N. 赤霉素处理后的中期Ⅱ；O. 赤霉素处理后的后期Ⅱ；P ~ S. 赤霉素处理后的末期Ⅱ；T. 赤霉素诱导的三分体。标尺 = 10 μm。 A, leptotene; B, zygotene; C, diakinesis; D, metaphase I; E, anaphase I; F, telophase I; G, prophase Ⅱ; H, metaphase Ⅱ; I, anaphase Ⅱ; J, telophase Ⅱ; K, tetrad; L, triad; M, prophase Ⅱ by GA3 treatment; N, metaphase Ⅱ by GA3 treatment; O, anaphase Ⅱ by GA3 treatment; P–S, telophase Ⅱ in treated PMCs by GA3 treatment; T, induced triad by GA3 treatment. Bar = 10 μm.Figure 3. Dynamic changes of microtubule cytoskeleton during meiosis of PMCs in P. bolleana第一次减数分裂结束后，花粉母细胞不进行细胞质分裂，立即进入第二次减数分裂。花粉母细胞发育至前期Ⅱ时，核周微管移向2个子细胞核内，与染色体相连（图3G）。随着微管的进一步聚合，2个子细胞核内分别出现纺锤体微管，连接着排列在赤道板上的染色体，此时花粉母细胞已发育至中期Ⅱ（图3H）。当花粉母细胞发育至后期Ⅱ时，子细胞内的纺锤体微管牵引着姐妹染色体单体移向两极（图3I），当姐妹染色体抵达子细胞的两极后，辐射状微管在核膜处扩散，延伸至整个细胞质，联系着4个子核，标志着第二次减数分裂结束（图3J）。在辐射状微管的驱动下，自细胞周沿发生向心胞质分裂，形成四分体（图3K），偶尔还可观察到三分体（图3L）。
与对照组相比，在第二次减数分裂时期赤霉素喷洒处理7次的处理组中，其微管骨架除末期Ⅱ呈现异常外（图3P ~ S），其他减数分裂时期的微管变化规律与对照组相似（图3M ~ O）。在新疆杨花粉母细胞第二次减数分裂时期赤霉素喷洒处理7次后，末期Ⅱ相邻子核之间的辐射状微管发生部分缺失（图3P ~ S），与对照组末期Ⅱ的微管骨架明显不同。末期Ⅱ辐射状微管的部分缺失可导致2个相邻子核间胞质分裂失败，从而产生核复原现象，形成三分体（图3T），进而发育成1个单倍性配子和1个2n配子。
3. 讨 论
新疆杨花粉母细胞减数分裂具有连续而不同步性，这可能是植物性母细胞减数分裂过程中的一种普遍现象。研究人员曾在银灰杨[25]、毛白杨[42]、通辽杨（Populus simonⅡ × P. nigra ‘Tongliao’）[43]、响叶杨[44]、美洲黑杨（Populus deltoides）[45]以及胡杨（Populus euphratica）等[46]物种的花粉母细胞减数分裂的研究中均发现不同步性。同一花枝不同部位花芽，其花粉母细胞减数分裂进程不同；同一花芽的不同小花内花粉母细胞减数分裂时期亦不相同；甚至同一小花内的不同花粉母细胞，其减数分裂进程也不相同。在恶劣环境条件下，这种复杂的花芽发育不同步性可以延长花期，有利于增强两性功能的有效作用时间，对其种群繁殖后代有利[47]。
利用化学诱变剂人工诱导植物花粉染色体加倍获得2n花粉时，诱变剂浓度、处理方式和处理时期是影响2n花粉得率的关键影响因素。李赟等[48]利用0.5%的秋水仙碱于细线末期至粗线期时对银白杨花粉母细胞进行7次注射处理，获得了高达82.33%的2n花粉。Zhou等[49]利用0.5%的秋水仙碱于粗线期时对银灰杨花粉母细胞进行11次注射处理，获得了30.3 %的2n花粉。本研究利用50、100、200 μmol/L赤霉素于花粉母细胞、第一次减数分裂以及第二次减数分裂时期对新疆杨雄花芽进行喷洒处理，不同处理组合均可获得一定比例的2n花粉，说明利用赤霉素处理诱导杨树2n花粉是可行的。其中，以50 μmol/L赤霉素于第二次减数分裂时期喷洒处理7次的最佳处理组合仅获得8.83%的2n花粉，原因主要包括2个方面：一是当花粉母细胞发育至第二次减数分裂前后，花芽一直被芽鳞包被，喷洒的赤霉素未能直接与花药接触而被蒸发，降低了赤霉素的处理效果；另一方面则是新疆杨同一花芽不同花药内花粉母细胞减数进程同步性差，影响了赤霉素的处理效果。
微管骨架在细胞分裂过程中发挥着重要作用。植物减数分裂的正常进行与微管骨架的有序分布和动态变化密切相关[47]。王君[18]对通辽杨的减数分裂免疫荧光分析研究发现：减数分裂末期Ⅱ细胞成膜体缺失导致胞质分裂失败，细胞内子核未能正常分离而共存于同1个小孢子内，形成二分体和三分体，进一步发育形成2n花粉。鲁敏[23]在响叶杨天然2n花粉发生机制的研究中发现：减数分裂后期Ⅱ纺锤体微管缺失导致姐妹染色单体不分离，末期Ⅱ成膜体微管缺失导致胞质分裂异常，二者是天然2n花粉产生的细胞学机制。张平冬等[50]在胡杨减数分裂过程中观察到平行纺锤体、三极纺锤体和融合纺锤体，这些异常的纺锤体与胡杨天然2n花粉发生密切相关。Liu等[33]对拟南芥施加赤霉素处理后研究发现：减数分裂末期Ⅱ相邻子核间的微管缺失，导致子核聚集在同1个小孢子中，发育形成2n花粉。本研究通过免疫荧光技术探索了赤霉素处理对新疆杨花粉母细胞微管骨架构建及分布的影响，发现赤霉素处理可导致减数分裂末期Ⅱ辐射状微管产生部分缺失，致使相邻子核间胞质分裂失败，从而产生核复原现象，形成三分体，进而发育成1个单倍性配子和1个2n配子。
针对杨树花药被芽鳞包被而导致赤霉素喷洒处理效果不理想的问题，下一阶段的工作应聚焦于赤霉素处理方式、药剂浓度以及处理时间等技术条件优化研究，旨在提高赤霉素诱导2n花粉的效率。
4. 结 论
在温室水培条件下，新疆杨花粉母细胞完成减数分裂历时约4 d，且存在明显的不同步性。利用浓度为50、100、200 μmol/L赤霉素于花粉母细胞、第一次减数分裂以及第二次减数分裂时期分别对新疆杨雄花芽进行喷洒处理，不同处理组合均可获得一定比例的2n花粉，说明利用赤霉素处理诱导新疆杨获得2n花粉是可行的。其中以50 μmol/L赤霉素于第二次减数分裂时期喷施处理7次为最佳处理组合，可以获得8.83%的2n花粉。在最佳处理组合条件下，新疆杨雄花芽经赤霉素喷洒处理可导致减数分裂末期Ⅱ辐射状微管产生部分缺失，致使相邻子核间胞质分裂失败，从而产生核复原现象，形成三分体，进而产生2n花粉。
-
![]()
图 1 红松叶片性状的主成分分析载荷图
Figure 1. Principal component analysis loading plot for leaf traits of Pinus koraiensis
![]()
图 2 紫椴叶片性状的主成分分析载荷图
Figure 2. Principal component analysis loading plot for leaf traits of Tilia amurensis
表 1 林分树种组成
Table 1 Composition of stand tree species
树种
Tree species每公顷胸高断面积/（m2·hm−2）
Basal area at breast height per hectare/(m2·ha−1)胸高断面积比例
Proportion of basal area at breast height/%水曲柳 Fraxinus mandshurica 7.860 24.281 紫椴 Tilia amurensis 4.296 13.272 五角槭 Acer pictum subsp. mono 4.116 12.715 红松 Pinus koraiensis 3.895 12.031 蒙古栎 Quercus mongolica 2.899 8.956 拧筋槭 Acer triflorum 1.987 6.138 白桦 Betula platyphylla 1.743 5.383 春榆 Ulmus davidiana var. japonica 1.375 4.249 白牛槭 Acer mandshuricum 0.778 2.403 枫桦 Betula costata 0.726 2.243 胡桃楸 Juglans mandshurica 0.653 2.017 香杨 Populus koreana 0.341 1.053 其他 Others 1.703 5.260 总计 Total 32.373 100.000 
下载: 导出CSV
表 2 样地内基础变量信息统计表
Table 2 Statistical table of basic variable information within the sample plots
变量 Variable 范围 Range 平均值 Mean 标准偏差 Standard deviation 海拔 Elevation (E)/m 468.57 ~ 486.49 477.24 4.20 坡度 Slope (Sl)/（°） 0.62 ~ 24.16 4.71 4.27 坡向 Aspect/（°） 0.89 ~ 359.26 124.18 85.05
下载: 导出CSV
表 3 试验树概况
Table 3 Description of test trees
树种
Tree species胸径
DBH (D)/cm树高
Tree height/m枝下高
Height under
branch/m红松 Pinus koraiensis 9.8 ± 3.4 8.6 ± 3.1 3.5 ± 1.3 紫椴 Tilia amurensis 8.3 ± 3.8 8.8 ± 3.7 3.1 ± 1.4
下载: 导出CSV
表 4 红松紫椴叶片性状
Table 4 Leaf traits of Pinus koraiensis and Tilia amurensis
叶片性状 Leaf trait 树种 Tree species 红松 P. koraiensis 紫椴 T. amurensis 光饱和净光合速率 Light saturation net photosynthetic rate/(μmol·m−2·s−1) 10.324 ± 2.672a 7.161 ± 1.774b 暗呼吸速率 Dark respiration rate/(μmol·m−2·s−1) 0.594 ± 0.250a 0.373 ± 0.099b 光饱和气孔导度 Light saturation stomatal conductance/(mol·m−2·s−1) 0.081 ± 0.026a 0.116 ± 0.043b 比叶面积 Specific leaf area/(cm2·g−1) 55.705 ± 5.433a 319.707 ± 48.634b 注：同行不同小写字母表示差异显著（P < 0.05）。Note: different lowercase letters in the same row indicate significant differences (P < 0.05).
下载: 导出CSV
表 5 红松、紫椴年均生长量的回归模型
Table 5 Regression models for annual mean growth of Pinus koraiensis and Tilia amurensis
树种
Tree species响应变量
Response variable解释变量
Explanatory variable系数 Coefficient 方程表现 Model performance 初始胸径
Initial DBH (D0)
INC地形
Topography叶性状
Leaf traitR2adj RMSE 红松
Pinus koraiensis年均胸径生长量
Annual mean DBH growthD0 + INC + Sl + SLA −0.171* −0.266** −0.178* 0.131 0.477 0.363 年均生物量生长量
Annual mean biomass growthD0 + INC + Sl + SLA 0.835*** −0.243* −0.206* 0.098 0.789 0.419 紫椴
Tilia amurensis年均胸径生长量
Annual mean DBH growthD0 + E + SLA −0.778*** −0.367** −0.159 0.746 0.425 年均生物量生长量
Annual mean biomass growthD0 + E 0.366* −0.536** 0.417 0.692 注：INC.领域拥挤指数；SLA.比叶面积。*代表P < 0.05，**代表P < 0.01，***代表P < 0.001。Notes: INC, neighborhood crowding index; SLA, specific leaf area. * represents P < 0.05, ** represents P < 0.01, and *** represents P < 0.001.
下载: 导出CSV
-
[1] Adler P B, Salguero-Gomez R, Compagnoni A, et al. Functional traits explain variation in plant life history strategies[J]. Proceedings of the National Academy of Sciences, 2014, 111(2): 740−745. doi: 10.1073/pnas.1315179111
[2] Poorter L, Castilho C V, Schietti J, et al. Can traits predict individual growth performance? A test in a hyperdiverse tropical forest[J]. New Phytologist, 2018, 219(1): 109−121. doi: 10.1111/nph.15206
[3] Ding Y, Zang R, Huang J, et al. Intraspecific trait variation and neighborhood competition drive community dynamics in an old-growth spruce forest in northwest China[J]. Science of The Total Environment, 2019, 678: 525−532. doi: 10.1016/j.scitotenv.2019.05.014
[4] van der Sande M T, Peña-Claros M, Ascarrunz N, et al. Abiotic and biotic drivers of biomass change in a Neotropical forest[J]. Journal of Ecology, 2017, 105(5): 1223−1234. doi: 10.1111/1365-2745.12756
[5] Guimarães Z T M, dos Santos V A H F, Nogueira W L P, et al. Leaf traits explaining the growth of tree species planted in a central Amazonian disturbed area[J]. Forest Ecology and Management, 2018, 430: 618−628. doi: 10.1016/j.foreco.2018.08.048
[6] Liu Z, Hikosaka K, Li F, et al. Plant size, environmental factors and functional traits jointly shape the stem radius growth rate in an evergreen coniferous species across ontogenetic stages[J]. Journal of Plant Ecology, 2021, 14(2): 257−269. doi: 10.1093/jpe/rtaa093
[7] Weintraub S R, Taylor P G, Porder S, et al. Topographic controls on soil nitrogen availability in a lowland tropical forest[J]. Ecology, 2015, 96(6): 1561−1574. doi: 10.1890/14-0834.1
[8] Fotis A T, Murphy S J, Ricart R D, et al. Above-ground biomass is driven by mass-ratio effects and stand structural attributes in a temperate deciduous forest[J]. Journal of Ecology, 2018, 106(2): 561−570. doi: 10.1111/1365-2745.12847
[9] Hao M, Zhang C, Zhao X, et al. Functional and phylogenetic diversity determine woody productivity in a temperate forest[J]. Ecology and Evolution, 2018, 8(5): 2395−2406. doi: 10.1002/ece3.3857
[10] Fortunel C, Lasky J R, Uriarte M, et al. Topography and neighborhood crowding can interact to shape species growth and distribution in a diverse Amazonian forest[J]. Ecology, 2018, 99(10): 2272−2283. doi: 10.1002/ecy.2441
[11] Hao M, Messier C, Geng Y, et al. Functional traits influence biomass and productivity through multiple mechanisms in a temperate secondary forest[J]. European Journal of Forest Research, 2020, 139(6): 959−968. doi: 10.1007/s10342-020-01298-0
[12] Yang J, Song X, Zambrano J, et al. Intraspecific variation in tree growth responses to neighbourhood composition and seasonal drought in a tropical forest[J]. Journal of Ecology, 2021, 109(1): 26−37. doi: 10.1111/1365-2745.13439
[13] Liu X, Swenson N G, Lin D, et al. Linking individual-level functional traits to tree growth in a subtropical forest[J]. Ecology, 2016, 97(9): 2396−2405. doi: 10.1002/ecy.1445
[14] Brown J H, Gillooly J F, Allen A P, et al. Toward a metabolic theory of ecology[J]. Ecology, 2004, 85(7): 1771−1789. doi: 10.1890/03-9000
[15] Niklas K J, Enquist B J. Invariant scaling relationships for interspecific plant biomass production rates and body size[J]. Proceedings of the National Academy of Sciences, 2001, 98(5): 2922−2927. doi: 10.1073/pnas.041590298
[16] Cheng D L, Li T, Zhong Q L, et al. Scaling relationship between tree respiration rates and biomass[J]. Biology Letters, 2010, 6(5): 715−717. doi: 10.1098/rsbl.2010.0070
[17] Paine C E T, Amissah L, Auge H, et al. Globally, functional traits are weak predictors of juvenile tree growth, and we do not know why[J]. Journal of Ecology, 2015, 103(4): 978−989. doi: 10.1111/1365-2745.12401
[18] Violle C, Navas M L, Vile D, et al. Let the concept of trait be functional![J]. Oikos, 2007, 116(5): 882−892. doi: 10.1111/j.0030-1299.2007.15559.x
[19] Chaturvedi R K, Raghubanshi A S, Singh J S. Leaf attributes and tree growth in a tropical dry forest: leaf attributes and tree growth[J]. Journal of Vegetation Science, 2011, 22(5): 917−931. doi: 10.1111/j.1654-1103.2011.01299.x
[20] Pérez-Harguindeguy N, Díaz S, Garnier E, et al. Corrigendum to: new handbook for standardised measurement of plant functional traits worldwide[J]. Australian Journal of Botany, 2016, 64(8): 715. doi: 10.1071/BT12225_CO
[21] Reich P B. The world-wide ‘fast-slow’ plant economics spectrum: a traits manifesto[J]. Journal of Ecology, 2014, 102(2): 275−301. doi: 10.1111/1365-2745.12211
[22] Katabuchi M, Wright S J, Swenson N G, et al. Contrasting outcomes of species- and community-level analyses of the temporal consistency of functional composition[J]. Ecology, 2017, 98(9): 2273−2280. doi: 10.1002/ecy.1952
[23] 李冠男, 张春雨, 赵秀海. 吉林蛟河针阔混交林幼苗动态及环境驱动因子[J]. 北京林业大学学报, 2021, 43(8): 41−49. Li G N, Zhang C Y, Zhao X H. Seedling dynamics and environmental driving factors of coniferous and broadleaved mixed forest in Jiaohe, Jilin Province of northeastern China[J]. Journal of Beijing Forestry University, 2021, 43(8): 41−49.
[24] 姚杰, 宋子龙, 张春雨, 等. 距离和密度制约对吉林蛟河阔叶红松林幼苗生长的影响[J]. 北京林业大学学报, 2019, 41(5): 108−117. Yao J, Song Z L, Zhang C Y, et al. Effects of distance and density dependence on seedling growth in a broadleaved Korean pine forest in Jiaohe of Jilin Province, northeastern China[J]. Journal of Beijing Forestry University, 2019, 41(5): 108−117.
[25] 郝珉辉, 代莹, 岳庆敏, 等. 阔叶红松林功能多样性与森林碳汇功能关系[J]. 北京林业大学学报, 2022, 44(10): 68−76. Hao M H, Dai Y, Yue Q M, et al. Relationship between functional diversity of broadleaved Korean pine forest and forest carbon sink function[J]. Journal of Beijing Forestry University, 2022, 44(10): 68−76.
[26] 胡乘风, 陈巧玲, 乔雪涛, 等. 阔叶红松林主要树种光合与光谱反射特性及初级生产力研究[J]. 北京林业大学学报, 2020, 42(5): 12−24. Hu C F, Chen Q L, Qiao X T, et al. Photosynthetic, spectral reflectance characteristics and primary productivity of main tree species in broadleaved Korean pine forest[J]. Journal of Beijing Forestry University, 2020, 42(5): 12−24.
[27] 胡乘风. 阔叶红松林冠层树种光合特性及光合生产潜力研究[D]. 北京: 北京林业大学, 2021. Hu C F. Study on photosynthetic characteristics and potential productivity of canopy tree species in broad-leaved Korean pine forest[D]. Beijing: Beijing Forestry University, 2021.
[28] Hardwick R C. Mathematical models in plant physiology[J]. Experimental Agriculture, 1977, 13(1): 112. doi: 10.1017/S0014479700007675
[29] Ball J T, Woodrow I E, Berry J A. A model predicting stomatal conductance and its contribution to the control of photosynthesis under different environmental conditions[M]//Biggins J. Progress in photosynthesis research. Dordrecht: Springer Netherlands, 1987: 221−224.
[30] 张林, 金冬梅. 针叶树种叶面积的主要测定方法[C]//第八届全国生物多样性保护与持续利用研讨会论文集. 北京: 气象出版社, 2008: 293−297. Zhang L, Jin D M. Principal methods of measuring leaf area of cornifer trees[C]//Procceedings of The 8th National Symposium on Biodiversity Conservation and Sustainable Use of Biodiversity in China. Beijing: China Meteorological Press, 2008: 293−297.
[31] Huston M A. Local processes and regional patterns: appropriate scales for understanding variation in the diversity of plants and animals[J]. Oikos, 1999, 86(3): 393. doi: 10.2307/3546645
[32] Chiang J M, Spasojevic M J, Muller L H C, et al. Functional composition drives ecosystem function through multiple mechanisms in a broadleaved subtropical forest[J]. Oecologia, 2016, 182(3): 829−840. doi: 10.1007/s00442-016-3717-z
[33] Li Y, Kröber W, Bruelheide H, et al. Crown and leaf traits as predictors of subtropical tree sapling growth rates[J]. Journal of Plant Ecology, 2017, 10(1): 136−145. doi: 10.1093/jpe/rtw041
[34] Hidaka A, Kitayama K. Divergent patterns of photosynthetic phosphorus-use efficiency versus nitrogen-use efficiency of tree leaves along nutrient-availability gradients[J]. Journal of Ecology, 2009, 97(5): 984−991. doi: 10.1111/j.1365-2745.2009.01540.x
[35] Coble A P, Fogel M L, Parker G G. Canopy gradients in leaf functional traits for species that differ in growth strategies and shade tolerance[J]. Tree Physiology, 2017, 37(10): 1415−1425. doi: 10.1093/treephys/tpx048
[36] Han T, Lu H, Ren H, et al. Are reproductive traits of dominant species associated with specific resource allocation strategies during forest succession in southern China?[J]. Ecological Indicators, 2019, 102: 538−546. doi: 10.1016/j.ecolind.2019.03.007
[37] Lusk C H, Warton D I. Global meta-analysis shows that relationships of leaf mass per area with species shade tolerance depend on leaf habit and ontogeny[J]. New Phytologist, 2007, 176(4): 764−774. doi: 10.1111/j.1469-8137.2007.02264.x
[38] Poorter H, Niinemets Ü, Ntagkas N, et al. A meta-analysis of plant responses to light intensity for 70 traits ranging from molecules to whole plant performance[J]. New Phytologist, 2019, 223(3): 1073−1105. doi: 10.1111/nph.15754
[39] Wang H, Zhang M, Nan H. Abiotic and biotic drivers of species diversity in understory layers of cold temperate coniferous forests in North China[J]. Journal of Forestry Research, 2019, 30(6): 2213−2225. doi: 10.1007/s11676-018-0795-2
[40] Lin D, Anderson-Teixeira K J, Lai J, et al. Traits of dominant tree species predict local scale variation in forest aboveground and topsoil carbon stocks[J]. Plant and Soil, 2016, 409(1−2): 435−446. doi: 10.1007/s11104-016-2976-0
[41] Asefa M, Song X, Cao M, et al. Temporal trait plasticity predicts the growth of tropical trees[J/OL]. Journal of Vegetation Science, 2021, 32(4): 292615[2023−11−09]. https://doi.org/10.1101/2020.09.11.292615.
[42] Fortunel C, Valencia R, Wright S J, et al. Functional trait differences influence neighbourhood interactions in a hyperdiverse Amazonian forest[J]. Ecology Letters, 2016, 19(9): 1062−1070. doi: 10.1111/ele.12642
[43] Zhang J, Huang S, He F. Half-century evidence from western Canada shows forest dynamics are primarily driven by competition followed by climate[J]. Proceedings of the National Academy of Sciences, 2015, 112(13): 4009−4014. doi: 10.1073/pnas.1420844112
[44] Clark J S, Bell D M, Hersh M H, et al. Climate change vulnerability of forest biodiversity: climate and competition tracking of demographic rates: forests, climate change, and competition[J]. Global Change Biology, 2011, 17(5): 1834−1849. doi: 10.1111/j.1365-2486.2010.02380.x
[45] Freckleton R P, Watkinson A R. Asymmetric competition between plant species[J]. Functional Ecology, 2001, 15(5): 615−623. doi: 10.1046/j.0269-8463.2001.00558.x
[46] Uriarte M, Swenson N G, Chazdon R L, et al. Trait similarity, shared ancestry and the structure of neighbourhood interactions in a subtropical wet forest: implications for community assembly: Traits, phylogeny, neighbourhood interactions[J]. Ecology Letters, 2010, 13(12): 1503−1514. doi: 10.1111/j.1461-0248.2010.01541.x
[47] Balachowski J A, Volaire F A. Implications of plant functional traits and drought survival strategies for ecological restoration[J]. J. James. Journal of Applied Ecology, 2018, 55(2): 631−640.
[48] Martínez-Garza C, Bongers F, Poorter L. Are functional traits good predictors of species performance in restoration plantings in tropical abandoned pastures?[J]. Forest Ecology and Management, 2013, 303: 35−45. doi: 10.1016/j.foreco.2013.03.046
[49] Kunstler G, Falster D, Coomes D A, et al. Plant functional traits have globally consistent effects on competition[J]. Nature, 2016, 529: 204−207. doi: 10.1038/nature16476
-
期刊类型引用(3)
1. 李婧，左欣欣，赵培伶，张陆玉，徐文鸾，张徐杨，倪德江，王明乐. 茶树高亲和硝酸盐转运蛋白家族基因NRT2的鉴定与表达. 应用与环境生物学报. 2022(01): 50-56 . 
百度学术
2. 杨利艳，高源，朱满喜，邓妍，王创云. 藜麦NRT2基因家族的鉴定及表达分析. 华北农学报. 2022(S1): 8-18 . 百度学术
3. 张云，赵艳菲，王雅平，牟彬，张嘉越，马浩然，韩玉珠. 延薯4号马铃薯对氮素的生理生化响应及转录组分析. 广东农业科学. 2021(02): 56-66 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 398
- HTML全文浏览量: 51
- PDF下载量: 57
- 被引次数: 4
