Effect of intercropping Alpinia katsumadai on soil fungal community within a teak (Tectona grandis) plantation
-
摘要:目的
柚木Tectona grandis 是热带珍贵树种之一，生长周期长，与药用植物套种是一种额外获得短期收益、达到“以短养长”的方法。因土壤微生物群落在土壤健康和质量方面有着重要作用，本文旨在探明柚木套种药用植物对土壤微生物群落的结构和功能影响，为柚木人工林科学经营和绿色发展提供理论参考。
方法以柚木林下套种草豆蔻Alpinia katsumadai 为研究对象，未套种草豆蔻的柚木为对照，通过采集两者根际土壤，利用高通量测序技术分析其真菌群落特征，从多样性指数、群落结构和组成、共现性网络来分析林药套种下土壤真菌群落的变化及差异。
结果柚木套种草豆蔻处理相比未套种的土壤真菌群落多样性无显著变化，多样性指数(Chao1、Shannon、Observed_species 和PD_whole_tree)差异不显著；但群落结构有显著不同，非度量多维尺度(NMDS)分析和主成分(PCA)分析显示柚木套种草豆蔻与未套种处理群落组成分离明显，线性判别分析(LEfSe) 显示两者富集最多的类群不同；套种处理土壤真菌网络结构更加集中和复杂(网络节点数和正向边比例上升，模块化指数增加)，子囊菌门和担子菌门分别为套种和未套种土壤真菌群落的中心响应类群。
结论柚木林下套种草豆蔻对土壤真菌群落结构有显著影响，套种草豆蔻提升了柚木根际真菌群落的网络复杂性和集中程度，形成了更紧密的土壤微生物环境，研究结果可从微生物角度为柚木人工林套种及提升生产力提供理论指导。
Abstract:ObjectiveTeak (Tectona grandis) is one of the precious tree species in the tropics, it has a long growth cycle. Intercropping with medicinal plants is an effective way to achieve additional short-term benefits and to achieve the goal of “raising the short with the long”. Since soil microbial community plays an important role in soil health and quality, the purpose of this study is to investigate the effects of intercropping teak with medicinal plants on the structure and function of soil microbial community, and provide theoretical references for the scientific management and green development of teak plantations.
MethodThis study investigated the effects of intercropping Alpinia katsumadai with teak in a forest plantation. A non-intercropped teak forest was used as a control. Soil samples were collected from the rhizospheres of both treatments and analyzed for fungal community characteristics using the high-throughput sequencing technology. The changes and differences in soil fungal communities under intercropping were analyzed based on diversity indices, community structure and composition, and co-occurrence networks.
ResultThe diversity of soil fungal communities in teak intercropped with A. katsumadai did not change significantly compared to non-intercropped forests. The diversity indices (Chao1, Shannon, Observed_species, and PD_whole_tree) did not differ significantly. However, the community structure was significantly different. Non-metric multidimensional scaling (NMDS) and principal component analysis (PCA) showed that the community composition of teak intercropped with A. katsumadai was clearly separated from that of non-intercropped forests. Linear discriminant analysis effect size (LEfSe) analysis showed that the most enriched taxa in the two treatments were also different. Soil fungal network structure in the intercropped treatment was more concentrated and complex (the number of network nodes, the proportion of positive edges and the modularity index increased). Ascomycota and Basidiomycota were the central response taxa for soil fungal communities in intercropped and non-intercropped forests, respectively.
ConclusionIntercropping with A. katsumadai significantly influences the soil fungal community structure in teak plantations. Intercropping A. katsumadai enhances the network complexity and centrality of the rhizosphere fungal community in teak, forming a more compact soil microbial environment. The research results can provide a theoretical guidance for teak plantation interplanting and improving productivity from the perspective of microorganisms.
-
Keywords:
- Intercropping of forest medicine /
- Fungus /
- Community structure /
- Teak /
- Co-occurrence network
-
荔枝、龙眼是我国南方优势名贵水果，主要分布于广东、广西、福建、海南、四川及云南等省区，被誉为“南国四大果品”之首。2019年全国荔枝、龙眼栽培面积分别高达54.1万和27.6万hm2，在农业经济中占有重要地位[1]。荔枝蒂蛀虫Conopomorphasinensis Bradley (Lepidoptera: Gracillariidae)是荔枝、龙眼栽培中频发的暴发性鳞翅目害虫，也是荔枝、龙眼果树害虫防控中的难点。首先，荔枝蒂蛀虫的初孵幼虫咬破卵壳后，钻入果实、花穗、嫩梢或叶肉等组织内为害，整个幼虫期都潜伏于寄主植物组织内部，无法防治；其次，该虫在果园世代重叠，蛹期参差不齐，成虫期长且高峰期不明显，唯有反复大量施用化学药剂才能保证成虫防控效果[2-3]。近年来，由杀卵剂介导的害虫卵期防控成为了荔枝蒂蛀虫防控的有效突破口。
昆虫依靠发达灵敏的感觉系统识别外界环境中的各种信息，进而完成取食、交配、生殖、躲避天敌、产卵等行为反应。气味结合蛋白(Odorant binding proteins, OBPs)和化学感受蛋白(Chemosensory proteins, CSPs)在昆虫感知外界化学信息过程中发挥着重要作用，其大量存在于化学感受器的淋巴液中，是昆虫化学信息感知的2类关键载体蛋白[4-5]。其中，CSPs是一类可溶性的小分子球状蛋白，典型特征是含有4个保守的半胱氨酸，能两两配对形成2个二硫键，用以保持CSPs蛋白结构的稳定性。CSPs在昆虫体内的许多部位都有表达，例如触角、性腺、翅膀、足、精卵巢等[6-7]。CSPs以广泛的分布及高度保守性被认为在昆虫对化学刺激的行为反应及多种生理生化过程中扮演着重要的角色。
有关昆虫CSPs研究绝大多数集中于成虫期。国家荔枝龙眼产业体系—虫害防控岗位科学家研究团队前期在荔枝蒂蛀虫卵转录组数据中鉴定了5个CSPs基因，为深入解析其结合特征及功能，本研究对荔枝蒂蛀虫卵5个CsinCSPs基因进行克隆，对其序列特征进行生物信息学分析，通过分子对接模拟预测荔枝、龙眼果园常用的几种杀虫剂与5个CsinCSPs的结合特征，利用qRT-PCR技术测定了5个CsinCSPs基因在除虫脲处理后的表达变化，以期为阐明荔枝蒂蛀虫卵期CSPs功能研究提供线索。
1. 材料与方法
1.1 供试昆虫
荔枝蒂蛀虫成虫采集自广东省农业科学院白云基地荔枝园，在温度(25±1)℃，湿度为(75±1)%、光周期为14 h光∶10 h暗的养虫室中，以φ为10%的蜂蜜水作为营养补充饲养，每天定时收集并更换产有荔枝蒂蛀虫卵的白色卫生卷纸，于养虫室内存放备用，方法参照文献[8]。
1.2 总RNA的提取与反转录
收集孵化1、2和3 d荔枝蒂蛀虫卵各200粒，分别用液氮研磨后，利用Trizol试剂提取总RNA，通过20 g/L琼脂糖凝胶电泳检测总RNA的完整性。采用PrimeScript RT Reagent试剂盒(Takara, Tokyo, Japan)以提取的总RNA为模板反转录得到cDNA第1链，−20 ℃保存备用。
1.3 序列克隆与分析
以荔枝蒂蛀虫卵转录组(数据待发表)测序结果为基础，从成功注释的unigene中搜索CSPs基因序列，将其对应的蛋白序列与NCBI中的nr数据库进行BLASTP，验证是否属于化学感受蛋白。使用Primer Premier 3.0软件根据目的基因片段设计用于5个CsinCSPs基因全长序列扩增和表达分析的引物(表1)。分别以2对特异性嵌套引物进行5′和3′RACE扩增。5′RACE反应程序为：94 ℃ 3 min；94 ℃ 30 s，65 ℃ 30 s，72 ℃ 1 min，33个循环；72 ℃ 10 min；4 ℃保存。3′RACE反应程序为：94 ℃ 3 min；94 ℃ 30 s，58 ℃ 30 s，72 ℃ 1 min，33个循环；72 ℃ 10 min；4 ℃保存。最后用DNAstar软件将所得序列片段拼接为完整的序列。
表 1 CsinCSPs克隆及实时荧光定量PCR所需引物Table 1. Primers for CsinCSPs cloning and qRT-PCR基因
Gene巢式PCR引物(5′→ 3′) Nested gene-specific primers 实时荧光定量PCR (5′→ 3′)
qRT-PCR5′扩增引物
Primers for 5′RACE3′扩增引物
Primers for 3′RACECsinCSP1 CAATATGGTCGTACTTGGTAG AAAGAGTGGAAACAGCTTCT CCGAAGGAGCTGAACTTAGAA AGACTACTGCAATTCGTGA AACAGCTTCTTGACAAATGG TCTCTGGAGTCTGTCAACTACT CsinCSP2 AGCACACGTTTGAAGTTCT AGCTGGTGGTCAGGAAAC CCTCATTACAACCGCTACGATAA TTGTAGTAAGCGAGAAGAAC AAGAGCGAGCCCAGATTCTC GAAGTTCTTCCCGTCCTTAGTG CsinCSP3 GCAACATCTCCTCGAAGT TCAAGCACAACACTCGATTG CACCACACACTGAGCCATTC CACATCGTATATGGTCAGGT AGTACGACCCTGAGGGTAAC GACGGCATCGAGGAACAATC CsinCSP4 GCACTTGATCTGTCTTTGG GTGTCCGCAATGCACACCAC TGGTGGCGTGCTGTCAT GTTGTTGTTGTTGTTGCG CGTACTCTATCCTACGTGC TTGTCATTGAGAGCGTCGTC CsinCSP5 CGCTCCATGTAGTACAAGT ACTTGTACTACATGGAGCG GGAACGGACTCTGGCTGAT CACAGCGATACAGGACAG AGTGGGCCAAGATTGTACG GACACGCTCCATGTAGTACAAG β-Actin AGATCTGGCACCACACCTACGAT
ACCGGTGGTACGAC蛋白相对分子质量及等电点预测使用ExPASy平台中的Compute pI/Mw ( http://expasy.org/tools/pi_tool.html)，信号肽预测使用ExPASy平台中的SignalP 4.1 ( http://www.Cbs.dtu.dk/services/SignalP/)，亚细胞定位预测使用LocTree2软件( https://rostlab.org/services/loctree3/)。将去除信号肽的氨基酸序列导入在线软件PSIPRED( http://bioinf.cs.ucl.ac.uk/psipred/)构建蛋白二级结构。采用MEGA 4.1软件中的邻位相连法(Neighbor-joining，NJ)构建系统进化树，Bootstrap为1000次。
1.4 3D结构预测及分子对接模拟
使用Phyre2软件构建蛋白三维结构( http://www.sbg.bio.ic.ac.uk/phyre2)。选取荔枝、龙眼果园最常用的除虫脲、吡虫啉、辛硫磷、毒死蜱、三唑磷、氯虫苯甲酰胺、高效氯氟氰菊酯和高效氯氰菊酯作为配体分子，从NCBI的Pubchem Compound数据库得到8种药剂的三维结构。使用AutoDock半柔性配体整体分子自动对接程序，将模建蛋白及配体分子进行分子对接，根据结合能大小讨论其相互作用。
1.5 药剂处理
前期以孵化后1 d的荔枝蒂蛀虫卵开展除虫脲毒力测定，其毒力方程为y=−0.49+1.32x，LC50 为2.37 mg/L[9]。处理组：将荔枝蒂蛀虫卵纸置于2.37 mg/L除虫脲(以丙酮为介质稀释)中浸泡处理10 s；对照组：将荔枝蒂蛀虫卵纸置于丙酮中浸泡处理10 s。然后分别于处理后6、24和48 h随机收集处理组和对照组虫卵各200粒，液氮处理后−80 ℃保存备用。样品总RNA的提取与反转录方法同“1.2”。
1.6 实时荧光定量PCR
根据5个CsinCSPs基因序列，用Primer Premier 3.0软件设计qRT-PCR 特异性引物。利用GoTaq qPCR Master Mix试剂盒(Promega，麦迪逊，美国)进行qRT-PCR。反应体系包括：10 μL GoTaq qPCR Master Mix，1 μL cDNA模板，0.5 μL正向引物，0.5 μL反向引物，8 μL DEPC水。反应程序：95 ℃预变性2 min；95 ℃变性15 s；60 ℃退火60 s；40个循环。以β-actin (GenBank accession No. KF598848)基因为内参基因计算相对表达量，数据分析使用比较Ct值方法分析(Schmittgen and Livak, 2008)。每次qRT-PCR需要3次技术重复，并且每个样本进行3次生物学重复。
2. 结果与分析
2.1 荔枝蒂蛀虫卵5个CsinCSPs基因的克隆与序列分析
以荔枝蒂蛀虫卵cDNA为模板进行PCR扩增，得到CsinCSP1~5的基因开放阅读框(Open reading frame, ORF)长度分别为438、378、360、336和318 bp，分别编码145、125、119、111和105个氨基酸(表2)。CsinCSP1~5蛋白相对分子质量分别为16 400、14 100、13 200、12 600和11 900，其中，CsinCSP3等电点为5.04，偏酸性，其余4个CsinCSPs等电点均大于9，偏碱性。CsinCSP1~5的N端分别包含由起始位置开始的16~22个信号肽序列，无跨膜结构。蛋白二级结构预测表明，5个CsinCSPs的保守半胱氨酸位点排布模式均为C1-X6-C2-X18-C3-X2-C4。α−螺旋是所有CsinCSPs的主要二级结构，其中CsinCSP1、CsinCSP2和CsinCSP3均含有6个α−螺旋，而CsinCSP4和CsinCSP5仅含有5个α−螺旋(图1)。亚细胞定位分析表明CsinCSP1~5均为非膜蛋白。
表 2 荔枝蒂蛀虫卵的CsinCSPs蛋白信息Table 2. Information of CsinCSPs in Conopomorpha sinensis eggs基因
Gene开放阅读框长度/ bp
Length of open
reading frame氨基酸数量
No. of amino acid相对分子质量
Relative molecular weight等电点
Isoionic
point亲水系数
Hydrophilic
coefficient信号肽数量
No. of signal
peptiedeCsinCSP1 438 145 16 400 9.06 −0.417 22 CsinCSP2 378 125 14 100 9.28 −0.240 19 CsinCSP3 360 119 13 200 5.04 −0.024 16 CsinCSP4 336 111 12 600 9.41 −0.523 16 CsinCSP5 318 105 11 900 9.58 −0.011 19 ![]()
图 1 荔枝蒂蛀虫卵5个CsinCSPs的二级结构分析Figure 1. Secondary structure analysis of five CsinCSPs from Conopomorpha sinensis eggs2.2 基于荔枝蒂蛀虫卵5个CsinCSPs氨基酸序列进化树的构建与同源性分析
通过与其他昆虫CSPs氨基酸序列进行比对，结果(图2)显示，与CsinCSP1~5相似性最高的同类蛋白分别为蚜科棉蚜Aphis gossypii CSP1、螟蛾科稻纵卷叶螟Cnaphalocrocis medinalis CSP1、螟蛾科二化螟Chilo suppressalis CSP2、凤蝶科柑橘凤蝶Papilio xuthus CSP2和枯叶蛾科云南松毛虫Dendrolimus houi CSP2。使用MEGA 4.1软件邻接法进行1000次重复计算后构建进化树(图2)。由图2可以看出，所有选中的昆虫CSPs形成2个明显的分支，CsinCSP1~3属于6个α−螺旋的类群，而CsinCSP4和CsinCSP5属于5个α−螺旋的类群。借助WebLogo软件对所有选中昆虫的保守区域进行motif可视化分析，结果表明4个半胱氨酸位点非常保守，其次第2个和第3个半胱氨酸位点中间的甘氨酸、脯氨酸和亮氨酸也有较高保守性(图3)。
![]()
图 2 CsinCSP1~5与其他昆虫CSPs的系统进化树分析Figure 2. Phylogenetic analysis of CsinCSP1−5 and CSPs in other insects![]()
图 3 CsinCSP1~5与其他昆虫CSPs的WebLogo比对分析“*”表示4个高度保守的半胱氨酸残基Figure 3. The WebLogo alignment of CsinCSP1-5 and CSPs in other insects“*” indicates four highly conservative cysteine residues2.3 荔枝蒂蛀虫卵5个CsinCSPs同源建模和结合功能预测
荔枝蒂蛀虫卵CsinCSP1~5蛋白的三维结构预测分析中，软件数据库比对建模筛选出的CsinCSP1~5结构模板为甘蓝夜蛾Mamestra brassicae CSP2的三维结构(PDB: d1n8va)。与其他昆虫的CSP一样，其编码蛋白主要由α−螺旋(α-helix)、β−折叠(β-strand)和卷曲环(loop)组成，这些结构形成一个与配体结合的空腔，空腔内部的氨基酸大部分为疏水性氨基酸，使得其内部可容纳疏水性的配体分子(图4)。
从NCBI的Pubchem Compound数据库得到荔枝、龙眼果园最常用的8种化学药剂的三维结构，然后利用AutoDock软件执行自动计算对接，以获得荔枝蒂蛀虫卵CsinCSPs与化学药剂分子的结合能。分子对接数据表明，5个CsinCSPs与所有配体对接相互作用能量值均为负值。其中，配体分子高效氯氟氰菊酯、高效氯氰菊酯与5个CsinCSPs结合能均为最小，结合亲和力最强；5个CsinCSPs与除虫脲、吡虫啉相互作用结合能均高于其他药剂的，表明CsinCSPs与除虫脲和吡虫啉的结合力最弱(表3)。
表 3 农药与荔枝蒂蛀虫卵CsinCSPs对接的结合能Table 3. The binding capacities of insecticides and CsinCSPs of Conopomorpha sinensis eggskJ·mol−1 配体分子
LigandCSP1 CSP2 CSP3 CSP4 CSP5 除虫脲 Diflubenzuron −8.07 −9.00 −8.45 −7.66 −7.83 吡虫啉 Imidacloprid −7.62 −8.57 −7.92 −8.22 −8.03 辛硫磷 Phoxim −8.41 −9.70 −9.43 −8.66 −8.45 毒死蜱 Chlorpyrifos −8.44 −9.32 −9.14 −8.36 −8.46 三唑磷 Triazophos −9.00 −9.84 −9.77 −8.57 −8.62 氯虫苯甲酰胺 Chlorantraniliprole −10.44 −10.86 −11.16 −9.95 −9.84 高效氯氟氰菊酯 λ-cyhalothrin −11.65 −12.59 −11.03 −10.99 −10.42 高效氯氰菊酯 β-cypermethrin −11.01 −11.57 −11.39 −10.18 −11.01 2.4 除虫脲处理对荔枝蒂蛀虫卵5个CsinCSPs基因表达量的影响
以亚致死剂量除虫脲处理荔枝蒂蛀虫1日龄卵，利用qRT-PCR检测虫卵5个CsinCSPs基因在处理后3个时间点的表达量，结果(图5)表明：与空白对照组相比，除虫脲处理后6、24、48 h，CsinCSP1、CsinCSP2和CsinCSP4的相对表达量变化不大。但是CsinCSP3在处理后24 h表达显著上升49.4倍，CsinCSP5在处理后6 h表达上升57.8倍，随后二者在处理后48 h表达则分别下降36.2和23.4倍，说明除虫脲能影响荔枝蒂蛀虫卵CsinCSP3和CsinCSP5的表达，但对其他3个CsinCSPs的表达无明显影响(图5)。
![]()
图 5 亚致死剂量的除虫脲对荔枝蒂蛀虫卵5个CsinCSPs基因表达量的影响Figure 5. Effects of diflubenzuron on expression of five CsinCSPs genes in Conopomorpha sinensis eggs at sublethal dose3. 讨论与结论
化学感受蛋白在昆虫感知外界信号分子的生化反应过程中具有重要生理学意义。本研究利用转录组测序和RACE-PCR的方法，从荔枝蒂蛀虫卵中克隆出5个CsinCSPs基因的全长序列，为首次从鳞翅目害虫卵中克隆同类蛋白序列。其中ORF最短和最长的CsinCSP分别为CsinCSP5和CsinCSP1，分别为318和438 bp，分别编码105和145个氨基酸，预测蛋白相对分子质量分别为11 900和16 400。普遍认为昆虫CSPs是一类相对分子质量约为10 000~12 000的酸性可溶性蛋白，编码100~115个氨基酸分子[10-11]。昆虫CSPs在进化上十分保守，在N端有高度保守区域，其序列多样性表现在蛋白序列的C端，有些CSPs则会在C端翻译提前终止而序列较短、相对分子质量更小，有些CSPs的C端有延伸而序列较长、相对分子质量更大，例如棉铃虫Helicoverpa armigera HarmiCSP12、HarmiCSP 14和HarmiCSP 25[12]；食蚜蝇Episyrphus balteatus EbalCSP3、EbalCSP 4和EbalCSP6[13]；柑橘木虱Diaphorina citri DcitCSP12[14]和褐飞虱Nilaparvata lugens NlugCSP3[15]等。本研究中CsinCSP1、CsinCSP2和CsinCSP3则属于C端有延伸的昆虫CSP类群。此外，在5个CsinCSPs的N端具有长度为16~22个氨基酸的信号肽序列，说明CsinCSPs具备分泌蛋白特征。亚细胞定位软件LocTree2预测结果显示，本研究中5个CsinCSPs主要为非膜蛋白，与上述推断一致。
5个荔枝蒂蛀虫卵CsinCSPs有保守的半胱氨酸位点，其位点排布模式均为C1-X6-C2-X18-C3-X2-C4，在C2和C3位点之间有多个高度保守的其他氨基酸位点，说明5个CsinCSPs基序的高度保守性。5个CsinCSPs在二级结构上相邻的半胱氨酸顺次连接形成2个二硫键，与报道的昆虫CSPs基本特征一致。在二级结构上CsinCSP1~3包括6个α−螺旋，而CsinCSP4和CsinCSP5有5个α−螺旋，在系统进化分析中证实了5个CsinCSPs属于不同CSP类群，利用在线软件预测CsinCSPs蛋白的三维结构也验证了5个CsinCSPs的α−螺旋数目不同这点。对功能蛋白而言，保守的基序和稳定的结构对其功能是非常重要的，因此由5个CsinCSPs二维及三维结构的区别可预测其生理功能的区别。与昆虫OBPs相比，CSPs在进化上更为保守，不同物种间序列相似度可高达40%~50%。昆虫CSPs序列进化的保守性也可能与其特征有关，比如它们对常见环境挥发性物质有较低结合力，而与通信有关的信号化学物质则有强结合力[6]。昆虫CSPs在三维结构上有5或6个α−螺旋，因其序列长度范围较广，三维结构多变，而被认为可结合更多种类的配体分子[7, 11]。在解析荔枝蒂蛀虫5个CsinCSPs三维结构基础上，本研究以荔枝、龙眼果园8种最常用药剂为配体，利用基于“锁−钥学说”的分子对接技术简单快速地模拟出配体与CsinCSPs蛋白之间的相互作用[16]。分子对接数据表明5个CsinCSPs与2个菊酯类杀虫剂分子的结合能力最强，其次为辛硫磷、毒死蜱、三唑磷和氯虫苯甲酰胺，与除虫脲和吡虫啉结合能力最弱。在实际生产应用中，除虫脲和吡虫啉可作为杀卵剂应用于防治荔枝蒂蛀虫[9, 17]。因此我们推测，荔枝蒂蛀虫卵CsinCSPs可能参与其他药剂的运输和降解，而无法包裹除虫脲和吡虫啉分子，最终导致除虫脲和吡虫啉对荔枝蒂蛀虫卵起到灭杀效果。但分子对接结果仅为模拟预测数据，后续需通过荧光竞争结合试验获得解离常数进一步确定其结合效力。
除了在昆虫嗅觉识别过程中发挥重要作用，CSPs还可能参与调节昆虫的生长发育、取食、感受机械刺激、免疫、信号传导等生理功能[11]。近年来研究发现，昆虫CSP与抗药性相关。对Bt抗性种群和Bt敏感种群的美洲条螟Diatraea saccharrlis幼虫进行基因表达水平检测发现，2个种群中的某些CSPs表达不同，表明这些CSPs可能参与了条螟的抗性反应[18]。亚致死剂量阿维菌素处理家蚕Bombyx mori成虫后，会导致BmorCSP3、BmorCSP 6和BmorCSP 7等多个CSP蛋白在多个组织中表达上调[19]。用闹羊花毒素III或苄氯菊酯处理小菜蛾Plutella xylostella成虫[20-21]、毒死蜱处理斜纹夜蛾Spodoptera litura幼虫[22]、噻虫嗪处理烟粉虱Bemisia tabaci成虫[23]、吡虫啉处理蜜蜂Apis cerana和Apis mellifera成虫[24-25]均得到同样的结果。此外，通过RNAi方法，以禾谷缢蚜Rhopalosiphum padi成虫为材料证实RpadCSP10与其对吡虫啉的抗性直接相关[26]；以赤拟谷盗Tribolium castaneum若虫为材料[27]，证明TcasCSP10与靶标害虫对敌敌畏和克百威的抗性相关。有学者推测：杀虫剂处理后昆虫CSPs表达变化，可能原因是成虫腹部表达的CSP蛋白可包裹有毒的杀虫剂分子，让其随粪便排出体外，对杀虫剂效力起到缓冲作用[9, 13]。但是同类研究在多种昆虫的幼虫期或成虫期均有所开展，在害虫卵期却鲜见文献报道。本研究利用害虫卵开展相关研究，结果显示亚致死剂量除虫脲处理荔枝蒂蛀虫卵后，CsinCSP3和CsinCSP5在处理后6和24 h的表达均上升40倍以上，随后二者在处理后48 h表达则下降20余倍，而另外3个CsinCSPs则变化不大，说明CsinCSP3和CsinCSP5可能参与荔枝蒂蛀虫卵对除虫脲的解毒过程。但是在分子对接分析中，CsinCSP2和CsinCSP3与除虫脲结合能力更强，CsinCSP5却与除虫脲结合能力较弱，说明分子对接结果不能完全预测CsinCSPs在昆虫化学药剂胁迫下解毒过程中的生理功能。
本研究以害虫卵为研究对象，对荔枝蒂蛀虫5个CsinCSPs基因进行了克隆及序列分析，并且在分子建模及分子对接的基础上，对除虫脲作用后虫卵中CsinCSPs基因表达情况进行了解析，结果证明在CsinCSPs可能参与荔枝蒂蛀虫卵在杀虫剂肋迫下的解毒过程，后续可通过RNAi或CRISPR基因编辑等技术在活体水平确定该类蛋白在害虫抗药性中的作用机制以及在昆虫卵期的相关功能。
-
![]()
图 1 柚木林下套种草豆蔻示意图
Figure 1. Diagram of intercropping Alpinia katsumadai within a Tectona grandis forest palntation
![]()
图 2 林药套种(T)与未套种(CK)根际土壤真菌多样性指数
Figure 2. Rhizosphere soil fungal diversity index of forest medicine intercropping (T) and non-intercropping (CK)
![]()
图 3 林药套种(T)与未套种(CK)根际土壤物种在门水平的组成
Figure 3. Rhizosphere soil species composition of forest medicine intercropping (T) and non-intercropping (CK) at phylum level
![]()
图 4 林药套种(T)与未套种(CK)担子菌门和子囊菌门相对丰度归一化指数
Figure 4. Normalized index of relative abundance of Basidiomycota and Ascomycota of forest medicine intercropping (T) and non-intercropping (CK)
![]()
图 5 林药套种(T)与未套种(CK)根际土壤非度量多维尺度法(NMDS) 分析和主成分分析(PCA)
Figure 5. Nonmetric multidimensional scaling (NMDS) analysis and principal component analysis (PCA) of rhizosphere soil of forest medicine intercropping (T) and non-intercropping (CK)
![]()
图 6 林药套种(T)与未套种(CK)根际土壤线性判别(LEfSe)分析
Figure 6. Linear discriminant analysis effect size (LEfSe) analysis of rhizosphere soil of forest medicine intercropping (T) and non-intercropping (CK)
![]()
图 7 林药套种(T)和未套种(CK)根际土壤真菌在门水平的共现性网络分析
Figure 7. Network analysis of rhizosphere soil fungi co-occurrence for forest medicine intercropping (T) and non-intercropping (CK) at phylum level
表 1 研究样地基本概况
Table 1 Basic information of the research plot
坡度/(°)
Slope海拔/m
Altitude树高/m
Height胸径或地径/cm
DBH or ground
diameter柚木
Tectona grandis0 61 12.6 33.5 草豆蔻
Alpinia katsumadai0 61 2.27 2.03 
下载: 导出CSV
表 2 林药套种(T)和未套种(CK)根际土壤真菌共现网络拓扑结构特征
Table 2 Topological characteristics of fungal co-occurrence networks in rhizosphere soil of forest medicine intercropping (T) and non-intercropping (CK)
处理
Treatment边
Edge节点
Note正向边比例/%
Positive edge proportion负向边比例/%
Negative edge proportion平均度
Average degree图密度
Graph density模块化
ModularityT 367 192 69.88 30.12 3.904 0.021 1.336 CK 332 188 48.77 51.23 3.458 0.018 0.830
下载: 导出CSV
-
[1] ZHOU X G, ZHANG J Y, RAHMAN M, et al. Interspecific plant interaction via root exudates structures the disease suppressiveness of rhizosphere microbiomes[J]. Molecular Plant, 2023, 16(5): 849-864. doi: 10.1016/j.molp.2023.03.009
[2] LIAN T X, MU Y H, JIN J, et al. Impact of intercropping on the coupling between soil microbial community structure, activity, and nutrient-use efficiencies[J]. PeerJ, 2019, 7: 6412. doi: 10.7717/peerj.6412
[3] ZHANG R Z, MENG L B, LI Y, et al. Yield and nutrient uptake dissected through complementarity and selection effects in the maize/soybean intercropping[J]. Food and Energy Security, 2021, 10(2): 282.
[4] QIANG W, HE L L, ZHANG Y, et al. Aboveground vegetation and soil physicochemical properties jointly drive the shift of soil microbial community during subalpine secondary succession in southwest China[J]. CATENA, 2021, 202: 105251. doi: 10.1016/j.catena.2021.105251
[5] WARD S E, BARDGETT R D, MCNAMARA N P, et al. Plant functional group identity influences short-term peatland ecosystem carbon flux: Evidence from a plant removal experiment[J]. Functional Ecology, 2009, 23(2): 454-462. doi: 10.1111/j.1365-2435.2008.01521.x
[6] DE DEYN G B, CORNELISSEN J, BARDGETT R D. Plant functional traits and soil carbon sequestration in contrasting biomes[J]. Ecology Letters, 2008, 11(5): 516-531. doi: 10.1111/j.1461-0248.2008.01164.x
[7] YU R P, LAMBERS H, CALLAWAY R M, et al. Belowground facilitation and trait matching: Two or three to tango?[J]. Trends in Plant Science, 2021, 26(12): 1227-1235. doi: 10.1016/j.tplants.2021.07.014
[8] CHEN Y, BONKOWSKI M, SHEN Y, et al. Root ethylene mediates rhizosphere microbial community reconstruction when chemically detecting cyanide produced by neighbouring plants[J]. Microbiome, 2020, 8(4): 40168.
[9] CHEN Y, SUN R B, SUN T T, et al. Organic amendments shift the phosphorus-correlated microbial co-occurrence pattern in the peanut rhizosphere network during long-term fertilization regimes[J]. Applied Soil Ecology, 2018, 124: 229-239. doi: 10.1016/j.apsoil.2017.11.023
[10] BANERJEE S, SCHLAEPPI K, VAN DER HEIJDEN M. Keystone taxa as drivers of microbiome structure and functioning[J]. Nature Reviews Microbiology, 2018, 16(9): 567-576. doi: 10.1038/s41579-018-0024-1
[11] DING K, ZHANG Y T, YRJÄLÄ K, et al. The introduction of Phoebe bournei into Cunninghamia lanceolata monoculture plantations increased microbial network complexity and shifted keystone taxa[J]. Forest Ecology and Management, 2022, 509: 120072. doi: 10.1016/j.foreco.2022.120072
[12] LIU Y, MA W, HE H, et al. Effects of sugarcane and soybean intercropping on the nitrogen-fixing bacterial community in the rhizosphere[J]. Frontiers in Microbiology, 2021, 12: 13349.
[13] LOREAU M, DE MAZANCOURT C. Biodiversity and ecosystem stability: A synthesis of underlying mechanisms[J]. Ecology Letters, 2013, 16(S1): 106-115. doi: 10.1111/ele.12073
[14] LI X G, DE BOER W, ZHANG Y A, et al. Suppression of soil-borne Fusariwn pathogens of peanut by intercropping with the medicinal herb Atractylodes lancea[J]. Soil Biology and Biochemistry, 2018, 116: 120-130. doi: 10.1016/j.soilbio.2017.09.029
[15] KUI L, CHEN B Z, CHEN J, et al. A comparative analysis on the structure and function of the Panax notoginseng rhizosphere microbiome[J]. Frontiers in Microbiology, 2021, 12: 673512. doi: 10.3389/fmicb.2021.673512
[16] 严毅, 何银忠, 王亚婷, 等. 云南海口林场中药林下种植模式初探[J]. 中国现代中药, 2016, 18(2): 173-177. [17] 黄桂华, 梁坤南, 付强, 等. 11年生柚木无性系遗传变异与优良无性系选择[J]. 东北林业大学学报, 2023, 51(8): 18-22. doi: 10.3969/j.issn.1000-5382.2023.08.003 [18] 赵威威, 周再知, 张青青, 等. 配比施肥对柚木早期生长及叶片养分含量的影响[J]. 西北农林科技大学学报(自然科学版), 2024, 52(1): 71-78. [19] LIU M M, YANG G, ZHOU W L, et al. Transcriptomic analysis reveals CBF-dependent and CBF-independent pathways under low-temperature stress in teak (Tectona grandis)[J]. Genes, 2023, 14(11): 2098.
[20] 付强, 黄桂华, 周强, 等. 柚木无性系早期选择年龄的研究[J]. 中南林业科技大学学报, 2022, 42(10): 30-38. [21] 陈天宇, 黄桂华, 王西洋, 等. 配比施肥对柚木无性系幼林生长的影响[J]. 中南林业科技大学学报, 2021, 41(4): 66-75. [22] 周强, 黄桂华, 张绍祥, 等. 肥料施肥对不同柚木无性系早期生长的影响[J]. 东北林业大学学报, 2022, 50(7): 6-10. doi: 10.3969/j.issn.1000-5382.2022.07.002 [23] 刘李花, 刘云, 尹加笔, 等. 柚木心材和边材的LC-MS代谢组学比较[J]. 林业工程学报, 2021, 6(3): 100-106. [24] YU Z, LIANG K N, WANG X B, et al. Alterations in arbuscular mycorrhizal community along a chronosequence of teak (Tectona grandis) plantations in tropical forests of China.[J]. Frontiers in Microbiology, 2021, 12: 737068. doi: 10.3389/fmicb.2021.737068
[25] 胡璇, 王丹, 于福来, 等. 草豆蔻的本草考证[J]. 中国实验方剂学杂志, 2020, 26(21): 210-219. [26] CHEN K Y, HU L, WANG C T, et al. Herbaceous plants influence bacterial communities, while shrubs influence fungal communities in subalpine coniferous forests[J]. Forest Ecology and Management, 2021, 500: 119656. doi: 10.1016/j.foreco.2021.119656
[27] LI T, LI Y Z, SHI Z, et al. Crop development has more influence on shaping rhizobacteria of wheat than tillage practice and crop rotation pattern in an arid agroecosystem[J]. Applied Soil Ecology, 2021, 165: 104016. doi: 10.1016/j.apsoil.2021.104016
[28] PROBER S M, LEFF J W, BATES S T, et al. Plant diversity predicts beta but not alpha diversity of soil microbes across grasslands worldwide[J]. Ecology Letters, 2015, 18(1): 85-95. doi: 10.1111/ele.12381
[29] BARBERÁN A, MCGUIRE K L, WOLF J A, et al. Relating belowground microbial composition to the taxonomic, phylogenetic, and functional trait distributions of trees in a tropical forest[J]. Ecology Letters, 2015, 18(12): 1397-1405. doi: 10.1111/ele.12536
[30] 陈爱玲, 林思祖, 陈世品, 等. 杉木连栽地上闽楠与杉木林分凋落物特征比较[J]. 福建林学院学报, 2006(4): 289-293. doi: 10.3969/j.issn.1001-389X.2006.04.001 [31] EDWARDS J, JOHNSON C, SANTOS-MEDELLÍN C, et al. Structure, variation, and assembly of the root-associated microbiomes of rice[J]. Proceedings of the National Academy of Sciences of the United States of America, 2015, 112(8): E911-E920.
[32] LIU F, HEWEZI T, LEBEIS S L, et al. Soil indigenous microbiome and plant genotypes cooperatively modify soybean rhizosphere microbiome assembly[J]. BMC Microbiology, 2019, 19(1): 201.
[33] WAGG C, SCHLAEPPI K, BANERJEE S, et al. Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning[J]. Nature Communications, 2019, 10: 4801.
[34] FAUST K, RAES J. Microbial interactions: From networks to models[J]. Nature Reviews Microbiology, 2012, 10(8): 538-550. doi: 10.1038/nrmicro2832
[35] MOUGI A, KONDOH M. Diversity of interaction types and ecological community stability[J]. Science, 2012, 337(6092): 349-351. doi: 10.1126/science.1220529
[36] YANG W, JING X Y, GUAN Y P, et al. Response of fungal communities and co-occurrence network patterns to compost amendment in black soil of northeast China[J]. Frontiers in Microbiology, 2019, 10: 1562. doi: 10.3389/fmicb.2019.01562
[37] FAN K, WEISENHORN P, GILBERT J A, et al. Wheat rhizosphere harbors a less complex and more stable microbial co-occurrence pattern than bulk soil[J]. Soil Biology and Biochemistry, 2018, 125: 251-260. doi: 10.1016/j.soilbio.2018.07.022
[38] YU Z, LIANG K N, HUANG G H, et al. Soil bacterial community shifts are driven by soil nutrient availability along a teak plantation chronosequence in tropical forests in China[J]. Biology-Basel, 2021, 10(12): 1329.
-
期刊类型引用(1)
1. 陈晨，保文星，陈旭，景永俊，李卫军. 改进YOLOv8的学生课堂行为识别算法：DMS-YOLOv8. 计算机工程与应用. 2024(24): 222-234 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 0
- 被引次数: 1
