Effects of applying compound fertilizer on yield, nutrient absorption and oligosaccharide accumulation of Morinda officinalis
-
摘要:目的
研究复合肥施用量对巴戟天产量、养分吸收和主要有效成分寡糖累积量的影响，为巴戟天的施肥管理提供理论依据。
方法以巴戟天为研究对象，通过盆栽试验，设置每株施用复合肥0、2、4、6、8和10 g (CK、T1、T2、T3、T4和T5)共6个处理，测定巴戟天在不同施肥处理下的产量、养分含量和寡糖含量，分析各指标间的相关性，并运用主成分分析方法进行综合评价。
结果施用复合肥可以显著提高巴戟天产量和N、P、K含量，并有利于提高寡糖的累积量。其中，单株施肥量为6 g处理的巴戟天产量和寡糖总累积量最高，与CK相比，分别提升114.47%和75.36%。施肥在一定程度上降低了巴戟天的寡糖含量，与CK相比，降幅最高达35.65%。巴戟天对复合肥中养分的利用效率表现为N > K > P，其中T3处理N、P、K累积量均高于其他处理。影响巴戟天寡糖累积量的最关键因素是产量，巴戟天产量与全株N/P/K含量、4种寡糖累积量呈显著或极显著正相关，但耐斯糖、1F−果呋喃糖基耐斯糖含量与产量、全株N/P/K含量呈显著或极显著负相关。主成分分析的综合得分由高到低为T3 > T5 > T4 > T2 > T1 > CK，T3处理最佳。
结论建议巴戟天的种植以N肥为主，适量补充K和P肥，适当施肥可显著提高巴戟天药材产量和寡糖累积量。
Abstract:ObjectiveTo study the effects of compound fertilizer application amount on yield, nutrient absorption and accumulation of main effective component oligosaccharides of Morinda officinalis, and provide a theoretical basis for the fertilization management of M. officinalis.
MethodUsing M. officinalis as the research object, through the pot experiment, a total of six treatments were set up to apply 0, 2, 4, 6, 8 and 10 g compound fertilizer (CK, T1, T2, T3, T4 and T5) per plant. The yield, nutrient content and oligosaccharide content of M. officinalis under different fertilization treatments were determined. Then we analyzed the correlation between the indicators, and comprehensively evaluated them by principal component analysis.
ResultThe application of compound fertilizer significantly improved the yield of M. officinalis, increased the N, P and K contents, and contributed to the increase in the accumulation of oligosaccharides. The highest yield and total oligosaccharide accumulation were achieved in the treatment of 6 g compound fertilizer per plant (T3), compared with no fertilizer (CK), with increases of 114.47% and 75.36% respectively. Fertilizer application reduced the oligosaccharide content of M. officinalis to some extent, with a maximum decrease of 35.65% compared to CK. The utilization efficiency of nutrient in compound fertilizer by M. officinalis showed N > K > P, and the accumulations of N, P and K in T3 treatment were higher than those of other treatments. The yield was the most critical factor affecting the accumulation of oligosaccharides in M. officinalis. The yield of M. officinalis was significantly or highly significantly positively correlated with the whole plant N, P, K contents and the accumulation of four oligosaccharides, but the contents of nystose and 1F-fructofuranosyl nystose were significantly or highly significantly negatively correlated with the yield and the whole plant N, P, K contents. The comprehensive score of principal component analysis from high to low was T3>T5>T4>T2>T1>CK, with T3 treatment being the best.
ConclusionIt is suggested that the cultivation of M. officinalis should focus on N fertilizer, and supplement K and P fertilizer appropriately. Proper fertilization can significantly increase the yield of M. officinalis and the accumulation of oligosaccharides.
-
Keywords:
- Morinda officinalis /
- Compound fertilizer /
- Yield /
- Oligosaccharide /
- Nutrient absorption
-
中国是世界上水禽养殖第一大国，肉鸭产量占世界总产量的70%左右，养殖量平均每年以8%~10%的速度增加。目前水禽养殖业已经成为我国解决“三农”问题的支柱产业，也是现代畜牧业的重要组成部分。虽然水禽业发展迅速，但也存在诸多制约因素，尤其是疫病多发、频发严重阻碍水禽业的健康发展。其中，由鸭疫里默氏杆菌Riemerella anatipestifer引起的鸭疫里默氏杆菌病已经成为影响水禽业经济最严重的细菌性传染病，发病率可达100%，死亡率高于75%，耐过的病鸭生长发育受阻，成为僵鸭，造成巨大的经济损失[1]。一直以来主要采用疫苗免疫与药物治疗相结合的方式防控该病，但是防控的效果并不理想。一方面，鸭疫里默氏杆菌有21个血清型且各血清型间缺乏交叉保护性，而目前市售的灭活疫苗尚未涵盖所有流行的血清型，当出现疫苗株血清型以外的鸭疫里默氏杆菌感染时，现有的灭活疫苗不能起到保护作用，导致疫苗应用受限；另一方面，虽然鸭疫里默氏杆菌对大多数抗生素或化学合成药物较敏感，但是由于长期连续使用或不规范使用药物导致耐药性普遍存在[2]，并且已经成为当前养殖业面临的最严峻问题。因此，本研究拟对广东地区鸭疫里默氏杆菌分离株进行血清型和耐药性调查及遗传进化关系研究，为开展针对性疫苗免疫预防及选用敏感药物治疗提供重要依据。
1. 材料与方法
1.1 材料
1.1.1 病料样本
病料于2015—2019年采自广东省清远、惠州、佛山等地区鸭场，采集疑似患鸭疫里默氏杆菌病病死鸭的脑、肝脏、心脏等组织，病料样本及分离鸭疫里默氏杆菌信息见表1，相同养殖场采集的样本分离的鸭疫里默氏杆菌记为1株。
表 1 采集样本与分离鸭疫里默氏杆菌信息Table 1. Collected samples and isolation information of Riemerella anatipestifer地区
District采集时间
Acquisition time养殖场数
Farm number样本数
Sample number分离株数
Isolate number代表性株数
Representative isolate number韶关 Shaoguan 2015−05—2019−04 19 69 15 5 河源 Heyuan 2015−07—2019−11 16 56 13 5 云浮 Yunfu 2015−09—2019−12 30 112 25 4 肇庆 Zhaoqing 2015−04—2019−10 17 63 12 4 清远 Qingyuan 2015−03—2019−12 23 91 19 5 佛山 Foshan 2015−05—2019−12 15 57 12 5 茂名 Maoming 2015−04—2019−11 20 80 16 5 惠州 Huizhou 2015−05—2019−11 23 80 17 5 汕头 Shantou 2015−03—2019−12 25 91 21 5 湛江 Zhanjiang 2015−05—2019−12 22 82 18 5 总计 Total 2015—2019 210 781 168 48 1.1.2 试剂与药品
改良胰蛋白胨大豆肉汤(TSB)培养基，含5%(φ)新生牛血清；胰酶大豆琼脂(TSA)培养基，含5%(φ)新生牛血清；麦康凯培养基按细菌学常规方法配制；革兰染色试剂和微量生化发酵管为广东环凯微生物科技有限公司产品；DNA marker DL 2000和Premix Taq为TaKaRa公司产品；药物标准品——氧氟沙星(Ofloxacin，OFX)、诺氟沙星(Norfloxacin，NOR)、盐酸四环素(Tetracycline hydrochloride，TCY)、头孢噻肟(Cefotaxime，CTX)、阿莫西林(Amoxicillin，AMX)、土霉素(Oxytetracycline，OXY)、盐酸金霉素(Chlortetracycline hydrochloride，CTE)、氨苄西林(Ampicillin，AMP)、庆大霉素(Gentamicin，GEN)、盐酸环丙沙星(Ciprofloxacin hydrochloride，CIP)、卡那霉素(Kanamycin，KAN)、大观霉素(Spectinomycin，STP)、磺胺嘧啶(Sulfadiazine，SDI)、磺胺二甲嘧啶(Sulfadimidine，SUL)和磺胺对甲氧嘧啶(Sulfametoxydiazine，SMD)为中国兽医药品监察所产品。
1.1.3 试验动物
1日龄商品代樱桃谷肉鸭，大约4000只，购自广州郊区某鸭场，未接种任何疫苗和抗血清，饲养于洁净环境中，至7~14日龄进行动物回归试验。
1.1.4 引物合成
引物参照文献[3]合成，以16S rRNA为扩增基因。引物( RA-F：5′-ACGTCATCCCACCTTCCTC-3′，RA-R：5′-GTTCAGACTAA GCGAAAG-3′ ) 由Invitrogen公司合成。
1.1.5 参考菌株
血清1型鸭疫里默氏杆菌参考株由广东省农科院动物卫生研究所寄生生物学研究室保存。
1.1.6 定型血清
1、2、3、4、5、6、7、8、9、10型鸭疫里默氏杆菌参考株定型血清，由广东省农科院动物卫生研究所寄生生物学研究室自制。
1.1.7 质控菌
大肠埃希菌ATCC25922购自中国兽医药品监察所。
1.2 方法
1.2.1 病原分离与鉴定
病原分离：将无菌采集的病死鸭的脑、肝脏、心脏组织，划线接种于TSA培养基平板和麦康凯培养基平板，于37 ℃厌氧培养24 h，挑取单个可疑菌落进行纯化培养。
革兰染色：挑取纯培养的单菌落进行革兰染色、镜检，观察细菌的形态和染色特性。
生化试验：将分离株的纯培养物接种于微量生化发酵管中，于37 ℃温箱中培养，观察结果。
PCR鉴定：以分离株的菌液为模板进行PCR扩增，PCR反应总体积为20 μL：Premix Taq 10 μL、上下游引物各1 μL、菌液1 μL、ddH2O 7 μL；反应条件：94 ℃预变性5 min，以94 ℃ 30 s、55 ℃ 30 s、72 ℃ 1 min进行30个循环，72 ℃延伸7 min；反应结束后，将PCR产物以10 g/L琼脂糖凝胶电泳进行检测，目的条带大小为325 bp。
动物致病性试验：将168株鸭疫里默氏杆菌分离株分别接种7~14日龄试验鸭，每株接种10只试验鸭，同时对照组设10只试验鸭，连续观察14 d；接种组腿部肌肉注射鸭疫里默氏杆菌分离株活菌1×108 CFU/mL，每只注射1 mL，对照组接种生理盐水，每只注射1 mL。隔离饲养，观察致病力，记录试验鸭发病及死亡情况；同时根据试验鸭死亡和发病情况在接种后1~5 d采集3~4只试验鸭的心脏和肝脏，分离鸭疫里默氏杆菌，并进行革兰染色、生化试验和PCR鉴定。
血清型鉴定：取洁净的载玻片分别滴加30 μL参考株阳性血清，再滴加30 μL分离株菌液，充分混匀，观察1~2 min，以出现清晰凝集者判为阳性，否则判为阴性。
1.2.2 最低抑菌浓度的测定
从168株鸭疫里默氏杆菌分离株中选择48株作为代表性菌株，进行最低抑菌浓度(Minimum inhibitory concentration，MIC)测定。代表性菌株选择依据：所选菌株均来自不同的规模化鸭场，每个鸭场养殖时间均在5年以上，每批肉鸭养殖规模均在20 000只以上，从韶关、河源、清远、佛山、茂名、惠州、汕头、湛江的规模化鸭场随机各选5株，从云浮和肇庆的规模化鸭场随机各选4株，总计48株。
抗菌药物贮存液的制备：将药物制成2 560 mg/L的原液，抗菌药物贮存液按照梯度稀释。
抗菌药物浓度范围：根据抗菌药物敏感性试验操作标准，药物浓度范围应包含耐药、中介和敏感分界点值，特殊情况除外。
接种物的制备：采用直接菌落悬液配制法，将培养18~24 h的鸭疫里默氏杆菌菌落调配成0.5麦氏比浊标准的悬液，再以TSB培养基按体积比1∶100稀释后备用；在15 min内接种完配制好的接种物，并取1份接种物在TSA平板上培养，以检查接种物的纯度。
抗菌药物稀释及菌液接种：以改良的试管两倍稀释法[4]测定鸭疫里默氏杆菌分离株的MIC，取14支灭菌试管为一排，除第1管加入1.6 mL TSB培养基外，其余每支试管均加入1 mL TSB培养基，于第1管中添加0.4 mL待测药物原液并混匀，吸取1 mL加入第2管，混匀后再吸取1 mL至第3管，如此连续倍比稀释至第12管，混匀后吸取1 mL弃掉，第13管为不加药物的阳性对照，第14管为不加药物和菌液、只加生理盐水的阴性对照；向第1~13管中各加入1 mL稀释菌液，此时第1~12管药物终质量浓度分别为256、128、64、32、16、8、4、2、1、0.5、0.25、0.125 mg/L。
培养：将接种后的试管置于37 ℃温箱中培养22~24 h。
结果判断：在阳性对照和阴性对照成立且质控菌株的MIC处于质控范围的情况下，肉眼可见澄清的最低药物浓度管即为测试菌株的MIC。
药物敏感性判定：根据标准[4-5]判断药物耐药、敏感和中介。
1.2.3 全基因组序列特征分析与核心基因组进化树的构建
基因组DNA文库的构建与序列测定：使用天根细菌基因组DNA提取试剂盒，按照操作说明，提取鸭疫里默氏杆菌分离株的基因组DNA；通过浓度测定和琼脂糖凝胶电泳，初步质检合格后，构建Illumina(250 bp)文库，利用Illumina Novaseq 6000平台进行测序，下机数据经质控和去除接头后运用SPAdes v3.6.2对原始数据进行拼接[6]，拼接后的数据使用Quast软件进行质控[7]。
序列特征分析。耐药基因分析：将拼接好的测序结果上传至CGE的ResFinder4.1数据库[8]以分析菌株所携带的耐药基因(ARGs)( https://cge.cbs.dtu.dk/services/ResFinder/)；ST分型：将拼接好的测序结果上传至多位点序列分型(Multi-locus sequence typing, MLST)数据库以分析菌株所属的ST型( https://cge.cbs.dtu.dk/services/MLST/)。
构建核心基因组进化树：从NCBI Genome下载数据库中全部共计59株鸭疫里默氏杆菌的全基因组序列，以标准菌株ATCC 11845(NCBI登录号：GCA_000252855.1)作为参考菌株，利用Parsnp软件比对全基因组序列并构建核心基因组进化树[9]，使用hierBAPS软件进行进化树分簇[10]，再利用FigTree v1.4.2和iTOL v4软件进行美化注释[11]。
2. 结果与分析
2.1 鸭疫里默氏杆菌分离与鉴定
分离株在TSA平板上长成光滑、湿润、半透明的奶油状菌落，在麦康凯平板上不生长；革兰染色呈阴性的单个或成双排列的短小杆菌；生化试验显示分离株均不发酵蔗糖、乳糖、木糖、葡萄糖、果糖、甘露糖、甘露醇、山梨醇，不产生硫化氢，不利用柠檬酸盐，甲基红试验、靛基质试验、硝酸盐还原试验阴性，氧化酶、过氧化物酶试验阳性；通过特异性PCR扩增，分离株均可扩增到约300 bp的条带，与目的条带大小一致(图1)。
![]()
图 1 鸭疫里默氏杆菌分离株PCR鉴定结果M：DNA标准DL2000；1、2：鸭疫里默氏杆菌分离株；3：阳性对照；4：阴性对照Figure 1. PCR identification results of Riemerella anatipestifer isolatesM: DNA marker DL2000; 1, 2: Riemerella anatipestifer isolates; 3: Positive control; 4: Negative control2.2 动物回归试验
所有鸭疫里默氏杆菌分离株对试验鸭均有不同程度的致病力，发病率达100%，致死率介于70%~100%，其中70.24%(118/168)的分离株的毒力较强，在接种后48 h内试验鸭全部死亡，剖检未见典型的鸭疫里默氏杆菌病病理变化，其余分离株接种试验鸭后病程较长，剖检可见典型的鸭疫里默氏杆菌病病理变化。所有试验鸭接种分离株后均表现鸭疫里默氏杆菌病临床症状，主要表现精神沉郁，采食、饮水减少或废绝，排绿色或黄绿色粪便，共济失调，头颈震颤，角弓反张及转圈运动等，严重者导致死亡。从接种鸭疫里默氏杆菌分离株试验鸭的肝脏或心脏中均可分离到接种菌，革兰染色均为阴性的短小杆菌，生化试验结果符合鸭疫里默氏杆菌特性，PCR鉴定结果与分离株初次分离鉴定结果一致。因此，根据细菌分离鉴定结果和动物回归试验结果确定从广东地区的鸭场共分离到168株鸭疫里默氏杆菌菌株。
2.3 血清型鉴定
鸭疫里默氏杆菌分离株血清型鉴定结果见表2，血清1、2、3、5、6、7、8、10型和未定型均有流行，其中高达54.17%(91/168)的分离株为1型，其次为2型，占27.97%(47/168)。因此，广东地区鸭疫里默氏杆菌分离株以1型为优势血清型。
表 2 各地区各血清型鸭疫里默氏杆菌分离株数量Table 2. Isolate number of each serotype for Riemerella anatipestifer isolates from different districts地区
District分离株数
Isolate number血清型 Serotype 1 2 3 4 5 6 7 8 9 10 未定型 Undefined 韶关 Shaoguan 15 7 2 2 0 0 1 0 1 0 0 2 河源 Heyuan 13 7 4 1 0 0 0 0 0 0 1 0 云浮 Yunfu 25 15 8 0 0 1 0 0 0 0 0 1 肇庆 Zhaoqing 12 6 3 1 0 0 0 0 1 0 1 0 清远 Qingyuan 19 10 6 1 0 1 0 0 0 0 0 1 佛山 Foshan 12 7 4 0 0 0 0 1 0 0 0 0 茂名 Maoming 16 8 5 1 0 0 0 0 1 0 1 0 惠州 Huizhou 17 12 3 0 0 0 0 1 0 0 1 0 汕头 Shantou 21 12 7 0 0 0 0 0 0 0 2 0 湛江 Zhanjiang 18 7 5 0 0 1 1 1 2 0 0 1 总计 Total 168 91 47 6 0 3 2 3 5 0 6 5 2.4 药物敏感性试验
药物对鸭疫里默氏杆菌代表性菌株的MIC测定结果及敏感性统计结果见表3。根据MIC分析各药物抑制50%和90%鸭疫里默氏杆菌菌株生长的MIC，即MIC50和MIC90，MIC50介于0.25~>256 mg/L，其中，头孢噻肟最低，仅0.25 mg/L，阿莫西林为1 mg/L，氨苄西林、土霉素、盐酸四环素、盐酸金霉素、氧氟沙星均为8 mg/L，盐酸环丙沙星为16 mg/L，诺氟沙星、大观霉素、庆大霉素分别为32、64、128 mg/L，磺胺嘧啶为256 mg/L，卡那霉素、磺胺二甲嘧啶、磺胺对甲氧嘧啶均为>256 mg/L；MIC90值介于8~>256 mg/L，其中，头孢噻肟最低，仅8 mg/L；土霉素和盐酸金霉素为16 mg/L；盐酸四环素、氧氟沙星、盐酸环丙沙星均为32 mg/L；诺氟沙星为64 mg/L；阿莫西林、氨苄西林、大观霉素均为128 mg/L；卡那霉素、庆大霉素、磺胺嘧啶、磺胺二甲嘧啶、磺胺对甲氧嘧啶较高，均为>256 mg/L。
表 3 48株鸭疫里默氏杆菌代表性分离株最低抑菌浓度(MIC)测定结果1)Table 3. The minimum inhibitory concentration (MIC) of 48 representative isolates of Riemerella anatipestifer药物类型
Drug type药物1)
DrugMIC/
(mg·L−1)MIC50/
(mg·L−1)MIC90/
(mg·L−1)占比/% Proportion 敏感
Sensitive中介
Intermediate耐药
Resistantβ−内酰胺类 AMX <0.125~256 1 128 75.00(36/48) 6.25(3/48) 18.75(9/48) β-lactams AMP <0.125~256 8 128 64.58(31/48) — 35.42(17/48) CTX <0.125~256 0.25 8 91.67(44/48) 4.17(2/48) 4.17(2/48) 四环素类 OXY <0.125~128 8 16 20.83(10/48) 14.58(7/48) 64.58(31/48) Tetracyclines TCY <0.125~128 8 32 18.75(9/48) 10.42(5/48) 70.83(34/48) CTE <0.125~32 8 16 18.75(9/48) 20.83(10/48) 60.42(29/48) 氨基糖苷类 GEN 0.25~>256 128 >256 8.33(4/48) — 91.67(44/48) Aminoglycosides KAN 2~>256 >256 >256 8.33(4/48) 2.08(1/48) 89.58(43/48) STP 16~256 64 128 45.83(22/48) 25.00(12/48) 29.17(14/48) 喹诺酮类 OFX 0.5~64 8 32 20.83(10/48) 6.25(3/48) 72.92(35/48) Quinolones NOR 0.25~128 32 64 14.58(7/48) 16.67(8/48) 68.75(33/48) CIP 0.5~64 16 32 14.58(7/48) 4.17(2/48) 81.25(39/48) 磺胺类 SDI 64~>256 256 >256 52.08(25/48) — 47.92(23/48) Sulfonamides SUL 32~>256 >256 >256 39.58(19/48) — 60.42(29/48) SMD 128~>256 >256 >256 33.33(16/48) — 66.67(32/48) 1)AMX：阿莫西林，AMP：氨苄西林，CTX：头孢噻肟，OXY：土霉素，TCY：盐酸四环素，CTE：盐酸金霉素，GEN：庆大霉素，KAN：卡那霉素，STP：大观霉素，OFX：氧氟沙星，NOR：诺氟沙星，CIP：盐酸环丙沙星，SDI：磺胺嘧啶，SUL：磺胺二甲嘧啶，SMD：磺胺对甲氧嘧啶；MIC50和MIC90为能抑制50%和90%细菌生长所需的MIC；“—”表示未检测出菌株
1) AMX: Amoxicillin, AMP: Ampicillin, CTX: Cefotaxime, OXY: Oxytetracycline, TCY: Tetracycline hydrochloride, CTE: Chlortetracycline hydrochloride, GEN: Gentamicin, KAN: Kanamycin, STP: Spectinomycin, OFX: Ofloxacin, NOR: Norfloxacin, CIP: Ciprofloxacin hydrochloride, SDI: Sulfadiazine, SUL: Sulfadimidine, SMD: Sulfametoxydiazine; MIC50 and MIC90 are MIC values that can inhibit 50% and 90% bacteria growth; “—” indicates no detectable isolate分析鸭疫里默氏杆菌代表性菌株对受试药物的敏感性，对庆大霉素的耐药率高达91.67%(44/48)；对卡那霉素和盐酸环丙沙星的耐药率分别为89.58%(43/48)和81.25%(39/48)；对土霉素、盐酸四环素、盐酸金霉素、氧氟沙星、诺氟沙星、磺胺二甲嘧啶、磺胺对甲氧嘧啶的耐药率介于60%~80%；对磺胺嘧啶和氨苄西林的耐药率介于30%~60%；对阿莫西林、头孢噻肟和大观霉素的耐药水平较低(30%以下)。
根据药敏试验结果，分析鸭疫里默氏杆菌代表性菌株的多重耐药性，结果如图2所示。48株鸭疫里默氏杆菌对5~12种药物耐药，其中9耐和10耐的菌株总数最多，占54.17%(26/48)，6耐和11耐、7耐与8耐的菌株数相同，分别为4株和5株，而5耐和12耐菌株数最少，各有2株。
![]()
图 2 48株鸭疫里默氏杆菌代表性分离株多重耐药性分布Figure 2. Results of multi-drug resistance of 48 Riemerella anatipestifer representative isolates统计并分析鸭疫里默氏杆菌分离株的耐药谱，48株鸭疫里默氏杆菌具有44种耐药谱型，其中40株为单一谱型，构成比为2.08%(1/48)，其余4种耐药谱型分别为2株所共有，构成比为4.16%(2/48) (表4)。
表 4 48株鸭疫里默氏杆菌代表性分离株耐药谱统计结果Table 4. Drug resistance spectrum of 48 representative isolates of Riemerella anatipestifer药物数量
Drug number耐药谱1)
Resistant spectrum菌株数
Isolate number占比/%
Proportion5 OFX-NOR-OXY-GEN-KAN
OFX-CTE-GEN-CIP-KAN1
12.08
2.086 OFX-NOR-TCY-GEN-CIP-KAN
OFX-NOR-GEN-CIP-KAN-SMD
NOR-GEN-CIP-KAN-SUL-SMD
CIP-KAN-STP-SUL-SMD-SDI1
1
1
12.08
2.08
2.08
2.087 OXY-TCY-CTE-GEN-KAN-SUL-NOR
OXY-GEN-KAN-SDI-SUL-SMD-CIP
OXY-TCY-CTE-GEN-KAN-OFX-CIP
OXY-TCY-CTE-GEN-KAN-NOR-CIP1
1
1
22.08
2.08
2.08
4.168 OFX-NOR-TCY-CTE-GEN-CIP-SUL-SMD
NOR-TCY-CTE-GEN-CIP-KAN-SDI-SMD
CIP-GEN-SDI-SUL-SMD-NOR-TCY-AMP
AMX-AMP-OXY-TCY-CTE-GEN-KAN-OFX
AMP-STP-SDI-SUL-SMD-GEN-KAN-OXY1
1
1
1
12.08
2.08
2.08
2.08
2.089 OFX-NOR-TCY-OXY-CTE-GEN-CIP-KAN-SUL
OFX-NOR-TCY-CIP-KAN-SDI-SUL-SMD-GEN
OFX-NOR-AMX-CTE-GEN-CIP-KAN-SUL-SMD
OXY-TCY-GEN-KAN-OFX-NOR-CIP-SDI-SUL
TCY-GEN-CIP-KAN-SDI-SMD-CTE-OFX-SUL
OXY-AMP-GEN-KAN-STP-SDI-SUL-SMD-NOR
OFX-NOR-OXY-GEN-CIP-KAN-SUL-SMD-SMD
OFX-TCY-OXY-GEN-KAN-STP-SDI-SUL-SMD
AMP-CTX-SDI-GEN-SMD-OFX-STP-KAN-CIP
AMX-AMP-OXY-TCY-STP-KAN-OFX-SDI-GEN
SMD-ERY-OFX-SDI-CIP-AMX-AMP-CTE-TCY
NOR-TCY-OXY-CTE-GEN-CIP-KAN-SUL-AMP2
1
1
2
1
1
1
1
1
1
1
14.16
2.08
2.08
4.16
2.08
2.08
2.08
2.08
2.08
2.08
2.08
2.0810 NOR-TCY-OXY-CTE-GEN-CIP-KAN-STP -SMD-SDI
AMP-OXY-TCY-CTE-GEN-KAN-OFX-NOR-CIP-SMD
OFX-NOR-TCY-CTE-GEN-CIP-KAN-SDI-SMD-SUL
OFX-NOR-TCY-OXY-CTE-GEN-CIP-KAN-SMD-SDI
OFX-NOR-AMX-AMP-GEN-CIP-KAN-SUL-SMD-CTE
OFX-NOR-TCY-AMX-OXY-CTE-GEN-CIP-KAN-SDI
OFX-NOR-TCY-OXY-CTE-GEN-CIP-KAN-SUL-SMD
TCY-AMP-GEN-CIP-KAN-STP-SDI-SUL-SMD-KAN
AMX-AMP-OXY-TCY-CTE-GEN-KAN-OFX-SDI-SUL
OFX-NOR-OXY-CTE-AMP-GEN-CIP-KAN-SUL-SMD
OFX-NOR-AMX-AMP-GEN-CIP-KAN-SDI-SMD-STP
TCY-OXY-GEN-KAN-STP-SDI-SUL-SMD-CTX-AMP1
1
1
1
1
1
1
1
1
1
1
12.08
2.08
2.08
2.08
2.08
2.08
2.08
2.08
2.08
2.08
2.08
2.0811 OFX-NOR-TCY-CTE-OXY-GEN-CIP-KAN-STP-SDI-SMD
OFX-NOR-TCY-OXY-CTE-AMP-GEN-CIP-KAN-STP-SUL
OFX-NOR-TCY-OXY-CTE-AMP-GEN-CIP-KAN-STP-SMD2
1
14.16
2.08
2.0812 AMX-OXY-TCY-CTE-STP-KAN-OFX-NOR-SUL-SMD-SMD-CIP
OFX-NOR-TCY-OXY-AMP-GEN-CIP-STP-SUL-SMD-SDI-CTE1
12.08
2.081) AMX：阿莫西林，AMP：氨苄西林，CTX：头孢噻肟，OXY：土霉素，TCY：盐酸四环素，CTE：盐酸金霉素，GEN：庆大霉素，KAN：卡那霉素，STP：大观霉素，OFX：氧氟沙星，NOR：诺氟沙星，CIP：盐酸环丙沙星，SDI：磺胺嘧啶，SUL：磺胺二甲嘧啶，SMD：磺胺对甲氧嘧啶
1) AMX: Amoxicillin, AMP: Ampicillin, CTX: Cefotaxime, OXY: Oxytetracycline, TCY: Tetracycline hydrochloride, CTE: Chlortetracycline hydrochloride, GEN: Gentamicin, KAN: Kanamycin, STP: Spectinomycin, OFX: Ofloxacin, NOR: Norfloxacin, CIP: Ciprofloxacin hydrochloride, SDI: Sulfadiazine, SUL: Sulfadimidine, SMD: Sulfametoxydiazine2.5 全基因组测序分析
由北京诺禾致源科技股份有限公司对48株鸭疫里默氏杆菌代表性菌株进行全基因组测序，经过Quast软件质检后，其中2株拼接质量不高，故仅有46株进行后续分析。共检出6种耐药基因，其中携带四环素类耐药基因tet(X)的菌株最多，阳性率高达82.60%(38/46)；大环内酯类耐药基因erm(F)阳性率为73.91%(34/46)，而大环内酯类另一耐药基因ere(D)的阳性率仅有19.57%(9/46)；β−内酰胺类耐药基因bla OXA-209的阳性率为54.35%(25/46)；氨基糖苷类耐药基因aadS和氟苯尼考耐药基因floR的阳性率相对较低，分别为34.78%(16/46)和30.43%(14/46)。耐药基因结果见表5。
表 5 46株鸭疫里默氏杆菌代表性分离株耐药基因携带情况Table 5. The resistance genes carried by 46 Riemerella anatipestifer representative isolates药物
Drug耐药基因
Resistant gene阳性率/%
Positive rate大环内酯类 erm(F) 73.91(34/46) Macrolides ere(D) 19.57(9/46) 四环素类
Tetracyclinestet(X) 82.60(38/46) 氟苯尼考
FlorfenicolfloR 30.43(14/46) 氨基糖苷类
AminoglycosidesaadS 34.78(16/46) β−内酰胺类 β-lactams bla OXA-209 54.35(25/46) 基于获得的46株鸭疫里默氏杆菌全基因组序列，进一步分析其MLST，其中18株ST分型成功，分别为ST38(4株)、ST13(2株)、ST24(2株)、ST35(2株)、ST44(2株)，其他ST型分别为ST2、ST16、ST21、ST34、ST42和ST43，均只有1株，其余28株为未知ST型；说明本研究测序菌株的ST型呈多样性，未出现优势ST型。
将本研究46株鸭疫里默氏杆菌菌株的全基因组测序结果与数据库中59株鸭疫里默氏杆菌菌株的全基因组数据进行比对，构建核心基因组进化树(图3)，数据库59株鸭疫里默氏杆菌菌株中，39株来源于中国，9株来源于英国，2株来源于印度，德国、俄罗斯和美国各有1株，另有6株未标注来源。总计105株鸭疫里默氏杆菌菌株分布在6个分簇群系中，本研究所有测序菌株与数据库中来自中国的菌株相似度较高，且集中分布在3个簇群中，其中Clade 1群系最大，包含67株，我国菌株占50.75%(34/67)，其次为Clade 3群系，包含17株，我国菌株占47.06%(8/17)，Clade 2 群系包含5株鸭疫里默氏杆菌，我国菌株数占60.00%(3/5)。数据库中大部分鸭疫里默氏杆菌菌株来自于中国，而欧洲来源的菌株占比较低，并且观察到来自世界各地的鸭疫里默氏杆菌菌株间存在地理分布差异。另外，数据库中我国菌株从1990年到2017年间均有收集，持续时间接近30年，时间跨度较大，但仍能与本研究测序菌株在进化树中聚集在一起，说明在我国鸭疫里默氏杆菌存在区域内长时间持续传播的特点。
![]()
图 3 105株鸭疫里默氏杆菌菌株核心基因组系统进化树Figure 3. Phylogenetic tree of core genomes for 105 isolates of Riemerella anatipestifer3. 讨论与结论
3.1 讨论
鸭疫里默氏杆菌病血清型复杂，在世界范围内流行，目前在国内多个地区均有鸭疫里默氏杆菌血清型流行情况的报道，广西百色分离株血清型比较单一，全部为1型[12]；山东分离株至少存在1、2、6、10和11型5种血清型，其中以1型和2型为优势血清型[13-14]；超过半数的江苏分离株为1型[15]；安徽分离株则多达7种血清型，但以1型和2型发生最广泛[16]；河南分离株以1型和2型多发，此外还流行10型和未定型[17]；血清1、2、3、4、8、10、15型在广东地区均有流行，以1型最多发[18-20]。与国内情况类似，国外也流行多种血清型，美国以1型为主，此外还流行1、2、5、11、13、15、19、21型[21]；1型和2型为英国多发血清型，同时5、9、13和15型也有流行[22]；泰国主要流行1、5、6、7、10和21型[23]；新加坡以1、5和10型这3种血清型为主[24]；1型和3型为丹麦的优势血清型[25]；韩国只有1、4和7型这3种血清型[26]；澳大利亚则流行1、6、8、9、13和14型[27]。本研究高达54.17%的分离株为血清1型，也是优势血清型。分析国内外鸭疫里默氏杆菌血清型流行情况发现多血清型流行已呈基本态势，其中1型为最普遍发生的血清型，也是重点防控的血清型。
在养殖过程中，由于长期应用抗菌药物治疗鸭疫里默氏杆菌病，导致广泛而严重的耐药性。调查显示山东分离株对氨基糖苷类药物高度耐药，并且具有多重耐药性[28]；安徽分离株对恩诺沙星、新霉素、安普霉素表现不同程度的耐药性[29]；贵州分离株则对氨基糖苷类和大环内酯类等药物表现高度耐药，多重耐药高达17耐[30]；同样，广东地区分离株耐药状况也比较严峻，朱元军等[31]报道分离株对丁胺卡那、卡那霉素、庆大霉素和链霉素高度耐药，多数分离株表现七重以上耐药。另外，耐药性还具有随药物使用时间延长进一步加剧的现象，从同一地区采集的分离株对红霉素和多黏菌素B的耐药率呈逐年升高的趋势[32-33]。本研究采样养殖场经常采用庆大霉素、卡那霉素、磺胺嘧啶、头孢类药物、喹诺酮类药物、环丙沙星、恩诺沙星、四环素、土霉素、多黏菌素B、林可霉素等药物进行治疗，从养殖场样本中分离的鸭疫里默氏杆菌对头孢噻肟、阿莫西林、大观霉素的耐药率低于30%，对土霉素、盐酸四环素、盐酸金霉素、庆大霉素、卡那霉素等药物均表现高度耐药。因此，总体上养殖场用药与耐药表型具有一定的对应性，在临床选择药物时，最好先进行药敏试验，选用抗菌效果较好的药物，同时也要结合本养殖场的用药史合理用药。
通过分析本研究测序菌株耐药表型与耐药基因的相关性，发现携带tet(X)基因的菌株对四环素类药物全部表现耐药，因此，tet(X)基因是介导四环素耐药的最重要基因[34]；氨基糖苷类耐药率高于耐药基因aadS的携带率，说明除耐药基因外还存在其他因素导致氨基糖苷类耐药，这一结果与蔡秀磊[35]研究结论一致，即氨基糖苷类耐药还与基因盒−整合子系统有关；bla OXA-209的阳性率与耐药率基本一致。由于本研究未对大环内酯类和氟苯尼考进行药物敏感性测试，因此，尚不清楚这2类药物耐药基因携带率与耐药表型的相关性。总体上，耐药表型与耐药基因具有一定的相关性，但是也要明确鸭疫里默氏杆菌分离株产生耐药性的原因比较复杂，除与耐药基因有关外，还与环境、药物等诸多因素有关。
通过系统进化树发现，来自世界各地的鸭疫里默氏杆菌菌株间存在地理分布差异，本研究的测序菌株主要分布在Clade 1簇群，且血清型全部为1型，由于数据库中其他菌株没有相关血清型信息，因此，具体血清型不明确，但可以推测，该簇群的分离株与血清1型密切相关。
鉴于本研究鸭疫里默氏杆菌血清型复杂及耐药性严重的现状，应首先针对当地或本养殖场主要流行的血清型选取相应的灭活疫苗进行免疫，当出现疫苗株以外的血清型感染时，则以当地或本养殖场的分离株制备灭活疫苗，同时选择敏感药物并合理用药辅助治疗，可以达到更有效的防治效果。
3.2 结论
本研究在对广东地区鸭疫里默氏杆菌分离株的血清型和耐药性初步了解的基础上，又结合全基因组测序数据对耐药基因、ST分型、遗传进化树做进一步分析，发现分离株的优势血清型为1型，耐药性严重，所携带的耐药基因与耐药表型具有一定的相关性，ST型呈多样性，与MLST 数据库中来自我国的菌株遗传背景相近，研究结果为鸭疫里默氏杆菌病疫苗免疫预防与药物治疗及掌握鸭疫里默氏杆菌遗传进化特征提供了依据。
-
![]()
图 1 复合肥不同施用量对巴戟天不同部位N、P、K含量的影响
各小图中，相同部位柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 1. Effects of different application amounts of compound fertilizer on contents of N, P and K in different parts of Morinda officinalis
In each figure, different lowercase letters on the columns of the same part indicate significant differences among different treatments (P<0.05, Duncan’s method)
![]()
图 2 复合肥不同施用量对巴戟天各部位N、P、K单株累积量的影响
各小图中，不同小写字母表示处理间地上部、根部养分累积量差异显著，不同大写字母表示处理间全株养分累积量差异显著(P<0.05，Duncan’s法)
Figure 2. Effects of different application amounts of compound fertilizer on the accumulation of N, P and K per plant in different parts of Morinda officinalis
In each figure, different lowercase letters on the columns indicate significant differences in aboveground and root nutrient accumulation among treatments, different capital letters indicate significant differences in whole plant nutrient accumulation among treatments (P < 0.05, Duncan’s method)
![]()
图 3 复合肥不同施用量对巴戟天4种寡糖含量的影响
相同寡糖柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 3. Effects of different application amounts of compound fertilizer on contents of four oligosaccharides in Morinda officinalis
Different lowercase letters on the columns of the same oligosaccharide indicate significant differences among treatments (P<0.05, Duncan’s method)
![]()
图 4 复合肥不同施用量对巴戟天4种寡糖单株累积量的影响
相同寡糖柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 4. Effects of different application amounts of compound fertilizer on accumulation of four oligosaccharides in Morinda officinalis
Different lowercase letters on the columns of the same oligosaccharide indicate significant differences among treatments (P<0.05, Duncan’s method)
![]()
图 5 巴戟天不同指标的相关性分析
“*”和“**”分别表示在P<0.05和P<0.01水平显著相关(Pearson法)
Figure 5. Correlation analyses of different indicators for Morinda officinalis
“*” and “**” indicate significant correlations at P < 0.05 and P < 0.01 levels, respectively (Pearson method)
![]()
图 6 复合肥不同施用量的主成分分析
PC1：第一主成分，PC2：第二主成分
Figure 6. Principal component analyses of different application amounts of compound fertilizer
PC1: The first principal component, PC2: The second principal component
表 1 复合肥不同施用量对巴戟天生长指标的影响1)
Table 1 Effects of different application amounts of compound fertilizer on growth indexes of Morinda officinalis
处理
Treatment地径/mm
Ground diameter总生物量/g
Total biomass根部产量/g
Root yieldCK 5.19±0.35a 6.12±0.31d 0.76±0.04c T1 5.25±0.14a 8.00±0.14c 0.89±0.05c T2 6.32±0.54a 9.65±0.61bc 1.11±0.08b T3 6.07±0.31a 13.96±0.78a 1.63±0.04a T4 6.22±0.44a 10.88±0.24b 1.22±0.08b T5 5.72±0.25a 10.40±0.82b 1.26±0.05b 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among different treatments (P<0.05, Duncan’s method)
下载: 导出CSV
表 2 复合肥不同施用量的肥料利用效率1)
Table 2 Fertilizer utilization efficiency of different application amounts of compound fertilizer
% 处理 Treatment N P K T1 36.97±1.37a 2.69±0.31ab 9.16±0.53b T2 27.14±4.41b 2.31±0.07b 8.04±2.05b T3 41.35±2.37a 3.23±0.20a 16.62±2.63a T4 20.97±1.29bc 2.11±0.17b 7.27±1.00b T5 18.38±1.20c 1.18±0.10c 6.11±0.56b 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s method)
下载: 导出CSV
表 3 主成分分析的载荷系数
Table 3 Loading coefficients for principal component analysis
评价指标
Evaluation
index第一主成分
The first
principal
component第二主成分
The second
principal
component总生物量 Total biomass 0.952 −0.083 产量 Yield 0.959 −0.154 全株N含量
Whole plant N content0.843 0.353 全株P含量
Whole plant P content0.708 0.475 全株K含量
Whole plant K content0.961 −0.141 蔗糖积累量
Sucrose accumulation0.946 −0.189 1−蔗果三糖积累量
1-Kestose accumulation0.935 −0.232 耐斯糖积累量
Nystose accumulation0.755 −0.372 1F−果呋喃糖基耐斯糖累积量
1F-Fructofuranosyl nystose accumulation0.618 0.655 特征值 Eigenvalue 6.682 1.057 方差贡献率/%
Contribution of variance74.246 11.748 累计贡献率/%
Cumulative contribution rate74.246 85.994
下载: 导出CSV
表 4 复合肥不同施用量的主成分分析综合得分
Table 4 Comprehensive score of principal component analysis of different application amounts of compound fertilizer
处理
Treatment主成分分值
Principal component value综合得分(F)
Comprehensive
score排序
RankingF1 F2 CK −3.92 −0.49 −2.96 6 T1 −2.28 0.27 −1.66 5 T2 −0.01 −0.86 −0.11 4 T3 3.48 −0.52 2.53 1 T4 1.07 1.17 0.93 3 T5 1.65 0.42 1.27 2
下载: 导出CSV
-
[1] 国家药典委员会. 中华人民共和国药典: 1部[M]. 北京: 中国医药科技出版社, 2020: 83. [2] YANG Z, HU J, ZHAO M. Isolation and quantitative determination of inulin-type oligosaccharides in roots of Morinda officinalis[J]. Carbohydrate Polymers, 2011, 83(4): 1997-2004. doi: 10.1016/j.carbpol.2010.11.006
[3] CHEN D, LI N, LIN L, et al. Confocal mirco-Raman spectroscopic analysis of the antioxidant protection mechanism of the oligosaccharides extracted from Morinda officinalis on human sperm DNA[J]. Journal of Ethnopharmacology, 2014, 153(1): 119-124. doi: 10.1016/j.jep.2014.01.021
[4] YANG X, CHEN D, CHEN T, et al. Oligosaccharides from Morinda officinalis slow the progress of aging mice by regulating the key microbiota-metabolite pairs[J]. Evidence-Based Complementary and Alternative Medicine, 2019, 2019: 9306834. doi: 10.1155/2019/9306834.
[5] ZHU J, PENG Q, XU Y, et al. Morinda officinalis oligosaccharides ameliorate depressive-like behaviors in poststroke rats through upregulating GLUT3 to improve synaptic activity[J]. FASEB Journal, 2020, 34(10): 13376-13395. doi: 10.1096/fj.201902546RR
[6] 葛阳, 康传志, 万修福, 等. 生产中氮肥施用及其对中药材产量和质量的影响[J]. 中国中药杂志, 2021, 46(8): 1883-1892. [7] 刘哲, 钟海蓉, 威则日沙, 等. 氮、磷、钾配施对华重楼根茎的产量和有效成分含量的影响[J]. 中草药, 2019, 50(24): 6103-6113. [8] 邹庆军, 汪涛, 郭巧生, 等. N、P、K调控对野菊药用部位性状及活性成分含量的影响[J]. 中草药, 2019, 50(2): 535-540. doi: 10.7501/j.issn.0253-2670.2019.02.036 [9] 苏文华, 张光飞, 王泽明, 等. 氮、磷和钾肥对灯盏花生长和有效成分积累的影响[J]. 中草药, 2009, 40(12): 1963-1966. [10] 苏现明, 王洪庆, 陈若芸, 等. 巴戟天属植物化学成分及药理活性研究进展[J]. 中药材, 2017, 40(4): 986-991. [11] LIANG J, LIANG J, HAO H, et al. The extracts of Morinda officinalis and its hairy roots attenuate dextran sodium sulfate-induced chronic ulcerative colitis in mice by regulating inflammation and lymphocyte apoptosis[J]. Frontiers in Immunology, 2017, 8: 905. doi: 10.3389/fimmu.2017.00905.
[12] 代文豪. 发酵法炮制巴戟天工艺优化及其活性成分变化研究[D]. 广州: 华南农业大学, 2018. [13] 李倩. 巴戟天炮制过程中化学成分变化规律的研究[D]. 广州: 广州中医药大学, 2015. [14] 朱露林, 沈燚, 易晓薇, 等. 不同产地及种植年限巴戟天中6种寡糖类成分的定量分析[J]. 中草药, 2021, 52(14): 4348-4352. [15] 冼丽铧, 吴道铭, 隆曼迪, 等. 土壤养分、种植年限和种植方式对巴戟天寡糖含量的影响[J]. 华南农业大学学报, 2021, 42(3): 75-85. [16] 邵玲, 李美映, 黎学荣, 等. 肇庆巴戟天两种种植模式的植物学性状及药效成分研究[J]. 中药材, 2019, 42(11): 2480-2485. [17] 隆曼迪. 种植环境及栽培管理措施对巴戟天生长和品质的影响[D]. 广州: 华南农业大学, 2019. [18] 吴礼树. 土壤肥料学[M]. 北京: 中国农业出版社, 2011. [19] 邵镪钎, 李丹, 蒋攀, 等. 氮、磷、钾配施对川明参产量和品质的影响[J]. 中草药, 2018, 49(16): 3926-3932. doi: 10.7501/j.issn.0253-2670.2018.16.029 [20] 梁琴, 陈兴福, 李瑶, 等. 化肥与有机肥配施对川芎产量的影响[J]. 中药材, 2015, 38(10): 2015-2020. [21] 陈晓玉, 贺超, 闫滨滨, 等. N、P、K配施对丹参产量和质量的影响[J]. 中草药, 2019, 50(3): 722-730. [22] 范巧佳, 刘灵, 郑顺林, 等. 春季施氮时期和数量对川芎产量、生物碱和阿魏酸的影响[J]. 四川农业大学学报, 2013, 31(2): 136-139. [23] 贾袭伟, 晋小军, 李慧, 等. 磷酸二氢钾施肥量对款冬生长及产量、质量的影响[J]. 中成药, 2022, 44(1): 314-317. [24] DIMKPA C O, FUGICE J, SINGH U, et al. Development of fertilizers for enhanced nitrogen use efficiency: Trends and perspectives[J]. Science of the Total Environment, 2020, 731: 139113. doi: 10.1016/j.scitotenv.2020.139113.
[25] 洪丕征, 刘世荣, 王晖, 等. 正常和环割条件下不同形态氮素添加对红椎幼苗光合特性的影响[J]. 生态学报, 2016, 36(14): 4485-4495. [26] 潘超美, 詹若挺, 丁平, 等. 不同有机肥对巴戟天生长及矿质养分吸收的影响[J]. 中药材, 2002, 25(10): 699-701. doi: 10.3321/j.issn:1001-4454.2002.10.002 [27] 漆小雪, 韦霄, 陈宗游, 等. 黄花蒿干物质的积累及青蒿素与N、P、K量的动态变化研究[J]. 中草药, 2011, 42(12): 2541-2544. [28] 董娟娥, 张康健, 梁宗锁. 植物次生代谢与调控[M]. 杨凌: 西北农林科技大学出版社, 2009. [29] 冯嘉仪, 谢姗宴, 吴道铭, 等. 氮磷钾配施对银杏果实和外种皮产量及品质的影响[J]. 生态学杂志, 2021, 40(6): 1650-1659. [30] 胡佳栋, 毛歌, 张志伟, 等. 不同施肥处理对党参产量和次生代谢物含量的影响研究[J]. 中国中药杂志, 2017, 42(15): 2946-2953. [31] 卢挺, 杨全, 唐晓敏, 等. 氮磷钾配比施肥对广金钱草产量及质量的影响[J]. 广西植物, 2014, 34(3): 426-430. [32] 吴向维. 基于巴戟天寡糖的喘可治注射剂质量评价研究[D]. 广州: 广州中医药大学, 2014. [33] 邓秋林, 杨正明, 陈雨, 等. 氮磷钾配施对瓦布贝母产量及总生物碱质量分数的影响[J]. 西北农业学报, 2019, 28(7): 1138-1146. [34] 鲁泽刚, 朱永全, 卢迎春, 等. 氮、磷、钾施用对灯盏花产量和主要提取物收获量的影响[J]. 植物科学学报, 2019, 37(1): 55-62. -
期刊类型引用(4)
1. 郭毫单，吴影，韦玉琼，曹力，白周亚，樊秋霞，彭楠，古绍彬. 凝结魏茨曼氏菌BC99应急能量棒研制及抗疲劳效应研究. 食品工业科技. 2025(01): 218-230 . 
百度学术
2. 陈倩玲，沙玉柱，刘秀，邵鹏阳，王翻兄，陈小伟，杨文鑫，谢转回，高敏，黄薇. 肠道微生物与线粒体互作调控动物脂肪沉积的研究进展. 畜牧兽医学报. 2024(06): 2293-2303 . 百度学术
3. 陶舒悦，梁万徽，汪杰，张云静，方玲，彭灿. 心通口服液治疗急性心肌梗死的代谢组学研究. 安徽中医药大学学报. 2023(06): 73-79 . 百度学术
4. 孟科，赵薇，郭晨浩，聂伟，陶毛孩，袁晓春，孙昊然，冯登侦. 不同品种绵羊肌内脂肪沉积相关miRNA的筛选与功能预测. 江苏农业学报. 2023(07): 1554-1566 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 175
- HTML全文浏览量: 27
- PDF下载量: 17
- 被引次数: 8
