1 材料与方法
1.1 试验设计
1.2 指标测定
1.2.1 生长性能指标测定
1.2.2 肝脏生理生化指标检测
1.2.3 肝脏组织切片分析
1.2.4 转录组测序及基因功能注释
1.2.5 转录组DEGs筛选
1.2.6 实时荧光定量PCR(RT-qPCR)分析
表1 本研究所用的引物序列Table 1 Primer sequences used in this study |
| 基因 Genes | 上游引物 Forward primer (5'—3') | 下游引物 Reverse primer (5'—3') | 产物大小 Product size/bp |
|---|---|---|---|
| EF1α | TGGGTGTTGGACAAACTGAA | CAACACCACCAGCAACAATC | 190 |
| LPL | AGTGTCAGTGAACCGTCGTAGC | CGCAGTTTCCTCCTAATCTCCT | 159 |
| LPIN1 | GCTCCTGAGCTGCACATAGAC | AGGTGGCGGCACTTTACTT | 100 |
| CYP24A1 | GTGAGGTGGCTCTGGTGGT | CCTTTGGAACAATGATGGTGAC | 193 |
| CYP2R1 | TTAGGTAAGCGGCACTGTCTG | GAATAAAGCCTTCGGGGAAC | 109 |
| POLA1 | CAACCAGTAAAAGCGGAAAAC | GAAGATGGAGGAAGGGGAAC | 115 |
| LIG1 | GGCTAGAGGCAGTGCAAGTG | TCGTGGGAAACGGAGTGA | 119 |
| PCNA | TTCTGATATGCAATGAAATCCC | ACCAAGGCTACTCCTCTGTCC | 185 |
| FOXM1 | TGTTTTCTTGCGGTCTGGTT | GAGTGACGTGAGCTTGGTGAT | 128 |
| AHR1 | GGGCAACCAGATTGCATCA | GCAACCACAGCCTAAAACAGA | 168 |
| KLF9 | CGTCTGTATCCCGTTGCG | TTGACTTTCCTCAGGTCCCTC | 199 |
| TGFβ1 | GCAAGACTCTGGACCTGGAG | TGCAGAGGGACAGTGTTCAG | 202 |
| TNFα | GCAGATCGCTGAGAACACAA | AGCTTCATTGTCGCCTTCAT | 226 |
| IL1β | CCACCTGGTCTTGAAGGAGA | CACCTGCTCGAAATCCTCAT | 172 |
EF1α:延长因子1α elongation factor 1α;LPL:脂蛋白脂肪酶 lipoprotein lipase;LPIN1:脂素1 lipin 1;CYP24A1:1α,25-二羟维生素D3-24-羟化酶 1α,25-dihydroxyvitamin D3-24-hydroxylase;CYP2R1:维生素D-25-羟化酶vitamin D-25-hydroxylase;POLA1:DNA聚合酶α催化亚基 DNA polymerase α catalytic subunit;LIG1:DNA连接酶1 DNA ligase 1;PCNA:增殖细胞核抗原 proliferating cell nuclear antigen;FOXM1:叉头蛋白M1 forkhead box protein M1;AHR1:芳香烃受体1 aryl hydrocarbon receptor 1;KLF9:Krueppel样因子9 Krueppel-like factor 9;TGFβ1:转化生长因子β1 transforming growth factor beta 1;TNFα:肿瘤坏死因子α tumor necrosis factor α;IL1β:白细胞介素1β interleukin 1 beta。 |
1.3 数据统计与分析
2 结果与分析
2.1 高投饲率对圆口铜鱼幼鱼生长性能的影响
表2 高投饲率对圆口铜鱼幼鱼生长性能的影响Table 2 Effects of high feeding rate on growth performance of juvenile largemouth bronze gudgeon |
| 项目 Items | 初始体重 IBW/g | 终末体重 FBW/g | 增重率 WGR/% | 特定生长率 SGR/(%/d) | 肥满度 CF/(g/cm3) | 存活率 SR/% |
|---|---|---|---|---|---|---|
| 对照组 Control group | 4.93±0.16 | 9.11±0.21a | 83.70±3.23a | 1.09±0.03a | 1.42±0.10a | 100±0 |
| 高投饲率组 High feeding rate group | 4.97±0.13 | 12.37±0.15b | 149.40±4.65b | 1.63±0.02b | 1.55±0.06b | 100±0 |
同列数据肩标无字母或相同字母者表示差异不显著(P>0.05),不同小写字母表示差异显著(P<0.05)。 | |
In the same column, values with no letter or the same letter superscripts mean no significant difference (P>0.05), while with different small letter superscripts mean significant difference (P<0.05). |
2.2 高投饲率对圆口铜鱼幼鱼肝脏生理生化指标的影响
表3 高投饲率对圆口铜鱼幼鱼肝脏生理生化指标的影响Table 3 Effects of high feeding rate on hepatic physiological and biochemical indicators of juvenile largemouth bronze gudgeon |
| 项目 Items | 对照组 Control group | 高投饲率组 High feeding rate group |
|---|---|---|
| 总蛋白 TP/(g/L) | 48.48±3.34 | 47.80±4.89 |
| 甘油三酯 TG/(μmol/g prot) | 200.19±19.48a | 448.10±47.77b |
| 总胆固醇 TC/(μmol/g prot) | 40.97±4.46a | 69.95±7.35b |
| 低密度脂蛋白胆固醇 LDL-C/(μmol/g prot) | 12.14±1.17 | 12.42±1.00 |
| 游离脂肪酸 NEFA/(μmol/g prot) | 68.98±9.63a | 133.37±20.40b |
| 谷丙转氨酶 ALT/(U/g prot) | 23.29±1.86a | 30.29±2.26b |
| 谷草转氨酶 AST/(U/g prot) | 33.79±3.13a | 45.09±4.38b |
| 碱性磷酸酶 AKP/(mg/g prot) | 24.63±2.51a | 31.45±2.65b |
| 总胆汁酸 TBA/(μmol/g prot) | 13.04±1.07b | 8.06±0.56a |
同行数据肩标无字母或相同字母者表示差异不显著(P>0.05),不同小写字母表示差异显著(P<0.05)。下表同。 | |
In the same row, values with no letter or the same letter superscripts mean no significant difference (P>0.05), while with different letter superscripts mean significant difference (P<0.05). The same as below. |
2.3 高投饲率对圆口铜鱼幼鱼肝脏抗氧化和免疫功能的影响
表4 高投饲率对圆口铜鱼幼鱼肝脏抗氧化和免疫指标的影响Table 4 Effects of high feeding rate on hepatic antioxidant and immune indicators of juvenile largemouth bronze gudgeon |
| 项目 Items | 对照组 Control group | 高投饲率组 High feeding rate group |
|---|---|---|
| 抗氧化指标 Antioxidant indicators | ||
| 活性氧 ROS/(U/mg prot) | 69.02±6.00a | 161.31±14.82b |
| 总抗氧化能力 T-AOC/(μmol/g prot) | 104.22±3.45a | 182.65±16.44b |
| 过氧化氢酶 CAT/(U/g prot) | 137.16±16.26b | 82.84±8.86a |
| 还原性谷胱甘肽 GSH/(μmol/g prot) | 12.14±1.12a | 27.30±3.63b |
| 超氧化物歧化酶 SOD/(U/mg prot) | 202.66±9.42b | 165.00±15.36a |
| 丙二醛 MDA/(μmol/g prot) | 11.09±9.63a | 29.59±2.94c |
| 免疫指标 Immune indicators | ||
| 免疫球蛋白G IgG/(mg/g prot) | 4.03±0.26a | 8.65±0.78b |
| 免疫球蛋白M IgM/(mg/g prot) | 0.55±0.04a | 0.85±0.05b |
| 免疫球蛋白Z IgZ/(mg/g prot) | 0.22±0.01 | 0.25±0.02 |
| 溶菌酶 LYS/(U/g prot) | 53.99±3.30a | 89.88±6.66b |
图1 高投饲率对圆口铜鱼幼鱼肝脏炎症因子基因表达的影响数据柱标注不同小写字母表示组间有差异显著(P<0.05)。 Fig.1 Effects of high feeding rate on expression of hepatic inflammatory cytokine genes of juvenile largemouth bronze gudgeon Data columns with different lowercase letters in the figure indicate significant difference between groups (P<0.05). |
2.4 圆口铜鱼幼鱼肝脏组织病理分析
图2 圆口铜鱼幼鱼肝脏组织病理分析A:肝脏组织病理2种表型(比例尺=50 μm),黑色箭头表示脂滴(LD);B:油红O染色脂滴相对面积分析,图中数据柱标注不同小写字母表示组间有差异显著(P<0.05);C:脂肪肝表型数量的统计。 Fig.2 Hepatic histopathological analysis of juvenile largemouth bronze gudgeon A: two phenotypes of hepatic histopathology (scale bar=50 μm), and the black arrow represents lipid drops (LD); B: relative area analysis of lipid drops with oil red O staining, data columns with different lowercase letters in the figure indicate significant difference between groups (P<0.05); C: statistics of the number of fatty liver phenotype. |
2.5 转录组DEGs KEGG通路富集分析
2.6 脂质代谢DEGs及差异表达转录因子分析
图5 脂质代谢相关差异表达基因的分析NL表示无明显异常肝脏,FL表示脂肪肝。NL represents no abnormality liver, and FL represents fatty liver. GPAT3:甘油-3-磷酸酰基转移酶3 glycerol-3-phosphate acyltransferase 3;LPIN1:脂素1 lipin 1;PLPP1:磷脂磷酸酶1 phospholipid phosphatase 1;FASN:脂肪酸合成酶 fatty acid synthase;ACC1:乙酰辅酶A羧化酶1 acetyl-CoA carboxylase 1;ACPM:线粒体酰基载体蛋白 mitochondrial acyl carrier protein; PNPLA2:含Patatin样磷脂酶域包含蛋白2 Patatin-like phospholipase domain-containing protein 2;MGLL:单甘油酯脂肪酶monotriglyceride lipase;LPL:脂蛋白脂肪酶 lipoprotein lipase;DGKB:二酰基甘油激酶β diacylglycerol kinase β;ACSL6:长链脂肪酸乙酰辅酶A连接酶6 long-chain fatty acid acetyl-CoA ligase 6。 Fig.5 Analysis of DEGs related to lipid metabolism |
表5 差异表达转录因子Table 5 Differential expression transcription factors |
| 缩写 Abbreviations | 基因名称 Gene names | log2(差异倍数) log2(FC) | 调控 Regulation | q值 q-value |
|---|---|---|---|---|
| 脂质代谢相关 Related to lipid metabolism | ||||
| AHR1 | 芳香烃受体1 Aryl hydrocarbon receptor 1 | 2.15 | 上调 | 9.80E-03 |
| NR0B2 | 核受体亚家族B0组成员2 Nuclear receptor subfamily 0 group B member 2 | 2.23 | 上调 | 3.49E-02 |
| NR1D1 | 核受体亚家族D1组成员1 Nuclear receptor subfamily 1 group D member 1 | -2.12 | 下调 | 2.35E-05 |
| KLF9 | Krueppel样因子9 Krueppel transcription factor 9 | -2.63 | 下调 | 8.76E-12 |
| KLF11 | Krueppel样因子11 Krueppel transcription factor 11 | -2.19 | 下调 | 2.01E-11 |
| 其他 Others | ||||
| DMRT2 | 双性和单抗3相关转录因子2 Double sex and mab3-related transcription factor 2 | 3.34 | 上调 | 4.71E-02 |
| FOXM1 | 叉头蛋白M1 Forkhead box protein M1 | -2.35 | 下调 | 3.56E-02 |
| GLIS1 | GLIS家族锌指1 GLIS family Zinc finger 1 | -4.01 | 下调 | 8.29E-04 |
| DNAJB11 | DNA J同源B家族11 DNA J homolog subfamily B member 11 | 2.14 | 上调 | 2.16E-05 |
| OTPB | 同源异型盒蛋白 Othopedia homeobox protein B | -3.51 | 下调 | 1.82E-07 |
| UHRF1 | 泛素蛋白连接酶E3 Ubiquitin-protein ligase E3 | 3.37 | 上调 | 5.62E-04 |
| SBSPON | 生长抑素b和血小板反应蛋白1型结构域含蛋白 Somatomedin-B and thrombospondin type-1 domain-containing protein | 2.34 | 上调 | 2.36E-02 |
| HLF | 肝性白血病因子 Hepatic leukemia factor | -2.56 | 下调 | 2.92E-17 |
| PPA | 无机焦磷酸酶 Inorganic pyrophosphatase | 2.46 | 上调 | 8.17E-08 |
| NFIL3 | 核因子白细胞介素3调节蛋白 Nuclear factor interleukin-3-regulated protein | 2.99 | 上调 | 8.56E-06 |
图6 脂质代谢相关转录因子与基因的相关性分析图中各基因和转录因子的中文和英文名称见图5和表5。图中*表示相关性显著(P<0.05)。 Fig.6 Correlation analysis of transcription factors and genes related to lipid metabolism The Chinese and English names of the genes and transcription factors in the figure are shown in Fig.5 and Table 5. * in the figure indicates significant correlation (P<0.05). |
2.7 转录组测序结果的qPCR验证
图7 10个差异表达基因的验证A:qPCR检测基因表达情况,图中各基因的中文和英文名称见表1。数据用平均值±标准误差表示。*表示差异显著(P<0.05),**表示差异极显著(P<0.01)。B:qPCR数据与转录组测序数据间的相关性。 Fig.7 Validation of 10 DEGs A: gene expression measured by qPCR, the Chinese and English names of the genes in the figure are shown in Table 1. Values are expressed as mean±SE. * mean significant difference (P<0.05), and ** mean extremely significant difference (P<0.01). B: the correlation between the qPCR data and RNA-seq data. |

京公网安备 11010802041926号