Classification and feature band extraction of diseased citrus plants based on UAV hyperspectral remote sensing
-
摘要:目的
结合传统与现代农业病虫害监测的优缺点，探索通过无人机高光谱遥感技术检测出患病的柑橘植株、通过人工田间调查方式判断其患病种类及患病程度的病虫害监测方法。
方法使用无人机获取原始高光谱图像，经过光谱预处理和特征工程后，采用连续投影算法提取对柑橘患病植株分类贡献值最大的特征波长组合，基于全波段使用BP神经网络和XgBoost算法、基于特征波段使用逻辑回归和支持向量机算法，建立分类模型。
结果基于全波段的BP神经网络和XgBoost算法的ROC曲线下面积(Area under curve，AUC)分别为0.883 0和0.912 0，分类准确率均超过95%；提取出698和762 nm的特征波长组合，基于特征波长使用逻辑回归和支持向量机算法建立的分类模型召回率分别达到了93.00%和96.00%。
结论基于特征波长建模在患病样本分类中表现出很高的准确率，证明了特征波长组合的有效性。本研究结果可为柑橘种植园的病虫害监测提供一定的数据和理论支撑。
Abstract:ObjectiveCombined with the advantages and disadvantages of traditional and modern agricultural pest monitoring, the method of monitoring pest and disease were discussed, which detected the diseased citrus plants by UAV hyperspectral remote sensing technology and judged the disease species and disease degree by artificial field investigation.
MethodThe original hyperspectral images were obtained by UAV. After spectral preprocessing and feature engineering, continuous projection algorithm was used to extract the feature wavelength combination which contributed the most to the classification of citrus diseased plants. Finally, the BP neural network and XgBoost algorithm were used based on the full band, and the logistic regression and support vector machine algorithm were used to establish the classification model based on the characteristic band.
ResultThe AUC scores of BP neural network and XgBoost were 0.8830 and 0.9120 respectively, and the accuracy rates of both methods were over 95%. The feature wavelength combination of 698 and 762 nm was extracted. Based on this characteristic band, the recall rates of logistic regression and support vector machine algorithm were 93.00% and 96.00% respectively.
ConclusionThe model based on characteristic band shows high accuracy in the classification of disease samples, which proves the effectiveness of characteristic wavelength combination. This result can provide some data and theoretical support for monitoring diseases and pests in citrus plantations.
-
中药大黄来源于蓼科植物大黄的根茎，最早记载于《神农百草经》，有药中“将军”之美称，被列为中药“四大金刚”之一，已经有
2000 多年的中医临床用药史[1]。中药大黄主要由大黄酸、大黄素、大黄酚、芦荟大黄素等游离蒽醌，番泻苷A/B等二蒽酮苷类结合型蒽醌，二苯乙烯，多糖，单宁等成分组成，具有泻下攻积、抗菌、抗感染、清热泻火、凉血解毒、逐瘀通经、保肝利胆等多种药理作用[2]。其中，蒽醌是大黄的主要活性成分，含量(w)约占其成分的3%~5%[3]。大黄酸(Rhein，RH)是大黄蒽醌的主要成分之一，结构式如图1所示，具有抗感染[4]、抗肿瘤[5]、抗糖尿病[6]、降脂[7]、抗氧化、抗菌、抗病毒、抗肾毒性等[8]多种药理作用。然而，大黄酸溶解度低、脂溶性差、毒性强、生物利用度差、易引起胃肠道不适，这些特性限制了其临床应用[9]。近年来，畜禽病原菌的抗生素耐药性逐渐变得严重，导致抗生素敏感性下降，研发抗生素增效剂、复方制剂等的防控策略逐渐成为应对畜禽耐药病原菌的研究热点。大黄酸对幽门螺杆菌、金黄色葡萄球菌等多种厌氧菌有强大的抑菌作用。对大黄酸进行结构改造，通过羟基、羧基取代成酯或酰胺等修饰，改善大黄酸理化性质，降低毒性，提高抗感染、抗癌、抗菌等活性，具有重要意义，可为畜禽病原菌耐药防控提供新思路，也为大黄酸类衍生物的进一步研究提供参考依据。目前，大黄酸衍生物双醋瑞因是大黄酸最具代表性的衍生物之一，可阻滞白细胞介素−1β(Interleukin-1β，IL-1β)的作用，减少破骨细胞的形成和抑制吸收因子的合成在临床上广泛应用于治疗骨关节炎[10]。双醋瑞因的上市为后续大黄酸衍生物的研究和开发提供了理论依据；但是，双醋瑞因存在严重的毒副作用，为了控制其重度腹泻及影响肝脏的风险，对其临床使用设置了限制性条件。双醋瑞因又名1, 8−二乙酰基−3−羧基蒽醌，用乙酰氯对大黄酸第1、8位羟基进行保护。有研究表明，第1、8位羟基对大黄酸毒性的影响较小[11]，因此，双醋瑞因对大黄酸毒性的改善并不明显。研究表明，大黄酸第3位羧基是大黄酸毒性的关键影响因素[11]，对第3位羧基进行修饰能提高大黄酸的跨膜能力，从而可能提高大黄酸的抗感染效果[12]。据报道，游离的脂肪酸在肠道炎症反应中具有潜在的应用价值，它可以与肠道细胞膜的特异性受体结合从而发挥免疫保护作用，降低肠道促炎因子水平，以治疗肠道炎症[13]。一方面，脂肪酸作为能量物质，提供了成年人每日所需能量的30%[14]；另一方面，脂肪酸参与生物膜的构成、受体激活、信号转导、神经系统发育和各种基因的调节[15]。脂肪酸衍生化可以使大黄酸具有更好的应用价值，保护大黄酸，避免其在体内快速代谢或降解，提高大黄酸稳定性[16]；还可以增强大黄酸分子的亲脂性和跨膜转运能力[17]；由于肿瘤组织对脂肪酸的吸收能力较强，也可以促进大黄酸分子对癌细胞的靶向积累。总之，通过脂肪酸对大黄酸进行衍生化具有广阔的应用前景。
本研究前期工作发现，大黄酸与脂肪酸具有协同抗感染作用，并且脂肪酸能在一定程度上降低大黄酸的毒性；但是，大黄酸和脂肪酸在体内代谢速率、生物利用度等因素影响下很难发挥协同作用。因此，本研究通过对大黄酸第3位羧基与羟基癸酸和癸醇进行酯化衍生化，制备出4种大黄酸衍生物，通过核磁共振氢谱(1H nuclear magnetic resonance，1H NMR)和傅里叶变换红外光谱(Fourier transforme infrared spectroscopy，FTIR)表征其结构，进一步探究其毒性、抗感染和抗氧化活性变化，筛选出具有应用前景的大黄酸衍生物，以期为大黄酸衍生物的结构设计及其生物活性等研究提供参考依据。
1. 材料与方法
1.1 仪器与材料
核磁共振波谱(德国Bruker公司)；磁力搅拌器(上海龙跃仪器设备有限公司)；流式细胞计数仪(美国贝克曼库尔特)；傅里叶变换红外光谱仪(美国Thermo Fisher Nicolet公司)；酶标定量测定仪(美国Beckman公司)；CO2细胞培养箱(新加坡艺思高科技有限公司)；RT-qPCR仪(美国Bio-Rad生命医学产品有限公司)；旋转蒸发仪N-1200B(日本EYELA东京理化)；小鼠单核巨噬细胞RAW 264.7(武汉普诺赛生命科技有限公司)；一氧化氮(NO)检测试剂盒(上海碧云天生物技术股份有限公司)；DMEM、特级胎牛血清、青链霉素(赛默飞生命科学产品与服务旗舰店)；总RNA极速提取试剂盒(广州晖鼎生物科技有限公司)；4−(二甲氨基)吡啶(4-Dimethylaminopyridine，DMAP)、N,N′−二环己基碳二亚胺(N, N′-Dicyclohexylcarbodiimide，DCC)(上海麦克林生化科技有限公司)；氘代二甲基亚砜(Dimethyl sulfoxide，DMSO)(上海源叶生物科技有限公司)；大黄酸(上海麦克林生化科技有限公司)；1−癸醇(1-Decanol，DA)、10−羟基癸酸(10-Hydroxydecanoic acid，10-HA)、10−羟基−2−癸烯酸(10-Hydroxy-2-decenoic acid，10-HDA)、5−羟基癸酸钠(5-Hydroxydecanoate sodium，5-HD)(上海毕得医药科技股份有限公司)结构式如图2所示；双氯荧光黄乙酸乙酯(2, 7-Dichlorofluorescin diacetate，DCFH-DA)(广州市齐云生物技术有限公司)。
![]()
图 2 癸醇(A)和羟基癸酸(B、C、D)的结构式Figure 2. Structural formulas of decanol (A) and hydroxydecanoic acid (B, C, D)1.2 试验与方法
1.2.1 大黄酸衍生物的制备与纯化
精密称取0.5 g(1.70 mmol)大黄酸，置于20 mL DMSO溶液中，于50 ℃、200 r/min搅拌1 h，完全溶解后，加入DA/10-HA/10-HDA/5-HD、DCC 3.40 mmol和DMAP 0.85 mmol，反应摩尔比按照大黄酸∶DA/10-HA/10-HDA/5-HD∶DCC∶DMAP=1.0∶1.0∶2.0∶0.5。在50 ℃条件下反应，每间隔1 h，用薄层层析法监测反应情况，确认产物生成，直至产物不再增加(展开剂为乙酸乙酯和石油醚，体积比为2∶3)，停止反应。大黄酸与羟基脂肪酸的反应结构式如图3所示。反应结束后待温度降至室温，将反应混合物用饱和碳酸氢钠溶液透析1 d (透析袋截留相对分子质量为500)，去离子水透析2 d，每隔12 h更换透析液，除去未反应的物质。随后，在−75 ℃、0.001 MPa条件下冷冻干燥48 h，得到冻干粗产物。将上述得到的粗产品用38~48 μm孔径的硅胶柱层析进行干法上样，使样品层与硅胶层之间的高度比为1∶20。洗脱液乙酸乙酯∶石油醚=2∶3(V/V)，直至监测到产物完全洗脱，获得最终产物。用旋转蒸发器除去溶剂，于50 ℃真空烘箱干燥24 h，得到黄色粉末固体产物。
![]()
图 3 大黄酸与羟基脂肪酸的反应结构式F-A为接枝脂肪酸或癸醇，R为F-A反应掉羟基的部分结构，R-FA为接枝后的大黄酸衍生物。Figure 3. Structural formula for the reaction of rhein with hydroxy fatty acidsF-A is grafted fatty acid or decanol, R is the part of the structure where F-A reacts to drop the hydroxyl group, R-FA is the grafted rhein acid derivative.1.2.2 大黄酸衍生物的抗感染活性测定
为验证大黄酸及其衍生物对RAW 264.7细胞存活率的影响，选用CCK-8试剂盒对RAW 264.7细胞存活率进行检测。分别用10~160 μmol/L的大黄酸及其衍生物处理RAW 264.7细胞24 h，10%(
$\varphi $ ) CCK-8试剂孵育0.5~4.0 h，酶标仪设定波长为450 nm，测定每孔光密度。为验证大黄酸衍生物的抗感染效果，构建脂多糖(Lipopolysaccharide，LPS)诱导的RAW 264.7细胞模型，测定大黄酸及其衍生物处理后RAW 264.7中NO释放量以及促炎因子TNF-α、IL-1β、IL-6表达量变化。用20 μmol/L的大黄酸、R-5HD、R-DA、R-10HA、R-10HDA在37 ℃、5%(
$\varphi $ ) CO2条件下预孵育8 h，加入LPS，终质量浓度为 1 μg/mL，在37 ℃、5%($\varphi $ ) CO2条件下孵育16 h(另外设置LPS阳性对照组，不加药物处理)。用细胞培养液上清液[DMEM+10% (φ) FBS溶液+1% (φ) 青链霉素溶液]稀释NaNO2标准品至0、1、2、5、10、20、40、60、100 μmol/L。按50 μL每孔，在96孔板中加入标准品及样品(每组设置3个复孔)，用酶标仪于540 nm波长处测定光密度，计算NO含量。提取RAW 264.7的RNA，反转录，利用RT-qPCR测定肿瘤坏死因子−α (Tumour necrosis factor-α，TNF-α)、IL-1β、白细胞介素−6 (Interleukin-6，IL-6)的表达量变化，具体步骤如下。首先，投入总量为1 μg的RNA，2 μL 4× gDNA wiper Mix，用RNase-free ddH2O定容至10 μL，充分混匀，于42 ℃离心2 min后，得到去DNA的反应液。然后，在上述反应液中加入4 μL 5× Hiscript Ⅲ qRT SuperMix，用RNase-free ddH2O定容至20 μL，得到反转录体系，充分混匀，短暂离心，经37 ℃ 15 min、85 ℃ 5 s的PCR反应后，得到cDNA。最后，按照下述反应程序进行RT-qPCR：95 ℃ 30 s，循环1次；95 ℃ 10 s，循环40次；60 ℃ 30 s，循环40次；72 ℃ 30 s，循环40次。RT-qPCR反应体系中所需的引物序列见表1。表 1 RT-qPCR引物Table 1. RT-qPCR primer引物名称
Primer name引物序列(5′→3′)
Primer sequence扩增序列
Amplified sequenceGAPDH-F ACCCAGAAGACTGTGGATGG GAPDH mRNA GAPDH-R CACATTGGGGGTAGGAACAC GAPDH mRNA TNF-ɑ-F CTCTTCAAGGGACAAGGCTG TNF-ɑ mRNA TNF-ɑ-R CGGACTCCGCAAAGTCTAAG TNF-ɑ mRNA IL-6-F CCGGAGAGGAGACTTCACAG IL-6 mRNA IL-6-R TCCACGATTTCCCAGAGAAC IL-6 mRNA IL-1β-F GACCTTCCAGGATGAGGACA IL-1β mRNA IL-1β-R AGGCCACAGGTATTTTGTCG IL-1β mRNA 1.2.3 大黄酸衍生物的抗氧化活性测定
为验证大黄酸及其衍生物对LPS诱导的RAW 264.7细胞中活性氧(Reactive oxygen species，ROS)产生的影响，用20 μmol/L的大黄酸以及大黄酸衍生物处理RAW 264.7细胞。ROS荧光探针DCFH-DA可穿透活细胞膜进入细胞内，并可被细胞内的ROS氧化，生成氧化乙啶；氧化乙啶可嵌入染色体DNA中，产生红色荧光。根据此原理可以评估细胞内ROS含量及其变化情况。用DCFH-DA染色，通过流式细胞计数仪检测大黄酸及其衍生物对RAW 264.7细胞中ROS产生的影响，评估大黄酸及其衍生物抗氧化水平。
1.2.4 统计学分析
采用Orgin Pro 9.0和Graph Pad 8进行统计分析，试验分组均至少有3次重复，采用t检验比较组别之间的差异性。若P>0.05，则无统计学意义；若P<0.05、P<0.01或P<0.001，则其有显著差异；若P<0.000 1，则其有极显著差异。
2. 结果与分析
2.1 大黄酸衍生物结构表征
为了表征大黄酸衍生物是否成功合成，用FTIR和1H NMR表征其结构。R-5HD的FTIR表征结果如图4A所示。大黄酸的羟基(—OH)峰位于3 420.40 cm−1，羰基(—C=O)峰位于1 648.98 cm−1；5-HD的羟基峰位于3 410.34 cm−1，羰基峰位于1 703.44 cm−1；R-5HD的羟基峰位于3 447.58 cm−1，羰基峰位于1 740.69 cm−1，新生成的酯键为1 275.67 cm−1。大黄酸蒽醌母核上的氢在1H NMR上分别对应a1(7.40×10−6)、b1/c1(7.70×10−6)、d1(7.82×10−6)、e1(8.09×10−6)，酚羟基氢对应h1/f1(11.87×10−6)，第3位羧基(—COOH)氢对应g1(13.79×10−6)；5-HD碳链上的氢对应b2(2.33×10−6)、c2(1.54×10−6)、d2(1.38×10−6)、e2(3.10×10−6)、g2(1.38×10−6)、h2/i2(1.25×10−6)、j2(1.28×10−6)、k2(0.88×10−6)，羟基氢对应f2(3.41×10−6)，羧基氢对应a2(11.87×10−6)。由大黄酸衍生物R-5HD的1H NMR表征(图4B)可以看出，大黄酸第3位羧基峰和5-HD碳链上的羟基峰消失，与5-HD对比，R-5HD中峰h3(4.02×10−6)发生明显的偏移。以上结果均表明R-5HD成功合成。
R-DA的FTIR表征结果如图4C所示，大黄酸的羟基峰位于3 420.40 cm−1，羰基峰位于1 648.98 cm−1；R-DA的羟基峰位于3 447.58 cm−1，羰基峰位于1 740.69 cm−1，新生成的酯键为1 275.67 cm−1。大黄酸蒽醌母核上的氢在1H NMR上分别对应a1(7.40×10−6)、b1/c1(7.70×10−6)、d1(7.82×10−6)、e1(8.09×10−6)，酚羟基的氢对应h1/f1(11.87×10−6)，第3位羧基的氢对应g1(13.79×10−6)；DA碳链上的氢对应b2/d2/e2/f2(1.26×10−6)、g2(1.30×10−6)、h2(1.43×10−6)、i2(1.50×10−6)、j2(4.35×10−6)，羟基氢对应k2(3.23×10−6)。由R-DA的1H NMR表征(图4D)可以看出，大黄酸第3位羧基峰和DA碳链上的羟基峰消失，与DA对比，R-DA中羟基相邻烷烃链上氢对应的h3(4.33×10−6)发生明显偏移。以上结果均表明R-DA成功合成。
R-10HA FTIR结构表征如图4E所示。大黄酸的羟基峰位于3 420.40 cm−1，羰基峰位于1 648.98 cm−1；10−HA的羟基峰位于3 410.34 cm−1，羰基峰位于
1703.44 cm−1；R-10HA的羟基峰位于3 447.58 cm−1，羰基峰位于1 740.69 cm−1，新生成的酯键为1 275.67 cm−1。大黄酸蒽醌母核上的氢在1H NMR上分别对应a1(7.40×10−6)、b1/c1(7.70×10−6)、d1(7.82×10−6)、e1(8.09×10−6)，酚羟基的氢对应h1/f1(11.87×10−6)，第3位羧基的氢对应g1(13.79×10−6)；10-HA碳链上的氢对应c2(1.54×10−6)、d2(1.33×10−6)、e2/f2(1.26×10−6)、g2(1.30×10−6)、h2(1.43×10−6)、i2(1.58×10−6)、j2(3.62×10−6)，羟基的氢对应k2(4.7×10−6)，羧基的氢对应a2(11.87×10−6)。由R-10HA的1H NMR 表征(图4F)可以看出，大黄酸第3位羧基峰和10-HA碳链上的羟基峰消失，与10-HA对比，R-10HA中峰h3(4.02×10−6)发生明显的偏移。以上结果均表明R-10HA成功合成。R-10HDA FTIR表征结果如图4G所示，大黄酸的羟基峰位于3 420.40 cm−1，羰基峰位于1 648.98 cm−1；10-HDA的羟基峰位于3 410.34 cm−1，羰基峰位于1 703.44 cm−1；R-10HDA的羟基峰位于3 447.58 cm−1，羰基峰位于1 740.69 cm−1，新生成的酯键为1 275.67 cm−1。大黄酸蒽醌母核上的氢在1H NMR上分别对应a1(7.40×10−6)、b1/c1(7.70×10−6)、d1(7.82×10−6)、e1(8.09×10−6)，酚羟基的氢对应h1/f1(11.87×10−6)，第3位羧基的氢对应g1(13.79×10−6)；10-HDA碳链上的氢对应b2(6.02×10−6)、c2(7.11×10−6)、d2(2.16×10−6)、e2(1.29×10−6)、f2(1.33×10−6)、g2(1.26×10−6)、h2(1.43×10−6)、i2(1.58×10−6)，羟基的氢对应k2(4.70×10−6)，羧基的氢对应a2(12.05×10−6)。由R-10HDA的1H NMR表征(图4H)可以看出，大黄酸第3位羧基峰和10-HDA碳链上的羟基峰消失，与10-HDA对比，R-10HDA中峰h3(4.02×10−6)发生明显偏移。以上结果均表明R-10HDA成功合成。
2.2 大黄酸衍生物的抗感染活性评估
2.2.1 大黄酸衍生物对RAW 264.7细胞的毒性
结果如图5所示，大黄酸浓度高于40 μmol/L时具有明显的毒性，RAW 264.7细胞存活率降低至60%；DA浓度为160 μmol/L时，RAW 264.7细胞存活率降低至60%。而R-DA、R-10HA、R-10HDA、R-5HD浓度为160 μmol/L时，RAW 264.7细胞存活率仍接近100%。由此可见，大黄酸衍生物显著降低大黄酸的细胞毒性。
![]()
图 5 不同浓度大黄酸及其衍生物对RAW 264.7细胞存活率的影响“*”“**”“***”分别表示相同处理浓度大黄酸衍生物组与大黄酸组在0.05、0.01和0.001水平差异显著(t检验)。Figure 5. Effect of different concentrations of rhein and its derivatives on survival rates of RAW 264.7 cells“*” “**” “***” indicate significant differences between rhein derivative group and rhein group at 0.05, 0.01 and 0.001 levels respectively in the same treatment concentration (t test).2.2.2 大黄酸衍生物的抗感染活性
大黄酸及其衍生物对NO释放量的影响如图6A所示。与空白对照组(CK)相比，阳性对照组(LPS模型组)NO释放量极显著增加(P<
0.0001 )。与LPS模型组相比，大黄酸、R-5HD、R-DA、R-10HA、R-10HDA均能够极显著抑制LPS诱导的RAW 264.7中NO的释放量，说明4种大黄酸衍生物均具有较强的抗感染活性，对第3位羧基进行结构修饰不会丧失大黄酸的抗感染活性。![]()
图 6 大黄酸及其衍生物的抗感染作用“*”“**”“***”“****”分别表示与LPS组(阳性对照组)在0.05、0.01、0.001和0.000 1水平差异显著(t检验)。Figure 6. Anti-inflammatory effects of rhein and its derivatives“*” “**” “***” “****” indicate significant differences differed from LPS group (positive control group) at 0.05, 0.01, 0.001 and 0.000 1 levels respectively (t test).大黄酸及其衍生物均能显著降低LPS诱导的RAW 264.7中促炎因子TNF-α、IL-1β、IL-6的mRNA相对表达量，结果如图6B~6D所示。与空白对照组(CK)相比，阳性对照组(LPS组)的TNF-α mRNA表达水平增加了47倍；与LPS组相比，大黄酸、R-5HD、R-DA、R-10HA、R-10HDA均显著降低了TNF-α mRNA表达水平(P<0.05)。LPS组的IL-1β mRNA表达水平比CK增加了8 813倍；与LPS组相比，大黄酸、R-5HD、R-DA、R-10HA、R-10HDA均显著降低了IL-1β mRNA表达水平。LPS组的IL-6表达水平比CK增加了4 885倍；与LPS组相比，大黄酸、R-5HD、R-DA、R-10HA、R-10HDA均显著降低了IL-6 mRNA表达水平。以上结果表明，大黄酸、R-5HD、R-DA、R-10HA和R-10HDA均能显著抑制炎症因子IL-1β、IL-6和TNF-α mRNA的表达水平，对大黄酸进行结构修饰使其抗感染效果更强，且R-DA的抗感染效果最强。
2.2.3 大黄酸衍生物的抗氧化活性
如图7所示，DCFH-DA探针孵育后，空白对照组(CK)的荧光强度为1 319.66；经过LPS刺激后细胞内ROS水平显著升高，荧光强度为60 968.00，比空白对照组增加了46倍(P<
0.0001 )。药物处理后大黄酸、R-5HD、R-DA、R-10HA以及R-10HDA组的荧光强度均显著下降，且R-5HD、R-DA、R-10HA、R-10HDA的抑制效果均优于大黄酸。以上结果表明，修饰后的R-5HD、R-DA、R-10HA、R-10HDA显著增强了大黄酸的抗氧化活性。![]()
图 7 大黄酸及其衍生物对LPS诱导RAW 264.7细胞释放ROS的影响“**”“***”“****”分别表示与LPS组(阳性对照组)在0.01、0.001和0.000 1水平差异显著(t检验)。Figure 7. Effect of rhein and its derivatives on LPS-induced ROS release from RAW 264.7 cells“**” “***” “****” indicate significant differences differed from LPS group (positive control group) at 0.01, 0.001 and 0.000 1 levels respectively (t test).3. 结论与讨论
本研究以大黄酸为原料，以羟基癸酸和癸醇为修饰物，经一步酯化反应合成了4种大黄酸衍生物(R-5HD、R-DA、R-10HA、R-10HDA)，通过FTIR以及1H NMR对合成的化合物进行结构鉴定，证实4种衍生物合成成功，经对大黄酸及其衍生物进行毒性研究以及抗感染、抗氧化活性研究，主要结论及其分析如下。
通过大黄酸及其衍生物对RAW 264.7细胞存活率影响的分析，经过羟基癸酸和癸醇修饰之后的大黄酸衍生物对RAW 264.7细胞的毒性降低，安全性增加。羟基癸酸和癸醇与药物分子的化学偶联可能会引起它们在体内的药效学/药代动力学变化，从而降低它们的毒性。经羟基癸酸和癸醇修饰之后的大黄酸衍生物具有更强的抗感染效果，能显著降低LPS诱导的NO释放水平和促炎因子TNF-α、IL-1β、IL-6的转录水平。综合来讲，在体外炎症模型下，R-DA的抗感染效果最为突出。
ROS作为LPS诱导炎症的重要中介，可引起DNA、蛋白质和脂质的损伤[18]。此外，过量的ROS产生和积累可诱导细胞损伤和氧化应激[19]。本研究中，LPS刺激后，ROS水平显著升高，经大黄酸以及4种衍生物孵育后，LPS诱导的RAW 264.7细胞中ROS水平显著降低，R-DA抑制ROS效果与大黄酸有显著差异，相比其他3种衍生物有更好的ROS抑制效果。有研究发现，大黄酸可以通过抑制诱导型一氧化氮合酶的表达起到抗感染作用，经修饰后的大黄酸衍生物作用后，细胞NO释放水平显著降低[20]，这与Ge等[21]研究结果一致。
在大黄酸构效关系研究中，大黄酸分子上第3位碳上的取代基是羧基，增大了大黄酸的极性，是溶解性的关键基团，因此，本研究采用酯化反应的方法，将第3位碳上的羧基转化为酯键，以制备具有良好脂溶性和抗感染活性的大黄酸衍生物。据调查，对蒽醌类药物的修饰改性主要是通过在羟基和羧基引入侧链基团增加药物的吸收，引入脂肪酸可以克服母体药物的脂溶性、吸收率、生物利用度等限制。此外，多数脂肪酸具有抗菌、抗感染、抗氧化活性，例如Omega-3和Omega-6等，在人体内发挥至关重要的作用。通过引入脂肪酸可以增强大黄酸的抗感染活性，降低毒性，同时也为机体的正常运转和代谢提供必需的化学物质。前期工作中经对大黄酸衍生物的抗菌活性进行初步验证，发现其能显著增强大黄酸对金黄色葡萄球菌和肠球菌的抗菌活性，相关研究后续将继续展开，为畜禽耐药性防控提供思路，也为大黄酸进一步的研究提供基础数据。
在体外抗感染过程中，R-DA可能通过抑制NF-κB途径减少LPS刺激的巨噬细胞中促炎细胞因子的产生，也可能是由于R-DA中羟基癸醇的脂溶性作用，与细胞膜上的蛋白受体靶向结合，增加细胞膜的通透性，使药物进入细胞内发挥作用，从而起到抗感染、抗氧化的作用。
本研究中，虽然制备的大黄酸衍生物具有更低的毒性以及更优的抗感染效果，但大黄酸衍生物的纯化方法、合成产率需要进一步优化和提高。因此，在本文的基础上，仍有许多工作值得深入探究，并且要对衍生物的抗感染机制做深入性研究，如衍生物的抗感染机制是否与大黄酸的作用机制一致，以及通过结构修饰后衍生物的生物利用度是否有所增加，以期为新药的开发利用提供参考价值。
-
![]()
图 1 试验区域及样本标注
粉、红、蓝、黄、白圆圈标记分别代表1、2、3、4级和患病未定级的黄龙病植株，三角形标记为缺素植株；没有标记的植株为完全健康植株
Figure 1. Experiment area and sample marking
The pink, red, blue, yellow and white circle markers represented the plants with Huanglongbing disease of grades 1, 2, 3, 4 and indefinite respectively, and the triangular markers represented the plants lacking in nutrients; The plants without markers were the complete healthy plants.
表 1 BP神经网络结构
Table 1 The structure of BP neural network
层(类型)
Layer (Type)输入节点数
Number of input node输出节点数
Number of ouput node激活函数
Activation function输入层(全连接) Input layer(Full connect) 125 32 ReLu 隐含层1(全连接) Hidden layer 1(Full connect) 32 32 ReLu 隐含层2(全连接) Hidden layer 1(Full connect) 32 16 ReLu 输出层(全连接) Output layer(Full connect) 16 2 Sigmoid 
下载: 导出CSV
表 2 BP神经网络和XgBoost测试结果(混淆矩阵)
Table 2 The test results of BP neural network and XgBoost (Confusion matrix)
模型 Model 项目 Item 预测为患病 Predicted illness 预测为健康 Predicted healthy 总计 Total BP神经网络
BP neural network真实为患病 True illness 78 22 100 真实为健康 True healthy 7 493 500 总计 Total 85 515 600 XgBoost 真实为患病 True illness 85 15 100 真实为健康 True healthy 13 487 500 总计 Total 98 502 600
下载: 导出CSV
表 3 特征波长下LR和SVM分类结果(混淆矩阵)
Table 3 The results of LR and SVM test bases of feature bands(Confusion matrix)
模型 Model 项目 Item 预测为患病 Predicted illness 预测为健康 Predicted healthy 总计 Total LR 真实为患病 True illness 93 7 100 真实为健康 True healthy 44 456 500 总计 Total 137 463 600 SVM 真实为患病 True illness 96 4 100 真实为健康 True healthy 61 439 500 总计 Total 157 443 600
下载: 导出CSV
表 4 4种模型评估得分
Table 4 The evaluation scores of four models
评价指标
Evaluating
indicator全波段 Full-wave band 特征波段 Feature band BP Net XgBoost LR SVM ACC 0.951 7 0.953 3 0.915 0 0.891 6 Recall 0.780 0 0.850 0 0.930 0 0.960 0 Precision 0.917 6 0.867 3 0.678 8 0.611 5 F1-score 0.843 8 0.858 5 0.784 9 0.747 1 AUC 0.883 0 0.912 0 0.921 0 0.919 0
下载: 导出CSV
-
[1] DENG X, LAN Y, HONG T, et al. Citrus greening detection using visible spectrum imaging and C-SVC[J]. Comput Electron Agr, 2016, 130: 177-183. doi: 10.1016/j.compag.2016.09.005
[2] 陈波, 姚林建. 光谱检测技术在柑橘黄龙病诊断中的研究进展[J]. 赣南师范大学学报, 2018, 39(6): 69-72. [3] 兰玉彬, 邓小玲, 曾国亮. 无人机农业遥感在农作物病虫草害诊断应用研究进展[J]. 智慧农业, 2019, 1(2): 1-19. [4] 纪景纯, 赵原, 邹晓娟, 等. 无人机遥感在农田信息监测中的应用进展[J]. 土壤学报, 2019, 56(4): 1-13. [5] 黄文江, 张竞成, 师越, 等. 作物病虫害遥感监测与预测研究进展[J]. 南京信息工程大学学报(自然科学版), 2018, 10(1): 30-43. [6] LI X, LEE W S, LI M, et al. Spectral difference analysis and airborne imaging classification for citrus greening infected trees[J]. Comput Electron Agr, 2012, 83: 32-46. doi: 10.1016/j.compag.2012.01.010
[7] KUMAR A, LEE W S, EHSANI R J, et al. Citrus greening disease detection using aerial hyperspectral and multispectral imaging techniques[J]. J Appl Remote Sens, 2012, 6(1): 63542. doi: 10.1117/1.JRS.6.063542
[8] 兰玉彬, 朱梓豪, 邓小玲, 等. 基于无人机高光谱遥感的柑橘黄龙病植株的监测与分类[J]. 农业工程学报, 2019, 35(3): 92-100. doi: 10.11975/j.issn.1002-6819.2019.03.012 [9] LAN Y, HUANG Z, DENG X. Comparison of machine learning methods for citrus greening detection on UAV multispectral images[J]. Comput Electron Agr, 2020, 171: 10524. doi: 10.1016/j.compag.2020.105234
[10] 李修华, 李民赞, LEE W S, et al. 柑桔黄龙病的可见−近红外光谱特征[J]. 光谱学与光谱分析, 2014, 34(6): 1553-1559. doi: 10.3964/j.issn.1000-0593(2014)06-1553-07 [11] 尚方信, 郭浩, 李钢, 等. 基于One-class SVM的噪声图像分割方法[J]. 计算机应用, 2019, 39(3): 874-881. doi: 10.11772/j.issn.1001-9081.2018071494 [12] TAX D M J, DUIN R P W. Support vector data description[J]. Mach Learn, 2004(54): 45-66.
[13] CHAWLA N V, BOWYER K W, HALL L O. SMOTE: Synthetic minority over-sampling technique[J]. J Artif Intell Res, 2002(16): 321-357.
[14] 成忠, 张立庆, 刘赫扬, 等. 连续投影算法及其在小麦近红外光谱波长选择中的应用[J]. 光谱学与光谱分析, 2010, 30(4): 949-952. doi: 10.3964/j.issn.1000-0593(2010)04-0949-04 [15] 高洪智, 卢启鹏, 丁海泉, 等. 基于连续投影算法的土壤总氮近红外特征波长的选取[J]. 光谱学与光谱分析, 2009, 29(11): 2951-2954. doi: 10.3964/j.issn.1000-0593(2009)11-2951-04 [16] 吴迪, 金春华, 何勇. 基于连续投影算法的光谱主成分组合优化方法研究[J]. 光谱学与光谱分析, 2009, 29(10): 2734-2737. doi: 10.3964/j.issn.1000-0593(2009)10-2734-04 [17] CHEN T, GUESTRIN C. XgBoost: A scalable tree boosting system[C]//The 22nd ACM SIGKDD International Conference. ACM, 2016: 785-794.
[18] 王术波, 韩宇, 陈建, 等. 基于深度学习的无人机遥感生态灌区杂草分类[J]. 排灌机械工程学报, 2018, 36(11): 1137-1141. [19] 高林, 杨贵军, 于海洋, 等. 基于无人机高光谱遥感的冬小麦叶面积指数反演[J]. 农业工程学报, 2016, 32(22): 113-120. doi: 10.11975/j.issn.1002-6819.2016.22.016
计量
- 文章访问数: 1130
- HTML全文浏览量: 4
- PDF下载量: 1555
