Effectiveness of spraying defoliation and ripening agent on four cotton cultivars using unmanned aerial vehicle
-
摘要:目的
研究棉花收获前植保无人飞机喷施脱叶催熟剂的效果。
方法使用植保无人飞机以22.5 L·hm−2的施药液量喷施采收前30 d的棉花‘金棉18#’、‘北泉9#’、‘植金13#’和‘合信47#’，以人工机械喷雾(施药液量450 L·hm−2)为对照，测试4种棉花冠层中的雾滴覆盖率；于药后4、7、14和21 d分别进行脱叶和吐絮效果调查，计算不同棉花品种的脱叶率和吐絮率；检测4种棉花的纤维品质和产量构成因子(铃质量、衣分)。
结果雾滴覆盖率在4种棉花品种中均呈现出上部冠层显著高于中下部冠层的规律，且植株越高，中下部冠层的雾滴覆盖率越低；除‘合信47#’人工喷雾处理的吐絮率没有达到机采棉采收标准外，其余处理的脱叶率和吐絮率均达到了机采棉的采收标准；4种棉花品种的纤维品质和产量构成因子(铃质量、衣分)均不受施药器械的影响。
结论植保无人飞机可作为棉花采收前喷施脱叶催熟剂的施药器械。
Abstract:ObjectiveTo evaluate the effects of spraying defoliation and ripening agent on cotton before harvest using unmanned aerial vehicle (UAV).
MethodFour cotton cultivars including ‘Jinmian 18#’, ‘Beiquan 9#’, ‘Zhijin 13#’ and ‘Hexin 47#’ were sprayed defoliation and ripening agent with the spraying volume of 22.5 L·hm−2 at 30 days before harvest by UAV. Aritificial spraying with an electric backpack sprayer was used as cantrol at the the spraying volume of 450 L·hm−2. Droplet coverage rates in the canopies of four cotton cultivars were tested. The defoliation rates and boll opening rates were investigated at the 4th, 7th, 14th and 21st days after the application. The fiber qualities and yield components (boll weight and lint percentage) of four cotton cultivars were detected.
ResultThe droplet coverage rates of four cotton cultivars showed a rule that coverage of the upper canopy was significantly higher than that of the middle and lower parts of the canopy. The higher the plant height was, the lower the droplet coverage rate of the middle and lower canopy was. Except for the boll opening rate of artificial spraying treatment on ‘Hexin 47#’, the defoliation rates and boll opening rates of other treatments met the requirements of machinary harvesting. The fiber quality and yield components (boll weight and lint percentage) of four cotton cultivars were not affected by the type of sprayers.
ConclusionThe UAV can be applied to spray defoliation and ripening agent before cotton harvest.
-
李斯特菌属Listeria是革兰阳性菌，隶属厚壁菌门Firmicutes芽孢杆菌纲Bacilli芽孢杆菌目Bacillales，该属包括17个种。狭义上的李斯特菌属包括6个种，即单增李斯特菌L. monocytogenes、绵羊李斯特菌L. ivanovii、威尔斯李斯特菌L. welshimeri、默氏李斯特菌L. murrayi、西尔李斯特菌L. seeligeri和英诺克李斯特菌L. innocua，公认的致病菌只有2种，即单增李斯特菌和绵羊李斯特菌[1-2]。其中单增李斯特菌致病性最强，可导致人与动物患李斯特菌病[3]。绵羊李斯特菌又称伊氏李斯特菌，仅能引起动物李斯特菌病[4]，但也有引起人患李斯特菌病的个别案例[2]。此外，非常罕见英诺克李斯特菌和西尔李斯特菌引起人李斯特病的报道，且二者不具有溶血活性，因此它们一般被认为是非致病菌。其中，英诺克李斯特菌常伴随单增李斯特菌出现，并被视为后者的指示菌[5]。李斯特菌属成员的中间宿主范围较广，常常生存于蔬菜、乳类和肉类制品中；适应性强，在2~42 ℃均能正常生长[6-7]。李斯特菌病的潜伏期一般小于24 h[8]，临床表现为脑膜炎、菌血症、眼内炎、颈淋巴结炎、胆囊炎、骨髓炎及关节炎等，易感人群为老人、孕妇、新生儿和免疫缺陷患者[9]。在动物中，由绵羊李斯特菌引起的李斯特菌病可表现为流产、败血症等，但不引起脑膜炎[10]。鉴于李斯特菌有危害性大、宿主谱广和适应性强的特点，很多国家都已经采取相关措施防控食品中李斯特菌的污染，一些发达国家如美国等已将李斯特菌病定为法定传染病，但目前我国仅在部分省市对李斯特菌病开展相关检测工作[11]。
目前对李斯特菌的研究多集中在单增李斯特菌上，对其他李斯特菌的研究较少。已知的单增李斯特菌毒力因子主要有溶血素基因hly编码的李斯特溶血素O(Listerriolysin O，LLO)、plcA编码的磷脂酰肌醇特异性磷脂酶C(Phospholipase C，PI-PLC)、plcB基因编码的磷脂酰胆碱磷脂酶(Phosphatidylcholine-specific phospholipase C，PC-PLC)、SmcL基因编码的神经磷脂酶C(Nerve phoshpolipase C，NPI-PLC)、肌动蛋白A (Actin A，Act A)、内化素(Internalins，Inl)和细胞壁水解酶(p60)[12]。其中LLO是重要的毒力因子，相对分子质量为60 000，是一种孔成形毒素，与细胞膜表面分子结合后可以直接插入细胞膜上，是李斯特菌逃逸初级吞噬和双层膜吞噬体所必须的毒力因子[13-14]。此外，LLO还可以调节巨噬细胞的细胞因子表达，从而诱导树突细胞凋亡。当李斯特菌借助细胞内吞侵入宿主细胞后被包裹在吞噬体中，随后李斯特菌分泌LLO和PI-PLC裂解吞噬体膜，从而逃逸进入宿主细胞浆中[15]。但此时，李斯特菌仍存在被宿主细胞自噬所清除的危险，而Act A在李斯特菌从宿主胞质中逃逸的过程中发挥着重要的作用[16]。Act A是一种表面蛋白，当李斯特菌逃出吞噬体进入细胞浆中时，Act A可以模仿Wiskott-Aldrich综合征蛋白(Wiskott-Aldrich syndrome protein，WASP)机制，调用宿主细胞自身的肌动蛋白，这些肌动蛋白在李斯特菌一端形成1条彗星状的尾巴，赋予李斯特菌在胞质中运动的能力，并向邻近的细胞迁移，实现细胞间侵染[17]。Act A还可以与宿主细胞中的Ena/VASP和Arp2/3复合体相结合，以此逃避宿主细胞的自噬识别系统[18]。
研究表明，英诺克李斯特菌与单增李斯特菌的关系密切，两者可能具有共同的毒力祖先，在进化过程中英诺克李斯特菌丢失了编码主要毒力因子的基因[19]。英诺克李斯特菌在基因组上与单增李斯特菌非常相似，例如两者都有300多个编码运输蛋白的基因，包括80多个碳转运蛋白，两者都含有糖酵解所必需的酶[20]。典型的英诺克李斯特菌已在进化过程中缺失了LIPI-1毒力岛，但Johnson等[21]和Moura等[22]研究发现某些天然非典型英诺克李斯特菌株却含有完整的LIPI-1毒力岛，这类菌株具有溶血活性，可以主动穿过肠上皮，传播到肝脏和脾脏，但危害程度远低于单增李斯特菌。由于英诺克李斯特菌与单增李斯特菌的亲缘关系最为接近，所以在一些食品的研究中，英诺克李斯特菌常常作为单增李斯特菌的替代菌种[23]。例如，Sorrentino等[24]发现在低pH条件下3−苯基乳酸可以抑制英诺克李斯特菌的生长。另外Tremonte等[25]发现，英诺克李斯特菌在面对酸胁迫时会显著上调普遍胁迫蛋白(Universal stress protein，USP)基因的表达，表明在抵御酸胁迫时USP基因发挥着至关重要的作用。
目前，李斯特菌引起的相关并发症尚无特效药，也鲜有其侵染后寄主基因组水平的详细反应情况的研究报道。由于典型的英诺克李斯特菌被认为是非致病菌[26]，人们对宿主与其互作的机制关注较少。鉴于英诺克李斯特菌与单增李斯特菌的亲缘关系最近，我们如果能够深入挖掘其侵染后宿主基因表达的变化情况，将为宿主调控、防御单增李斯特菌的危害提供对比与参考。利用英诺克李斯特菌可以侵染果蝇Drosophila melanogaster S2细胞的特性，我们在英诺克李斯特菌侵染果蝇S2细胞3 h后，对果蝇S2细胞进行了转录组测序。
1. 材料与方法
1.1 细胞与菌株的培养
英诺克李斯特菌(ATCC，33090，由上海海洋大学提供)用李斯特菌培养基(HB4160，海博生物技术有限公司)活化和培养。果蝇S2细胞(由中国科学院上海植物生理生态研究所提供)在含有体积分数为5%胎牛血清(FBS500-S，AusGeneX，澳大利亚)的Insect SIM SF培养基(MSF1-1，信和生物公司)中培养。
1.2 英诺克李斯特菌侵染果蝇S2细胞
果蝇S2细胞铺板培养18 h后，用50 nmol/L水溶性胆固醇(134428，BIOBERRY)孵育30 min。
试验组处理：胆固醇孵育30 min后加入英诺克李斯特菌悬浮液(约20个细菌)在28 ℃条件下孵育1.5 h。然后用pH 7.4的PBS缓冲液洗涤果蝇S2细胞，在含100 µmol/L CuSO4和10 µg/mL庆大霉素的Insect SIM SF培养基基中继续培养果蝇S2细胞3 h。对照组处理：不加英诺克李斯特菌悬浮液，其他操作同试验组。培养结束后，收集试验组和对照组果蝇S2细胞进行后续试验。
1.3 RNA的抽提及转录组测序
用TRIzol试剂(9108-1，Takara，日本)抽提感染英诺克李斯特菌后的果蝇S2细胞中的总RNA。使用2100 Bioanalyze(Agilent Technologies，美国)定量总RNA模板。通过凝胶电泳确定每个样品中总RNA的完整性。根据浓度和质量值，分别选取试验组与对照组3个最好的RNA样品用于文库制备和Illumina Novaseq 6000测序。原始配对的末端读数修整后由SeqPrep( https://github.com/jstjohn/SeqPrep)和Sickle( https://github.com/najoshi/sickle)使用默认参数进行质量控制。然后使用TopHat软件( http://tophat.cbcb.umd.edu/，版本2.1.1)以定向模式将原始的读数与参考基因组比对。根据每百万条映射片段中每千个碱基的转录本所包含的片段(Fragments per kilobase of transcript per million fragments mapped，FPKM)计算每个转录本的表达水平。使用RSEM软件( http://deweylab.biostat.wisc.edu/rsem/)定量基因丰度。R统计软件包EdgeR( http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html)用于基因差异表达分析。此外，使用Goatools( https://github.com/tanghaibao/Goatools)和KOBAS 2.1.1 ( http://kobas.cbi.pku.edu.cn/download.php)进行GO和KEGG分析。对照组和试验组样品均设置3个生物学重复。
1.4 实时定量PCR检测
对差异基因进行qPCR验证：用TRIzol试剂分别提取对照组与试验组果蝇S2细胞的总RNA，使用PrimeScriptTM RT试剂盒(RR047A，Takara，日本)反转录生成cDNA，进行qPCR检测。以果蝇的rp49基因作为内参基因。对照组和试验组样品均设置3个生物学重复。
利用qPCR对英诺克李斯特菌的增殖情况进行检测：英诺克李斯特菌侵染果蝇S2细胞3 h后，收集果蝇S2细胞，用基因组提取试剂盒(D3129-01，Magen)提取总基因组，以此为模板，以果蝇rp49基因作为内参基因，用英诺克李斯特菌的凋亡抑制蛋白(Inhibitor of apoptosis protein，IAP)基因特异性引物进行qPCR检测英诺克李斯特菌的增殖情况[27]。
1.5 数据分析
采用Microsoft Office Excel2010软件绘图。试验组和对照组的差异显著性通过t检验进行统计分析。
按以下公式计算每个基因的差异表达倍数：
$$ {\text{差异倍数}}{\rm{ = }}\frac{{\text{试验组基因表达量}}}{{\text{对照组基因表达量}}}{\text{。}} $$ 2. 结果与分析
2.1 英诺克李斯特菌在果蝇S2细胞内的增殖情况
利用qPCR对英诺克李斯特菌的增殖情况进行检测，结果表明，英诺克李斯特菌的IAP基因拷贝数剧烈上升，果蝇S2细胞中侵染了大量的英诺克李斯特菌(图1)。因此，我们在英诺克李斯特菌侵染果蝇S2细胞3 h后收集果蝇S2细胞提取总RNA用于转录组测序以及后续的qPCR验证。
![]()
图 1 英诺克李斯特菌侵染果蝇S2细胞3 h后英诺克李斯特菌增殖情况图中数据为平均值±标准差；“***”表示差异达到0.001的显著水平(t检验)Figure 1. Proliferation of Listeria innocua after its infection in Drosophila S2 cells for three hoursDatum in the figure is $\overline x $ ± SD; “***” indicates significant difference at the level of 0.001 (t test)2.2 果蝇S2细胞转录组测序数据质控与评估
对照组和试验组分别平均产生了59 125 670和60 779 885条原始数据，错误率均在0.025 6%以下(表1)。果蝇S2细胞被英诺克李斯特菌侵染前基因组的GC质量分数在50%左右，被侵染3 h后基因组GC质量分数变化不大(表1)。对照组和试验组的3次生物学重复之间的相关性系数均在0.994 9以上，表明各组之间重复性良好，可进行下一步分析。
表 1 果蝇S2细胞转录组测序数据统计表Table 1. Statistics of transcriptome sequencing data for Drosophila S2 cells分组
Grouping样本编号
Sample ID原始序列读数
No. of raw reads有效序列读数
No. of clean reads错误率/%
Error rateGC质量分数/%
GC content对照组
Control group1 59 639 278 59 073 172 0.025 6 51.16 2 61 599 522 61 100 672 0.024 3 51.54 3 56 138 212 55 717 464 0.024 2 51.67 试验组
Test group4 64 092 782 63 648 024 0.024 0 51.43 5 63 143 600 62 708 658 0.024 1 51.56 6 55 103 274 54 732 148 0.024 3 51.40 2.3 果蝇S2细胞基因表达差异
英诺克李斯特菌侵染宿主果蝇S2细胞3 h后，韦恩分析表明对照组与试验组共同表达的基因有7 808个，仅在试验组表达的基因有98个，仅在对照组表达的基因有187个。
与对照组相比，试验组共有54个基因呈现显著性差异(上/下调倍数≥2)，其中21个基因上调，33个基因下调，变化幅度最大的是DmCR46081，上调了232倍，其次是DmDef(DmDefensin)、DmDro(DmDrosomycin)、DmDpt A(DmDptericin A)和DmDpt B(DmDptericin B)等抗菌肽基因。DmCR46081基因编码rRNA前体加工蛋白，提示果蝇S2细胞被英诺克李斯特菌侵染后某些基因的转录和翻译水平大幅度提升(表2)。
表 2 试验组与对照组果蝇S2细胞基因表达的差异倍数Table 2. Fold change in gene expression of Drosophila S2 cells between test group and control group基因
Gene基因描述
Gene description差异倍数
Fold change基因
Gene基因描述
Gene description差异倍数
Fold changeDmCR46081 CR46081 232.000 DmMuc30E Mucin 30E 0.473 DmCG34330 Uncharacterized protein, transcript variant A 14.649 DmWhamy WHAMM and JMY related, transcript variant C 0.464 DmDef Defensin 9.805 DmGgt-1 Gamma-glutamyl transpeptidase, transcript variant B 0.463 DmDro Drosocin, transcript variant A 7.629 DmCG15528 Uncharacterized protein 0.454 DmDpt A Diptericin A 7.447 DmCG9279 Uncharacterized protein, transcript variant C 0.447 DmDpt B Diptericin B 5.825 DmCG7848 Uncharacterized protein 0.432 DmEdin Elevated during infection 5.272 DmMuc68Ca Mucin 68Ca 0.423 DmMtk Metchnikowin 4.888 DmCG8219 Uncharacterized protein 0.414 DmCec A2 Cecropin A2 4.174 DmCG1246 Uncharacterized protein, transcript variant E 0.403 DmAtt D Attacin D 3.878 DmTsp42Ep Tetraspanin 42Ep, transcript variant A 0.390 DmAtt A Attacin A 3.820 DmAOX2 Aldehyde oxidase 2 0.384 DmAtt B Attacin B, transcript variant A 3.558 DmCG4409 Uncharacterized protein 0.384 DmPGRP-SD Peptidoglycan recognition protein SD 3.380 DmSmydA-7 SET and MYND domain containing, arthropod-
specific, member 7, transcript variant A0.346 DmPGRP-SA Peptidoglycan recognition protein SA,
Transcript variant B3.282 DmCR32636 CR32636 0.322 DmCec B Cecropin B 3.209 DmCd Cardinal 0.314 DmCG31274 Uncharacterized protein 3.131 DmVha100-4 Vacuolar H+ ATPase 100 kD subunit 4 0.303 DmTsf1 Transferrin 1, transcript variant A 2.963 DmCR45215 CR45215 0.294 DmPsf1 Psf1 2.436 DmCG17140 Uncharacterized protein, transcript variant B 0.280 DmCG43175 Uncharacterized protein, transcript variant A 2.419 DmGr59f Gustatory receptor 59f 0.278 DmObp99a Odorant-binding protein 99a, transcript variant B 2.151 DmCG34124 Uncharacterized protein 0.276 DmSp212 Serine-peptidase 212 2.080 DmCtr1B Copper transporter 1B 0.251 DmNdl Nudel 0.495 DmCG13992 Uncharacterized protein, transcript variant C 0.250 DmCep290 Cep290, transcript variant A 0.483 DmCG10514 Uncharacterized protein 0.237 DmCyp4d14 Cyp4d14, transcript variant A 0.481 DmVha68-3 Vacuolar H+ ATPase 68 kD subunit 3 0.235 DmCc2d2a Coiled-coil and C2 domain containing 2A 0.478 DmDhd Deadhead, transcript variant A 0.231 DmCG34056 Uncharacterized protein 0.475 DmCG12105 Uncharacterized protein, transcript variant B 0.230 DmMtg Mind the gap, transcript variant G 0.473 DmCG42326 Uncharacterized protein, transcript variant D 0.149 2.4 果蝇S2细胞差异表达基因的功能注释分析
GO注释分析显示，21个上调的基因共涉及10个GO条目，可分为生物进程、细胞组分和分子功能3大类，其中富集基因数目最多的GO条目分别是胞外区(9个)、应激反应(5个)、代谢进程(4个)、免疫系统进程(4个)、多细胞生物进程(3个)、催化活性(3个)、结合(3个)、细胞进程(1个)、细胞成分组织或生物合成(1个)、胞外区组件(1个)(表3)。
表 3 英诺克李斯特菌侵染后果蝇S2细胞上调基因GO注释Table 3. GO annotation of upregulated genes in Drosophila S2 cells after Listeria innocua infection分类 Category GO条目 GO term 所含基因 Contained gene 生物进程
Biological process应激反应 Response to stimulus DmDef, DmPGRP-SD, DmPGRP-SA, DmCec A2, DmCec B 细胞成分组织或生物合成
Cellular component organization or biogenesisDmCR46081 免疫系统进程 Immune system process DmPGRP-SD, DmPGRP-SA, DmCec A2, DmCec B 代谢进程 Metabolic process DmSp212, DmCR46081, DmPGRP-SD, DmPGRP-SA 细胞进程 Cellular process DmCR46081 多细胞生物进程 Multi-organism process DmPGRP-SD, DmCec A2, DmCec B 细胞组分
Cellular component胞外区 Extracellular region DmDpt B, DmDef, DmTsf1, DmAtt B, DmAtt A, DmPGRP-SD, DmAtt D, DmCec A2, DmCec B 胞外区组件 Extracellular region part DmDef 分子功能
Molecular function结合 Binding DmObp99a, DmPGRP-SD, DmPGRP-SA 催化活性 Catalytic activity DmSp212, DmPGRP-SD, DmPGRP-SA 下调的33个基因共涉及22个GO条目，其中富集基因数最多的GO条目分别是催化活性(11个)、代谢进程(10个)、细胞膜(10个)、单组织进程(8个)、细胞膜组件(8个)、结合(8个)、细胞进程(7个)、定位(6个)、细胞(6个)、细胞组件(6个)、细胞器(4个)、细胞器组件(3个)、大分子复合物(3个)、转录活性(3个)、排毒(1个)、生物进程调控(1个)、多细胞组织进程(1个)、生物调控(1个)、应激反应(1个)、胞外区(1个)、抗氧化活性(1个)、电子转运活性(1个)(表4)。
表 4 英诺克李斯特菌侵染后果蝇S2细胞下调基因GO注释Table 4. GO annotation of downregulated genes in Drosophila S2 cells after Listeria innocua infection分类 Category GO条目 GO term 所含基因 Contained gene 生物进程
Biological process多细胞组织进程 Multicellular organismal process DmGr59f 生物进程调控 Regulation of biological process DmDhd 排毒 Detoxification DmCd 应激反应 Response to stimulus DmCd 生物调控 Biological regulation DmDhd 定位 Localization DmVha68-3, DmCG9279, DmCG17140, DmCG8219, DmCtr1B, DmVha100-4 单组织进程 Single-organism process DmVha68-3, DmCG9279, DmCG7848, DmAOX2, DmCd, DmDhd, DmVha100-4, DmCyp4d14 代谢进程 Metabolic process DmCG15528, DmVha68-3, DmNdl, DmCG7848, DmAOX2, DmMtg, DmCd, DmDhd, DmGgt-1, DmCyp4d14 细胞进程 Cellular process DmCG15528, DmVha68-3, DmCG9279, DmCG7848, DmCd, DmDhd, DmGgt-1 细胞组分
Cellular component细胞 Cell DmCG9279, DmCG17140, DmCG8219, DmGr59f, DmDhd, DmVha100-4 细胞器 Organelle DmCG9279, DmCG17140, DmDhd, DmVha100-4 细胞器组件 Organelle part DmCG9279, DmCG17140, DmVha100-4 细胞组件 Cell part DmCG9279, DmCG17140, DmCG8219, DmGr59f, DmDhd, DmVha100-4 细胞膜 Membrane DmVha68-3, DmNdl, DmCG34056, DmCG1246, DmCG17140,
DmTsp42Ep, DmCd, DmGr59f, DmCtr1B, DmVha100-4细胞膜组件 Membrane part DmVha68-3, DmNdl, DmCG1246, DmTsp42Ep, DmCd, DmGr59f, DmCtr1B, DmVha100-4 胞外区 Extracellular region DmMtg 大分子复合物 Macromolecular complex DmVha68-3, DmCG9279, DmVha100-4 分子功能
Molecular function电子转运活性 Electron carrier activity DmAOX2 转录活性 Transporter activity DmVha68-3, DmCtr1B, DmVha100-4 结合 Binding DmVha68-3, DmCG9279, DmCG7848, DmAOX2, DmMtg, DmCd, DmCG8219, DmCyp4d14 催化活性 Catalytic activity DmCG15528, DmVha68-3, DmNdl, DmCG7848, DmAOX2, DmCG34056, DmCd, DmDhd, DmGgt-1, DmCG10514, DmCyp4d14 抗氧化活性 Antioxidant activity DmCd 结果显示，上调基因多集中在胞外区，下调基因多集中在代谢进程、细胞膜和催化活性，说明被英诺克李斯特菌侵染后果蝇S2细胞的某些代谢反应下降。
2.5 果蝇S2细胞差异表达基因的富集分析
2.5.1 上调基因KEGG富集
KEGG富集分析表明，显著上调的21个基因只注释到1个KEGG 通路，即Toll/Imd通路(P<0.001)，包括4个果蝇抗菌肽基因DmDpt A、DmDef、DmCec A2(Cecropin A2)和DmCec B(Cecropin B)。
2.5.2 下调基因KEGG富集
KEGG富集分析结果(表5)表明，下调的33个基因注释到23个KEGG 通路。
表 5 英诺克李斯特菌侵染后果蝇S2细胞下调基因KEGG富集表Table 5. KEGG enrichment of downregulated genes in Drosophila S2 cells after Listeria innocua infection通路 Pathway P 所含基因 Contained gene 霍乱弧菌感染 Vibrio cholerae infection 0.033 DmVha68-3, DmVha100-4 类风湿关节炎 Rheumatoid arthritis 0.034 DmVha68-3, DmVha100-4 幽门螺杆菌感染中的上皮细胞信号传导
Epithelial cell signaling in Helicobacter pylori infection0.036 DmVha68-3, DmVha100-4 集合管酸分泌 Collecting duct acid secretion 0.037 DmVha68-3, DmVha100-4 突触小泡循环 Synaptic vesicle cycle 0.042 DmVha68-3, DmVha100-4 咖啡因代谢 Caffeine metabolism 0.043 DmAOX2 牛磺酸和牛磺酸的代谢 Taurine and hypotaurine metabolism 0.049 DmGgt-1 吞噬体 Phagosome 0.049 DmVha68-3, DmVha100-4 花生四烯酸代谢 Arachidonic acid metabolism 0.072 DmGgt-1 氧化磷酸化 Oxidative phosphorylation 0.078 DmVha68-3, DmVha100-4 矿物质吸收 Mineral absorption 0.082 DmCtr1B 粘蛋白型 O−聚糖的生物合成 Mucin type O-glycan biosynthesis 0.089 DmCG34056 药物代谢−其他酶 Drug metabolism - other enzymes 0.155 DmAOX2 NOD样受体信号通路 NOD-like receptor signaling pathway 0.158 DmDhd 丙酮酸代谢 Pyruvate metabolism 0.163 DmCG7848 铂耐药 Platinum drug resistance 0.182 DmCtr1B 过氧化物酶体 Peroxisome 0.197 DmAOX2 结核 Tuberculosis 0.198 DmVha100-4 谷胱甘肽代谢 Glutathione metabolism 0.204 DmGgt-1 流体剪切应力与动脉粥样硬化 Fluid shear stress and atherosclerosis 0.209 DmDhd mTOR信号通路 mTOR signaling pathway 0.216 DmVha68-3 溶酶体 Lysosome 0.247 DmVha100-4 嘌呤代谢 Purine metabolism 0.262 DmAOX2 2.6 免疫相关通路中的基因分析
英诺克李斯特菌侵染果蝇S2细胞后的转录组测序数据经GO注释和KEGG富集分析显示，显著上调的基因主要与免疫相关。从24个KEGG通路中筛选出10个与免疫相关的KEGG 通路，分别是Toll/Imd通路、类风湿关节炎、吞噬体、霍乱弧菌感染、幽门螺杆菌感染中的上皮细胞信号传导、NOD样受体信号通路、铂耐药、结核、流体剪切应力与动脉粥样硬化、mTOR信号通路，共涉及到8个基因，其中上调幅度最大的是DmDef，上调了9.805倍，下调幅度最大的是DmDhd，下调了4.329倍。
2.7 果蝇S2细胞转录组中抗菌肽的qPCR验证
英诺克李斯特菌侵染果蝇S2细胞后，变化最为显著的基因类群为抗菌肽。转录组分析显示，英诺克李斯特菌侵染果蝇S2细胞3 h后共10个抗菌肽基因显著上调，包括DmDef、DmDro、DmDpt A、DmDpt B、DmCec A2、DmCec B、DmMtk(DmMetchnikowin)、DmAtt D(DmAttacin D)、DmAtt A(DmAttacin A)、DmAtt B(DmAttacin B)，其中上调幅度最大的是DmDef(表2)。
我们选取9个抗菌肽基因，利用qPCR进行了验证，其中上调幅度最大的是DmMtk，上调8.180倍，DmAtt A 上调7.533倍，DmDro 上调7.204倍， DmDef 上调4.569倍(图2)，与转录组测序结果的趋势基本一致，说明英诺克李斯特菌侵染果蝇S2细胞后，抗菌肽基因受到显著的诱导表达。
![]()
图 2 免疫相关基因表达的qPCR验证结果Figure 2. qPCR verification results for expression of immune related genes3. 讨论与结论
已知的果蝇免疫信号通路有3条，分别是Toll通路、Imd通路和JAK/STAT通路。果蝇虽然没有获得性免疫，但是依然能对不同病原菌的入侵做出不同的先天免疫反应，革兰阳性菌和真菌所引起的免疫反应主要由Toll通路介导，而革兰阴性菌和细胞壁含DAP型肽聚糖的阳性菌引起的免疫反应主要由Imd通路介导[28-29]。Toll通路或Imd通路被激活后可诱导下游抗菌肽的表达[30]。已知的果蝇抗菌肽主要有7种，分别是Attacin、Diptericin、Drosocin、Defensin、Cecropin、Metchnikowin和Drosomycin。抗菌肽的表达由Toll和Imd通路调控，若Toll通路和Imd通路缺失，将不能诱导抗菌肽的表达[31-33]。
英诺克李斯特菌属于革兰阳性菌，其细胞壁肽聚糖成分主要为DAP型肽聚糖[34]，DAP型肽聚糖可以激活果蝇的Imd通路。已有研究表明，果蝇的PGRP-SD作为胞外受体不参与Toll通路，而是参与Imd通路的激活[35]。本研究发现DmPGRP-SD出现显著上调，提示英诺克李斯特菌侵染果蝇S2细胞后Imd通路可能被激活。抗菌肽Attacin、Cecropin和Diptericin的表达主要受Imd通路的调控[36-37]，本研究结果显示果蝇的抗菌肽基因DmAtt A、DmAtt B、DmAtt D 和 DmCec A2 和 DmCec B均出现显著上调，显示了Imd通路的激活。而抗菌肽Drosomycin、Defensin和Metchnikowin的表达主要受Toll通路调控[38]，本研究结果中显示果蝇的抗菌肽基因DmDro、DmDef和DmMtk也出现了显著上调，提示英诺克李斯特菌可能同时引起了Imd和Toll通路的激活，但具体是Toll通路还是Imd通路起主要作用，则需要进一步的研究。抗菌肽被认为具有广谱性，但Hanson等[39]研究发现抗菌肽针对某些病原微生物具有特异性，Imd通路调控的抗菌肽对单增李斯特菌并不起作用[34]。Hanson等[39]研究还发现果蝇的spz基因缺失后对英诺克李斯特菌表现出极高的敏感性，提示在抵御英诺克李斯特菌的过程中，Toll通路调控的抗菌肽起重要作用。鉴于上述研究结果和单增李斯特菌与英诺克李斯特菌相近的亲缘关系，猜测宿主在抵御英诺克李斯特菌和单增李斯特菌侵染的过程中Toll通路调控的抗菌肽可能起到关键作用，但Toll通路的激活机制需要进一步深入研究。
英诺克李斯特菌作为一种胞内寄生菌，通过内吞作用侵入宿主细胞后会被包裹在吞噬体中，英诺克李斯特菌在宿主细胞中的增殖需要它从宿主的吞噬体中逃逸至细胞质中，我们注意到液泡型H+离子腺苷三磷酸酶(Vacuolar H+-ATPase，V-ATPase)基因DmVha68-3和DmVha100-4均出现了显著下调，KEGG富集结果显示它们富集到吞噬体等通路中。DmVha68-3编码果蝇V-ATPase-A亚基，DmVha100-4编码V-ATPase-a亚基，为V-ATPase复合体形成的关键亚基[40]，V-ATPase广泛存在于细胞器膜上，如内体、溶酶体及吞噬体膜，可以将H+泵入这些细胞器中以维持这些细胞器的相对酸性内环境，促进细胞器对内吞或是细胞自噬包裹的底物或是入侵的微生物进行降解[41-42]。单增李斯特菌的毒力因子LLO正常发挥作用时需要酸性环境，在pH 5.6时LLO的溶血活性最高，当pH超过7.0时LLO的活性基本丧失[43]。典型的英诺李斯特菌缺失毒力岛，而其侵染果蝇S2细胞后上述2个DmV-ATPase基因的下调是否跟毒力岛缺失有关目前尚不清楚，值得我们今后深入思考。此外，是否可以通过调控宿主细胞的V-ATPase活性从而抑制单增李斯特菌毒力因子LLO的活性对其进行防控，也是有益的思考。
-
![]()
图 3 不同棉花冠层部位的雾滴覆盖率
各图中柱子上方不同小写字母表示差异显著(P<0.05, Duncan’s 法)
Figure 3. Droplet coverage rates of different canopy layers in four cotton cultivars
Different lowercase letters on the columns in the same figure represent significant difference (P<0.05, Duncan’s test)
![]()
图 4 4种棉花品种冠层雾滴总覆盖率
Figure 4. Total droplet coverage rates for canopies of four cotton cultivars
![]()
图 5 喷药后不同调查时间4种棉花品种的脱叶率
“*”表示不同处理间差异显著(P<0.05，t检验)
Figure 5. Defoliation rates of four cotton cultivars at different investigation time after application
“*” represent different significance between different treatment (P<0.05，t test)
![]()
图 6 喷药后不同调查时间4种棉花品种的吐絮率
Figure 6. Boll opening rate of four cotton cultivars at different investigation time after application
表 1 供试棉花品种的基本信息
Table 1 Basic information of tested cotton cultivars
棉花品种
Cultivar播种日期
Sowing
date株高/cm
Plant
height药前吐絮率/%
Boll opening rate before treatment金棉18#
Jinmian 18#2018−04−29 81 33.2 北泉9#
Beiquan 9#2018−04−29 72 35.4 植金13#
Zhijin 13#2018−05−02 70 34.7 合信47#
Hexin 47#2018−05−02 68 31.5 
下载: 导出CSV
表 2 不同施药处理下4种棉花的纤维品质和产量构成因子1)
Table 2 The fiber quality and lint cotton yield components of four cotton cultivars under different spraying treatments
品种
Cultivar施药处理
Spraying treatment纤维长度/mm
Fiber length整齐度指数/%
Fiber uniformity断裂比强度/
(cN·tex−1)
Fiber tenacity马克隆值
Micronaire伸长率/%
Elongation单铃质量/g
Boll weight衣分/%
Lint percentage金棉18#
Jinmian 18#无人飞机
UAV27.93±0.64a 84.10±1.05a 27.20±0.46a 5.33±0.35a 6.77±0.06a 5.41±0.08a 47.52±0.49a 人工喷雾
Backpack sprayer28.27±0.90a 85.10±1.10a 27.28±0.43a 5.40±0.26a 6.78±0.03a 5.37±0.07a 48.38±1.10a 北泉9#
Beiquan 9#无人飞机
UAV29.52±0.17a 86.14±0.50a 32.23±0.48a 4.66±0.05a 6.81±0.04a 5.88±0.49a 42.90±0.73a 人工喷雾
Backpack sprayer29.53±0.21a 85.63±0.78a 31.70±0.70a 4.57±0.23a 6.80±0.00a 5.88±0.11a 43.33±1.52a 植金13#
Zhijin 13#无人飞机
UAV27.82±0.99a 85.17±0.76a 26.93±1.37a 4.74±0.13a 6.72±0.03a 5.61±0.10a 47.58±0.62a 人工喷雾
Backpack sprayer27.80±0.95a 85.13±0.81a 26.70±1.28a 4.70±0.10a 6.73±0.06a 5.59±0.08a 47.21±0.52a 合信47#
Hexin 47#无人飞机
UAV28.58±0.77a 83.97±1.05a 28.38±1.46a 5.03±0.39a 6.73±0.06a 5.32±0.34a 46.80±1.33a 人工喷雾
Backpack sprayer28.28±0.65a 83.33±2.02a 28.52±1.26a 4.82±0.59a 6.75±0.05a 5.16±0.47a 45.77±1.24a 1)表中数据为平均值±标准差；同一品种不同施药处理不同的小写字母表示差异显著(P < 0.05， t检验)
1)Data in the table are mean±standard deviation；Different lowercase letters of the same cultivar with different sprayer treatment indicate significant difference (P < 0.05， t test)
下载: 导出CSV
-
[1] ZHANG T Z, QIAN N, ZHU X F, et al. Variations and transmission of QTL alleles for yield and fiber qualities in upland cotton cultivars developed in China[J]. PLoS One, 2013, 8(2): e57220. doi: https://doi.org/10.1371/journal.pone.0057220.
[2] DEGUINE J P, FERRON P, RUSSELL D. Sustainable pest management for cotton production: A review[J]. Agronomy for Sustainable Development, 2008, 28(1): 113-137. doi: 10.1051/agro:2007042
[3] 路战远, 咸丰, 张建中, 等. 棉花机械化采收技术技术规程[J]. 现代农业科技, 2017(11): 57. doi: 10.3969/j.issn.1007-5739.2017.11.033 [4] 李新裕, 陈玉娟, 闫志顺. 棉花脱叶技术研究[J]. 中国棉花, 2000, 27(7): 14-15. doi: 10.3969/j.issn.1000-632X.2000.07.006 [5] 秦维彩, 薛新宇, 崔龙飞, 等. 棉花脱叶剂施药机喷雾参数优化与试验[J]. 中国农机化学报, 2017, 38(4): 25-32. [6] 耿涛, 戴路, 许占伟. 棉花高密度种植群体结构的研究[J]. 新疆农业科学, 2003, 40(5): 269-272. doi: 10.3969/j.issn.1001-4330.2003.05.004 [7] 王刚, 刘辉, 赵海, 等. 新疆兵团棉花机械采收存在的问题及对策[J]. 中国棉花, 2011, 38(9): 37-38. doi: 10.3969/j.issn.1000-632X.2011.09.014 [8] 蒙艳华, 兰玉彬, 梁自静, 等. 无人飞机施药液量对棉花脱叶效果的影响[J]. 中国棉花, 2019, 46(6): 10-15. doi: 10.11963/1000-632X.myhlyb.20190611 [9] 娄尚易, 薛新宇, 顾伟, 等. 农用植保无人飞机的研究现状及趋势[J]. 农机化研究, 2017, 39(12): 1-6. [10] MENG Y H, LAN Y B, MEI G Y, et al. Effect of aerial spray adjuvant applying on the efficiency of small unmanned aerial vehicle on wheat aphids control[J]. International Journal of Agricultural and Biological Engineering, 2018, 11(5): 46-53. doi: 10.25165/j.ijabe.20181105.4298
[11] QIN W C, QIU B J, XUE X Y, et al. Droplet deposition and control effect of insecticides sprayed with an unmanned aerial vehicle against plant hoppers[J]. Crop Protection, 2016, 85: 79-88. doi: 10.1016/j.cropro.2016.03.018
[12] MENG Y H, SONG J L, LAN Y B, et al. Harvest aids efficacy applied by unmanned aerial vehicles on cotton crop[J]. Industrial Crops and Products, 2019, 140. doi: https://doi.org/10.1016/j.indcrop.2019.111645.
[13] XIN F, ZHAO J, ZHOU Y T, et al. Effects of dosage and spraying volume on cotton defoliants efficacy: A case study based on application of unmanned aerial vehicles[J]. Agronomy, 2018, 8(6). doi: https://doi.org/10.3390/agronomy8060085.
[14] MENG Y H, HAN Y X, LIANG Z J, et al. Harvest-aid application strategy in different cotton planting densities by unmanned aerial vehicle[J]. International Journal of Precision Agricultural Aviation, 2019, 2(1): 31-41. doi: 10.33440/j.ijpaa.20190201.0027
[15] 胡红岩, 任相亮, 姜伟丽, 等. 植保无人飞机飞行高度与方向对棉田雾滴沉积分布的影响[J]. 中国棉花, 2017, 44(12): 16-19. [16] 胡红岩, 任相亮, 马小艳, 等. 无人飞机喷施与人工喷施棉花脱叶剂效果对比[J]. 中国棉花, 2018, 45(7): 13-15. [17] 中华人民共和国农业农村部农药鉴定所. 农药田间药效试验准则(二)第134部分: 棉花生长调节剂试验, GB/T 17980.134—2004[S]. 北京: 中国标准出版社, 2004. [18] DU M W, LI Y, TIAN X L, et al. The phytotoxin coronatine induces abscission-related gene expression and boll ripening during defoliation of cotton[J]. PLoS One, 9(5): e97652. doi: https://doi.org/10.1371/journal.pone.0097652.
[19] DU M W, REN X M, TIAN X L, et al. Evaluation of harvest aid chemicals for the cotton-winter wheat double cropping system[J]. Journal of Integrative Agriculture, 2013, 12(2): 273-282. doi: 10.1016/S2095-3119(13)60226-9
[20] 辛明华, 李小飞, 韩迎春, 等. 新疆南疆不同棉花品种的生长特征及其品质分析[J]. 新疆农业科学, 2020, 57(4): 729-739. [21] 何勇, 吴剑坚, 方慧, 等. 植保无人飞机雾滴沉积效果研究综述[J]. 浙江大学学报(农业与生命科学版), 2018, 44(4): 392-398. [22] MENG Y H, SU J, SONG J L, et al. Experimental evaluation of UAV spraying for peach trees of different shapes: Effects of operational parameters on droplet distribution[J]. Computers and Electronics in Agriculture, 2020, 170. doi: https://doi.org/10.1016/j.compag.2020.105282.
[23] 张金龙, 董合林, 陈国栋, 等. 不同熟性棉花品种棉铃空间分布及产量品质形成的差异[J]. 西北农业学报, 2017, 2: 234-241. doi: 10.7606/j.issn.1004-1389.2017.02.011 [24] 张强, 赵冰梅, 朱玉永, 等. 不同植保施药机械喷施棉花脱叶剂效果评价[J]. 新疆农垦科技, 2020, 7: 29-32. doi: 10.3969/j.issn.1001-361X.2020.10.023 [25] 陈秀玲, 陈振, 张玮涛, 等. 北疆地区不同棉花品种机采性能差异研究与综合评价[J/OL]. 南京农业大学学报, http://kns.cnki.net/kcms/detail/32.1148.s.20200327.1600.004.html. -
期刊类型引用(2)
1. 邓健浩，刘力，舒婵，吴坤，吴坤钟，姚兵华，李康，田铃，黄志君. 单增李斯特菌、格氏李斯特菌或威尔斯李斯特菌侵染果蝇S2细胞后的转录分析. 华南农业大学学报. 2025(01): 12-24 . 
本站查看
2. 徐佳卉，刘华伟，孙小桐，唐章晨，黄敏，赵萍. 家蚕核转录因子NF-κB对不同细胞系中免疫相关基因表达的调控分析. 生物工程学报. 2024(10): 3649-3665 . 百度学术
其他类型引用(0)
计量
- 文章访问数:
- HTML全文浏览量: 0
- PDF下载量:
- 被引次数: 2
