吉富罗非鱼胃肠道几丁质酶的克隆、组织分布和纯化

2021-07-27 09:48黄金凤冯礼燊李文笙孙彩云
广东农业科学 2021年6期
关键词:罗非鱼克隆鱼类

黄金凤,李 波,白 莹,冯礼燊,李文笙,孙彩云

(中山大学生命科学学院/有害生物控制与资源利用国家重点实验室/广东省水生经济动物良种繁育重点实验室/广东省重要经济鱼类健康养殖工程技术研究中心,广东 广州 510006)

【研究意义】几丁质是由β-1,4 糖苷键连接N-乙酰氨基葡萄糖的一种线性多糖,是自然界中唯一存在的阳离子型可食用纤维。几丁质贮量十分丰富,仅海洋生物每年所合成的几丁质已超过10 亿t,是自然界中贮量仅次于纤维素的第二大天然多糖[1]。几丁质多存在于真菌细胞壁、昆虫消化道和多数节肢动物的外骨骼中,因此研究关于降解几丁质的酶具有重要意义。几丁质酶(EC 3.2.1.14)是一种糖苷水解酶,催化N-乙酰氨基葡萄糖从几丁质上脱去,从而降解几丁质。几丁质酶根据其催化结构域的氨基酸序列相似性主要分为糖基水解酶(GH)家族18、19、23 和48,该酶在自然界分布广泛,在细菌到高等植物和动物等生物体中均有分布,植物的几丁质酶属于GH19,鱼类的几丁质酶均为GH18,GH23 和GH48 分别在细菌和真菌中分布[2-4]。几丁质酶具有多种生理作用,如促进甲壳类和昆虫类蜕皮,降解鱼类胃肠道里的虾蟹壳帮助消化,在植物里则参与抵抗含有几丁质的致病菌[4],在鱼皮肤粘液中抵抗寄生虫等[5]。

【前人研究进展】鱼类中的几丁质酶最早在海鲷的胃肠道中被发现[6]。随着基因组测序的快速发展,已研究获得斑马鱼、日本比目鱼等多种鱼类的几丁质酶序列[7]。肉食性海水鱼类的天然食谱中含有丰富的几丁质,如含有高结晶度α-几丁质的虾蟹和浮游动物,含有低结晶度β-几丁质的乌贼等[4,8]。研究表明,鱼类消化道组织具有几丁质降解活性,其中胃中的几丁质酶活性最高,降解几丁质寡糖的活性比其他组织高,如白姑鱼、日本鲭、斑鱾、琉球棘角鱼和梭鱼等。而降解单糖的β-N-乙酰氨基葡萄糖酶的活性在各组织中广泛分布,而非局限于消化道。不同的几丁质降解酶的分布很可能取决于该鱼类品种的饮食习惯和栖息地位置[4,9]。已有研究表明,从鱼胃分离获得的几丁质酶约38~62 ku,与植物几丁质酶大小(24.5~29 ku)有明显差异。在实验室条件下,鱼胃的丁质酶降解几丁质低聚糖的最适pH 值为2.0~5.0,降解几丁质多聚糖的最适pH 值为1.5~5.5。这表明鱼类胃中的几丁质酶适应胃组织的酸性环境来降解几丁质底物[4]。

【本研究切入点】迄今为止,有关罗非鱼几丁质酶的研究报道仅有1 篇,它揭示了罗非鱼血清、胃和肠中存在几丁质酶降解活性,能够降解4-甲基伞形酮-β-D-N,N-二乙酰壳三糖〔4MU-(GlcNAc)2〕或4-甲基伞形酮-β-D-N,N,N-三乙酰壳三糖〔4MU-(GlcNAc)3〕,其中血清中的几丁质酶具有最高的比活性[10]。但目前仍未见关于罗非鱼几丁质酶序列及其功能的报道。【拟解决的关键问题】本研究克隆了吉富罗非鱼3 种预测的酸性几丁质酶,并对克隆结果进行了进化树分析,同时检测几丁质酶在罗非鱼中的组织分布,以期获得相关功能信息。提取罗非鱼胃组织蛋白并纯化其中的几丁质酶,测定酶活,比较罗非鱼和其他鱼类几丁质酶的活性,为罗非鱼几丁质酶的深入研究奠定基础。

1 材料与方法

1.1 试验材料

本研究所用的吉富罗非鱼(70~100 g)约为30 尾,均为雄性,购自广东省国家罗非鱼良种场。试验前,供试罗非鱼在12 h ∶12 h 黑暗光照循环条件下,室温循环淡水中驯化1 周,每天9:30 和16:30进行饱食投喂。采样前24 h内对罗非鱼停食,在第2 天早上进行取样。

主要试剂:总RNA 提取试剂Trizol Reagent、反转录试剂盒M-MLV First Strand cDNA Synthesis Kit 购自英潍捷基(上海)贸易有限公司;聚合酶Blend taq plus 和连接试剂Ligation High,购自东洋纺(上海)生物科技有限公司;DNase Ⅰ 购自美国纽英伦生物技术公司;E.Z.N.A® Gel Extraction Kit、E.Z.N.A® Plasmid Extraction Kit,购自Omega Bio-Tek 公司;SYBR® Green Real-time PCR Master Mix 试剂盒、100 bp DNA Ladder,购自广州东盛生物科技有限公司;DEPC、X-gal 和IPTG 购自生工生物工程(上海)公司。引物由英潍捷基(上海)贸易有限公司合成。所用大肠杆菌DH5α 和pUCM-T载体购自英潍捷基(上海)贸易有限公司,感受态细菌则由本实验室制备。其他试剂均为国产分析纯。

主要仪器:凝胶成像分析系统和CFX Connect荧光定量PCR 仪购于伯乐生命医学产品有限公司,VeritiTM PCR 仪、高速冷冻离心机5810R 和移液枪购于德国艾本德股份有限公司,核酸蛋白测定仪NanoDrop2000 购于美国应用生物系统公司,三洋制冰机购于日本SANYO 公司,微波炉购于广东美的集团有限公司。

1.2 试验方法

1.2.1 罗非鱼tChit1a、tChi3和tChit基因cDNA 的克隆 根据NCBI 的尼罗罗非鱼基因组中3 种预测的酸性几丁质酶序列tChit1a(XM_039611991.1)、tChi3(XM_005464320.4)和tChit(XM_003459030.3),采用Primer 5 软件设计特异性引物(表1)。

用丁香酚将罗非鱼麻醉,断头处死,取出胃、肠用DEPC 水清洗,按照说明书用Trizol Reagent 进行总RNA 提取,用DNase Ⅰ去除基因组DNA,用反转录试剂盒反转录合成第一链cDNA,并以此为模板按照聚合酶Blend tap plus所要求的反应条件进行PCR 扩增,所用引物分别 为tChi1a-ORF-F/R、tChi3-ORF-F/R 和tChit-ORF-F/R(表1)。PCR 反应程序均为:94 ℃ 3 min;94 ℃ 30 s、56 ℃ 30 s、72 ℃ 45 s,40 个循环;72 ℃ 7 min;4 ℃ 保存。对PCR 产物进行胶回收,加入pUCM-T 载体和Ligation High 室温孵育连接,然后通过热激法转化至大肠杆菌DH5α 中,阳性克隆经菌落PCR 鉴定后,送至生工生物工程(上海)公司进行测序。

表 1 罗非鱼3 种几丁质酶tChit1a、tChi3 和tChit 的ORFs 克隆引物及组织分布研究引物Table 1 Primers used in cloning of ORFs and tissue distribution study of three chitinases tChit1a,tChi3 and tChit in tilapia

1.2.2 罗非鱼tChit1a、tChi3和tChit基因序列的生物信息学分析 利用SignalP 5.0 寻找几丁质酶tChit1a、tChi3和tChit基因的信号肽(http://www.cbs.dtu.dk/services/SignalP/)[11];利用ClustalW2软件进行氨基酸序列同源性分析(http://www.ebi.ac.uk/Tools/msa/clustalo/)[12]。通过NCBI 在线网站进行氨基酸序列相似度分析(https://blast.ncbi.nlm.nih.gov/Blast.cgi)[13];采 用MEGA6.0 软 件Neighbor-Joining 法进行进化树分析[14]。

1.2.3 罗非鱼tChit1a、tChi3和tChit基因的mRNA 组织分布 用丁香酚麻醉罗非鱼,断头处死,迅速取出端脑、中脑、小脑、延髓、脊髓、下丘、垂体、鳃、心脏、肾、脾、脂肪、肝、胃、前肠、中肠、后肠、白肌、红肌和性腺等组织,放入RNase-free 的1.5 mL 离心管中,立即放入液氮,之后置于-80 ℃冰箱中贮藏。按照1.2.1 克隆步骤进行总RNA 提取和反转录,以获得的第一链cDNA 为模板,按照SYBR® Green 说明书进行实时定量PCR,所用引物为tChit1a-qt-F/R、tChi3-qt-F/R 和tChit-qt-F/R(表1)。PCR 反应程序为:95 ℃ 5 min;95 ℃ 15 s、56 ℃ 15 s、72 ℃ 30 s,40 个循环;每个循环结束后在82 ℃收集荧光15 s;95 ℃ 10 s,然后以0.5 ℃/2 s 的递增速度,65~95 ℃下进行温度梯度退火及收集荧光值。所测的值以18S rRNA 为内参进行标准化,2-ΔΔCt法分析数据。

1.2.4 罗非鱼胃几丁质酶的纯化 用丁香酚麻醉罗非鱼,断头处死,取出其胃组织进行清洗,剪下内壁粘膜层置于PBS 缓冲液中破碎(70 Hz/ 120 s),随后12 000 g 离心20 min 后收集上清用于纯化。用几丁质树脂制备重力柱,根据说明书清洗柱子,随后将样品过几丁质柱2 次,缓冲液清洗杂蛋白,然后加入洗脱液,关闭出样口,4 ℃过夜。第2 天依次收集滤液,并进行蛋白浓度测定和Western blot 检测。

1.2.5 SDS-PAGE及Western blot验证 制备12% SDS 聚丙烯酰胺凝胶,上样后80 V 电泳30 min,120 V 电泳60 min。电泳结束,其中一块PAGE胶用于G250 染色,另一块PAGE 胶则按照正极-海绵-滤纸-NC 膜-PAGE 胶-滤纸-海绵-负极的顺序夹紧,放在电泳槽中,装置在冰浴中恒流200 mA 电泳60 min;电转完的NC 膜用TBST缓冲液漂洗后浸入5%脱脂奶粉中封闭30 min;用1 ∶4 000 的斜带石斑鱼1 型几丁质酶抗体稀释液孵育过夜,最后放入Goat anti-rabbit IgG按1∶5 000 稀释液中孵育1 h,最后滴加ECL 荧光混合液到膜上,在化学发光仪中显色并拍照记录。

1.2.6 酶活测定 采用4-甲基伞形酮(4-MU)荧光底物测定酶活。酶活定义为:在37 ℃下,每分钟释放1 μmol 4-MU 所需要的酶量,用U 表示。比活力用每克蛋白所含的酶活力单位(U/mg)数表示。1 μg 待测样品、不同pH 的PBS 缓冲液和底物0.5 mmol/L 的4-甲基伞形酮-β-D-N,N-二乙酰壳三糖〔4MU-(GlcNAc)2〕或4-甲基伞形酮-β-D-N,N,N-三乙酰壳三糖〔4MU-(GlcNAc)3〕,在37 ℃下孵育70 min,然后用 0.5 mol/L Na2CO3缓冲液终止反应。用pH 7.4 的PBS 缓冲液配制梯度浓度0、0.1、0.2、0.4、0.6、0.8、1 nmol 的4-MU 标准样品,而待测样品是1 μg 纯化产物与2.5 nmol 的4MU-(GlcNAc)2/3,标准样品和待测样品的反应终体积均为100 μL,均在37 ℃、500 r/min 反应70 min,最后在每孔加入 100 μL 的 0.5 mol/L 的Na2CO3溶液终止反应,用酶标仪在激发波长370 nm、发射波长448 nm 下的荧光强度测定。

2 结果与分析

2.1 罗非鱼3 种几丁质酶的cDNA 及其编码的氨基酸序列

从吉富罗非鱼胃中克隆得到tChit1a的ORF序列,全长1 362 bp,编码453 个氨基酸(AA);Expasy 在线软件预测蛋白大小为50.03 ku,等电点为6.14;SignalP 在线软件预测结果显示,前19 AA 为信号肽,紧接着为糖苷酶18 家族催化域343 AA,然后为43 AA 的连接区域,最后为48 AA 的几丁质结合域(图1)。

从吉富罗非鱼前肠和中肠中克隆得到tChi3的ORF 序 列,全 长1 362 bp,编 码453 AA;Expasy 在线软件预测蛋白大小为50.03 ku,等电点为5.39;SignalP 在线软件预测结果显示,前19 AA 为信号肽,紧接着为糖苷酶18 家族催化域343 AA,然后是43 AA 的连接区域,最后为48 AA 的几丁质结合域(图1)。

图1 罗非鱼3 种几丁质酶氨基酸序列比较Fig.1 Comparison of amino acid sequences of three chitinases of tilapia

从吉富罗非鱼中肠中克隆得到tChit的ORF序列,大小为1 422 bp,编码473AA;Expasy 在线软件预测蛋白大小为52.44 ku,等电点为5.68;SignalP 在线软件预测结果显示,前20 AA 为信号肽,紧接着为糖苷酶18 家族催化域343 AA,然后为49 AA 的连接区域,最后为61 AA 的几丁质结合域(图1)。

比对吉富罗非鱼3 种几丁质酶序列,结果显示tChit1a 与tChi3 的氨基酸序列相似度为83.66%,而tChit 与tChit1a、tChi3 的相似度则较低,分别为49.89%和50.11%。从结构(图1)来看,3 种几丁质酶结构相似,在序列N 端具有信号肽和第18 家族糖苷催化域,C 端具有几丁质结合域,中间由连接域连接。其中,tChit1a 和tChi3 结构相似,几丁质酶催化活性位点序列一致,而tChit的连接域长于tChit1a 和tChi3 的连接域,有较多连续的苏氨酸(Thr),且在催化活性位点有2 个氨基酸不一样。

2.2 罗非鱼3 种几丁质酶的进行树分析

本研究供试罗非鱼与其他物种的几丁质酶进化树分析结果(图2)显示,罗非鱼3 种几丁质酶既不属于哺乳动物的酸性几丁质酶和几丁三糖苷酶,也不属于鱼类酸性几丁质酶1 和2,可归类为鱼类三型几丁质酶,即非酸性几丁质酶。

图2 罗非鱼与其他物种的几丁质酶进化树分析Fig.2 Phylogenetic tree analysis of chitinases among tilapia and other species

2.3 3 种几丁质酶在罗非鱼不同组织中的表达

3 种几丁质酶在罗非鱼不同组织中的表达结果显示,tChit1a 在中肠的表达量最高,其次是前肠(图3A);而tChi3 在前肠的表达量最高,其次是中肠(图3B);tChit 主要在中肠后部表达(图3C)。可见,罗非鱼3 种几丁质酶均主要在肠道中表达,在其他组织中表达量非常低。

图3 几丁质酶tChit1a、tChi3 和tChit mRNA 在罗非鱼不同组织中的分布Fig.3 Distribution of tChit1a,tChi3 and tChit in different tissues of tilapia

2.4 罗非鱼胃几丁质酶的纯化和酶活测定

利用几丁质树脂纯化罗非鱼胃几丁质酶,SDS-PAGE 胶在40 ku 附近有两条明显条带,50 ku 附近没有明显条带(图4A)。本研究所克隆获得的3 种几丁质酶tChit1a、tChi3 和tChit 预测的成熟肽大小分别为48.01、48.11、50.425 ku,但Western blot 验证未在预测大小附近发现明显条带(图4B)。用实验室前期制备的斜带石斑鱼1型几丁质酶多克隆抗体检测,同样没有检测到同源性高的条带。

图4 罗非鱼胃几丁质酶的纯化结果Fig.4 Purification of chitinases extracted from the stomach of tilapia

用4-MU法检测纯化产物的几丁质酶活性,结果如图5所示,纯化产物的几丁质降解酶活性为弱酸性,在pH 4~7之间活性较高,最适pH值为5.0。1 μg纯化产物与2.5 nmol的4MU-(GlcNAc)2/3反应70 min,在最适pH值下,降解4MU-(GlcNAc)2的酶活为1.73 U/g,降解4MU-(GlcNAc)3的酶活为4.89 U/g。

3 讨论

本研究从吉富罗非鱼中克隆获得3 种几丁质 酶tChit1a、tChi3和tChit的ORF,这3 种 几丁质酶主要分布于中肠和前肠。tChit1a、tChi3和tChit 结构相似,均属于18 家族几丁质酶,预测的成熟肽大小分别为48.01、48.11、50.425 ku。其中tChit1a 和tChi3 的全长氨基酸序列相似性较高,达83.66%,催化活性位点序列一致,而tChit 在该位点与它们则有两个氨基酸差异,表明tChit1a 和tChi3 亲缘关系更近。进化树分析表明,3 种几丁质酶不属于鱼类酸性几丁质酶-1(AFCase-1)和鱼类酸性几丁质酶-2(AFCase-2),而与鱼类几丁质酶-3(FCase-3)亲缘关系更近。其中tChit 与金枪鱼几丁质酶3、比目鱼几丁质酶3 和虹鳟几丁质酶聚成一起,而tChit1a 和tChi3 则聚成另一分支,说明它们亲缘关系较近。为更有效降解食物中的几丁质,硬骨鱼鱼类胃中表达两种不同的几丁质酶,即AFCase-1 和AFCase-2,它们具有不同A:12% SDS-PAGE;B:Western blot,一抗:斜带石斑鱼1 型几丁质酶多克隆抗体(1 ∶4 000)M:Marker;1:纯化样品;箭头指示克隆获得的几丁质酶预测大小位置

A:12% SDS-PAGE;B:Western blot by using polyclonal anti-Chitinase 1 ofEpinephelus lanceolatus(1 ∶4 000)M:Marker;1:Purified product;Arrowhead showed the predicted protein of tChit1a,tChi3 and tChit的降解模式,AFCase-1 优先水解GlcNAc 链的非还原端的第二糖苷键,而AFCase-2 优先水解第三糖苷键[15]。但关于FCase-3 的研究较少。在褐菖鲉肝脏和肾脏表达的一种几丁质酶属于FCase-3[16]。有关FCase-3 的降解特性还没有报道。此外,罗非鱼3 种几丁质酶的等电点均在5~6.5 间,略小于7.0,说明这3 种酶均偏中性,推测它们发挥催化作用的环境pH 也偏中性,具体功能需要进一步探索。此外,为继续探索罗非鱼酸性几丁质酶,本研究分析了来源于NCBI 的尼罗罗非鱼基因组注释的另外12 种罗非鱼几丁质酶,包 括XM_003440467.5、XM_003459039.4、X M_0 1 3 2 6 6 3 5 4.2、X M_0 1 9 3 4 9 7 5 6.2、X M_0 2 5 9 0 1 8 9 2.1、X M_0 2 5 9 0 1 8 9 1.1、X M_0 2 5 9 0 1 8 9 0.1、X M_0 1 9 3 4 9 7 5 0.2、X M_0 1 9 3 4 9 7 4 9.2、X M_0 0 5 4 4 8 5 6 1.3、XM_019355318.2 和XM_003459030.3,进化分析结果显示,它们不属于AFCase-1 和AFCase-2,均归类为FCase-3。

研究表明,tChit、tChit1a 和tChi3 主要在罗非鱼前肠和中肠表达,在胃中的表达量相对较低。其他物种研究表明,大菱鲆的3 种几丁质酶在肠道高表达,其次是脑、皮肤和脾脏[17];而日本沙丁鱼几丁质酶的高表达组织是性腺[15];褐菖鲉几丁质酶主要在肝和肾脏中表达[16]。不同的几丁质酶同工酶可能发挥不同的作用,包括免疫、消化和蜕皮等[4,18]。本研究中,tChit、tChit1a 和tChi3 几丁质酶均在肠道高表达,因此,推测这3种几丁质酶可能主要将上游进行了部分分解的几丁质寡糖再进一步水解,从而帮助机体消化吸收几丁质。

为了解罗非鱼是否存在高活性的酸性几丁质酶,本研究通过纯化胃蛋白以期获得酸性几丁质酶。已有多种具有活性的几丁质酶在鱼类胃中分离出来,且有多种同工酶,例如,在白姑鱼的胃组织中纯化获得两种大小分别为42、56 ku 的几丁质酶[19-20];在褐菖鲉胃中发现3 种大小分别为46、52、56 ku 的几丁质酶[21];在日本沙丁鱼胃中纯化得到大小为45、56 ku 的几丁质酶[15]。在六线鱼胃组织中纯化出47、51、62 ku 等3 种几丁质酶[22]。这些纯化获得的几丁质酶具有较高的酶活性,能够降解pNP-(GlcNAc)n,推测几丁质酶活性与这些鱼类栖息地和食物来源有关,它们均以海洋富含几丁质的贝类、甲壳类或浮游生物为食,胃组织中高活性的几丁质酶有助于它们消化食物中的几丁质。几丁质酶存在于许多鱼类的胃中,但并不局限于胃组织,在其他组织中同样发现几丁质酶的表达和酶活性,例如在褐菖鲉肾脏中发现几丁质酶mRNA 的表达[16],在白姑鱼肾脏检测到几丁质酶活性[23],在日本鲭鱼幽门盲囊中检测到几丁质酶mRNA[23]。但还未有研究展示胃以外组织的纯化几丁质酶结果。

在前期研究中,我们成功纯化获得鞍带石斑鱼和赤点石斑鱼胃组织中的几丁质酶,并通过了斜带石斑鱼1 型几丁质酶多克隆抗体的Western blot 验证。本研究中,我们首先从吉富罗非鱼的胃组织中提取总蛋白,通过几丁质树脂亲和层析法分离纯化几丁质酶,SDS-PAGE 结果显示,我们只获得约40 ku 的两条条带,且未能用多克隆抗体检测到这两条条带。我们使用了较DNS法[24]更为灵敏的4-MU 法分析纯化产物酶活性,结果显示,在最适pH 5 的条件下,纯化产物降解4MU-(GlcNAc)2的酶活为1.73 U/g,降 解4MU-(GlcNAc)3的酶活为4.89 U/g。以4-硝基 苯基-β-D-乙酰氨基葡萄糖〔pNP-(GlcNAc)n〕为底物,从六线鱼胃组织纯化获得3 种几丁质酶的最佳酶活分别为3.04、1.25、2.37 U/mg[22],从白姑鱼胃中纯化获得的两种几丁质酶降解pNP-(GlcNAc)3分别为0.109、2.07 U/mg[15]。本研究的罗非鱼纯化产物明显低于这两种鱼的胃几丁质酶。而检测真鲷、海鲈鱼、六线鱼、太平洋鳕鱼和黄尾鰤等胃组织粗酶液的几丁质酶活性在0.24~1.59 U/g 之间[9],与本研究获得的纯化产物几丁质酶活性接近。推测罗非鱼胃中的几丁质酶蛋白表达量较低,通过几丁质树脂亲和层析法并不能很好地纯化几丁质酶,后续可能需要用多种纯化方法相结合,如在现基础上再结合阴/阳离子交换层析方法等,以进一步分离天然的几丁质酶,探究它们对不同形式完整几丁质底物的降解特性及反应机制。几丁质在海洋贮藏量非常丰富[25],其中包括大量虾蟹壳。在处理虾蟹壳废料方面,除了鱼类几丁质酶潜在的降解虾蟹壳几丁质的作用外,细菌几丁质降解酶也是重要的潜在候选者。对虾酱中筛选分离出来的细菌进行全基因组测序,灰色链霉菌CS1801 包含103 个糖苷水解酶家族基因,体现了几丁质的降解潜力[26]。在海鲜工业废物(虾壳倾倒区)分离出的地衣芽孢杆菌(SSCL-10)的最高比活度为3.4 U/mL。对菌株进行化学诱变并通过紫外线照射来筛选潜在的几丁质分解菌群,诱变模型已成功使突变型产生的几丁质酶产量提高了5~6 倍[27]。

4 结论

本研究发现,不同于以贝类、甲壳类或者浮游生物为食的鱼类,罗非鱼3 种几丁质酶tChit1a、tChi3和tChit主要在前肠或中肠中高表达,而在胃中的表达量很低。推测主要原因为罗非鱼是栖息在水体中下层的鱼类,是以植物性饵料为主的杂食性鱼类,其食物中的虾蟹壳所占比例较小,这种食性影响到罗非鱼几丁质酶的表达。因此,这3 种在罗非鱼肠道不同位置表达的几丁质酶可能不同程度地参与罗非鱼的消化吸收,但降解几丁质可能不是其主要功能。

猜你喜欢
罗非鱼克隆鱼类
茂名罗非鱼产业链上新举措
克隆狼
罗非鱼生态健康养殖研究进展
美国消费鳕鱼超过罗非鱼
罗非鱼上半年行情飘红,价格达“最佳”水平!今年能否水涨船高?
鱼类运动会
奇妙的古代动物 泥盆纪的鱼类
属于“我们”
属于“我们”
Cloning Pets克隆宠物