Effects of straw mulching on nutrient release and nutrient balance in potato soil
-
摘要:目的
探索马铃薯种植过程中稻草覆盖还田模式下稻草N、P、K释放和土壤N、P、K养分平衡的特征。
方法采用盆栽试验的方法，对稻草覆盖还田后的腐解速率及养分释放状况进行研究。
结果稻草覆盖还田82 d后，其C、P含量变化较小，K含量降低93.96%，而N含量却提高了49.22%；覆盖稻草可为土壤提供大量养分且以K居多，每覆盖还田100 kg稻草可为马铃薯生长提供N 0.42 kg、P 0.14 kg(P2O5)和K 1.70 kg(K2O)；稻草腐解率为70.87%；各养分累积释放率表现为：K>C>P>N，C、N和P累积释放率分别为71.89%、56.85%和68.75%，K累积释放率最高，为98.22%。覆盖稻草和施用钾肥对土壤K平衡影响较大，不施钾肥和稻草土壤的K亏缺达163.80 kg/hm2，覆盖稻草或施用钾肥均可有效降低土壤钾素亏缺量，稻草与钾肥配合施用可使土壤K收支平衡，甚至盈余。
结论在稻草覆盖还田情况下，可适当减少钾肥用量，以降低生产成本，缓解我国钾肥资源匮乏的困境。
Abstract:ObjectiveTo study the characteristics of nitrogen(N), phosphorus(P) and potassium (K) nutrient release and soil N, P, K nutrient balance under straw mulching pattern during potato planting process.
MethodWe studied the decomposition rate and nutrient release status of mulching straw using pot experiment.
ResultStraw carbon (C) and P contents changed slightly after mulching for 82 days, while K content decreased by 93.96%, and N content increased by 49.22%. The mulching straw provided a large amount of nutrients (especially K) with 100 kg mulching straw contributing 0.42 kg N, 0.14 kg P (P2O5), and 1.70 kg K (K2O) to soil used for potato growth. Straw decomposition rate was 70.87%. The nutrient accumulative release rates were performed as: K (98.22%) > C (71.89%) > P (68.75%) > N (56.85%). Straw mulching and K fertilization largely affected soil K balance. The K deficient amount was 163.80 kg/hm 2 without straw mulching and K fertilizer application. K deficient amount was effectively reduced by straw mulching and K fertilization, and even K had surplus with combined application of straw and K fertilizer.
ConclusionThe straw mulching pattern can reduce K fertilizer application amount and production cost, and relieve the dilemma of K fertilizer resource deficiency in China.
-
Keywords:
- potato /
- nutrient release /
- nutrient balance /
- straw mulching /
- decomposition
-
转录因子能够调控基因的转录，在发育过程中发挥着重要的作用，并且转录因子的突变和许多疾病的发生有着紧密的联系[1]。E26转录因子(E26 transcription-specific，ETS)家族是转录因子家族之一，它们在细胞的增殖、分化、迁移以及细胞的互作过程中起着重要的作用，而根据序列的相似性，又将ETS家族分成几个亚家族，在哺乳动物中发现了将近20多个亚家族，其中包括ETS、ELG、TEL、PEA3等[2]。PEA3亚家族由3个成员组成：ETV1、ETV4和ETV5[3]。ETV5也称ETS相关分子(ETS related molecule，ERM)，在成年哺乳动物的组织中有着广泛的表达，尤其在睾丸[4]和卵巢[5]中表达量较高。
精原干细胞是雄性动物产生精子的基础，精原干细胞的研究对解决雄性不育症以及畜禽育种有着重要的推动作用，精原干细胞移植是近年来比较难以攻克的研究问题，其中如何制备内源性精原细胞消融的受体就是一个棘手的问题，传统方法常将白消安注射到雄性受体的睾丸中使其内源性的精原细胞消融[6]，但是此方法对动物的伤害很大，并且雄性受体睾丸的精原细胞也无法做到完全消融，不可避免地会有一定程度的恢复，对后续试验的影响很大。随着基因编辑技术的发展，研究者常利用基因修饰的手段敲除或失活某种在雄性哺乳动物睾丸精子发生过程中起关键作用的基因，进而获得雄性睾丸内源性精原细胞消融理想的移植受体。Chen等[7]利用Cre-LoxP技术靶向敲除了小鼠的ETV5基因，制备了ETV5 纯合敲除（ETV5-/-）小鼠，并发现敲除小鼠内源性精原细胞发生了消融且呈现出“唯支持细胞综合征”的表型，成功地制备了内源性精原细胞消融的移植受体。同时研究发现，ETV5-/-雄性小鼠对雌性小鼠失去了性欲，精子也没有受精能力，无论是自然交配、人工授精还是体外受精都不能使野生型(WT)的雌性小鼠受孕。而且，有报道称ETV5-/-雌性小鼠也是不孕的[8]，这是因为它们的卵巢出现了结构上的缺陷，完全无法排卵，和野生型公鼠交配后不能形成阴道栓。ETV5敲除小鼠的成功制备实现了内源性精原细胞完全消融的雄性受体，为精原干细胞的研究提供了重要的材料。
然而，ETV5敲除小鼠不仅生殖系统的发育受到阻碍，在其他身体发育机能上也有缺陷。Schlesser等[9]靶向干扰小鼠的ETV5基因，和野生型相比ETV5纯合突变小鼠的体质量减少12%~31%。Jamsai等[10]利用ENU诱变将小鼠的ETV5基因产生错义突变，纯合ETV5突变的小鼠除了出现“唯支持细胞综合征”外，还产生了肾脏不对称、多趾以及发育不良等症状。Zhang等[11]为研究精原干细胞的移植，利用CRISPR/Cas9技术制备了ETV5敲除小鼠，试验结果和Schlesser等[9]以及Jamsai等[10]的研究结果相似，敲除ETV5基因的小鼠除了出现内源性精原细胞的消融，也出现了生长发育缓慢以及多趾的现象。
为了探究ETV5敲除小鼠生长发育缓慢的分子机制，我们通过对ETV5纯合敲除公鼠和野生型公鼠的肌肉组织进行了转录组测序，获得了2组小鼠肌肉组织的转录本数据，并且通过生物信息学分析筛选了ETV5−/−公鼠和野生型公鼠肌肉组织的差异表达基因，以及敲除ETV5影响了小鼠其他哪些功能基因及信号通路。
1. 材料与方法
1.1 ETV5敲除小鼠的制备及基因型的检测
ETV5敲除小鼠的制备是通过CRISPR/Cas9技术对C57BL/6小鼠的ETV5基因(GenBank accession number: NM_023794.2；Ensembl: ENSMUSG00000013089) 进行靶向敲除得到的(由苏州赛业生物科技有限公司完成)。设计了ETV5敲除小鼠的sgRNA，包括sgRNA1（GAACGGCCATTGTCGGTGGCAGG）和sgRNA2（CTTCTATGCTAATAACGGGTGGG）共获得49只F0代小鼠，1周龄时剪小鼠脚趾抽提DNA，通过PCR检测筛选出3只F0代ETV5基因杂合敲除小鼠(2公1母)。将杂合子公鼠和母鼠进行合笼繁殖，留下它们的公鼠后代，采集这些公鼠的尾巴，抽提组织DNA，并进行基因型的PCR鉴定，鉴定小鼠基因型的PCR引物包括上游引物（Primer-F: 5′-CAACTGGTGCCCTTCCCAGTCT-3′）、中间引物（Primer-M: 5′-TTAGACATTAGGCAGAGCCAGTGATG-3′）、下游引物（Primer-R: 5′-GCCGCTCTTAAACCTGTTCATTCG-3′）。
1.2 小鼠的体质量测量
留下基因型检测结果为野生型(WT)的公鼠(C57BL/6小鼠)和纯合敲除(ETV5-/-)的公鼠，待小鼠6周龄时，选择WT公鼠和ETV5-/-公鼠各3只，依次放在电子秤上称量6只小鼠的体质量并记录。
1.3 样品收集和准备
1.3.1 肌肉样本的采集
称量过体质量的6只小鼠，分别进行颈椎脱臼法处死，将小鼠解剖采取小鼠后腿部位的肌肉，分别放入标记好的冻存管中，立即置于液氮中。
1.3.2 RNA 提取与检测
采用Trizol法提取小鼠肌肉组织的总RNA，用琼脂糖凝胶电泳检测分析样品RNA的完整性，用Nanodrop 2000检测RNA浓度和纯度，用Agilent 2100 bioanalyzer精确检测RNA 完整性。
1.4 文库的构建及测序
转录组测序由北京诺禾致源科技股份有限公司完成，将2组小鼠肌肉的总RNA分别构建测序文库。首先通过磁珠法Oligo(dT)对带有PolyA尾的mRNA进行富集，随后将其分割成短片段mRNA，用随机引物进行扩增，逆转录合成第1条cDNA 链，随后采用RNaseH、DNA Polymerase I、dNTPs为原料降解RNA并且合成第2条cDNA链。使用AMPure XP beads 纯化双链产物，最后将DNA末端修复成平末端，加PolyA尾、接头，最终获得文库。文库质量合格后进行Illumina测序。
1.5 数据分析
1.5.1 数据质控
将测序片段的图像数据转化为序列数据(Reads)，对原始数据进行过滤，去除带接头的、无法确定碱基信息的、低质量的Reads。同时，对Clean data进行Q20、Q30和GC含量计算，供后续分析。
1.5.2 RNA比对及基因表达水平的测定
获得基础数据后，直接从基因组网站下载参考基因组和基因模型注释文件，使用HISAT2构建参考基因组的索引，并将Clean data与参照基因组进行比对，HISAT2可以基于基因模型注释文件生成拼接连接的数据库。利用featureCounts软件计算映射到每个基因的读数，再根据基因的长度计算每个基因的FPKM(FPKM指每百万碱基对测序的转录本序列片段的每千碱基片段的预期数量)，并计算映射到该基因的读数。
1.5.3 差异表达分析
使用DESeq2 软件进行2个组合之间差异表达分析[12](每个组3个生物学重复)。差异表达基因(Differentially expressed gene,DEG)是通过DESeq2 选择阈值为Padj<0.05，|log2(差异倍数)|≥0筛选出来的基因(Padj是Benjamini和Hochberg法矫正后的P值)。通过clusterProfiler软件使用Gene Ontology数据库对差异表达基因进行GO富集分析。利用KEGG(Kyoto encyclopedia of genes and genomes)数据库找出差异表达基因显著富集的信号通路[13]，用于分析分子水平的信息，使用clusterProfiler软件分析KEGG 通路中差异表达基因的统计富集。
1.5.4 实时荧光定量PCR验证
从差异表达基因中挑选了一部分进行了实时荧光定量PCR(qPCR)验证。反应体系为10 μL：5 μL的2×SYBR Green MasterMix，1 μL cDNA模板，0.3 μL 上、下游引物和3.4 μL ddH2O。反应条件为：95 ℃10 min；95 ℃15 s，60 ℃15 s，72 ℃10 s，共40个循环；72 ℃10 min。每个基因做3个重复，以β-actin作为内参基因，用2−ΔΔCt方法计算基因相对表达量。
2. 结果与分析
2.1 小鼠的基因型鉴定
如图1所示，如果小鼠是野生型小鼠，经过PCR能根据引物Primer-M和Primer-R得到460 bp长度的片段，如果小鼠是ETV5纯合敲除的小鼠，则扩增出Primer-F到Primer-R之间的1段710 bp的序列，同理，如果是杂合敲除的小鼠，则2种片段都有。因此，可以根据PCR产物的大小判断小鼠的基因型。
![]()
图 1 ETV5敲除小鼠基因型鉴定原理图Figure 1. Identification principle of ETV5 knockout mouse genotype选取部分后代小鼠中的公鼠，提取组织DNA后，用PCR法检测基因型，结果如图2所示，泳道2、3、4只扩增到1个条带(460 bp)，说明所代表的这3只小鼠是野生型，泳道5、6、7只扩增到1个条带(710 bp)，说明所代表的这3只小鼠是纯合敲除，泳道8、9、10扩增到2个条带(710和460 bp)，说明所代表的这3只小鼠是杂合敲除。
![]()
图 2 ETV5基因敲除小鼠基因型鉴定1： Trans DNA Marker II，2~4：野生型，5~7：ETV5纯合敲除，8~10：ETV5杂合敲除Figure 2. Genotype identification ofETV5 knockout mice1: Trans DNA Marker II, 2−4: Wild type, 5−7: ETV5 homozygous knockout, 8−10: ETV5 hybrid knockoutETV5纯合敲除小鼠的PCR产物测序结果如图3所示，将测序结果和NCBI上的小鼠ETV5基因序列进行比对，结果显示，从22253855~22259641缺失了5607个碱基。
![]()
图 3 ETV5纯合敲除小鼠测序结果和NCBI序列的比对Figure 3. Comparison of sequencing results of ETV5 homozygous knockout mice and NCBI sequence2.2 2组小鼠的体质量比较
在2组公鼠42日龄时，对它们的体质量进行测量，得到结果如图4所示，WT组的体质量显著高于ETV5-/-组的体质量。
![]()
图 4 ETV5基因敲除与野生型小鼠在42日龄时的体质量“*”表示处理间差异显著(t检验)Figure 4. The body weight of the ETV5 knockout mice and wild type mice at 42 days of age“*” indicates significant difference between treatments(t test)2.3 参考基因组比对结果
根据比对结果，分别统计Reads在基因组外显子区域(Exon)、内含子区域(Intron)以及基因间区(Intergenic)所占的比例。WT和ETV5−/−组的比对结果基本一致，大约95%的Reads比对到外显子区域。
2.4 差异表达基因
我们筛选出了574个差异表达基因，其中上调基因292个，下调基因282个。对这574个差异表达基因进行染色体定位，结果如图5所示，发现差异表达基因在11号染色体最多，其次是5号染色体。图6的热图显示筛选出的差异表达基因的上、下调水平。部分差异表达基因见表1。
表 1 部分差异表达基因Table 1. Partial DEGs上调基因
Up-regulated
geneFC1) 下调基因
Down-regulated
geneFC1) Chrna9 7.7151 Acaa2 −0.9966 Cdsn 6.9906 Ccdc85a −1.0183 Adrb3 6.8571 Homer2 −1.0468 Retn 4.8669 Crispld2 −1.0508 Slc7a10 4.2691 Rcan2 −1.0931 Fasn 3.9777 Rhobtb1 −1.1039 Acaa1b 3.8316 Khdrbs3 −1.1346 Cidec 3.5500 Spon1 −1.1809 Katnal2 3.0976 Nudt18 −1.2179 Elovl6 3.0859 Rps6ka1 −1.2612 Spon2 3.0029 Gdap1 −1.2711 Adig 2.8624 Nceh1 −1.3692 Rbp4 2.8034 Map2k6 −1.4616 Eda2r 2.5549 Grb14 −1.5417 Agpat2 2.4480 Dbp −1.7583 Smoc1 2.4085 Plin5 −1.7990 Igf2 2.2913 Hpn −1.9161 Ankrd1 2.1497 Smox −2.0291 Meg3 2.0777 Fabp3 −2.0493 B3galt2 2.0766 Mylk4 −2.0812 Slc1a3 2.0576 Mchr1 −2.1755 Perp 2.0052 Fabp3-ps1 −2.1760 Peg3 1.9896 Cacna2d4 −2.2257 Igfn1 1.9846 Amd2 −2.6730 Dapk1 1.8872 Bdh1 −2.9090 Pde3b 1.7084 Wisp2 −2.9357 Tbc1d31 1.6355 Amd-ps4 −2.9575 Fst 1.3919 Amd1 −3.0079 Thbs1 1.3821 Gapdh −3.7085 Il20rb 1.3807 Chrna2 −4.8778 Cd28 1.3191 Fam19a4 −5.1516 Apoe 0.9978 Etv5 −6.8862 1) FC表示log2 (差异倍数)
1) FC indicates log2 (Fold change)![]()
图 6 小鼠肌肉组织中的差异表达基因聚类热图红色代表上调，绿色代表下调Figure 6. Clustering heat map of DEGs in mouse muscle tissueRed color represents up-regulation, and green color represents down-regulation2.5 GO富集分析
为了进一步解析筛选出来的差异表达基因的功能，我们进行了GO富集分析，结果如图7所示。富集在生物学过程的最多，其中在生物过程中富集程度前5的都与脂肪代谢相关，包括：脂肪酸代谢、脂质氧化、脂质分解代谢、脂肪细胞分化、脂肪酸氧化。显著富集的细胞组成相对较少，只有5个，分别是：线粒体内膜、细胞器内膜、脂滴、细胞外基质、受体复合体。前5个显著富集的分子功能分别是：辅酶结合、生长因子结合、硫化合物结合、脂肪−酰基−辅酶A结合、胰岛素样生长因子结合。
![]()
图 7 ETV5基因敲除与野生型小鼠肌肉组织差异表达基因GO富集分析Figure 7. Enrichment analysis of DEGs betweenETV5 knockout mice and wild type mice muscle2.6 KEGG富集分析
我们还对差异表达基因进行了KEGG富集分析，图8结果表明有20条KEGG通路得到显著富集，其中差异表达基因主要富集到代谢通路中，例如：脂肪酸代谢、cAMP信号通路、PPAR信号通路、AMPK信号通路、不饱和脂肪酸生物合成、脂肪酸降解和丙酮酸代谢等。
![]()
图 8 ETV5基因敲除与野生型小鼠肌肉组织差异表达基因的KEGG富集分析Figure 8. KEGG enrichment analysis of DEGs between ETV5 knockout mice and wild type mice muscle2.7 qPCR验证分析
将分析得到的部分差异表达基因进行了qPCR验证，结果显示与测序结果一致(图9)，敲除ETV5基因影响了部分其他基因的表达，其中腺苷甲硫氨酸脱羧酶1(Adenosylmethionine decarboxylase 1，Amd1)[14]的cDNA编码肌源性碱性螺旋−环−螺旋(Basic helix-loop-helix，bHLH)因子的一个家族成员，此外，Amd1基因在成人体壁肌肉中有明显表达，但在心脏和其他非肌肉组织中未见表达，肌源性bHLH因子参与了肌肉发育和形成。另外，胆碱能受体烟碱α2(Cholinergic receptor nicotinic alpha 2，Chrna2)对脂肪的形成和代谢有重要作用，Jun等[15]敲除了小鼠的Chrna2基因，发现由于小鼠体内脂肪细胞自主调节异常，导致了小鼠在寒冷环境中产热和代谢功能的障碍，揭示了Chrna2通路在米色脂肪生物发生和能量稳态中的生物学意义。除此之外，差异表达基因中还有Hpn[16]，Wisp2[17]，Cxcl14[18]，Mylk4[19]，Plin5[20]，Mchr1[21]等一些基因影响了ETV5敲除小鼠的肌肉发育、脂肪累积和体质量变化。
![]()
图 9 qPCR检验差异表达基因在ETV5基因敲除和野生型小鼠肌肉组织的相对表达量柱子上方的“*”“**”“***”分别表示2组差异达到0.05、0.01和0.001的显著水平(t检验)Figure 9. Detecting the relative expression of DEGs in muscle tissue ofETV5 knockout and wild type mice by qPCR“*”“**”“***”on bars indicate significant difference between two groups at the levels of 0.05, 0.01 and 0.001, respectively (t test)3. 讨论与结论
Wang等[22]前期通过基因编辑技术制备了ETV5基因敲除的小鼠模型用来研究该基因与精原干细胞发育的关系，但是基因编辑带来的影响并不是单一的，虽然敲除小鼠的ETV5基因的确得到了内源性精原细胞消融的“唯支持细胞综合征”表型，但是同时也影响了小鼠的整体发育和生长，而产生这些表型的分子机制很少有人报道。为此，我们通过转录组测序得到了ETV5敲除小鼠和野生型小鼠的转录本，分析两者之间存在的差异表达基因，通过生物信息学的方法分析得到了敲除ETV5基因可能对小鼠产生影响的分子机制。结果表明，敲除ETV5后的小鼠与野生型小鼠相比差异表达基因大部分为编码蛋白的基因，差异表达基因中的部分基因对动物的生长发育、肌肉发育以及脂肪累积等有不同的影响，从而影响了ETV5敲除小鼠的肌肉发育和脂肪累积，出现体质量下降以及发育不良的现象。
GO富集分析的结果显示，富集前5的生物学过程都与脂肪代谢相关，分别是：脂肪酸代谢、脂质氧化、脂质分解代谢、脂肪细胞分化、脂肪酸氧化。这表明，敲除ETV5基因对小鼠的脂肪代谢影响很大，可能导致了小鼠储存脂肪的障碍。KEGG富集分析的结果显示，所有差异表达基因主要富集在脂肪酸代谢、cAMP信号通路、PPAR信号通路、AMPK信号通路等方面。PPAR通路可以调节脂质代谢、脂肪形成、维持代谢的稳态以及炎症基因的表达。PPAR在不同的物种中找到了3种亚型，PPARα通过调节参与肝脏和骨骼肌脂质代谢基因的表达以在细胞脂质中发挥作用[23]，PPARβ/δ参与脂质分解代谢、葡萄糖稳态、炎症、生存、增殖、分化以及哺乳动物皮肤、骨骼和肝脏的再生[24]，PPARγ可调节糖和脂代谢、内皮功能和炎症[25]。测序结果中的Plin5、Fabp3等部分差异表达基因在PPAR通路中起作用。GO和KEGG分析的结果表明，敲除ETV5基因可能会影响小鼠的脂肪代谢、能量代谢以及整体的发育。
本研究对ETV5敲除小鼠和野生型小鼠的肌肉组织进行了转录组测序分析，以探究敲除ETV5基因对小鼠肌肉形成以及生长发育相关基因和信号通路的影响。结果表明，我们筛选出了574个差异表达基因，这些差异表达基因大多与脂肪代谢以及生长发育相关，其中发现Amd1可能影响了ETV5敲除小鼠的肌肉发育，Chrna2可能影响了ETV5敲除小鼠的脂肪代谢。GO和KEGG分析显示差异表达基因主要富集在脂肪代谢、能量代谢的通路，这也能对应我们观察到ETV5敲除小鼠发育不良的表型。这些结果为敲除ETV5基因导致小鼠发育的缓慢和延迟提供了解释，以及为研究ETV5基因的体内多样性功能提供了参考。
-
![]()
图 1 稻草钾肥配施对稻草C、N、P和K含量的影响
各图中，相同指标柱子上方的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
Figure 1. Influences of combined application of straw and K fertilizer on C, N, P, K contents of rice straw
In each figure, different lowercase letters on bars of the same indicator indicated significant difference(P<0.05, Duncan’ s test)
![]()
图 2 稻草钾肥配施对稻草腐解率及C、N、P和K累积释放率的影响
相同指标柱子上方的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
Figure 2. Influences of combined application of paddy-rice straw and K fertilizer on paddy-rice decomposition rates and accumulative release rates of C, N, P, K
Different lowercase letters on bars of the same indicator indicated significant difference(P<0.05, Duncan’ s test)
表 1 盆栽试验各处理的施肥和覆盖稻草量
Table 1 Amounts of fertilizer application and straw mulching in pot experiments
kg·hm−2 处理
TreatmentN P2O5 K2O 稻草
Paddy-rice strawK0S0(CK) 160.7 36.9 0 0 K0S2 160.7 36.9 0 6 000 K1S2 160.7 36.9 99.8 6 000 K2S2 160.7 36.9 199.7 6 000 K3S2 160.7 36.9 299.5 6 000 K2S0 160.7 36.9 199.7 0 K2S1 160.7 36.9 199.7 3 000 K2S3 160.7 36.9 199.7 9 000 K1S1 160.7 36.9 99.8 3 000 
下载: 导出CSV
表 2 马铃薯收获时稻草养分累积释放量1)
Table 2 Accumulative release amounts of nutrients at harvesting period of potato
kg·hm−2 处理 Treatment C N P K K0S2 1 528.69±19.86b 25.76±0.14b 3.55±0.02b 102.32±0.01b K1S2 1 508.22±140.85b 26.60±3.96b 3.65±0.31b 102.35±0.39b K2S2 1 573.74±109.86b 27.39±2.30b 3.58±0.19b 102.35±0.23b K3S2 1 456.05±75.99b 24.39±1.96b 3.15±0.28b 101.47±0.48b K2S1 764.81±19.93c 12.23±0.69c 1.70±0.07c 51.07±0.12c K2S3 2 221.43±135.76a 36.58±2.90a 4.94±0.28a 152.74±0.46a K1S1 723.99±17.39c 11.13±0.64c 1.71±0.03c 51.03±0.02c 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 3 稻草钾肥配施对土壤N平衡的影响1)
Table 3 Influences of combined application of paddy-rice straw and K fertilizer on soil N balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学N
Chemical N稻草N
Straw NK0S0 160.65 0 89.20a 71.45±1.34e 1.80 80.19 K0S2 160.65 43.35 74.43b 129.57±3.62b 2.76 175.42 K1S2 160.65 43.35 88.63a 115.37±7.56bc 2.33 133.34 K2S2 160.65 43.35 73.63ab 130.37±4.12b 2.78 178.77 K3S2 160.65 43.35 77.41ab 126.59±9.37b 2.72 171.59 K2S0 160.65 0 83.97ab 76.68±0.29e 1.91 91.33 K2S1 160.65 21.68 82.01ab 100.31±1.54d 2.23 122.47 K2S3 160.65 65.03 75.40ab 150.28±3.78a 3.01 200.91 K1S1 160.65 21.68 80.42ab 101.90±2.54cd 2.27 127.15 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 4 稻草钾肥配施对土壤P平衡的影响1)
Table 4 Influences of combined application of paddy-rice straw and K fertilizer on soil phosphorus balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学P
Chemical P稻草P
Straw PK0S0 16.11 0 13.53a 2.58±0.07d 1.19 19.09 K0S2 16.11 5.0 12.77a 8.33±1.33ab 1.69 50.60 K1S2 16.11 5.0 15.01a 6.1±1.18bc 1.42 42.32 K2S2 16.11 5.0 11.47a 9.64±1.76a 1.95 59.84 K3S2 16.11 5.0 12.77a 8.34±1.59ab 1.71 49.99 K2S0 16.11 0 12.87a 3.24±0.78cd 1.26 26.08 K2S1 16.11 2.5 13.07a 5.54±0.45bcd 1.43 42.74 K2S3 16.11 7.5 12.98a 10.63±0.48a 1.82 82.34 K1S1 16.11 2.5 12.48a 6.13±0.77bc 1.50 50.21 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 5 稻草钾肥配施对土壤K平衡的影响1)
Table 5 Influences of combined application of paddy-rice straw and K fertilizer on soil potassium balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学K
Chemical K稻草K
Straw KK0S0 0 0 163.80b −163.80±2.36f 0 −100.00 K0S2 0 86.60 177.01ab −90.40±13.68e 0.49 −50.64 K1S2 83.19 86.60 214.60a −44.81±16.21d 0.80 −20.26 K2S2 166.38 86.60 178.28ab 74.70±23.28b 1.45 45.39 K3S2 249.56 86.60 179.59ab 156.58±28.92a 1.95 94.73 K2S0 166.38 0 186.66ab −20.28±10.75cd 0.90 −10.47 K2S1 166.38 43.30 207.44ab 2.24±5.09c 1.01 1.16 K2S3 166.38 129.90 202.50ab 93.78±4.32b 1.46 46.40 K1S1 83.19 43.30 183.02ab −56.53±13.26de 0.70 −30.41 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
-
[1] 屈冬玉. 马铃薯产业应在我国农村区域经济发展中发挥更大作用[C]//中国马铃薯作物学会马铃薯专业委员会. 马铃薯产业与农村区域发展. 北京: 中国作物学会, 2013: 1-3. [2] 谢开云, 屈冬玉, 金黎平, 等. 中国马铃薯生产与世界先进国家的比较[J]. 世界农业, 2008, 29(5): 35-38. doi: 10.3969/j.issn.1002-4433.2008.05.012 [3] 米健, 罗其友, 高明杰. 南方冬作区马铃薯发展趋势、区域格局与增产潜力[J]. 中国农业资源与区划, 2012, 33(3): 73-79. [4] 屈冬玉, 金黎平, 谢开云. 中国马铃薯产业10年回顾[M]. 北京: 中国农业科学技术出版社, 2010: 66-69. [5] 高媛, 韦艳萍, 樊明寿. 马铃薯的养分需求[J]. 中国马铃薯, 2011, 25(3): 182-187. doi: 10.3969/j.issn.1672-3635.2011.03.018 [6] MULDER A. Potato diseases-diseases, pests and defects[M]. Holland: Aardappelwereld B V and NIVAP, 2005: 209-221.
[7] GEORGE M S, LU G, ZHOU W. Genotypic variation for potassium uptake and utilization efficiency in sweet potato (Ipomoea batatas L.)[J]. Field Crop Res, 2002, 77(1): 7-15. doi: 10.1016/S0378-4290(02)00043-6
[8] 白由路. 高价格下我国钾肥的应变策略[J]. 中国土壤与肥料, 2009, 45(3): 1-4. doi: 10.3969/j.issn.1673-6257.2009.03.001 [9] 陈杨, 樊明寿, 李斐, 等. 氮素营养诊断技术的发展及其在马铃薯生产中的应用[J]. 中国农学通报, 2009, 25(3): 66-71. [10] 张新明, 张洪秀, 李水源, 等. 恩平市典型种植户冬作马铃薯施肥状况调查分析[J]. 安徽农业科学, 2011, 39(36): 22286-22288. doi: 10.3969/j.issn.0517-6611.2011.36.034 [11] LAL R. World crop residues production and implications of its use as a biofuel[J]. Environ Int, 2005, 31(4): 575-584. doi: 10.1016/j.envint.2004.09.005
[12] 毕于运, 王红彦, 王道龙, 等. 中国稻草资源量估算及其开发利用[J]. 中国农学通报, 2011, 27(15): 137-143. [13] 李国学, 张福锁. 固体废物堆肥化和有机复混肥的生产[M]. 北京: 化学工业出版社, 2000: 1-20. [14] 曹国良, 张小曳, 王亚强, 等. 中国区域农田秸秆露天焚烧排放量的估算[J]. 科学通报, 2007, 52(15): 1826-1831. doi: 10.3321/j.issn:0023-074x.2007.15.017 [15] 谭德水, 金继运, 黄绍文, 等. 施钾和秸秆还田对栗钙土区土壤养分及小麦产量的长期效应研究[J]. 干旱地区农业研究, 2009, 27(2): 194-198. [16] 张静, 温晓霞, 廖允成, 等. 不同玉米秸秆还田量对土壤肥力及冬小麦产量的影响[J]. 植物营养与肥料学报, 2010, 16(3): 612-619. doi: 10.11674/zwyf.2010.0314 [17] CARUSO G, CARPUTO D, CONTI S, et al. Effect of mulching and plant density on out-of-season organic potato growth, yield and quality[J]. Adv Hortic Sci, 2013, 27(3): 115-121.
[18] DVORAK P, KUCHTOVA P, TOMASEK J. Response of surface mulching of potato (Solanum tuberosum) on SPAD value, colorado potato beetle and tuber yield[J]. Int J Agric Biol, 2013, 15(4): 798-800.
[19] 农业部办公厅. 农业部关于印发《测土配方施肥技术规范(2011年修订版)》的通知: 农农发[2011]3号[A].(2011-09-22)[2018-11-02]. http://jiuban.moa.gov.cn/zwllm/tzgg/tz/201109/t20110922_2293389.htm. [20] 鲍士旦. 土壤融化分析[M]. 3版. 北京: 中国农业出版社, 2000. [21] 王志勇, 白由路, 杨俐苹, 等. 低土壤肥力下施钾和秸秆还田对作物产量及土壤钾素平衡的影响[J]. 植物营养与肥料学报, 2012, 18(4): 900-906. [22] TANG Q, ZHANG C. Data processing system (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research[J]. Insect Sci, 2013, 20(2): 254-260. doi: 10.1111/j.1744-7917.2012.01519.x
[23] 闫超. 水稻秸秆还田腐解规律及对土壤养分和酶活性的影响[D]. 哈尔滨: 东北农业大学, 2012. [24] 李世清, 李生秀. 有机物料在维持土壤微生物体氮库中的作用[J]. 生态学报, 2001, 21(1): 136-142. doi: 10.3321/j.issn:1001-9332.2001.01.043 [25] ABRO S A, TIAN X H, WANG X D, et al. Decomposition characteristics of maize (Zea mays L.) straw with different carbon to nitrogen (C/N) ratios under various moisture regimes[J]. Afr J Biotechnol, 2011, 10(50): 10149. doi: 10.5897/AJB10.2261
[26] SUN B, WANG X, WANG F, et al. Assessing the relative effects of geographic location and soil type on microbial communities associated with straw decomposition[J]. Appl Environ Microbiol, 2013, 79(11): 3327-3335. doi: 10.1128/AEM.00083-13
[27] LU P, LIN Y H, YANG Z Q, et al. Effects of application of corn straw on soil microbial community structure during the maize growing season[J]. J Basic Microbiol, 2015, 55(1): 22-32. doi: 10.1002/jobm.201300744
[28] 李逢雨, 孙锡发, 冯文强, 等. 麦秆、油菜秆还田腐解速率及养分释放规律研究[J]. 植物营养与肥料学报, 2009, 15(2): 374-380. doi: 10.3321/j.issn:1008-505X.2009.02.018 [29] 戴志刚, 鲁剑巍, 李小坤, 等. 不同作物还田秸秆的养分释放特征试验[J]. 农业工程学报, 2010, 26(6): 272-276. doi: 10.3969/j.issn.1002-6819.2010.06.047 [30] 张晋京, 窦森. 玉米秸秆分解期间胡敏酸、富里酸动态变化的研究[J]. 土壤通报, 2005(1): 134-136. doi: 10.3321/j.issn:0564-3945.2005.01.037 [31] 张亮. 关中麦玉轮作区施氮对秸秆还田小麦产量和秸秆养分释放的影响[D]. 杨凌: 西北农林科技大学, 2012. [32] 鲁如坤, 时正元, 施建平. 我国南方6省农田养分平衡现状评价和动态变化研究[J]. 中国农业科学, 2000, 33(2): 63-67. doi: 10.3321/j.issn:0578-1752.2000.02.010 [33] 鲁如坤, 刘鸿翔, 闻大中, 等. 我国典型地区农业生态系统养分循环和平衡研究:Ⅳ:农田养分平衡的评价方法和原则[J]. 土壤通报, 1996, 27(5): 197-199. [34] 谭德水, 金继运, 黄绍文, 等. 不同种植制度下长期施钾与秸秆还田对作物产量和土壤钾素的影响[J]. 中国农业科学, 2007, 40(1): 133-139. doi: 10.3321/j.issn:0578-1752.2007.01.018
计量
- 文章访问数: 1393
- HTML全文浏览量: 3
- PDF下载量: 1346
