中国农业科学 ›› 2025, Vol. 58 ›› Issue (18): 3583-3597.doi: 10.3864/j.issn.0578-1752.2025.18.002
李明1(), 程宇坤1,2, 白斌3, 雷斌4, 耿洪伟1,2
收稿日期:
2025-03-12
接受日期:
2025-05-06
出版日期:
2025-09-18
发布日期:
2025-09-18
通信作者:
联系方式:
李明,E-mail:1304712256@qq.com。
基金资助:
LI Ming1(), CHENG YuKun1,2, BAI Bin3, LEI Bin4, GENG HongWei1,2
Received:
2025-03-12
Accepted:
2025-05-06
Published:
2025-09-18
Online:
2025-09-18
摘要:
【目的】穗部性状是影响小麦产量的重要因素,通过对小麦穗部性状进行全基因组关联分析,旨在发掘控制小麦穗部性状显著位点,为小麦穗部性状的遗传改良研究提供理论参考。【方法】以261份冬小麦品种(系)为材料,测定小麦穗部表型性状,结合小麦90K SNP芯片,采用固定随机循环概率模型(Farm CPU)进行穗部性状全基因组关联分析,并对筛选出稳定且显著的位点进行单倍型分析。【结果】在3个环境下,11个穗部性状均表现出广泛的表型变异,变异系数为3.63—64.29,各穗部性状遗传力为0.42—0.84,且基因型、环境、基因型×环境间均呈现出极显著差异。通过全基因组关联分析,共检测到171个与11个性状显著关联的位点(P<0.001),其中,20个关联位点在2个及以上的环境中被检测到,分别与穗长(3个)、穗下节长(7个)、不育小穗数(1个)、可育小穗数(2个)、总小穗数(2个)、穗粒数(1个)、穗粒重(2个)和千粒重(2个)等8个穗部性状关联,表型贡献率为0.95%—18.54%。在7B染色体上检测到1个与穗粒重和穗粒数显著关联的多效位点Ra_c10072_677,贡献率为2.62%—6.16%。筛选在2个以上环境条件下均能检测到与穗下节长显著关联的wsnp_Ex_rep_c69639_68590556(可解释遗传变异的5.94%)标记,进行单倍型分析。共鉴定出3种单倍型Hap1、Hap2和Hap3,分布频率分别为77.40%、13.70%和8.80%。结合表型分析,含单倍型Hap3(30.58 cm)的261份冬小麦品种(系)穗下节长平均值显著高于(P<0.001)含Hap1(28.67 cm)和Hap2(27.49 cm)品种(系)的穗下节长。冬小麦单倍型分布频率呈显著差异,其中,在北部冬麦区Hap1单倍型的品种(系)占比较大,在黄淮冬麦区Hap2单倍型的品种(系)占比较大,而Hap3单倍型的品种(系)在所有冬麦区无较高出现频率。对3个环境下检测到稳定遗传的显著关联位点进行候选基因挖掘,筛选到4个与小麦穗部相关的候选基因。这些基因与MYB转录因子和F-box结构域等有关,可作为影响穗部性状的重要候选基因。【结论】小麦穗部性状在不同基因型间差异显著,在2个及以上环境中检测到20个稳定存在的关联位点,在7B染色体上鉴定到与穗下节长显著相关的3个不同单倍型,并筛选出4个与穗部性状相关的候选基因。
李明, 程宇坤, 白斌, 雷斌, 耿洪伟. 冬小麦穗部性状GWAS分析及优异单倍型筛选[J]. 中国农业科学, 2025, 58(18): 3583-3597.
LI Ming, CHENG YuKun, BAI Bin, LEI Bin, GENG HongWei. Genome-Wide Association Study on Spike Architecture Traits and Elite Haplotype Mining in Winter Wheat[J]. Scientia Agricultura Sinica, 2025, 58(18): 3583-3597.
表1
261份品种(系)小麦穗部性状统计分析"
性状 Trait | 环境 Environment | 范围 Range | 平均值±标准差 Mean±SD | 变异系数 CV (%) | 偏度 Skewness | 峰度 Kurtosis | 均方值Mean square | 遗传力 H2 | ||
---|---|---|---|---|---|---|---|---|---|---|
基因型 Genotype (G) | 环境 Environment (E) | G×E互作 G×E | ||||||||
穗长 SL (cm) | E1 | 5.61—11.14 | 8.39±1.03 | 12.23 | 0.12 | -0.01 | 4.80*** | 33.77*** | 1.02*** | 0.81 |
E2 | 6.08—11.97 | 8.66±1.03 | 11.93 | 0.54 | 0.65 | |||||
E3 | 5.97—12.90 | 8.9±1.14 | 12.80 | 0.58 | 0.61 | |||||
BLUP | 6.66—11.14 | 8.65±0.77 | 8.86 | 0.40 | 0.55 | |||||
穗下节长 PL (cm) | E1 | 20.09—45.14 | 29.48±4.74 | 16.80 | 0.43 | -0.23 | 97.01*** | 5072.93*** | 16.39*** | 0.84 |
E2 | 15.58—40.36 | 25.21±4.42 | 17.54 | 0.47 | 0.14 | |||||
E3 | 20.25—49.29 | 31.28±4.78 | 15.29 | 0.49 | 0.39 | |||||
BLUP | 22.80—36.51 | 28.66±2.88 | 10.04 | 0.38 | -0.34 | |||||
不育小穗数 SSN | E1 | 0.00—4.13 | 1.49±0.76 | 51.01 | 0.59 | 0.45 | 365.77*** | 51.27*** | 395.34*** | 0.55 |
E2 | 0.00—4.60 | 1.05±0.67 | 64.29 | 1.33 | 3.69 | |||||
E3 | 0.00—4.10 | 1.24±0.66 | 53.08 | 0.58 | 0.58 | |||||
BLUP | 0.65—1.98 | 1.26±0.27 | 21.32 | 0.29 | -0.32 | |||||
可育小穗数 FSN | E1 | 12.73—24.77 | 17.14±1.67 | 9.75 | 0.39 | 1.08 | 8.62*** | 88.72*** | 3.46*** | 0.66 |
E2 | 13.80—22.30 | 17.76±1.54 | 8.69 | 0.31 | -0.32 | |||||
E3 | 12.80—22.80 | 17.92±1.61 | 8.99 | 0.16 | 0.01 | |||||
BLUP | 15.78—20.06 | 17.61±0.82 | 4.66 | 0.22 | -0.16 | |||||
小穗结实率 SF (%) | E1 | 0.75—1.00 | 0.92±0.04 | 4.46 | -0.60 | 0.76 | 0.00*** | 0.08*** | 0.00*** | 0.55 |
E2 | 0.81—1.00 | 0.94±0.03 | 3.63 | -1.00 | 1.69 | |||||
E3 | 0.80—1.00 | 0.93±0.04 | 3.78 | -0.53 | 0.23 | |||||
BLUP | 0.89—0.96 | 0.93±0.01 | 1.53 | -0.24 | -0.22 | |||||
总小穗数 TSN | E1 | 15.17—25.13 | 18.63±1.56 | 8.38 | 0.22 | 0.41 | 8.31*** | 39.44*** | 3.14*** | 0.69 |
E2 | 15.80—24.00 | 18.81±1.6 | 8.48 | 0.36 | -0.25 | |||||
E3 | 14.50—24.40 | 19.17±1.52 | 7.93 | 0.15 | 0.17 | |||||
BLUP | 16.75—21.07 | 18.87±0.87 | 4.60 | -0.03 | -0.31 | |||||
穗粒数 GNS | E1 | 17.00—67.47 | 40.99±7.07 | 17.25 | 0.43 | 1.39 | 158.71*** | 3458.58*** | 78.16*** | 0.58 |
E2 | 29.70—70.70 | 45.32±7.34 | 16.19 | 0.59 | 0.28 | |||||
E3 | 26.60—68.80 | 45.57±7.33 | 16.08 | 0.33 | 0.06 | |||||
BLUP | 35.33—51.80 | 43.97±2.98 | 6.77 | 0.22 | -0.02 | |||||
穗粒重 GWPS (g) | E1 | 0.81—3.05 | 1.63±0.31 | 19.27 | 0.74 | 1.64 | 0.37*** | 25.83*** | 0.23*** | 0.53 |
E2 | 1.04—3.47 | 2.07±0.40* | 19.25 | 0.76 | 0.72 | |||||
E3 | 0.82—2.85 | 1.78±0.39 | 22.08 | 0.27 | -0.30 | |||||
BLUP | 1.53—2.16 | 1.83±0.11 | 6.20 | 0.26 | 0.12 | |||||
千粒重 TGW (g) | E1 | 27.31—51.72 | 39.3±4.63 | 11.77 | -0.47 | -0.04 | 91.76*** | 7452.80*** | 33.25*** | 0.72 |
E2 | 33.32—60.15 | 45.74±4.89 | 10.69 | 0.06 | -0.30 | |||||
E3 | 18.21—52.64 | 39.1±5.81 | 14.87 | -0.34 | -0.08 | |||||
BLUP | 34.04—47.27 | 41.43±2.71 | 6.54 | -0.42 | -0.35 | |||||
小穗着生密度 SD | E1 | 1.55—3.58 | 2.26±0.3 | 13.30 | 1.02 | 2.25 | 0.21*** | 0.80*** | 0.08*** | 0.68 |
E2 | 1.69—3.08 | 2.21±0.22 | 9.86 | 0.37 | 0.79 | |||||
E3 | 1.52—3.09 | 2.19±0.23 | 10.33 | 0.48 | 1.33 | |||||
BLUP | 1.86—2.66 | 2.22±0.14 | 6.26 | 0.38 | 0.49 | |||||
每小穗粒数 GNSL | E1 | 1.03—3.32 | 2.24±0.3 | 13.53 | -0.11 | 0.91 | 0.26*** | 3.85*** | 0.19*** | 0.42 |
E2 | 1.53—3.43 | 2.41±0.34 | 14.29 | 0.22 | 0.00 | |||||
E3 | 1.52—3.54 | 2.37±0.33 | 13.84 | 0.28 | 0.18 | |||||
BLUP | 2.09—2.54 | 2.34±0.08 | 3.52 | -0.11 | -0.23 |
表2
与小麦穗部性状显著关联的稳定位点"
性状 Trait | 标记 Marker | 染色体 Chr. | 位置 Position (Mb) | P值 P-value | 贡献率 R2(%) | 环境 Environment |
---|---|---|---|---|---|---|
穗长SL | wsnp_Ex_c19094_28015035 | 2B | 784.30 | 1.50E-04—4.06E-04 | 18.00—18.54 | E1/BLUP |
Tdurum_contig62138_385 | 2D | 45.80 | 7.59E-05—3.20E-04 | 5.57—5.89 | E2/BLUP | |
wsnp_Ex_c26128_35374652 | 3B | 440.89 | 9.17E-05—3.28E-04 | 5.64—7.15 | E1/E2/BLUP | |
穗下节长PL | Kukri_rep_c87640_135 | 3A | 573.99 | 6.43E-05—2.03E-04 | 5.81—9.23 | E1/E2 |
RAC875_c929_170 | 5B | 462.15 | 6.43E-05—8.20E-04 | 1.53—5.94 | E1/E2 | |
Tdurum_contig42655_1727 | 6B | 12.67 | 5.42E-07—3.09E-04 | 1.60—5.94 | E3/BLUP | |
Jagger_c3477_441 | 6D | 20.85 | 3.34E-04—3.51E-04 | 3.28—5.94 | E1/E2 | |
Kukri_rep_c70199_506 | 7A | 709.41 | 5.47E-04—8.35E-04 | 5.01—5.64 | E1/E2 | |
wsnp_Ex_rep_c69639_68590556 | 7B | 198.94 | 3.86E-04—9.04E-04 | 0.95—5.94 | E1/E2/BLUP | |
wsnp_Ex_c17176_25816766 | 7B | 700.39 | 3.48E-04—5.68E-04 | 5.94—9.32 | E1/E2 | |
不育小穗数SSN | Tdurum_contig42729_242 | 6A | 35.64 | 6.71E-04—7.92E-04 | 1.45—2.62 | E2/BLUP |
可育小穗数FSN | Kukri_c7218_1145 | 3A | 504.80 | 5.11E-04—5.21E-04 | 3.61—3.66 | E1/BLUP |
wsnp_Ex_c18318_27140346 | 4B | 12.53 | 5.90E-04—9.38E-04 | 1.05—3.55 | E1/BLUP | |
总小穗数TSN | wsnp_Ex_c1924_3630519 | 7A | 582.64 | 3.49E-04—6.84E-04 | 1.60—4.89 | E2/BLUP |
Excalibur_c87219_204 | 7B | 747.47 | 2.01E-04—3.95E-04 | 1.00—4.12 | E1/BLUP | |
穗粒数GNS | Ra_c10072_677 | 7B | 295.33 | 2.65E-04—6.77E-04 | 2.62—3.84 | E1/BLUP |
穗粒重GWPS | Excalibur_rep_c68296_2217 | 4A | 327.75 | 5.96E-04—8.19E-04 | 5.13—9.31 | E1/BLUP |
Ra_c10072_677 | 7B | 295.33 | 9.75E-05—2.60E-04 | 2.85—6.16 | E2/BLUP | |
千粒重TGW | BS00076889_51 | 1B | 22.15 | 3.82E-05—5.87E-04 | 1.16—4.61 | E1/BLUP |
RFL_Contig5418_347 | 3B | 567.25 | 3.63E-04—4.64E-04 | 1.09—4.30 | E1/BLUP |
[1] |
苏培森, 隋超, 郭尚敬, 臧娇娇, 满全财. 过氧化物酶基因TaPRX-2A调控小麦穗部性状的功能分析. 麦类作物学报, 2023, 43(3): 296-301.
|
|
|
[2] |
魏艳丽, 王彬龙, 李瑞国, 蒋会利, 张安静. 大穗小麦穗部性状的遗传分析. 麦类作物学报, 2015, 35(10): 1366-1371.
|
|
|
[3] |
doi: 10.3389/fpls.2018.01584 pmid: 30429867 |
[4] |
|
[5] |
|
[6] |
马艳明, 娄鸿耀, 王威, 孙娜, 颜国荣, 张胜军, 刘杰, 倪中福, 徐麟. 新疆冬小麦籽粒品质性状遗传差异与关联分析. 作物学报, 2024, 50(6): 1394-1405.
|
doi: 10.3724/SP.J.1006.2024.31030 |
|
[7] |
doi: 10.1186/s12870-019-1781-3 pmid: 31035920 |
[8] |
董继梓, 陈林渠, 郭浩儒, 张梦宇, 刘志霄, 韩磊, 田赵飒爽, 徐宁浩, 郭庆杰, 黄振洁, 等. 小麦穗长主效QTL—qSl-2D的遗传和育种选择效应解析. 中国农业科学, 2023, 56(20): 3917-3930. doi: 10.3864/j.issn.0578-1752.2023.20.001.
|
|
|
[9] |
丁浦洋, 周界光, 赵聪豪, 唐华苹, 牟杨, 唐力为, 邓梅, 魏育明, 兰秀锦, 马建. 小麦小穗数调控基因WAPO1的单倍型、遗传效应、地理分布及育种利用分析. 作物学报, 2022, 48(9): 2196-2209.
doi: 10.3724/SP.J.1006.2022.11078 |
|
|
[10] |
寇程, 高欣, 李立群, 李扬, 王中华, 李学军. 小麦粒重基因TaGW2-6A等位变异的组成分析及育种选择. 作物学报, 2015, 41(11): 1640-1647.
doi: 10.3724/SP.J.1006.2015.01640 |
|
|
[11] |
张颖, 石婷瑞, 曹瑞, 潘文秋, 宋卫宁, 王利, 聂小军. ICARDA引进-小麦苗期抗旱性的全基因组关联分析. 中国农业科学, 2024, 57(9): 1658-1681. doi: 10.3864/j.issn.0578-1752.2024.09.004.
|
|
|
[12] |
严勇亮, 时晓磊, 张金波, 耿洪伟, 肖菁, 路子峰, 倪中福, 丛花. 春小麦籽粒主要品质性状的全基因组关联分析. 中国农业科学, 2021, 54(19): 4033-4047. doi: 10.3864/j.issn.0578-1752.2021.19.001.
|
|
|
[13] |
常利芳, 李欣, 郭慧娟, 乔麟轶, 张树伟, 陈芳, 畅志坚, 张晓军. 小偃麦衍生系表型遗传多样性分析及综合评价. 草业学报, 2022, 31(11): 61-74.
doi: 10.11686/cyxb2022004 |
|
|
[14] |
|
[15] |
刘佩瑶, 冉莉萍, 杨佳庆, 王海博, 熊飞, 余徐润. 小麦穗形态建成和生理特征及外界影响因素的研究进展. 作物杂志, 2025(1): 1-9.
|
|
|
[16] |
刘阳娜, 刘丽华, 张风廷, 张立平, 苑少华, 张明明, 李宏博, 庞斌双, 赵昌平. 利用永久F2群体定位小麦穗部性状相关的QTL. 麦类作物学报, 2021, 41(5): 525-531.
|
|
|
[17] |
陈彩锦, 张尚沛, 师尚礼, 曾燕霞, 吴娟, 尚继红, 高婷, 张蓉. 基于GGE双标图对苜蓿品种丰产性和稳定性综合评价. 草地学报, 2021, 29(5): 912-918.
doi: 10.11733/j.issn.1007-0435.2021.05.007 |
|
|
[18] |
刘易科, 朱展望, 陈泠, 邹娟, 佟汉文, 朱光, 何伟杰, 张宇庆, 高春保. 基于SNP标记揭示我国小麦品种(系)的遗传多样性. 作物学报, 2020, 46(2): 307-314.
|
|
|
[19] |
张帅, 张希兰, 张娜, 赵明辉, 乔文臣, 孙丽静, 李辉, 傅晓艺, 何明琦, 纪军, 李俊明. 黄淮麦区北片小麦穗部相关性状的全基因组关联分析. 华北农学报, 2020, 35(S1): 31-39.
doi: 10.7668/hbnxb.20190900 |
doi: 10.7668/hbnxb.20190900 |
|
[20] |
王亚飞, 李世景, 徐萍, 张正斌, 景蕊莲. 黄淮和长江中下游冬麦区小麦品种(系)农艺性状及其聚类分析. 中国生态农业学报(中英文), 2020, 28(3): 395-404.
|
|
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
周淼平, 张鹏, 杨学明, 张平平, 宋桂成, 何漪. 扬麦158/西风重组自交系群体千粒重QTL的初步定位. 江苏农业学报, 2025, 41(1): 1-8.
|
|
|
[31] |
|
[32] |
陈玲玲, 刘亭萱, 谷勇哲, 宋健, 王俊, 邱丽娟. 大豆叶柄夹角相关基因GmILPA1单倍型分析. 植物遗传资源学报, 2021, 22(6): 1698-1707.
doi: 10.13430/j.cnki.jpgr. 20210419003 |
|
|
[33] |
朱新开, 郭文善, 李春燕, 封超年, 彭永欣. 小麦株高及其构成指数与产量及品质的相关性. 麦类作物学报, 2009, 29(6): 1034-1038.
|
|
|
[34] |
黄杰, 周璐琪, 葛昌斌, 王君, 曹燕燕, 宋丹阳, 廖平安. 小麦穗下节性状与灌浆特性、穗部相关因素的关系. 江苏农业科学, 2024, 52(10): 104-111.
|
|
|
[35] |
朱治, 李龙, 李超男, 毛新国, 郝晨阳, 朱婷, 王景一, 常建忠, 景蕊莲. 小麦转录因子TaMYB5-3B与株高和千粒重相关. 作物学报, 2023, 49(4): 906-916.
doi: 10.3724/SP.J.1006.2023.21029 |
|
|
[36] |
韩静然, 王长有, 赵宁娟, 吉万全, 赵毓, 刘雪利. 冬小麦新种质N0238D矮秆性状的遗传分析. 西北农业学报, 2010, 19(6): 60-63.
|
|
|
[37] |
于月华, 耿洪伟, 王莉萍, 高文伟, 陈全家, 张巨松, 曲延英. 小麦TaSIM基因结构及表达分析. 西北农业学报, 2012, 21(11): 49-52.
|
|
|
[38] |
doi: 10.1146/annurev-arplant-080720-093057 pmid: 34143645 |
[39] |
doi: 10.1186/s12870-023-04535-x pmid: 37904113 |
[40] |
|
[1] | 蒲丽霞, 张佳芮, 叶建萍, 黄秀兰, 樊高琼, 杨洪坤. 二氢赤霉素与秸秆覆盖对旱地小麦分蘖成穗与产量的影响[J]. 中国农业科学, 2025, 58(9): 1735-1748. |
[2] | 李云丽, 刁邓超, 刘雅睿, 孙玉晨, 孟祥宇, 邬陈芳, 汪妤, 吴建辉, 李春莲, 曾庆东, 韩德俊, 郑炜君. 小麦苗期耐热性全基因组关联分析[J]. 中国农业科学, 2025, 58(9): 1663-1683. |
[3] | 吴郁, 曲翔汝, 杨丹, 伍芩, 陈国跃, 江千涛, 魏育明, 许强. 广泛非靶向代谢组学解析小麦抗条锈病反应中叶绿体代谢产物[J]. 中国农业科学, 2025, 58(7): 1333-1343. |
[4] | 尹波, 于爱忠, 王鹏飞, 杨学慧, 王玉珑, 尚永盼, 张冬玲, 刘亚龙, 李悦, 王凤. 绿肥还田结合氮肥减施对干旱灌区麦田土壤水热特征及小麦产量的影响[J]. 中国农业科学, 2025, 58(7): 1366-1380. |
[5] | 潘丽媛, 王永军, 李海军, 侯富, 李菁, 李丽丽, 孙苏阳. 基于转录组和WGCNA筛选小麦籽粒蛋白质累积相关调控基因[J]. 中国农业科学, 2025, 58(6): 1065-1082. |
[6] | 唐宇, 雷毕欣, 王传伟, 严轩韬, 王浩, 郑杰, 张文静, 马尚宇, 黄正来, 樊永惠. 彩色小麦花青素积累对花后高温胁迫的响应机理[J]. 中国农业科学, 2025, 58(6): 1083-1101. |
[7] | 张洪程, 邢志鹏, 张瑞宏, 单翔, 奚小波, 程爽, 翁文安, 胡群, 崔培媛, 魏海燕. 小麦整合化无人化丰产栽培的特征与技术途径[J]. 中国农业科学, 2025, 58(5): 864-876. |
[8] | 张玲, 曹磊, 蔡诚, 颜欣怡, 项博才, 艾佳, 詹新杨, 宋有洪, 朱玉磊. 自然老化条件下冬小麦种子活力及生理指标变化[J]. 中国农业科学, 2025, 58(5): 877-889. |
[9] | 佘文婷, 孙蕊卿, 党海燕, 李文虎, 张丰, 田怡, 徐隽峰, 丁玉兰, 王朝辉. 我国主要麦区农户小麦籽粒硫的含量、分布及影响因素[J]. 中国农业科学, 2025, 58(5): 956-974. |
[10] | 刁邓超, 李云丽, 孟祥宇, 季颂涵, 孙玉晨, 马学红, 李杰, 冯永佳, 李春莲, 吴建辉, 曾庆东, 韩德俊, $\boxed{\hbox{王长发}}$, 郑炜君. 小麦TaGRAS34-5A的克隆及耐热功能研究[J]. 中国农业科学, 2025, 58(4): 617-634. |
[11] | 史帆, 李文广, 易树生, 杨娜, 陈玉萌, 郑伟, 张雪辰, 李紫燕, 翟丙年. 有机无机肥配施下旱地麦田土壤有机碳组分含量的变化特征[J]. 中国农业科学, 2025, 58(4): 719-732. |
[12] | 牟树佳, 董莉霞, 李广, 燕振刚, 逯玉兰. 基于鲸鱼算法的旱地春小麦农田土壤N2O排放模型参数优化[J]. 中国农业科学, 2025, 58(3): 537-547. |
[13] | 罗一诺, 李艳霏, 李文虎, 张思琦, 牟文燕, 黄宁, 孙蕊卿, 丁玉兰, 佘文婷, 宋文斌, 李小涵, 石美, 王朝辉. 我国新育成小麦品种(系)籽粒不同部位铁含量及影响因素[J]. 中国农业科学, 2025, 58(3): 416-430. |
[14] | 仇海龙, 李盼, 张殿凯, 樊志龙, 胡发龙, 陈桂平, 范虹, 何蔚, 殷文, 赵连豪. 西北绿洲灌区麦后复种绿肥对减量施氮春小麦生长及产量的补偿效应[J]. 中国农业科学, 2025, 58(3): 443-459. |
[15] | 张斯佳, 杨杰, 赵帅, 李莉威, 王贵彦. 华北平原多样化作物与小麦-玉米轮作对土壤质量的影响[J]. 中国农业科学, 2025, 58(2): 238-251. |
|