郑玉冲 ,张琳琦 ,刘彬彬**

(1.中国科学院遗传与发育生物学研究所农业资源研究中心/河北省土壤生态学重点实验室/中国科学院农业水资源重点实验室 石家庄 050022;2.中国科学院大学 北京 100049)

微生物与植物间存在密切且复杂的相互作用,对植物生长具有重要意义[1-2]。植物根区微生物种类繁多,它们与植物共同进化,其多样性和功能是影响植物生长和健康的重要因素[3-4]。植物根区微生物的组成和功能受多种因素影响,如土壤环境、农田管理措施、植物基因型以及生长阶段等[5-6]。

植物基因型是影响根区微生物组成的重要因素。不同基因型植物可以通过根系分泌物选择特定微生物从而影响自身生长发育,如玉米(Zea mays)分泌黄酮类化合物驱动根际富集了草酸杆菌(Oxalobacteraceae),进而促进玉米氮素吸收[7]。Kwak 等[8]选择了对青枯菌(Ralstonia solanacearum)具有抗性的番茄(Solanum lycopersicum)品种以及易受青枯菌侵染的番茄品种,通过16S rRNA 基因测序发现两个番茄品种的根际土壤微生物存在差异,宏基因组学数据表明抗性品种根际土壤中存在更加丰富的黄杆菌(Flavobacteriaceae)基因组。有研究发现抗尖孢镰刀(Fusarium oxysporum)品种的菜豆(Phaseolus vulgaris)根际土壤细菌中假单胞菌科(Pseudomonadaceae)、芽孢杆菌科(Bacillaceae)和噬纤维菌科(Cytophagaceae)等微生物的丰度高于易感品种,网络分析表明,抗尖孢镰刀菌菜豆品种的根际土壤细菌网络比易感品种更复杂、联系更紧密[9]。不同氮素耐受性的苹果(Malus domestica)砧木根际土壤细菌群落结构差异显著,其中假黄单胞菌属(Pseudoxanthomonas)、芽孢杆菌属(Bacillus)和不动杆菌属(Acinetobacter)的相对丰度与苹果砧木氮含量呈正相关,并且在低氮耐受性的苹果砧木根际土壤中相对丰度较高[10]。Mahoney 等[11]对9 个冬小麦(Triticum aestivum)品种的根际土壤微生物组进行研究,结果表明有24个操作分类单元(Operational Taxonomic Unit,OTU)在小麦品种间存在显著差异,并且与品种相关的OTUs可能对宿主植物具有有益作用。Kavamura 等[12]研究发现高杆小麦品种根际土壤细菌群落的放线菌门、拟杆菌门(Bacteroidota)和变形菌门的相对丰度较高,而半矮秆小麦品种的疣微菌门(Verrucomicrobia)、浮霉菌门(Planctomycetes)和酸杆菌门(Acidobacteria)的相对丰度较高。因此,植物基因型对根区微生物组成有显著影响,深入研究作物与根区微生物之间复杂的相互作用,为促进作物生产提供了新途径。

氮素是作物生长发育过程中所需的关键营养元素[13],对作物的生命活动以及产量具有重要意义。化肥的使用对促进粮食生产至关重要,但在农业生产过程中为了追求作物产量,过量的氮输入也给生态环境带来了负面影响,如土壤酸化以及温室气体排放等[14-15]。同时,土壤氮素水平通常会驱动土壤微生物多样性、结构和功能的变化[16-17]。研究表明,过量施氮显著降低了玉米根际土壤和非根际土壤中固氮细菌群落的香农指数,此外,还降低了伯克氏菌属(Burkholderia)的相对丰度,提高了鞘脂菌属(Sphingobium)的相对丰度[18]。缺氮胁迫使得小麦根际土壤微生物群落的Alpha 多样性增加,糖酵解、过氧化物代谢和磷酸肌醇代谢等过程增强[19]。然而,关于不同氮素水平下作物品种影响根区微生物群落的相关研究较少。因此,研究不同作物品种微生物组在不同氮素水平下的差异对于实现农业的绿色发展至关重要。

小麦是我国重要的粮食作物,其生产对于保障国家粮食安全具有重要意义。‘科农9204’是高产氮高效小麦品种,童依平等[20]分析了包括‘科农9204’和‘京411’在内的12 个冬小麦品种的氮素利用效率,结果发现,‘科农9204’在高氮和低氮条件下均具有较高的以生物量为基础和以籽粒产量为基础的氮素利用效率。Shi 等[21]从基因组水平和转录组水平阐述了‘科农9204’氮素高效吸收和利用的分子机制,该研究发现‘科农9204’相比‘京411’在生殖发育阶段具有较高的氮素利用效率和发达的根系,并且营养吸收相关基因持续高表达,表明‘科农9204’氮高效的机理是受多方面调控的结果。此外,小麦品种氮素利用效率受施氮水平影响较大[22],因此应在不同氮素水平对小麦氮素吸收利用能力进行综合评估。

宿主品种和土壤氮素水平会影响根区微生物群落的组成,但小麦在不同生长阶段下根区微生物组的品种特异性差异以及对土壤氮素水平的响应机制仍不明确。本研究选择‘科农9204’‘科农2011’‘京411’和‘百农207’4 个小麦品种,分别在2 个氮素水平进行种植,于分蘖期、拔节期和灌浆期采集根际土和根样品。采用16S rRNA 基因高通量测序技术,分析不同小麦品种间根区微生物群落组成和多样性,同时分析小麦生理指标与根际土壤微生物的关系,以期获得不同品种小麦根区微生物的差异,为筛选对小麦生长有益的微生物提供依据。

1 材料和方法

1.1 试验设计

本研究田间试验在中国科学院栾城农业生态系统试验站(37°53′N,114°41′E,海拔50 m)进行。田间试验设置2 个氮素水平,高氮水平尿素施用量为300 kg(N)·hm-2,低氮水平不施肥。分别在分蘖期、拔节期和灌浆期进行灌溉,每次灌溉水量约90 mm。氮肥条件为主区,品种为副区,播种方式为机械条播,每小区6 行,行长8 m,高氮和低氮间过渡区为2 m,两个施氮区各3 次重复,共计24 个小区。供试小麦品种分别为: ‘科农9204’ (KN9204,分蘖力强、株型紧凑、氮吸收利用效率较高)、‘科农2011’ (KN2011,抗倒伏、耐旱)、‘京411’ (J411,抗寒性强、适应性广)和‘百农207’ (B207,根系活力强、抗寒性好)。播种时间为2021 年10 月中旬,取样时间分别为2021 年11 月底(分蘖期)、2022 年4 月初(拔节期)和2022 年5 月中旬(灌浆期)。

1.2 地上部指标测定

在每个取样期,每个品种选择3 株小麦进行地上部指标测定。剪下小麦地上部分装入信封,置于烘箱在105 ℃下杀青30 min,然后于65 ℃条件下将小麦植株样品烘干至恒重,称重记录,计算地上部生物量(ADW)。将小麦干物质样品粉碎,采用硫酸-双氧水消解-凯氏定氮法测定全氮含量(ANC)。地上部氮素吸收量(ANA)为地上部生物量和全氮含量的乘积[23]。

1.3 样品采集

在每个取样期,每个小区选择5 处进行样品采集,每处选择5 株小麦将根系挖出,混合为该小区的代表样本,3 次重复。紧附在小麦根部的土壤作为根际土,首先抖落小麦根部附着的松散土壤,将小麦根部放入装有20 mL 无菌PBS 缓冲液(8.1 mmol·L-1Na2HPO4,1.9 mmol·L-1NaH2PO4,145 mmol·L-1NaCl,pH 7.2～7.4)的50 mL 离心管中,再于摇床中180 r·min-1震荡30 min。用无菌镊子取出根,静置后的沉淀为根际土壤样品,收集于2 mL 离心管内。将小麦根部用灭菌水冲洗3 次,置于含PBS 缓冲液的离心管中,超声震荡5 次,每次30 s,两次间歇30 s,无菌水清洗后于2%次氯酸钠中浸泡5 min,无菌水清洗3 次,再用70%的乙醇浸泡1 min。最后用无菌水清洗3 次,使用无菌滤纸吸干根表面水分,将根样品在无菌研钵中用液氮研磨后,用作根内生菌样品DNA 的提取。

1.4 DNA 提取和高通量测序

所有土壤和根系样品的DNA 均使用FastDNA Spin Kit for Soil 试剂盒(MP Biomedicals,USA)提取,具体方法参照说明书。DNA 的浓度和纯度用Nano-Drop One (ThermoFisher Scientific,USA)超微量分光光度计测定。将提取的DNA 送至上海派森诺生物科技有限公司进行16S rRNA 基因高通量测序,PCR扩增引物为799F (AACMGGATTAGATACCCKG)和1193R (ACGTCATCCCCACCTTCC)[24],对16S rRNA基因的V5-V7 区进行扩增,采用测序平台Illumina NovaSeq-PE250 测序。测序数据通过QIIME2 (Quantitative Insights Into Microbial Ecology)进行分析[25],对原始数据进行双端序列拼接、去除引物、去噪并生成特征表,用 Silva138 数据库进行微生物物种注释[26],在各分类水平上进行分类学分析、细菌群落Alpha多样性分析、Beta 多样性分析、差异物种分析以及相关性分析。

1.5 统计分析

采用R 语言4.1.2 的Duncan test 方法分析不同小麦品种间地上部指标以及细菌群落Alpha 多样性差异的显著性。采用R 语言中的vegan 包和ggplot2包基于Bray-Curtis 距离完成主坐标分析(PCoA)。采用R 语言中的vegan包的Adonis 功能分析不同小麦品种对细菌群落结构的影响。采用线性判别分析(LEfSe)筛选不同品种小麦根际土壤微生物的差异物种,其中分蘖期低氮水平下的筛选标准为LDA≥4且P＜0.05,其他条件下的筛选标准均为LDA≥2 且P＜0.05。采用R 语言中stats 包的Spearman 相关性分析方法分析微生物和小麦地上部指标的相关性。

2 结果与分析

2.1 不同小麦品种地上部性状分析

图1 为不同品种小麦的地上部生物量、地上部全氮含量和地上部氮素吸收量。分蘖期高氮水平下,‘科农9204’地上部生物量为46 g·m-2,比‘科农2011’和‘京411’分别高44%和60% (P＜0.05)。灌浆期高氮水平下,‘科农9204’的地上部生物量为4954 g·m-2,比‘百农207’高108% (P＜0.05)。灌浆期低氮水平下,‘科农9204’的地上部生物量为1633 g·m-2,比‘科农2011’和‘京411’分别高53%和68% (P＜0.05)。

图1 不同小麦品种地上部性状对比分析Fig.1 Comparative analysis of aboveground part traits of different wheat cultivars

在分蘖期、拔节期和灌浆期,高氮水平下4 个小麦品种的地上部全氮含量无显著差异。在分蘖期,低氮水平下‘科农2011’的地上部全氮含量分别比‘科农9204’‘京411’和‘百农207’高18%、8%和8% (P＜0.05)。

分蘖期高氮水平下,‘科农9204’的地上部氮素吸收量为1.62 g·m-2,比‘科农2011’和‘京411’分别高32%和51% (P＜0.05)。拔节期高氮水平下,‘科农9204’的地上部氮素吸收量为14 g·m-2,比‘百农207’高40%(P＜0.05)。灌浆期高氮水平下,‘科农9204’的地上部氮素吸收量为77 g·m-2,比‘百农207’高119% (P＜0.05)。灌浆期低氮水平下,‘科农9204’的地上部氮素吸收量为21 g·m-2,比‘科农2011’‘京411’和‘百农207’分别高132%、90%和88% (P＜0.05)。

2.2 不同小麦品种根区细菌群落分析

2.2.1 多样性分析

在分蘖期,高氮水平下‘京411’根际土壤细菌群落的Alpha 多样性显著低于其他3 个品种(P＜0.05),低氮水平下显著低于‘科农9204’和‘科农2011’(P＜0.05)。在拔节期,低氮水平下‘科农9204’根际土壤细菌群落的Alpha 多样性显著高于‘科农2011’和‘百农207’ (P＜0.05)。在分蘖期,低氮水平下‘科农9204’根内生细菌群落的Alpha 多样性显著高于‘京411’和 ‘百农207’ (P＜0.05) (表1)。

表1 不同品种小麦根际土壤和根内生细菌群落的Alpha 多样性Table 1 Alpha diversity of rhizosphere soil and endosphere bacterial community in different wheat cultivars

为探究小麦品种对根区细菌群落Beta 多样性的影响,基于Bray-Curtis 距离在ASV (amplicon sequence variants)水平对不同小麦品种根区细菌群落进行主坐标分析和Adonis 分析。结果发现,拔节期高氮水平下,小麦品种对根际土壤群落组成无显著影响,在其他条件下,小麦品种间根际土壤细菌群落组成差异显著(P＜0.05) (图2)。拔节期低氮水平下,小麦品种对根内生细菌群落组成无显著影响,在其他条件下,小麦品种间根内生细菌群落组成差异显著(P＜0.05)(图3)。

图2 不同施氮水平下不同小麦品种根际土壤细菌群落的主坐标分析(PCoA)Fig.2 Principal coordinate analysis (PCoA) of rhizosphere soil bacterial communities of different wheat cultivars at different nitrogen levels

图3 不同施氮水平下不同小麦品种根内生细菌群落的主坐标分析(PCoA)Fig.3 Principal coordinate analysis (PCoA) of root endosphere bacterial communities of different wheat cultivars at different nitrogen levels

2.2.2 物种组成分析

不同小麦品种根际土壤细菌群落门水平的组成如图4a。变形菌门和放线菌门是小麦根际土壤细菌群落的优势菌门。分蘖期时,小麦根际土壤中的变形菌门相对丰度显著高于拔节期和灌浆期(P＜0.05)。在分蘖期,低氮水平下‘科农9204’和‘科农2011’根际土壤变形菌门的相对丰度显著低于‘京411’和‘百农207’ (P＜0.05)。放线菌门和变形菌门是小麦根内生细菌群落的优势菌门(图4b)。高氮水平下,拔节期时小麦根内拟杆菌门相对丰度显著高于分蘖期和灌浆期(P＜0.05),拔节期和灌浆期时小麦根内绿弯菌门(Chloroflexi)的相对丰度显著高于分蘖期(P＜0.05)。在分蘖期,高氮水平下‘科农2011’根内放线菌门的相对丰度显著低于‘科农9204’‘京411’和‘百农207’(P＜0.05),低氮水平下‘京411’和‘百农207’根内变形菌门的相对丰度显著低于‘科农9204’和‘科农2011’(P＜0.05)。

图4 不同施氮水平下不同小麦品种根际土壤细菌(a)和根内生细菌(b)群落门水平相对丰度Fig.4 Relative abundance of rhizosphere soil (a) and root endosphere (b) bacterial community at the phylum level of different wheat cultivars at different nitrogen levels

2.2.3 根际土壤微生物差异分析

通过LEfSe 分析筛选不同小麦品种根际土壤微生物群落中的差异物种。在分蘖期高氮水下,相比其他3 个品种,‘科农9204’的根际土壤细菌群落中显著富集的菌群包括寡养单胞菌属(Stenotrophomonas)、生丝微菌科(Hyphomicrobiaceae)、福格斯氏菌属(Vogesella)、色杆菌科(Chromobacteriaceae)、代尔夫特菌属(Delftia)和Ferrovibrio(图5a)。在分蘖期,低氮水平下(因差异物种较多,图中仅展示LDA≥4 的微生物),‘科农9204’ 的根际土壤细菌群落中显著富集的菌群为Alpha 变形菌纲(Alphaproteobacteria)(图5b)。

图5 不同施氮水平下分蘖期不同小麦品种根际土壤微生物差异分析Fig.5 Analysis on the difference in rhizosphere soil microorganisms among different wheat cultivars at tillering stage at different nitrogen levels

在拔节期高氮水平下,微杆菌科(Microbacteriaceae)对4 个品种根际土壤细菌群落组间差异影响最大,并且在‘科农9204’ 根际土壤细菌群落中显著富集。此外,‘科农9204’ 的根际土壤细菌群落还显著富集四折叠球菌属(Quadrisphaera)、动孢菌目(Kineosporiales)、动孢囊菌科(Kineosporiaceae)和Rubellimicrobium(图6a)。在拔节期低氮水平下,对4 个品种根际土壤细菌群落组间影响最大的差异物种为Alpha 变形菌纲,并且在‘科农9204’ 的根际土壤细菌群落中显著富集。此外,根瘤菌目、根瘤菌科(Rhizobiaceae)、黄色杆菌科(Xanthobacteraceae)和芽单胞菌属也在‘科农9204’中显著富集(图6b)。

图6 不同施氮水平下拔节期不同品种小麦根际土壤微生物差异分析Fig.6 Analysis on the difference in rhizosphere soil microorganisms among different wheat cultivars at jointing stage at different nitrogen levels

在灌浆期高氮水平下,‘科农9204’ 根际土壤细菌群落中的弗兰克氏菌目和全噬菌纲(Holophagae)显著高于其他3 个品种(图7a)。在灌浆期低氮水平下,‘科农9204’ 的根际土壤细菌群落中显著富集的菌群包括黄色杆菌科和链孢囊菌属(Streptosporangium) (图7b)。

图7 不同施氮水平下灌浆期不同品种小麦根际土壤微生物差异分析Fig.7 Analysis on the difference in rhizosphere soil microorganisms among different wheat cultivars at filling stage at different nitrogen levels

2.3 小麦地上部性状和根际土壤微生物相关性分析

通过对小麦根际土壤细菌群落中相对丰度前20的菌属和地上部相关指标进行相关性分析。结果表明,根际土壤细菌群落与地上部性状存在显著相关性。节杆菌属、德沃斯氏菌属(Devosia)、土壤红杆菌属(Solirubrobacter)、红色杆菌属、链霉菌属、原小单孢菌属(Promicromonospora)、壤霉菌属(Agromyces)、叶杆菌属(Phyllobacterium)、Gaiella、类诺卡氏菌属和鞘氨醇单胞菌属(Sphingomonas)的相对丰度与小麦地上部氮素吸收量和地上部生物量呈显著正相关(P＜0.05),而马赛菌属、砂单胞菌属、假单胞菌属和黄杆菌属与全氮含量呈显著正相关(P＜0.05)(图8)。

图8 小麦根际土壤微生物与地上部性状的相关性分析Fig.8 Correlation analysis between rhizosphere soil microorganisms and traits of wheat aboveground part

3 讨论与结论

植物微生物组主要包括根际、叶际和内生微生物,对于植物的养分吸收、生长发育以及抵御生物和非生物胁迫至关重要[3,8,27-29]。影响植物微生物组成的因素包括人为干扰因素(如耕作)[30]、土壤因素(如理化性质)[31]以及宿主基因型[32]等。Zhang 等[33]研究发现水稻(Oryza sativa)亚种间根系微生物群落的Alpha 多样性存在显著差异,其中氮利用效率较高的籼稻品种多样性要显著高于粳稻,说明籼稻的根系相比粳稻招募了更多种类的微生物,并且含有更多氮代谢相关的菌属。在本研究中,灌浆期低氮水平下,‘科农9204’的地上部氮素吸收量显著高于其他3 个品种。不同小麦品种影响了根区微生物群落多样性。低氮水平下,‘科农9204’的根际土壤细菌多样性在分蘖期显著高于‘京411’,在拔节期显著高于‘科农2011’,在分蘖期和拔节期均显著高于‘百农207’。此外,‘科农9204’的根内细菌群落多样性在分蘖期低氮水平下显著高于‘京411’和‘百农207’。较高的微生物多样性可能意味着生态过程的多样性,这对于土壤健康以及植物生产力至关重要[34]。研究发现,根区微生物群落的组成会受植物基因型影响[35],这与本研究结果一致。小麦品种影响了根际土壤细菌群落组成,这可能是因为不同品种小麦根系分泌物的组成和数量不同,而植物可以通过改变根系分泌物驱动根际过程以适应环境[36-37]。

通过对4 个小麦品种根际土壤微生物进行差异分析,发现在‘科农9204’的根际土壤中显著富集了根瘤菌目、弗兰克氏菌目和芽单胞菌属等微生物物种。其中根瘤菌目中包含固氮细菌,如根瘤菌属和慢生根瘤菌属[42]。弗兰克氏菌作为植物根际促生菌可与非豆科植物共生并实现生物固氮,研究表明,弗兰克氏菌提高了植株生物量、总氮含量和叶绿素含量,显著促进植物生长,并增强了植物对非生物胁迫的抵抗力[43-44]。芽单胞菌属是农业生态系统中的主要群体之一,与土壤微生物群落的稳定性相关[45]。本研究发现‘科农9204’的根际土壤中显著富集的微生物包含Alpha 变形菌纲,该纲中的微生物具有重要生物学作用,包括固氮、氨氧化和光合作用等[46]。根际土壤微生物伴随植物的整个生育周期,和植物性状联系十分密切,其多样性和功能在作物生长和养分吸收中发挥至关重要的作用[47]。相关性分析表明节杆菌属、红色杆菌属、链霉菌属和类诺卡氏菌属与氮素吸收量和地上部生物量呈显著正相关。其中节杆菌属是重要的植物根际促生菌[48],链霉菌属、红色杆菌属和类诺卡氏菌属丰度的增加可能影响根际土壤的氮代谢过程[49]。

本研究揭示了不同品种小麦根区微生物群落组成差异以及与小麦氮素吸收量相关的微生物物种类群。这些结果为利用微生物区系组成调控作物生产和养分利用提供了依据。植物基因型会影响植物微生物群落的组成,但宿主基因如何驱动特征微生物装配并发挥功能仍需要进一步探索,后续研究可将宿主遗传学和微生物组联系起来,如利用全基因组关联分析(GWAS)研究植物与微生物组的互作机制。此外,后续可以进一步通过宏基因组学测序对不同品种小麦根区微生物群落的功能进行探究。