1 材料与方法
1.1 试验材料
1.2 试验设计
表1 饲料组成及营养水平(干物质基础)Table 1 Composition and nutrient levels of diets (DM basis) |
| 项目 Items | 芦根粉添加水平 Rhizoma phragmitis powder supplemental levels/% | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| 0(对照Control) | 0.5 | 1.0 | 3.0 | ||||||
| 原料Ingredients | |||||||||
| 鱼粉Fish meal | 5.00 | 5.00 | 5.00 | 5.00 | |||||
| 豆粕Soybean meal | 23.00 | 23.00 | 23.00 | 23.00 | |||||
| 菜籽粕Rapeseed meal | 16.00 | 16.00 | 16.00 | 16.00 | |||||
| 棉籽粕Cottonseed meal | 10.00 | 10.00 | 10.00 | 10.00 | |||||
| 小麦Wheat | 16.00 | 16.00 | 16.00 | 16.00 | |||||
| 米糠粕Rice bran meal | 16.00 | 16.00 | 16.00 | 16.00 | |||||
| 膨润土Bentonite clay | 2.00 | 2.00 | 2.00 | 2.00 | |||||
| 微晶纤维素Microcrystalline cellulose | 3.00 | 2.50 | 2.00 | ||||||
| 芦根粉Rhizoma phragmitis powder | 0.50 | 1.00 | 3.00 | ||||||
| 大豆油Soybean oil | 3.60 | 3.60 | 3.60 | 3.60 | |||||
| 氯化钠NaCl | 0.40 | 0.40 | 0.40 | 0.40 | |||||
| 磷酸二氢钙Ca(H2PO4)2 | 2.00 | 2.00 | 2.00 | 2.00 | |||||
| 羧甲基纤维素Carboxymethyl cellulose | 1.85 | 1.85 | 1.85 | 1.85 | |||||
| 氯化胆碱Choline chloride | 0.15 | 0.15 | 0.15 | 0.15 | |||||
| 预混料Premix1) | 1.00 | 1.00 | 1.00 | 1.00 | |||||
| 合计Total | 100.00 | 100.00 | 100.00 | 100.00 | |||||
| 营养水平Nutrient levels2) | |||||||||
| 粗蛋白质CP | 33.75 | 33.29 | 33.58 | 33.56 | |||||
| 粗脂肪EE | 5.92 | 5.96 | 5.87 | 5.81 | |||||
| 粗灰分Ash | 10.35 | 10.34 | 10.27 | 10.40 | |||||
| 总能GE/(MJ/kg) | 18.74 | 18.67 | 18.73 | 18.72 | |||||
1)预混料为每千克饲料提供 The premix provided the following per kg of diets:VA 4 797 IU,VD 1 000 IU,VE 50 IU,VK 10 mg,VB1 20 mg,VB2 20 mg,VB6 20 mg,VB12 0.02 mg,叶酸 folic acid 5 mg,泛酸 pantothenic acid 50 mg,肌醇 inositol 100 mg,VB3 100 mg,生物素 biotin 0.1 mg,VC 100 mg,KCl 200 mg,KI (1%) 60 mg,CoCl2·6H2O (1%) 50 mg,CuSO4·5H2O 30 mg,FeSO4·H2O 400 mg,ZnSO4·H2O 400 mg,MnSO4·H2O 150 mg,Na2SeO3·5H2O (1%) 65 mg,MgSO4·H2O 2 000 mg,沸石粉 zeolite powder 3 700.2 mg,小麦粉 wheat meal 2 468 mg。 2)营养水平为根据AOAC(1995)实测值。Nutrient levels were measured values based on AOAC (1995). |
1.3 饲养管理
1.4 样品采集
1.5 测定指标
1.5.1 生长性能
1.5.2 血浆生化指标
1.5.3 肝脏抗氧化指标
1.5.4 肝脏和脾脏组织结构
1.5.5 肝脏转录组分析
1.5.6 实时荧光定量PCR(qRT-PCR)验证转录组测序结果
表2 实时荧光定量PCR引物序列Table 2 Primer sequences for qRT-PCR |
| 基因 Genes | 引物序列 Primer sequences (5'—3') | GenBank登录号 GenBank accession numbers |
|---|---|---|
| 核因子-κB抑制因子α NFKBIA | F:CTGGTGGATGACTGTGGT R:TGCTCCTGTGCGTTTAT | KU554445.1 |
| 肿瘤坏死因子受体超家族成员14 TNFRSF14 | F:CAGTGCCGATTTGAG R:TCTGATGTCTTTGTCCC | XM_048209044.1 |
| 干扰素α/β受体 IFNAR | F:TCGGCGTTTCAATCC R:GCTGCTCAGAGGGTCA | XM_048192472.1 |
| 补体9 C9 | F:GGCGGGACAATCTCA R:TTCCGTTTGGTACACCTC | XM_048188403.1 |
| 丝裂原活化蛋白激酶14 MAPK14 | F:TTCTCGCCCGCTACCA R:TCCCGTCATCTCATCGTCAG | XM_048209822.1 |
| 核因子-κB NF-κB | F:AACAGAGCCGAGGGTA R:AGATGTCCGCCATAGAA | XM_048204723.1 |
| 三联体结构域蛋白23 TRIM23 | F:AATGCGGTGAAGGTAAAG R:ACGGGTCAAGCAGTCG | XM_048174744.1 |
| 含整合素αFG-GAP重复蛋白1 ITFG1 | F:CAGAGTGAAAGGGTGGA R:TTGATTCGGTGGTGAG | XM_048184707.1 |
| 干扰素调节因子1 IRF1 | F:TTCAGGGAGGGATTTACCA R:CGCCGCAGCCTTTGTT | XM_048166669.1 |
| NOD样受体家族CARD结构域蛋白3 NLRC3 | F:CTCTGCCGTTTCGTG R:TGTAAATCTGCCTCCCT | XM_048155547.1 |
| 白细胞介素-10受体 IL-10R | F:GAAGCCTCCCAACG R:AACCATTTCAGCCTTTA | XM_048189105.1 |
| 热休克蛋白70 HSP70 | F:CGCAGTTATCACAGTTCCA R:AAAGTCCTCGCCACCC | XM_048208288.1 |
| β-肌动蛋白 β-actin | F:CGTGCTGTTTTCCCTTCCATT R:CAATACCGTGCTCAAAGGATACTT | AY170122.2 |
1.6 数据统计分析
2 结果与分析
2.1 饲料中添加芦根粉对团头鲂生长性能的影响
表3 饲料中添加芦根粉对团头鲂生长性能的影响Table 3 Effects of dietary supplementation with Rhizoma phragmitis powder on growth performance of Megalobrama amblycephala |
| 项目 Items | 芦根粉添加水平 Rhizoma phragmitis powder supplemental levels/% | P值 P-value | |||
|---|---|---|---|---|---|
| 0(对照Control) | 0.5 | 1.0 | 3.0 | ||
| 增重率Weight gain rate/% | 74.18±0.11 | 73.08±8.36 | 71.77±8.18 | 68.22±10.81 | 0.886 |
| 摄食量Feed intake/(g/尾) | 93.37±0.44 | 93.61±0.10 | 93.45±0.27 | 93.52±0.16 | 0.804 |
2.2 饲料中添加芦根粉对团头鲂血浆生化指标的影响
表4 饲料中添加芦根粉对团头鲂血浆生化指标的影响Table 4 Effects of dietary supplementation with Rhizoma phragmitis powder on plasma biochemical indices of Megalobrama amblycephala |
| 项目 Items | 芦根粉添加水平 Rhizoma phragmitis powder supplemental levels/% | P值 P-value | |||
|---|---|---|---|---|---|
| 0(对照Control) | 0.5 | 1.0 | 3.0 | ||
| 丙氨酸氨基转移酶ALT/(U/L) | 60.37±7.75a | 33.07±9.91b | 36.30±2.61b | 35.43±9.68b | 0.010 |
| 天冬氨酸氨基转移酶AST/(U/L) | 290.50±13.16a | 171.87±28.25b | 171.10±21.96b | 162.73±16.85b | <0.001 |
| 总蛋白TP/(g/L) | 42.97±6.94a | 36.37±2.26ab | 42.07±0.12a | 31.30±2.66b | 0.020 |
| 白蛋白ALB/(g/L) | 16.43±2.89 | 15.87±2.07 | 17.40±0.50 | 12.93±1.20 | 0.090 |
| 球蛋白GLB/(g/L) | 26.53±4.24a | 20.50±2.29bc | 24.67±0.60ab | 18.37±1.48c | 0.015 |
| 尿素氮UN/(mmol/L) | 1.05±0.22 | 0.91±0.18 | 0.87±0.07 | 0.69±0.07 | 0.070 |
| 尿酸UA/(μmol/L) | 71.13±21.67 | 93.30±66.42 | 70.93±37.84 | 147.07±24.11 | 0.161 |
| 肌酐CREA/(μmol/L) | 3.53±1.24 | 2.60±0.26 | 4.03±0.06 | 3.80±1.13 | 0.255 |
2.3 饲料中添加芦根粉对团头鲂肝脏抗氧化指标的影响
图1 饲料中添加芦根粉对团头鲂肝脏抗氧化指标的影响*表示与对照组相比差异显著(P<0.05),**表示与对照组相比差异极显著(P<0.01)。 Fig.1 Effects of dietary supplementation with Rhizoma phragmitis powder on liver antioxidant indices of Megalobrama amblycephala (n=3) * indicated significant difference compared with the control group (P<0.05), and ** indicated extremely significant difference compared with the control group (P<0.01). |
2.4 饲料中添加芦根粉对团头鲂肝脏和脾脏组织结构的影响
图2 饲料中添加芦根粉对团头鲂肝脏和脾脏组织结构的影响A和E为对照组肝脏,B和F为0.5%芦根粉添加组肝脏,C和G为1.0%芦根粉添加组肝脏,D和H为3.0%芦根粉添加组肝脏;I和M为对照组脾脏,J和N为0.5%芦根粉添加组脾脏,K和O为1.0%芦根粉添加组脾脏,L和P为3.0%芦根粉添加组脾脏。△表示淋巴细胞,⇨表示脂肪细胞,☆表示肝细胞空泡化,➡表示含铁血黄素沉积。 Fig.2 Effects of dietary supplementation with Rhizoma phragmitis powder on tissue structure of liver and spleen of Megalobrama amblycephala (400×) A and E were liver in the control group, B and F were liver in 0.5% Rhizoma phragmitis powder supplemental group, C and G were liver in 1.0% Rhizoma phragmitis powder supplemental group, and D and H were liver in 3.0% Rhizoma phragmitis powder supplemental group. I and M were spleen in the control group, J and N were spleen in 0.5% Rhizoma phragmitis powder supplemental group, K and O were spleen in 1.0% Rhizoma phragmitis powder supplemental group, and L and P were spleen in 3.0% Rhizoma phragmitis powder supplemental group. △ indicated lymphocyte cells, ⇨ indicated lipocytes, ☆ indicated hepatocyte vacuolations, and ➡ indicated hemosiderin deposition. |
2.5 饲料中添加芦根粉对团头鲂肝脏转录组的影响
2.5.1 转录组测序数据统计
表5 团头鲂肝脏转录组测序数据统计表Table 5 Statistical table of liver transcriptome sequencing data of Megalobrama amblycephala |
| 芦根粉添加水平 Rhizoma phragmitis powder supplemental levels/% | 净读数 Clean reads | GC含量 GC content/% | Q30/% | 比对读数 Mapped reads | 比对率 Mapped ratio/% | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
0(对照Control) | 58 806 774 | 46.39 | 95.28 | 53 825 045 | 91.53 | ||||||
| 42 065 638 | 45.67 | 95.47 | 38 654 725 | 91.89 | |||||||
| 41 536 680 | 46.68 | 96.10 | 38 907 421 | 93.67 | |||||||
0.5 | 44 349 524 | 45.38 | 94.97 | 40 612 604 | 91.57 | ||||||
| 41 078 448 | 45.20 | 94.59 | 37 404 598 | 91.06 | |||||||
| 44 442 736 | 45.08 | 95.01 | 40 212 493 | 90.48 | |||||||
1.0 | 53 450 278 | 46.25 | 93.93 | 48 308 990 | 90.38 | ||||||
| 57 460 816 | 45.98 | 94.66 | 52 291 567 | 91.00 | |||||||
| 49 421 026 | 46.10 | 95.06 | 45 131 786 | 91.32 | |||||||
3.0 | 38 636 264 | 47.06 | 95.77 | 36 224 666 | 93.76 | ||||||
| 42 793 226 | 45.65 | 94.61 | 39 257 040 | 91.74 | |||||||
| 42 523 356 | 45.78 | 95.09 | 38 523 714 | 90.59 |
2.5.2 差异表达基因(differentially expressed genes,DEGs)鉴定
2.5.3 DEGs GO功能和KEGG通路富集分析
图4 团头鲂肝脏转录组差异表达基因GO功能富集分析柱状图A、B和C分别表示0.5%、1.0%和3.0%芦根粉添加组与对照组进行比较。A, B and C represented 0.5%, 1.0% and 3.0% Rhizoma phragmitis powder supplemental groups compared with the control group, respectively. GO_term:GO条目;up:上调 up-regulated;down:下调 down-regulated;Num of Genes:基因数量;biological process:生物过程;cellular component:细胞组分;molecular function:分子功能;cellular process:细胞过程;metabolic process:代谢过程;biological regulation:生物调节;response to stimulus:应激反应;developmental process:发育过程;signaling:信号传递;multicellular organismal process:多细胞生物过程;localization:定位;multi-organism process:多生物体过程;immune system process:免疫系统过程;reproduction:生殖;reproductive process:生殖过程;interspecies interaction between organisms:生物种间相互作用;biological adhesion:生物黏附;locomotion:运动;growth:生长;rhythmic process:节律过程;behavior:行为;detoxification:解毒作用;pigmentation:色素沉着;biomineralization:生物矿化;cellular anatomical entity:细胞解剖实体;intracellular:细胞内;protein-containing complex:含蛋白质复合物;binding:结合;catalytic activity:催化活性;transcription regulator activity:转录调节活性;molecular function regulator:分子功能调节因子;transporter activity:转运体活性;structural molecule activity:结构分子活性;molecular transducer activity:分子转导活性;antioxidant activity:抗氧化活性;translation regulator activity:翻译调节活性;molecular carrier activity:分子载体活性;cargo receptor activity:货物受体活性;protein tag:蛋白质标签。 Fig.4 Histograms of GO functional enrichment analysis of DEGs in liver transcriptome of Megalobrama amblycephala |
图5 团头鲂肝脏转录组差异表达基因GO功能富集分析气泡图A、B和C分别表示0.5%、1.0%和3.0%芦根粉添加组与对照组进行比较。A, B and C represented 0.5%, 1.0% and 3.0% Rhizoma phragmitis powder supplemental groups compared with the control group, respectively. Cellular Component:细胞组分;Molecular Function:分子功能;Biological Process:生物过程;qvalue:Q值;gene_number:基因数量;Diff:差异性 difference;down:下调 down-regulated;up:上调 up-regulated;up&down:上调和下调 up-regulated and down-regulated;Rich factor:富集因子。 限于篇幅,仅注释主要的条目名称 due to space limitation, only main term names were annotated。intracellular membrane-bounded organelle:细胞内膜结合细胞器:cytoplasm:细胞质;membrane-bounded organelle:膜结合细胞器;nucleus:细胞核;DNA-binding transcription factor activity:DNA结合转录因子活性;sequence-specific DNA binding:序列特异性DNA结合。 Fig.5 Bubble diagrams of GO functional enrichment analysis of DEGs in liver transcriptome of Megalobrama amblycephala |
图6 团头鲂肝脏转录组差异表达基因KEGG通路富集分析A、C和E分别表示0.5%、1.0%和3.0%芦根粉添加组与对照组进行比较的柱状图,B、D和F分别表示0.5%、1.0%和3.0%芦根粉添加组与对照组进行比较的气泡图。A, C and E represented histograms of 0.5%, 1.0% and 3.0% Rhizoma phragmitis powder supplemental groups compared with the control group, and B, D and F represented bubble diagrams of 0.5%, 1.0% and 3.0% Rhizoma phragmitis powder supplemental groups compared with the control group, respectively. Cellular Processes:细胞过程;Environmental Information Processing:环境信息处理;Genetic Information Processing:遗传信息处理;Human Diseases:人类疾病;Metabolism:代谢;Organismal Systems:生物体系统;Annotated Genes:注释基因;Statistics of Pathway Enrichment:通路富集统计;qvalue:Q值;gene_number:基因数量;Diff:差异性 difference;down:下调 down-regulated;up:上调 up-regulated;up&down:上调和下调 up-regulated and down-regulated;Rich factor:富集因子。 限于篇幅,仅注释主要的通路名称 due to space limitation, only main pathway names were annotated。Apoptosis:细胞凋亡;Apelin signaling pathway:Apelin信号通路;Herpes simplex virus l infection:单纯疱疹病毒1型感染;Oxidative phosphorylation:氧化磷酸化;Insulin signaling pathway:胰岛素信号通路;MAPK signaling pathway:丝裂原活化蛋白激酶信号通路;Salmonella infection:沙门氏菌感染;Arachidonic acid metabolism:花生四烯酸代谢;TGF-beta signaling pathway:转化生长因子-β信号通路;Autophagy-animal:细胞自噬-动物;mTOR signaling pathway:哺乳动物雷帕霉素靶蛋白信号通路;adipocytokine signaling pathway:脂肪细胞因子信号通路。 Fig.6 KEGG pathway enrichment analysis of DEGs in liver transcriptome of Megalobrama amblycephala |
表6 免疫相关差异表达基因KEGG通路富集分析Table 6 KEGG pathway enrichment analysis of DEGs related to immunity |
| 通路编号 Pathway ID | KEGG通路名称 KEGG pathway names | A | B | C | |||
|---|---|---|---|---|---|---|---|
| up&down | up&down | up&down | |||||
| ko04620 | Toll样受体信号通路Toll-like receptor signaling pathway | 5 | 5 | 4 | |||
| ko04625 | C型凝集素受体信号通路C-type lectin receptor signaling pathway | 7 | 8 | 7 | |||
| ko04622 | 视黄酸诱导基因-Ⅰ样受体通路RIG-Ⅰ-like receptor signaling pathway | 2 | 5 | 3 | |||
| ko04621 | NOD样受体信号通路NOD-like receptor signaling pathway | 7 | 6 | 5 | |||
| ko04010 | 丝裂原活化蛋白激酶信号通路MAPK signaling pathway | 8 | 13 | 10 | |||
| ko04140 | 细胞自噬Autophagy | 8 | 4 | 10 | |||
| ko04142 | 溶酶体Lysosome | 4 | 4 | 3 | |||
| ko04210 | 细胞凋亡Apoptosis | 15 | 12 | 8 | |||
A、B和C分别表示0.5%、1.0%和3.0%芦根粉添加组与对照组进行比较。 up&down:上调和下调 up-regulated and down-regulated。 | |
A, B and C represented 0.5%, 1.0% and 3.0% Rhizoma phragmitis powder supplemental groups compared with the control group, respectively. |
