Transcriptome analysis of Drosophila S2 cells after infection of Listeria monocytogenes, L. grayi or L. welshimeri
-
摘要:目的
利用果蝇作为宿主比较单增李斯特菌Listeria monocytogenes、格氏李斯特菌L. grayi以及威尔斯李斯特菌L. welshimeri侵染后触发果蝇S2细胞基因表达模式的异同，为探究宿主的先天免疫防御机制以及不同种李斯特菌的致病差异机制提供一定的研究基础。
方法单增李斯特菌、格氏李斯特菌以及威尔斯李斯特菌分别侵染果蝇S2细胞3 h后，利用转录组测序技术检测果蝇S2细胞中的mRNA表达谱，并利用生物信息学方法筛选出差异表达基因(DEGs)，然后对这些基因进行GO注释和KEGG分析，最后采用qRT-PCR技术验证5个Toll和Imd通路相关的抗菌肽基因的转录水平。
结果单增李斯特菌、格氏李斯特菌和威尔斯李斯特菌侵染3 h后果蝇S2细胞中共有18个DEGs，特有DEGs分别为104、28和33个。3种李斯特菌侵染组共有的DEGs被注释到代谢过程、细胞过程、细胞器、对刺激反应和免疫系统过程等共有的GO条目；此外，单增李斯特菌侵染组特有的DEGs被注释到生物黏附、参与化学突触传递的突触前过程、生物过程的负调节、核酸结合转录因子活性和信号转导器活性等特有GO条目；格氏李斯特菌侵染组特有的DEGs被注释到多个生物体过程、细胞外区域等特有GO条目；威尔斯李斯特菌侵染组特有的DEGs则被注释到发育过程及膜包裹腔等特有GO条目。KEGG分析结果显示，3种李斯特菌侵染组的DEGs均富集到了Toll/Imd信号通路。qRT-PCR 验证5个Toll和Imd通路相关的抗菌肽基因转录水平，定量结果与测序结果基本吻合。
结论3种不同李斯特菌侵染后触发S2细胞的转录表达模式具有一定的相似性，但同时也存在不同之处，这些差异存在的信号通路可能与3种李斯特菌在S2细胞中的致病力强弱存在相关性。
Abstract:ObjectiveTo compare the gene expression patterns in S2 cells after infection of Listeria monocytogenes, L. grayi, and L. welshimeri using Drosophila as a host, and provide a research basis for exploring the innate immune defense mechanisms of the host and the differential pathogenic mechanisms of different Listeria strains.
MethodAfter infection of L. monocytogenes, L. grayi or L. welshimeri for 3 h, the mRNA expression profiles in Drosophila S2 cells were detected using transcriptome sequencing technology. Analysis of differentially expressed genes (DEGs) was performed using bioinformatics tools of GO annotation and KEGG analysis. Finally, the mRNA levels of five antimicrobial peptide genes associated with the Toll and Imd signaling pathways were validated using qRT-PCR technology.
ResultAfter infection of L. monocytogenes, L. gray or L. welshimeri, there were 18 DEGs in common, as well as 104, 28, and 33 unique DEGs in the Drosophila S2 cells respectively. The DEGs commonly existed in three Listeria infected groups were annotated in GO terms including metabolic process, cellular process, organelle, response to stimulus, immune system process and etc. The unique DEGs existed in L. monocytogenes infected group were annotated in GO terms such as biological adhesion, presynaptic process involved in chemical synaptic transmission, negative regulation of biological process, nucleic acid binding transcription factor activity, signal transducer activity, and etc. The unique DEGs in L. grayi infected group were annotated in GO terms such as multi-organism process and extracellular regions, while the unique DEGs in L. welshimeri infected group were annotated in GO terms such as developmental processes and membrane-enclosed lumen. The KEGG analysis results showed that the DEGs of three Listeria infected groups were enriched in the Toll/Imd signaling pathways. The quantitative results of qRT-PCR validation of five antimicrobial peptide genes related to Toll and Imd signaling pathways were consistent with the transcriptomic results.
ConclusionThe expression patterns in S2 cells affected by the infection of these three Listeria species exhibited some similarities but also differences. The variation of the signaling pathways may be related to the virulence of these Listeria species in S2 cells.
-
广西地处亚热带地区，光热充沛，11月份晚稻收获后农田进入空窗期，利用冬闲田发展马铃薯产业空间大。然而广西冬季雨水偏少，灌水成为制约广西冬种马铃薯产业发展的条件之一。马铃薯实际生产中盲目灌水和过量施肥现象普遍存在，而滴灌施肥可以根据作物需水需肥规律和土壤水分养分状况精确控制灌水量、施肥量和灌水施肥时间，将水分养分直接供应到根区，实现作物“按需灌水施肥”，从而提高作物产量和水分养分利用效率[1]，同时滴灌施肥也影响土壤碳组分，因此，研究合适的滴灌施肥模式将为调控土壤碳库提供新的途径。目前常用土壤可溶性有机碳、易氧化有机碳和微生物量碳、碳库管理指数等表征土壤碳库，而水肥管理会影响土壤碳库和酶活性。有研究表明，长期合理施肥显著提高土壤有机碳、易氧化有机碳、可溶性有机碳、微生物量碳含量及碳库管理指数[2]，与传统施肥相比，滴灌施肥增加各层次土壤易氧化有机碳和可溶性有机碳含量[3]。其他研究也发现，滴灌施肥对提高土壤易氧化有机碳有积极的作用[4-6]。土壤水分含量影响土壤有机碳矿化速率和外界有机碳分解速率[7]，从而使土壤有机碳的含量发生变化。俞华林等[8]发现，适量灌水会增加土壤有机碳含量，但少量或过量灌水降低土壤有机碳矿化速率。当土壤水分过量时，土壤透气性和土壤微生物生长环境变差，土壤中有机碳不易被土壤中的微生物分解，而外界的碳源则易被微生物降解腐烂成有机物质，原有的有机碳不会分解且外源有机碳增加，从而使土壤有机碳含量增加[9]。水肥管理也会影响土壤酶如蔗糖酶、纤维素酶和过氧化氢酶等酶活性，滴灌施肥有利于提高土壤中酶活性[10]，而土壤酶活性会影响土壤碳组分。研究发现，各种形态有机碳组分与土壤蔗糖酶和纤维素酶活性均呈显著的正相关关系[11]。
近年来国内外学者较多关注滴灌施肥对马铃薯生长、产量、品质和水分利用效率的影响，而滴灌灌水量和滴灌施肥比例协同作用对种植马铃薯土壤碳库管理指数的影响研究较少，且土壤酶活性如何影响土壤有机碳组分和碳库管理指数有待深入研究。因此，在南宁市防雨棚内开展不同滴灌灌水量和滴灌施肥比例的田间试验，通过测定马铃薯收获后土壤有机碳及活性组分和酶活性，计算土壤碳库管理指数(Carbon pool management index，CPMI)，分析土壤有机碳组分和碳库管理指数与酶活性之间的关系，以获得种植马铃薯土壤有机碳库调控的水肥管理模式，并揭示土壤酶活性对土壤有机碳组分和碳库管理指数的影响。
1. 材料与方法
1.1 试验地点与材料
田间试验在南宁市广西大学校内移动防雨棚中进行，该移动棚通风、透光，可以保障作物生长期间自然光照和温度，通过电控传感器在降雨时遮盖，非降雨时移开。供试土壤为赤红壤，pH6.60(水土质量比2.5∶1.0，pH计法)，有机质10.6 g·kg−1(重铬酸钾容量法–外加热法)，全氮0.99 g·kg−1(半微量开氏法)，碱解氮53.6 mg·kg−1(NaOH碱解扩散法)，速效磷68.7 mg·kg−1[0.05 mol·L−1 HCl–0.025 mol·L−1 H2SO4浸提，比色法]，速效钾217.9 mg·kg−1(1 mol·L−1中性NH4OAc浸提，火焰光度法)，田间持水量30.5%(环刀法)，容重1.4 g·cm−3(室内环刀法)[12]。供试马铃薯品种为费乌瑞它。
1.2 试验方法
依据马铃薯在不同时期的需水规律及前人研究结果[13-14]，试验设高、低2种滴灌灌水量，其中，高灌水量：苗期、块茎形成期、块茎膨大期和淀粉积累期土壤含水量分别保持在田间持水量的60%~70%、70%~80%、75%~85%和50%~60%；低灌水量：苗期、块茎形成期、块茎膨大期和淀粉积累期土壤含水量分别保持在田间持水量的50%~60%、60%~70%、70%~80%和40%~50%。设3种滴灌施肥比例，即NK100-0：N、K肥以100%作基肥土施，不追肥；NK70-30：N、K肥以70%作基肥土施，30%作滴灌追肥(苗期7.5%，块茎形成期15%，块茎膨大期7.5%)；NK50-50：N、K肥以50%作基肥土施，50%作滴灌追肥(苗期12.5%，块茎形成期25%，块茎膨大期12.5%)。试验共设6个处理，具体如表1所示，每个处理重复3次，共18个小区，每小区面积8.64 m2(3.6 m×2.4 m)。
表 1 田间试验处理及N、K肥的基、追肥比例Table 1. Treatments for field experiment and radio of base fertilizer and topdressing for N,K fertilizer处理
Treatment滴灌灌水量
Drip irrigation amount滴灌施肥比例
Fertigation ratio基肥/%
Base fertilizer追肥 Topdressing/% 苗期
Seedling stage块茎形成期
Tuber formation stage块茎膨大期
Tuber expansion stageT1 高灌水量
High irrigation amountNK100-0 100 0 0 0 T2 NK70-30 70 7.5 15 7.5 T3 NK50-50 50 12.5 25 12.5 T4 低灌水量
Low irrigation amountNK100-0 100 0 0 0 T5 NK70-30 70 7.5 15 7.5 T6 NK50-50 50 12.5 25 12.5 各小区均施用化学肥料N 150 kg·hm−2、P2O5 90 kg·hm−2和K2O 300 kg·hm−2，以及堆沤后牛粪15 t·hm−2。氮肥用尿素[w(N)为 46.4%]，磷肥用钙镁磷肥[w(P2O5)为18.0%]，钾肥用硫酸钾[w(K2O)为52.0%]。牛粪中养分：w(有机质)为14.3%、w(N)为0.76%、w(P2O5)为0.85%、w(K2O)为0.59%。牛粪和钙镁磷肥全部做基肥土施。灌溉方式采用地表滴灌，滴头流量一致，滴头设在马铃薯植株两侧，用水表计量灌水。N肥和K肥按上述施肥方式施用，事先按设计要求配好肥料溶液，通过滴灌带进行灌溉施肥，灌溉方法采用交替滴灌。
1.3 田间试验及管理
于2017年11月4日将沤熟牛粪施入试验小区，11月5日翻地，11月10日将部分尿素、钙镁磷肥以及硫酸钾作为基肥土施。11月11日切马铃薯块茎，每个种薯块茎留2~3个芽眼，用质量分数为0.5％的高锰酸钾溶液和丁硫克百威水溶液浸泡拌种后晾干，11月14日播种，12月4日移栽或补齐未发芽位置的马铃薯苗。用TRIME-PICO-IPH TDR水分测定仪(德国IMKO)测定土壤含水量，确保土壤含水量在试验设定范围内。12月6日施苗肥，12月11日进行第一次中耕培土。12月20日施块茎形成肥，12月25日进行第2次培土(培土到植株附近，芽块顶部到垄背顶部达到15~20 cm左右，做成梯形垄)。2018年1月4日，施块茎膨大肥，2月8日喷农药(棉铃虫核型多角体病毒，预防马铃薯晚疫病)，试验于2018年3月5日收获马铃薯。
1.4 土壤样品采集与测定
于3月6日(马铃薯收获后次日)用5点法在马铃薯相邻植株中间采集0~20 cm耕作层土壤，将土样混匀，迅速运回实验室，部分新鲜土样过孔径2 mm筛网，除去根系、砂石等后，保存于4 ℃冰箱，直接用于土壤有机碳组分和酶活性的测定。剩余土样风干后过0.149 mm筛后进行土壤总有机碳含量的测定。
土壤总有机碳(Total organic carbon，TOC)含量用高温外加热重铬酸钾氧化–容量法测定[12]；活性有机碳(Labile organic carbon，LOC)含量用浓度为333 mmol·L−1的高锰酸钾溶液氧化土样，并于565 nm下通过测定光密度得到[12]；微生物量碳(Microbial biomass carbon，MBC)和可溶性有机碳(Dissolved organic carbon，，DOC)含量分别用三氯甲烷熏蒸和不用三氯甲烷熏蒸后，用浓度为0.5 mol·L−1硫酸钾溶液提取，采用高温外加热重铬酸钾氧化–容量法测定[12]。
土壤蔗糖酶活性用3，5–二硝基水杨酸溶液比色法测定，其活性以1 g干土1 d生成葡萄糖的质量(mg)表示；纤维素酶活性也用3，5–二硝基水杨酸溶液比色法测定，以1 g干土3 d生成葡萄糖的质量(mg)表示1个活性单位(U)；过氧化氢酶活性用高锰酸钾滴定法测定，其活性以1 g干土消耗浓度为0.02 mol·L−1的KMnO4溶液体积(mL)表示，3种酶活性测定的具体操作步骤见《土壤酶及其研究法》[15]。
土壤碳库指数(Carbon pool index，CPI)和碳库管理指数的计算参照杜爱林等[16]的方法进行。
1.5 数据统计分析
试验数据采用Excel 2016和SPSS 24.0软件进行分析。方差分析包括滴灌灌水量和滴灌施肥比例主效应，以及它们之间的交互效应。用Duncan’s法对不同处理进行多重比较。用Pearson相关系数表示土壤总有机碳及其组分和碳库管理指数与酶活性之间的相关性。
2. 结果与分析
2.1 不同水肥处理对土壤有机碳及其组分的影响
由表2方差分析可知，滴灌灌水量和滴灌施肥。比例对土壤总有机碳(TOC)影响显著(P<0.05)。土壤TOC质量分数在5.46~7.12 g·kg−1之间。多重比较结果显示，相同滴灌施肥比例下，高灌水量土壤TOC含量显著高于低灌水量土壤。在高灌水量下，NK50-50施肥处理土壤TOC含量分别比NK100-0和NK70-30处理提高15.2%和7.1%。在低灌水量下，NK50-50施肥处理土壤TOC含量比NK100-0和NK70-30处理提高12.6%和9.8%。
表 2 不同处理对土壤有机碳及其组分的影响1)Table 2. Effects of different treatments on soil organic carbon and its components处理
Treatment滴灌灌水量
Drip irrigation amount滴灌施肥比例
Fertigation ratiow/(g·kg−1) w/(mg·kg−1) 总有机碳
Total organic carbon(TOC)活性有机碳
Labile organic carbon(LOC)可溶性有机碳
Dissolved organic carbon(DOC)微生物量碳
Microbial biomass carbon(MBC)T1 高灌水量
High irrigation amountNK100-0 6.18±0.15bc 0.44±0.03b 323.0±57.0ab 374.8±25.3ab T2 NK70-30 6.65±0.24ab 0.49±0.01b 369.5±27.5a 384.8±20.3a T3 NK50-50 7.12±0.24a 0.55±0.02a 328.7±14.8ab 370.6±3.1b T4 低灌水量
Low irrigation amountNK100-0 5.46±0.15d 0.43±0.01b 189.5±49.8b 325.8±8.5b T5 NK70-30 5.60±0.16cd 0.44±0.01b 241.3±93.5ab 343.0±9.6ab T6 NK50-50 6.15±0.18bc 0.47±0.01b 215.6±7.6ab 324.1±18.7b 显著性检验
(P值)
Significance test
(P value)滴灌灌水量 Drip irrigation amount 0.004 0.008 0.011 0.005 滴灌施肥比例 Fertigation ratio 0.001 0.003 0.626 0.567 滴灌灌水量×滴灌施肥比例
Drip irrigation amount × Fertigation ratio0.674 0.125 0.979 0.975 1) 同列数据后的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among treatments (P<0.05，Duncan’s test)滴灌灌水量和滴灌施肥比例对土壤活性有机碳(LOC)影响显著(P<0.05)(表2)。土壤LOC质量分数介于0.43~0.55 g·kg−1之间。NK50-50下，高灌水量土壤LOC含量显著高于低灌水量土壤。高灌水量下，NK50-50处理土壤LOC含量较NK100-0增加25.0%，且差异显著，而在低灌水量下，不同滴灌施肥比例土壤LOC含量之间的差异并不显著。
滴灌灌水量对于土壤可溶性有机碳(DOC)影响显著(P<0.05)(表2)。土壤DOC质量分数介于189.5~369.5 mg·kg−1之间。相同滴灌施肥比例下，高灌水量土壤DOC含量与低灌水量土壤之间的差异不显著，相同滴灌灌水量下，不同滴灌施肥比例土壤DOC含量之间的差异也不显著。低灌水量下，NK70-30土壤DOC含量比NK100-0高27.3%。
滴灌灌水量对土壤微生物量碳(MBC)影响显著(P<0.05)(表2)。土壤MBC质量分数在324.1~384.8 mg·kg−1之间。相同滴灌施肥比例下，高灌水量土壤MBC含量与低灌水量土壤MBC含量之间的差异不显著；相同滴灌灌水量下，不同滴灌施肥比例土壤MBC含量之间的差异也不显著。
此外，滴灌灌水量和滴灌施肥比例之间的交互作用对土壤TOC、LOC、DOC和MBC含量的影响均不显著(P>0.05)。T3处理土壤TOC和LOC含量相对较高，而T2处理土壤DOC和MBC含量相对较高。在相同滴灌施肥比例下，高灌水量土壤有机碳及其组分较低灌水量土壤高。
2.2 不同水肥处理对土壤酶活性的影响
由表3方差分析可知，滴灌灌水量对土壤蔗糖酶活性影响显著(P<0.05)，但滴灌施肥比例和滴灌灌水量×滴灌施肥比例对土壤蔗糖酶活性的影响并不显著(P>0.05)。多重比较结果显示，NK100-0和NK50-50下，高灌水量土壤蔗糖酶活性较相应低灌水量土壤分别提高18.9%和18.2%，但差异不显著。土壤蔗糖酶活性以T3处理较高。
表 3 不同处理对土壤酶活性的影响1)Table 3. Effects of different treatments on soil enzyme activity处理
Treatment滴灌灌水量
Drip irrigation amount滴灌施肥比例
Fertigation ratio蔗糖酶活性/(mg·g−1·d−1)
Sucrase activity纤维素酶活性/U
Cellulase activity过氧化氢酶活性/(mL·g−1)
Catalase activityT1 高灌水量
High irrigation amountNK100-0 7.17±0.36ab 0.73±0.06a 0.45±0.03a T2 NK70-30 7.29±0.14a 0.75±0.04a 0.47±0.03a T3 NK50-50 7.39±0.24a 0.75±0.03a 0.46±0.02a T4 低灌水量
Low irrigation amountNK100-0 6.03±0.56b 0.64±0.06a 0.39±0.06a T5 NK70-30 6.30±0.18ab 0.67±0.06a 0.45±0.04a T6 NK50-50 6.25±0.44ab 0.66±0.03a 0.44±0.04a 显著性检验
(P值)
Significance
Test
(Pvalue)滴灌灌水量 Drip irrigation amount 0.003 0.062 0.311 滴灌施肥比例 Fertigation ratio 0.799 0.906 0.602 滴灌灌水量×滴灌施肥比例
Drip irrigation amount × Fertigation ratio0.969 0.999 0.873 1)同列数据后的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among treatments (P<0.05，Duncan’s test)滴灌灌水量、滴灌施肥比例以及滴灌灌水量×滴灌施肥比例对土壤纤维素酶和过氧化氢酶活性的影响均不显著(P>0.05)(表3)。各处理土壤纤维素酶和过氧化氢酶活性之间的差异不显著。
2.3 不同水肥处理对土壤碳库管理指数的影响
由表4方差分析可知，滴灌灌水量对土壤碳库指数影响显著(P<0.05)，但对碳库管理指数(CPMI)影响不显著(P>0.05)。滴灌施肥比例对土壤CPI和CPMI影响均不显著(P>0.05)。滴灌灌水量×滴灌施肥比例对土壤CPI和CPMI均有显著影响(P<0.05)。
在相同滴灌施肥比例下，高灌水量土壤CPI和CPMI均高于低灌水量土壤。高灌水量下，NK50-50施肥处理土壤的CPI和CPMI比NK100-0分别提高15.1%和25.8%；低灌水量下，NK50-50施肥处理土壤的CPI和CPMI比NK100-0分别提高12.6%和8.4%。土壤CPI和CPMI以T3处理最高。
2.4 土壤有机碳及其组分和碳库管理指数与酶活性之间的相互关系
土壤有机碳及其组分和碳库管理指数与酶活性之间的相关性分析结果如表5所示。土壤TOC、DOC、MBC和CPI均与蔗糖酶活性之间呈显著正相关(相关系数分别为0.61，0.48，0.46和0.60)，而土壤碳库指数与其他2种酶活性之间的相关性均不显著。
表 4 不同处理对土壤碳库管理指数的影响1)Table 4. Effects of different treatments on soil carbon pool management index处理
Treatment滴灌灌水量
Drip irrigation amount滴灌施肥比例
Fertigation ratio碳库指数
Carbon pool index
(CPI)碳库管理指数
Carbon pool management index
(CPMI)T1 高灌水量
High irrigation amountNK100-0 1.26±0.03bc 121.65±7.57b T2 NK70-30 1.35±0.05ab 134.36±4.23b T3 NK50-50 1.45±0.05a 153.04±5.71a T4 低灌水量
Low irrigation amountNK100-0 1.11±0.03d 120.08±4.93b T5 NK70-30 1.14±0.03cd 122.43±4.23b T6 NK50-50 1.25±0.04bc 130.19±2.63b 显著性检验
(P值)
Significance test
(P value)滴灌灌水量 Drip irrigation amount 0.001 0.111 滴灌施肥比例 Fertigation ratio 0.113 0.194 滴灌灌水量×滴灌施肥比例
Drip irrigation amount × Fertigation ratio0.000 0.001 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1)Different lowercase letters in the same column indicate significant differences among treatments (P<0.05，Duncan’s test)表 5 土壤有机碳及其组分含量和碳库管理指数与酶活性的相关性分析1)Table 5. Correlation analyses of soil organic carbon and fraction contents and carbon pool management index with enzyme activity指标
Index蔗糖酶
Sucrase纤维素酶
Cellulase过氧化氢酶
Catalase总有机碳 Total organic carbon (TOC) 0.61** 0.24 0.33 活性有机碳 Labile organic carbon (LOC) 0.29 0.23 0.14 可溶性有机碳 Dissolved organic carbon (DOC) 0.48* 0.02 0.29 微生物量碳 Microbial biomass carbon (MBC) 0.46* 0.29 0.03 碳库指数 Carbon pool index (CPI) 0.60** 0.24 0.31 碳库管理指数 Carbon pool management index (CPMI) 0.23 0.24 0.20 1)“*”和“**”分别表示达0.05和0.01水平的显著相关(n=3，Pearson法)
1)“*” and “**” indicate significant correlations at 0.05 and 0.01 levels, respectively(n=3, Pearson method)3. 讨论与结论
本研究表明，在相同滴灌灌水量下，与NK100-0相比，NK50-50和NK70-30滴灌施肥下的土壤TOC、LOC和DOC含量都有所提高。NK100-0处理土壤总有机碳及其组分含量等都较低，原因是该处理的肥料全部用作基肥施入土壤，后期养分供应不足，而且部分N肥易通过挥发或反硝化损失，影响N肥施用效果。而NK50-50和NK70-30交替灌溉追施N、K肥使两侧根区土壤处于交替干燥和湿润状态，在提供作物所需水分和养分的同时，使根区土壤处于良好的通气状态，为土壤微生物提供了有益的生存条件，故交替滴灌施肥比例的增加有利于土壤有机碳组分的增加[17]；再加上在马铃薯成熟期化学N、K肥配施能够促进作物根系生长，通过增加地下生物量来提高土壤有机质含量，进而有助于有机碳及其组分的增加[18]。
本研究表明，滴灌灌水量对于土壤有机碳及其组分的影响都达到显著水平。在相同滴灌施肥比例下，高灌水量土壤TOC、LOC、DOC和MBC含量都高于低灌水量土壤。相关研究发现，土壤含水量从土壤含水量<50%变成50%~100%时，土壤微生物活性通常会受到抑制，使土壤有机碳矿化分解缓慢，进而使土壤有机碳及其组分增加[19]。
土壤碳库管理指数作为反映土壤碳素动态变化灵敏而有效的指标，与土壤有效碳的关系密切，可反映和评估土壤碳素动态变化[20]。土壤碳库管理指数可用于衡量土壤质量，CPMI值越大，表明土壤质量越好[21]。本研究表明，在相同的灌水量下，NK50-50施肥处理土壤的CPI和CPMI均高于NK100-0，说明提高滴灌施肥比例会增加土壤CPMI，这与滕秋梅等[22]和张鹏等[23]的研究结果一致。说明适量N、K肥的加入可促进植物生长，增强土壤养分循环功能。究其原因，可能是N、K肥施入后主要提高的是LOC含量，导致碳库管理指数较高。凋落物和根系分泌物转化为有机质时，一部分有机质活化后为植物生长提供养分，一部分有机质转化为惰性碳库固存下来，这2个比例维持在一定范围内[24]。
土壤酶在土壤养分周转及土壤功能稳定中有重要作用。在影响土壤酶活性因子中，土壤水分对酶活性的影响具有异质性。本研究表明，土壤蔗糖酶活性在高灌水量下较高，说明灌水量的增加会提高土壤蔗糖酶活性，这与田幼华等[25]、高丽敏等[26]研究结果一致，但与万忠梅等[27]的研究结果相反，这可能是由于不同作物的需水量不同。而本研究结果可能是因为土壤水分的增加，加快了微生物胞外酶和底物的运输速率，可为酶促反应提供良好的反应环境，进而蔗糖酶和纤维素酶活性得到提高[28]。但滴灌施肥比例对土壤蔗糖酶、纤维素酶、过氧化氢酶活性的影响不显著，这与大多数研究结果并不相同。这可能是由于本研究是在相同的施肥量下，滴灌施肥比例对各种酶活性的影响较小；而大多数研究是通过设置不同的施肥梯度实现的。
蔗糖酶对蔗糖分解的催化作用具有专一性，能将土壤中蔗糖分子分解成果糖和葡萄糖，为土壤微生物提供营养物质，促进土壤有机碳积累与分解转化，从而直接或者间接地影响有机碳矿化过程[29]。本研究表明，土壤总有机碳与蔗糖酶活性呈极显著正相关，以往研究也有相似的结果[30]，说明土壤蔗糖酶活性影响土壤有机碳的积累。本研究发现，土壤有机碳组分与纤维素酶和过氧化氢酶活性之间的关系不显著，然而，马瑞萍等[11]对黄土高原不同植物群落土壤团聚体中有机碳和酶活性研究表明，土壤纤维素酶活性与各种组分有机碳之间的关系均呈显著正相关。张英英[31]研究发现，不同耕作措施下甘肃旱地农田0~30 cm土层土壤活性有机碳与纤维素酶和过氧化物酶活性之间的关系呈显著正相关，与本试验结果不同，可能是试验条件和土壤类型不同的原因所致。
综上所述，在高灌水量(苗期、块茎形成期、块茎膨大期和淀粉积累期土壤含水量分别保持在田间持水量的60%~70%、70%~80%、75%~85%和50%~60%)和NK50-50施肥处理(N、K肥以50%作基肥土施，50%作滴灌追肥)下土壤总有机碳及其组分、蔗糖酶活性和碳库管理指数较高，因此，高灌水量和N、K肥基、追肥比50∶50处理为广西冬种马铃薯种植土壤有机碳库调控的水肥耦合模式。此外，土壤TOC、DOC、MBC含量和CPI均与蔗糖酶活性呈显著正相关，说明土壤蔗糖酶活性会影响土壤有机碳及其组分。
-
![]()
图 1 单增李斯特菌、格氏李斯特菌和威尔斯李斯特菌侵染果蝇S2细胞后的增殖情况
各图中，柱子上方的不同小写字母表示差异显著（P<0.05，LSD法）。
Figure 1. The proliferation of Listeria monocytogenes, L. grayi and L. welshimeri after their infection in Drosophila S2 cells
In each figure, different lowercase letters on bars indicate significant differences (P<0.05, LSD method).
![]()
图 2 单增李斯特菌(LM)、格氏李斯特菌(LG)和威尔斯李斯特菌(LW)侵染3 h后果蝇S2细胞中差异表达的基因数量
Figure 2. Number of DEGs in Drosophila S2 cells after Listeria monocytogenes (LM), L. grayi (LG) and L. welshimeri (LW) infection for 3 h
![]()
图 3 单增李斯特菌侵染后果蝇S2细胞特有DEGs的GO注释分析
Figure 3. GO annotation of unique DEGs in Drosophila S2 cells after Listeria monocytogenes infection
![]()
图 4 格氏李斯特菌侵染后果蝇S2细胞特有DEGs的GO注释分析
Figure 4. GO annotation of unique DEGs in Drosophila S2 cells after Listeria grayi infection
![]()
图 5 威尔斯李斯特菌侵染后果蝇S2细胞特有DEGs的GO注释分析
Figure 5. GO annotation of unique DEGs in Drosophila S2 cells after Listeria welshimeri infection
![]()
图 6 单增李斯特菌、格氏李斯特菌以及威尔斯李斯特菌侵染后果蝇S2细胞共有DEGs的GO注释分析
Figure 6. GO annotation of common DEGs in Drosophila S2 cells after Listeria monocytogenes, L. grayi and L. welshimeri infection
![]()
图 7 单增李斯特菌侵染3 h后果蝇S2细胞差异表达抗菌肽基因的 qRT-PCR 验证
Figure 7. qRT-PCR identification of differentially expressed antibacterial peptide genes in Drosophila S2 cells after Listeria monocytogenes infection for 3 h
![]()
图 8 格氏李斯特菌侵染3 h后果蝇S2细胞差异表达抗菌肽基因的 qRT-PCR 验证
Figure 8. qRT-PCR identification of differentially expressed antibacterial peptide genes in Drosophila S2 cells after Listeria grayi infection for 3 h
![]()
图 9 威尔斯李斯特菌侵染3 h后果蝇S2细胞差异表达抗菌肽基因的 qRT-PCR 验证
Figure 9. qRT-PCR identification of differentially expressed antibacterial peptide genes in Drosophila S2 cells after Listeria welshimeri infection for 3 h
表 1 qRT-PCR检测基因引物
Table 1 Primers for qRT-PCR detection
基因
Gene正向引物序列(5′→3′)
Forward primer sequence反向引物序列(5′→3′)
Reverse primer sequence产物长度/bp
Product lengthPutative oxidoreductase GGATAAAGGCATTATTGAACGT CAACTACTGGGTGCGGGAT 227 Species-specific
Oxidoreductase-codingAACTCCGTCCATCCATCACC TCGTCCGAGGCTAGGAATAAG 173 ScrA ACCACCAAGACAAGCAGTTAGA CATGTAGAAATGATTGCCCAAA 237 DmRp49 GACAGTATCTGATGCCCAACA CTTCTTGGAGGAGACGCCGT 170 DmAttD AATAACGCCGCTCTGAATG ATTGAGTCCTCCGCCAAA 75 DmDptA ACTGAATGGAGGATATGG AATCTGTAGGTGTAGGTG 90 DmDef ATGAAGTTCTTCGTTCTC CATGATCCTCTGGAATTG 100 DmDptB TTCACCGCTAGTCTTCTA GGGATCAGGATAGGGATAA 75 DmDro TCCTGCTGCTTGCTTGCG ATGGGTGGTGGCCTGATGG 163 DmDrs ATGATGCAGATCAAGTAC AGGGACCCTTGTATCTTC 109 DmMtk TACATCAGTGCTGGCAGAG TTGGTTGGTTAGGATTGAAGG 83 
下载: 导出CSV
表 2 单增李斯特菌侵染后果蝇S2细胞特有上调基因的KEGG富集分析
Table 2 KEGG enrichment of unique up-regulated genes in Drosophila S2 cells after Listeria monocytogenes infection
通路
PathwayP 所含基因
Contained gene背腹轴形成 Dorso-ventral axis formation 0.019 aos 铂类药物抗药性 Platinum drug resistance 0.041 GstE9 化学诱发癌变 Chemical carcinogenesis 0.042 GstE9 药物代谢−细胞色素 P450 Drug metabolism-cytochrome P450 0.042 GstE9 细胞色素 P450代谢异物质 Metabolism of xenobiotics by cytochrome P450 0.042 GstE9 谷胱甘肽代谢 Glutathione metabolism 0.049 GstE9
下载: 导出CSV
表 3 单增李斯特菌侵染后果蝇S2细胞特有下调基因的KEGG富集分析
Table 3 KEGG enrichment of unique down-regulated genes in Drosophila S2 cells after Listeria monocytogenes infection
通路
PathwayP 所含基因
Contained gene苯丙氨酸代谢 Phenylalanine metabolism 0.042 CG11796 泛醌和其他萜类醌的生物合成 Ubiquinone and other terpenoid-quinone biosynthesis 0.047 CG11796 胰液分泌 Pancreatic secretion 0.050 iotaTry, RyR
下载: 导出CSV
表 4 格氏李斯特菌侵染后果蝇S2细胞特有上调基因的KEGG富集分析
Table 4 KEGG enrichment of unique up-regulated genes in Drosophila S2 cells after Listeria grayi infection
通路
PathwayP 所含基因
Contained geneToll和Imd信号传导途径 Toll and Imd signaling pathway 0.000 Def, Drs, CecC, CecB, CecA1, pirk 补体和凝血级联反应 Complement and coagulation cascades 0.003 Hml 细胞外基质−受体相互作用 ECM-receptor interaction 0.016 Hml
下载: 导出CSV
表 5 格氏李斯特菌侵染后果蝇S2细胞特有下调基因的KEGG富集分析
Table 5 KEGG enrichment of unique down-regulated genes in Drosophila S2 cells after Listeria grayi infection
通路
PathwayP 所含基因
Contained gene2型糖尿病 Type II diabetes mellitus 0.010 CG7069 癌细胞中的中心碳代谢 Central carbon metabolism in cancer 0.012 CG7069 丙酮酸代谢 Pyruvate metabolism 0.015 CG7069 糖酵解/糖异生 Glycolysis/gluconeogenesis 0.017 CG7069 胰高血糖素信号传导途径 Glucagon signaling pathway 0.021 CG7069 病毒致癌 Viral carcinogenesis 0.043 CG7069 嘌呤代谢 Purine metabolism 0.046 CG7069
下载: 导出CSV
表 6 威尔斯李斯特菌侵染后果蝇S2细胞特有下调基因的KEGG富集分析
Table 6 KEGG enrichment of unique down-regulated genes in Drosophila S2 cells after Listeria welshimeri infection
通路
PathwayP 所含基因
Contained gene抗原处理与呈递 Antigen processing and presentation 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 军团菌病 Legionellosis 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 甲型流感 Influenza A 0.000 Hsp70Bbb, deltaTry, Hsp70Bc, Hsp68, Hsp70Bb 弓形虫感染 Toxoplasmosis 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 寿命调控途径−多个物种 Longevity regulating pathway-multiple species 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 雌激素信号传导途径 Estrogen signaling pathway 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 麻疹 Measles 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb MAPK信号传导途径 MAPK signaling pathway 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 剪接体 Spliceosome 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 内吞作用 Endocytosis 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb 内质网蛋白加工 Protein processing in endoplasmic reticulum 0.000 Hsp70Bbb, Hsp68, Hsp70Bc, Hsp70Bb
下载: 导出CSV
表 7 单增李斯特菌、格氏李斯特菌以及威尔斯李斯特菌侵染后果蝇S2细胞共有下调基因的KEGG富集分析
Table 7 KEGG enrichment of common down-regulated genes in Drosophila S2 cells after Listeria monocytogenes, L. grayi and L. welshimeri infection
通路
PathwayP 所含基因
Contained gene咖啡因代谢 Caffeine metabolism 0.007 AOX2 维生素B6代谢 Vitamin B6 metabolism 0.007 CG14212 牛磺酸和低牛磺酸代谢 Taurine and hypotaurine metabolism 0.010 Ggt-1 花生四烯酸代谢 Arachidonic acid metabolism 0.020 Ggt-1 矿物质吸收 Mineral absorption 0.024 Ctr1B 集合管酸分泌 Collecting duct acid secretion 0.037 Vha100-4 类风湿性关节炎 Rheumatoid arthritis 0.050 Vha100-4
下载: 导出CSV
-
[1] COLLINS M D, WALLBANKS S, LANE D J, et al. Phylogenetic analysis of the genus Listeria based on reverse transcriptase sequencing of 16S rRNA[J]. International Journal of Systematic Bacteriology, 1991, 41(2): 240-246. doi: 10.1099/00207713-41-2-240
[2] 何冬梅, 邓峰, 赖蔚苳, 等. 单核细胞增生李斯特菌生物学研究进展[J]. 华南预防医学, 2006, 32(6): 26-29. [3] 柴文琴. 李斯特菌属不同种型在模拟食品加工环境相关胁迫下的存活率与耐药性研究[D]. 杭州: 浙江大学, 2020. [4] KAPTCHOUANG TCHATCHOUANG C D, FRI J, DE SANTI M, et al. Listeriosis outbreak in South Africa: A comparative analysis with previously reported cases worldwide[J]. Microorganisms, 2020, 8(1): 135. doi: 10.3390/microorganisms8010135
[5] CHARLIER C, PERRODEAU É, LECLERCQ A, et al. Clinical features and prognostic factors of listeriosis: The MONALISA national prospective cohort study[J]. The Lancet Infectious Diseases, 2017, 17(5): 510-519. doi: 10.1016/S1473-3099(16)30521-7
[6] 中华人民共和国国家卫生健康委员会, 国家市场监督管理总局. 食品安全国家标准 预包装食品中致病菌限量: GB 29921—2021[S]. 北京: 中国标准出版社, 2021. [7] 中华人民共和国家卫生和计划生育委员会. 食品安全国家标准 食品中致病菌限量: GB 29921—2013[S]. 北京: 中国标准出版社, 2014. [8] VÁZQUEZ-BOLAND J A, KUHN M, BERCHE P, et al. Listeria pathogenesis and molecular virulence determinants[J]. Journal of Personalized Medicine, 2001, 14(3): 584-640.
[9] 李钊, 刘阳泰, 李卓思, 等. 单增李斯特菌毒力因子及调控机制研究进展[J]. 食品与发酵工业, 2024, 50(11): 327-335. [10] IRETON K, MORTUZA R, GYANWALI G C, et al. Role of internalin proteins in the pathogenesis of Listeria monocytogenes[J]. Molecular Microbiology, 2021, 116(6): 1407-1419. doi: 10.1111/mmi.14836
[11] ZHOU H, COVENEY A P, WU M, et al. Activation of both TLR and NOD signaling confers host innate immunity-mediated protection against microbial infection[J]. Frontiers in Immunology, 2019, 9: 3082. doi: 10.3389/fimmu.2018.03082
[12] MITCHELL G, CHENG M I, CHEN C, et al. Listeria monocytogenes triggers noncanonical autophagy upon phagocytosis, but avoids subsequent growth-restricting xenophagy[J]. Proceedings of the National Academy of Sciences of the United States of America, 2017, 115(2): E210-E217.
[13] CHEN G Y, PENSINGER D A, SAUER J D. Listeria monocytogenes cytosolic metabolism promotes replication, survival, and evasion of innate immunity[J]. Cellular Microbiology, 2017, 19(10): e12762-e12771. doi: 10.1111/cmi.12762
[14] GABALLA A, GUARIGLIA-OROPEZA V, WIEDMANN M, et al. Cross talk between SigB and PrfA in Listeria monocytogenes facilitates transitions between extra- and intracellular environments[J]. Microbiology and Molecular Biology Reviews, 2019, 83(4): e00034-19.
[15] BIERNE H, HAMON M, COSSART P. Epigenetics and bacterial infections[J]. Cold Spring Harbor Perspectives in Medicine, 2012, 2(12): a010272
[16] MINÁROVITS J. Microbe-induced epigenetic alterations in host cells: The coming era of patho-epigenetics of microbial infections: A review[J]. Acta Microbiologica et Immunologica Hungarica, 2009, 56(1): 1-19. doi: 10.1556/AMicr.56.2009.1.1
[17] HAMON M A, COSSART P. Histone modifications and chromatin remodeling during bacterial infections[J]. Cell Host & Microbe, 2008, 4(2): 100-109.
[18] COSTA A C, PINHEIRO J, REIS S A, et al. Listeria monocytogenes interferes with host cell mitosis through its virulence factors InlC and ActA[J]. Toxins, 2020, 12(6): 411-424. doi: 10.3390/toxins12060411
[19] 袁江北. 定量蛋白质组学揭示免疫细胞对李斯特菌感染的防御机制[D]. 重庆: 重庆大学, 2020. [20] IMLER J L. Overview of Drosophila immunity: A historical perspective[J]. Developmental & Comparative Immunology, 2014, 42(1): 3-15.
[21] GOTTAR M, GOBERT V, MICHEL T, et al. The Drosophila immune response against Gram-negative bacteria is mediated by a peptidoglycan recognition protein[J]. Nature, 2002, 416: 640-644. doi: 10.1038/nature734
[22] SCHMID M W, NG E Y W, LAMPIDIS R, et al. Evolutionary history of the genus Listeria and its virulence genes[J]. Systematic and Applied Microbiology, 2005, 28(1): 1-18. doi: 10.1016/j.syapm.2004.09.005
[23] PIZARRO-CERDA J, COSSART P. Listeria monocytogenes: Cell biology of invasion and intracellular growth[J]. Microbiology Spectrum, 2018, 6(6). doi: 10.1128/microbiolspec.GPP3-0013-2018.
[24] FREITAG N E, PORT G C, MINER M D. Listeria monocytogenes: From saprophyte to intracellular pathogen[J]. Nature Reviews Microbiology, 2009, 7: 623-628. doi: 10.1038/nrmicro2171
[25] 刘甜, 罗开珺. 果蝇Toll和IMD信号通路中的功能结构域[J]. 环境昆虫学报, 2011, 33(3): 388-395. doi: 10.3969/j.issn.1674-0858.2011.03.015 [26] LEULIER F, RODRIGUEZ A, KHUSH R S, et al. The Drosophila caspase Dredd is required to resist gram-negative bacterial infection[J]. EMBO Reports, 2000, 1(4): 353-358. doi: 10.1093/embo-reports/kvd073
[27] BALDWIN D N, VANCHINATHAN V, BROWN P O, et al. A gene-expression program reflecting the innate immune response of cultured intestinal epithelial cells to infection by Listeria monocytogenes[J]. Genome Biology, 2002, 4(1): R2. doi: 10.1186/gb-2002-4-1-r2
[28] 吴坤钟, 李荣松, 曹阳, 等. 英诺克李斯特菌侵染果蝇S2细胞后的转录组分析[J]. 华南农业大学学报, 2020, 41(5): 17-26. doi: 10.7671/j.issn.1001-411X.202001023 [29] 王旭. 基于比较基因组学和转录组学技术揭示单增李斯特菌毒力因子的研究[D]. 上海: 上海海洋大学, 2016. [30] PIZARRO-CERDÁ J, COSSART P. Subversion of cellular functions by Listeria monocytogenes[J]. Journal of Pathology, 2006, 208(2): 215-223. doi: 10.1002/path.1888
[31] KRYPOTOU E, SCORTTI M, GRUNDSTRÖM C, et al. Control of bacterial virulence through the peptide signature of the habitat[J]. Cell Reports, 2019, 26(7): 1815-1827. e5.
[32] SAISAWANG C, WONGSANTICHON J, KETTERMAN A J. A preliminary characterization of the cytosolic glutathione transferase proteome from Drosophila melanogaster[J]. Biochemical Journal, 2012, 442(1): 181-190. doi: 10.1042/BJ20111747
[33] XIE A, ODATE S, CHANDRAMOULY G, et al. H2AX post-translational modifications in the ionizing radiation response and homologous recombination[J]. Cell Cycle (Georgetown, Tex), 2010, 9(17): 3602-3610. doi: 10.4161/cc.9.17.12884
[34] WU M H, ZHANG X D, WEI W, et al. CRISPR/Cas9 mediated genetic resource for unknown kinase and phosphatase genes in Drosophila[J]. Scientific Reports, 2020, 10: 7383. doi: 10.1038/s41598-020-64253-4
[35] SCHERFER C, KARLSSON C, LOSEVA O, et al. Isolation and characterization of hemolymph clotting factors in Drosophila melanogaster by a pullout method[J]. Current Biology, 2004, 14(7): 625-629. doi: 10.1016/j.cub.2004.03.030
[36] 单颖. 单核细胞增生李斯特菌: 无菌斑马鱼感染模型及Mmp-9在抗细菌感染中的作用机制[D]. 杭州: 浙江大学, 2016. -
期刊类型引用(3)
1. 贾新蕾，黄增朝，杨林狄，吕静，李妍萍，简纪常，黄郁葱. 不同培养温度的鱼源海豚链球菌转录组分析. 热带生物学报. 2024(01): 109-121 . 
百度学术
2. 杨林狄，贾新蕾，黄增朝，吕静，梁华芳，黄郁葱. 银鼓鱼海豚链球菌的分离、鉴定及毒力基因检测. 大连海洋大学学报. 2023(02): 233-241 . 百度学术
3. 徐伟，施慧，汪玮，张鼎元，许文军，柴学军. 小黄鱼海豚链球菌的分离鉴定. 中国预防兽医学报. 2022(07): 725-730 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 117
- HTML全文浏览量: 30
- PDF下载量: 54
- 被引次数: 5
