Effects of different photoperiods on gonadal development of Micropterus salmoides
-
摘要:
探究光照时长对大口黑鲈 (Micropterus salmoides) 性腺发育的影响,以解决反季节繁殖中因光照不足导致的亲鱼产卵量少、出苗率低等问题。选取450尾体质量为 (420.8±30.5) g、性腺发育至III期的大口黑鲈,随机分为3组 (6L18D、12L12D和18L6D组),每天分别在6、12和18 h的光照时长下养殖,其他时间为黑暗,每组设3个平行,包括雌、雄个体共50尾,实验周期为60 d。通过测定和分析大口黑鲈体质量、肝体比 (Hepatosomatic index, HSI)、脏体比 (Viscerasomatic index, VSI)、性腺指数 (Gonadsomatic index, GSI)、性类固醇激素水平,以及性腺组织学切片观察和性腺发育相关基因的表达分析,评估光周期对其性腺发育和生长性能的影响。结果显示,18L6D组大口黑鲈体质量虽高于12L12D和6L18D组,但差异不显著 (p>0.05),其HSI、VSI和GSI均显著高于12L12D和6L18D组 (p<0.05)。性腺组织学切片显示,18L6D和12L12D组的精巢与卵巢大多发育至V期,6L18D组多处于IV期,少数发育至V期。18L6D与12L12D组雄鱼睾酮 (Testosterone, T) 浓度差异较小 (p>0.05),但均显著高于6L18D组 (p<0.05)。18L6D组雌鱼雌二醇 (Estradiol, E2) 浓度显著高于其他2组 (p<0.05)。18L6D组雄鱼精巢中dmrt1和gsdf基因与雌鱼卵巢中foxl2和cyp19a1a基因的表达量均显著高于其他2组 (p<0.05)。综上,适宜范围内增加光照时长有助于促进大口黑鲈性腺发育。
Abstract:To investigate the effects of photoperiod on gonadal development in largemouth bass (Micropterus salmoides), aiming to address issues such as reduced spawning and low hatching rates in broodstock caused by insufficient light during off-season reproduction. We randomly divided 450 individuals [Body mass of (420.8±30.5) g and gonadal development at Stage III]into three groups (Group 6L18D, 12L12D, 18L6D), and cultured them under photoperiods of 6, 12 and 18 h per day respectively; for the rest of the time the fish were was in darkness. Each group had three replicates with a total of 50 individualsincluding both males and females, and the experiment lasted 60 d. We evaluated the effects of photoperiod on their gonadaldevelopment and growth performance by measuring the body mass (Wt), hepatosomatic index (HSI), viscerasomatic index (VSI), gonadsomatic index (GSI), levels of sex steroid hormones, as well as observing the histological sections of gonads and analyzing the expression of genes related to gonadal development. The results show that although the Wt value in 18L6D group was higher than that in 12L12D and 6L18D groups without significant differences (p>0.05), its HSI, VSI and GSI values were significantly higher (p<0.05). The histological sections of gonads show that most of the testes and ovaries in 18L6D and 12L12D groups developed to Stage V, while those in 6L18D group were mostly to Stage IV, and a few developed to Stage V. The difference in the testosterone (T) concentration of male between 18L6D and 12L12D group was insignificant (p>0.05), but both were significantly higher than that in 6L18D group. The estradiol (E2) concentration of female in 18L6D group was significantly higher than that in the other two groups. The expression levels of dmrt1 and gsdf genes in the testes of male, and foxl2 and cyp19a1a genes in the ovaries of female in 18L6D group were significantly higher than those in the other two groups (p<0.05). In conclusion,increasing photoperiod within a suitable range helps promote the gonadal development of M. salmoides.
-
Keywords:
- Micropterus salmoides /
- Photoperiod /
- Gonadal development /
- Growth performance
-
琼胶是由海洋藻类生物合成的一种天然高分子多糖[1-2]。琼胶经过水解后生成聚合度为2~20的琼胶寡糖,不仅溶解性好、生物利用度高,还具有改善水产品储藏性能[3]、抗氧化[4]、抗肿瘤[5]、抗炎症[6]、降血糖[7-8]、改善肠道菌群结构[9-10]以及美白保湿[11]等多项理化和生理功能,在功能性食品、化妆品、药品等多个领域均具有很高的经济价值和广泛的应用前景[12]。
近年我国对琼胶的开发存在规模较小、附加值低、易污染环境等不足,开发功能性琼胶寡糖等高附加值的琼脂衍生物对高效利用海洋资源、增加经济效益、减轻化学污染等意义重大。目前工业上主要利用酸水解[13]和酶解[14]两种方式制备琼胶寡糖,其中利用琼胶酶生物催化水解具有反应条件温和、催化效率高、底物特异性强等优势[15]。琼胶酶根据水解模式和产物的不同可分为α-琼胶酶 (EC 3.2.1.158) 和β-琼胶酶 (EC 3.2.1.81)[16]。α-琼胶酶可特异性水解α-1,3-糖苷键,生成以3,6-内醚-α-L-半乳糖为还原性末端的琼寡糖,β-琼胶酶可特异性水解β-1,4-糖苷键,生成以β-D-半乳糖为还原性末端的新琼寡糖。β-琼胶酶因来源丰富、稳定性好得以广泛应用,但仍存在菌株产酶性状不稳定、产酶量和催化活性低、产物专一性差等问题,难以满足工业化生产需求。因此,寻找水解活力更高、性状稳定的琼胶酶基因,并在基因工程菌中实现异源表达,成为生产琼胶酶最经济、高效的方法。
目前,已有部分细菌来源的琼胶酶实现了在大肠杆菌 (Escherichia coli)[17-18]、枯草芽孢杆菌(Bacillus subtilis)[19]、短芽孢杆菌 (B. pumilus)[20]、酿酒酵母 (Saccharomyces cerevisiae)[21]和毕赤酵母 (Pichia pastoris)[22]中的表达,其中大肠杆菌因具有生长与表达速度快、易于基因组修饰、表达量高等优势得到最广泛的应用。现有β-琼胶酶在大肠杆菌中的表达量大多高于原始菌株,但绝大多数是胞内分泌表达,不利于后续的分离、纯化和应用,且表达系统的选择多具有随机性,尚未见有结合生物信息学分析等手段指导进行异源表达的报道。本研究利用生物数据库资源挖掘技术和生物信息学分析手段,获得一个来源于海洋细菌——淡黄色噬琼胶菌 (Agarivorans gilvus) WH0801的β-琼胶酶 (β-AGA酶) 基因。通过分子生物学手段和发酵调控策略实现了其在大肠杆菌BL21(DE3)中的高效胞外表达,以期为新琼寡糖资源的高效制备与充分利用提供新途径。
1. 材料与方法
1.1 实验材料
大肠杆菌JM109、大肠杆菌BL21(DE3)、质粒pET-20b(+) 保藏于本实验室;pUCm-T simple质粒、限制性内切酶 Nco Ⅰ和 Xho Ⅰ、Primer STAR GXL DNA聚合酶、T4 DNA 连接酶、碱性磷酸酶均购自TaKaRa (大连)公司;琼脂糖、质粒小量提取试剂盒、氨苄青霉素等试剂均购自生工生物工程(上海)股份有限公司;分子级酵母粉和胰蛋白胨均购自Oxoid (英国)公司;大豆蛋白胨、酵母提取物 H07014、酵母提取物 H07002均购自上海统园食品技术有限公司;蛋白胨(鱼粉)、牛肉浸膏、酵母浸膏及其他常规培养基均购自国药集团;所用化学试剂均为分析纯[23]。
1.2 实验方法
1.2.1 β-琼胶酶的生物信息学分析
分别以来源于嗜琼胶菌Agarivorans sp. QM38 (NCBI登录号:ABM90422.1)[24]、弧菌Vibrio sp. PO-303 (NCBI登录号:BAG71428.1)[25]和火色杆菌Flammeovirga sp. OC4 (NCBI登录号:AJW82061.1)[26]的β-琼胶酶氨基酸序列为探针,利用美国国家生物信息中心(NCBI) 的BLASTp功能 (https://blast.ncbi.nlm.nih.gov/),设置E值为1.0×10−10,长度覆盖率90%以上且同一性大于30%,进行同源性蛋白分子筛选[27],利用MEGA软件中的邻接 (Neighbor-Joining) 法构建系统发生树[28]。
利用在线生物信息学预测工具分别对β-琼胶酶的基本理化性质、分泌特性、结构特征、功能潜力等进行预测分析[29],利用RobeTTAFold在线服务器 (https://robetta.bakerlab.org/) 进行蛋白质三级结构分析和预测。
1.2.2 β-AGA酶基因的克隆
以淡黄色噬琼胶菌WH0801基因组DNA为模板,根据其中的β-AGA酶序列分别设计含Nco I和Xho I酶切位点的两端引物,通过PCR扩增得到β-aga基因片段[23]。引物设计如下,正向 (P1):5'GCGATGGCCATGGCCACATTTACTAAAAGCAAAATCGCAACCGTTCTT 3' (下划线为Nco I酶切位点));反向 (P2):5'GGTGGTGCTCGAGTTTTTTGTAACGCAGATTATATAGATCACGGTTGAA 3' (下划线为Xho I酶切位点)。
引物由生工生物工程(上海)股份有限公司合成。PCR反应体系及扩增条件参照李兆丰 [30]进行。将得到的β-aga基因克隆到质粒pUCm-T simple,获得克隆载体pUCm-T simple/β-aga,并转化大肠杆菌JM109。提取含100 μg·mL−1氨苄青霉素平板筛选所得阳性菌株的质粒进行酶切电泳和测序验证,序列测定由生工生物工程(上海)股份有限公司完成。
1.2.3 表达载体pET-20b(+)/β-aga的构建和转化
采用限制性内切酶Nco I和Xho I对克隆载体pUCm-T simple/β-aga双酶切,经碱性磷酸酶消化后回收目的片段,再用T4 DNA连接酶与经同样双酶切处理的pET-20b(+)质粒片段连接后转化大肠杆菌JM109,提取质粒并测序鉴定,获得含正确序列的表达载体pET-20b(+)/β-aga,将其转化表达宿主大肠杆菌BL21(DE3),获得基因工程菌E. coli BL21(DE3)(pET-20b(+)/β-aga)。
1.2.4 重组β-AGA酶的生产
种子培养和发酵培养参照李兆丰[30]的方法进行,发酵结束后收集上清液即为重组β-AGA酶粗酶液。
1.2.5 发酵条件优化
以体积分数为4%的接种量将种子培养液接至50 mL的发酵培养基中,在不同的发酵初始培养基、诱导剂添加量及添加时间、发酵温度、发酵时间、发酵pH、碳源和氮源条件下进行发酵,测定发酵上清液中重组β-AGA酶活力,以确定最佳的发酵条件,每次优化的结果用于之后的实验。
1.2.6 正交试验
在单因素试验结果的基础上,选择3个影响最大的因素进行正交试验,测定发酵液中重组β-AGA酶活力,获得发酵产酶的最佳条件。
1.2.7 β-AGA酶水解活力的测定
以质量分数为0.25%的琼脂糖溶液为底物,β-AGA酶的水解活力参照Liu等[31]和Dong等[32]采用3,5-二硝基水杨酸 (DNS) 法测定。
1.2.8 数据处理
实验结果用3次平行实验的“平均值±标准偏差 (
$\overline { X}\pm { \rm {SD}} $ ) ”表示。采用GraphPad Prism 9.2软件分析数据和作图;显著性分析 (P<0.05) 通过SPSS 20.0软件的单因素方差分析 (One-way ANOVA) 中的Turkey (T)程序进行。2. 结果
2.1 β-琼胶酶的生物信息学分析
2.1.1 β-琼胶酶系统进化树的构建
来源于嗜琼胶菌Agarivorans sp. QM38、弧菌Vibrio sp. PO-303和火色杆菌Flammeovirga sp. OC4的β-琼胶酶已被报道具有优良的环境适应力和水解活力,因此以这3个氨基酸序列为探针,利用NCBI数据库中的BLASTp功能进行搜索,保留同源性大于30%的序列,共得到1 039条序列。对筛选出的序列构建系统进化树 (图1),可根据其菌属来源分为11大类,考虑到可获得性和潜在工业价值,最终选取其中与模板序列亲缘关系较近的来源于淡黄色噬琼胶菌的β-琼胶酶(β-AGA酶,NCBI登录号:AQT38174.1) 作为后续研究对象。
2.1.2 β-AGA酶理化性质预测
利用ProtParam预测了β-AGA酶的理化性质。β-AGA酶的氨基酸数量为955个,氨基酸组成见表1,其中含量较多的是丙氨酸(82个,8.6%)、天冬氨酸(81个,8.5%);其中带负电荷的氨基酸残基有135个,带正电荷的氨基酸残基有84个;仅含有1个半胱氨酸,推测可能具有较少的二硫键结构。β-AGA酶的分子式为C4767H7184N1238O1484S21,相对分子量为106.254 kD,理论等电点为 4.66,不稳定指数为24.45,总亲水性平均值为−0.428,表明β-AGA酶为酸性、稳定的亲水性蛋白。
表 1 β-AGA酶的氨基酸组成Table 1. Amino acid composition of β-AGAase氨基酸
Amino acid数量
Amount占比
Proportion/%丙氨酸 Alanine 82 8.6 精氨酸 Arginine 27 2.8 天冬酰胺 Aspartic acid 51 5.3 天冬氨酸 Asparagine 81 8.5 半胱氨酸 Cysteine 1 0.1 谷氨酰胺 Glutamine 38 4.0 谷氨酸 Glutamic acid 54 5.7 甘氨酸 Glycine 73 7.6 组氨酸 Histidine 15 1.6 异亮氨酸 Isoleucine 40 4.2 亮氨酸 Leucine 70 7.3 赖氨酸 Lysine 57 6.0 蛋氨酸 Methionine 20 2.1 苯丙氨酸 Phenylalanine 48 5.0 脯氨酸 Proline 40 4.2 丝氨酸 Serine 74 7.7 苏氨酸 Threonine 56 5.9 色氨酸 Tryptophane 26 2.7 酪氨酸 Tyrosine 39 4.1 缬氨酸 Valine 63 6.6 2.1.3 β-AGA酶二级结构分析
利用SOPMA预测了β-AGA酶可能包含的二级结构 (图2),其中303个氨基酸参与构成α-螺旋,30个氨基酸参与构成β-转角,170个氨基酸参与构成延伸链,452个氨基酸参与构成无规则卷曲。
2.1.4 β-AGA酶三维结构预测
采用RobeTTAFold对β-AGA酶结构进行在线预测,得到5个预测结构 (图3),置信度均为0.8 (满分为1),比对各模型的误差估算图,其中模型2具有相对较小的误差,可以作为基础进行进一步结构解析。
2.1.5 β-AGA酶亲水性/疏水性分析
利用Prot Scale分析了β-AGA酶的亲水性/疏水性 (图4)。β-AGA酶蛋白的最低分值(−3.100)峰区位于第537位脯氨酸,代表其亲水性最强;最高分值(2.422)峰区位于第16位苯丙氨酸,代表其疏水性最强。分析可得,β-AGA酶蛋白的亲水区多于疏水区,总亲水平均值为−0.400,是一种亲水性蛋白。
2.1.6 β-AGA酶亚细胞定位分析
利用PSORTb version 3.0分析了β-AGA酶的亚细胞定位,结果显示无法准确预测其亚细胞定位。同时利用Loc Tree3分析其亚细胞定位 (图5),结果表明β-AGA酶有一定可能定位在细胞周质中,亚细胞定位分数为36,精确度为89%。
2.1.7 β-AGA酶信号肽分析
利用Signal P 5.0 Server分析了β-AGA酶的信号肽 (图6),可知β-AGA酶的S-mean和Y-max的平均值 (D) 小于临界值,代表不具有信号肽序列,推测β-AGA酶为非分泌性蛋白。
2.1.8 β-AGA酶功能结构域分析
利用NCBI的保守域数据库 (CDD) 进行功能结构域分析 (图7)。β-AGA酶在N端48—211、240—424位氨基酸处含有2个Agarase_CBM结构域,该CBM样域在结构上与某些CBM家族非常相似,其中的一个Loop参与了催化位点通道顶部的形成,在β-琼胶酶的外切水解模式中起到了识别底物的作用。这些结构域与琼脂糖链的非还原末端特异性结合,并将附加的催化模块引导至受琼胶酶水解的区域,帮助β-AGA酶捕捉底物,促进参与1,4-β-D-半乳糖苷键水解的催化模块活性。
2.2 β-AGA酶大肠杆菌表达系统的构建
β-AGA酶在大肠杆菌中表达载体pET-20b(+)/β-aga的构建见图8。通过双酶切法构建获得分泌型表达载体pET-20b(+)/β-aga。测序验证序列正确性后将表达载体转化宿主大肠杆菌BL21(DE3),获得E.coli BL21(DE3)(pET-20b(+)/β-aga),实现了β-AGA酶在大肠杆菌中的胞外分泌表达。经测定,发酵所得重组β-AGA酶的初始酶活为2.79 U·mL−1,比酶活为121.30 U·mg−1。
2.3 重组β-AGA酶的发酵优化
2.3.1 种子生长曲线与种龄的确定
重组β-AGA酶的种子生长曲线见图9-a,菌体在培养前3 h处于延滞期,3~10 h处于对数生长期,10 h后进入稳定期。选择种龄4~9 h的β-AGA酶种子按体积分数为4%的接种量接种,以发酵时间24 h、发酵温度25 ℃的条件进行培养,离心后测定上清液中的胞外酶活力。对数中期 (≥5 h) 的种龄对重组β-AGA酶活力影响不显著,综合考虑时间和经济成本,选择种龄5 h的种子培养液进行发酵 (图9-b)。
2.3.2 发酵初始培养基的选择
发酵初始培养基对酶活的影响结果见图10。TB培养基对重组β-AGA酶的胞外发酵生产效果最好,胞外酶活力最高,达到3.20 U·mL−1。分析其原因主要为TB培养基营养丰富,更适合菌体生长,且酵母粉含量更高,能够促进重组酶释放到大肠杆菌细胞的周质空间并形成高渗透压便于胞外产酶;同时,TB培养基的pH缓冲能力较强,能够有效维护重组菌的生长。
2.3.3 诱导剂对重组β-AGA酶生产的影响
诱导剂浓度和添加时间也会在一定程度上影响重组β-AGA酶的生产(图11-a)。当IPTG诱导浓度介于0.005~0.050 mmol·L−1时酶活较稳定,IPTG浓度为0.025 mmol·L−1时酶活最高 (2.589 U·mL−1)。诱导剂添加时间对重组β-AGA酶生产的影响见图11-b,添加诱导剂前培养时间超过0.5 h,酶活均较稳定,培养时间为1 h时酶活最高 (3.96 U·mL−1)。
2.3.4 发酵时间、温度和pH的选择
发酵时间对重组β-AGA酶生产的影响见图12-a,发酵32~64 h其胞外酶活均较稳定,发酵第48小时时酶活达到最高。
发酵温度优化结果见图12-b,重组β-AGA酶活力在20~30 ℃处于较高水平,25 ℃下发酵48 h酶活最高,达到4.37 U·mL−1;超过30 ℃时酶活显著下降。
发酵初始培养基的pH会影响微生物的代谢、发酵及分泌,在pH 7~8.5酶活保持稳定 (图12-c),综合考虑环境和经济因素,选择7.0作为重组β-AGA酶发酵初始pH。
2.3.5 碳源和氮源对重组β-AGA酶生产的影响
发酵出发培养基优化结果表明,在9种常见的培养基中,TB培养基最有利于重组β-AGA酶的分泌。以此为基础,对其碳源组成进行优化以达到降低成本、提高效率的目的。以去除TB培养基中的碳源为空白,选取并添加与该培养基碳含量相等的果糖、麦芽糖、蔗糖等9种不同的碳源,结果表明碳源的种类和添加量对重组β-AGA酶表达分泌存在较大影响 (图13)。除了葡萄糖,其他碳源在重组β-AGA酶发酵中均能取得比甘油更好的效果,其中果糖作为碳源酶活达到最高 (6.28 U·mL−1)。进一步优化果糖的添加量,当果糖添加量为6 g·L−1时,胞外产酶活力最高 (7.45 U·mL−1)。
氮源是微生物合成核酸、蛋白质等的重要原料。氮源分为无机氮和有机氮,有机氮更适合大肠杆菌重组菌的发酵生产。以去除TB培养基中的氮源为空白,选取并添加与TB发酵培养基氮含量相等的不同种类的蛋白胨、酵母提取物,发现氮源的种类和添加量对重组β-AGA酶发酵产酶的影响存在明显差异 (表2)。酵母提取物作为氮源时的表达量整体明显优于蛋白胨。其中,使β-AGA酶发酵量达到最高的氮源种类和浓度为添加量30 g·L−1的酵母提取物H07014,此时胞外产酶活力达到15.33 U·mL−1。
表 2 有机氮源对重组β-AGA酶发酵的影响Table 2. Effect of organic nitrogen source on fermentation of recombinant β-AGAaseU·mL−1 氮源种类
Nitrogen source添加量 Additive amount/(g·L−1) 0 6 12 18 24 30 ① 1.37±0.03a 2.14±0.02b 2.94±0.03c 3.79±0.04d 4.48±0.03e 5.83±0.02f ② 1.37±0.03a 3.03±0.02b 4.26±0.04c 6.22±0.03d 8.04±0.04e 9.21±0.09f ③ 1.37±0.03a 3.28±0.02b 4.35±0.02c 6.19±0.03d 8.47±0.06e 8.64±0.07e ④ 1.37±0.03a 3.65±0.03b 4.80±0.02c 6.73±0.04d 8.51±0.05e 9.58±0.07f Ⅰ 1.37±0.03a 4.10±0.02b 6.67±0.05c 8.92±0.06d 10.44±0.05e 11.87±0.08f Ⅱ 1.37±0.03a 4.76±0.02b 7.59±0.04c 10.15±0.05d 12.48±0.06e 15.33±0.10f Ⅲ 1.37±0.03a 4.28±0.02b 7.02±0.04c 9.58±0.04d 11.57±0.07e 13.25±0.08f Ⅳ 1.37±0.03a 3.82±0.01b 5.04±0.03c 6.96±0.03d 8.81±0.04e 9.94±0.03f 注:数值为 3 次平行实验的平均值,同一行中不同上标字母表示显著性差异(P<0.05)。
①. 分子级蛋白胨;②. 蛋白胨(鱼粉);③. 牛肉浸膏;④. 大豆蛋白胨;Ⅰ. 分子级酵母提取物;Ⅱ. 酵母提取物 H07014;Ⅲ. 酵母提取物 H07002;Ⅳ. 酵母浸膏。 Note: Each value represents the average value of three independent measurements, and those with different letters within the same row are significantly different (P<0.05).
①. Molecular level peptone; ②. Peptone (fish meal); ③. Beef extract; ④. Soya peptone; I. Molecular level yeast extract; II. Yeast extract H07014; III. Yeast extract H07002; IV. Yeast extract.2.3.6 正交试验
单因素试验结果表明发酵温度、碳源添加量、IPTG添加时间3个因素对重组菌株产酶影响最大,选择这3个因素做三因素三水平正交试验,正交试验设计见表3,试验结果见表4。最佳组合为A2B2C3,即碳源添加量为6 g·L−1,发酵温度为25 ℃,IPTG添加时间为2.0 h,在此条件下发酵48 h后所得粗酶液酶活为16.72 U·mL−1。
表 3 正交试验因素和水平表Table 3. Factors and levels in orthogonal array design水平
LevelA:碳源添加量
Carbon source
additive amount/
(g·L−1)B:发酵温度
Fermentation
temperature/
℃C:IPTG添加时间
IPTG additive
time/h1 5 20 1.0 2 6 25 1.5 3 7 30 2.0 表 4 发酵条件优化正交试验结果Table 4. Orthogonal array design layout and experimental results试验号
Test No.A:碳源
添加量
Carbon source
additive amount/
(g·L−1)B:发酵
温度
Fermentation
temperature/
℃C:IPTG
添加时间
IPTG additive
time/h酶活
Enzyme
activity/
(U·mL−1)1 1 (5) 1 (20) 1 (1.0) 13.99±0.04e 2 1 2 (25) 2 (1.5) 14.15±0.05f 3 1 3 (30) 3 (2.0) 13.02±0.07b 4 2 (6) 1 2 15.63±0.06g 5 2 2 3 16.72±0.06h 6 2 3 1 13.39±0.04c 7 3 (7) 1 3 13.64±0.05d 8 3 2 1 14.18±0.03f 9 3 3 2 12.67±0.02a K1 13.72 14.42 13.85 K2 15.25 15.02 14.15 K3 13.50 13.03 14.46 R 1.75 1.99 0.61 注:同一列中不同上标字母表示显著性差异 (P<0.05)。 Note: Different superscript letters within the same column indicate significant difference (P<0.05). 3. 讨论
对琼胶酶进行异源表达可以突破天然酶的分泌限制,利用基因工程菌大幅提高其产量,可实现优良酶制剂的高效制备。目前国内外对琼胶酶在大肠杆菌中的克隆表达主要定位于细胞质或周质空间,存在酶产量有限、制备工艺复杂、分离纯化困难等瓶颈[33],限制了其工业化应用。本研究通过生物信息学分析,发现来源于淡黄色噬琼胶菌的β-AGA酶氨基酸序列中不包含信号肽,因此引入pET-20b(+)质粒中的信号肽进行表达载体的构建,成功实现了重组β-AGA酶的胞外分泌表达。
发酵过程中菌体的生长状态和产酶量可通过多种因素进行系统调控。菌体在生长的不同阶段具备不同的生理特征,处于对数生长期的种子因其活性高、繁殖快、适应性强,更适合用于接种发酵。由于宿主菌E.coli BL21(DE3)中存在由lac I基因调控的T7 RNA聚合酶基因激活机制[34],因此在发酵过程中通常需要添加诱导剂IPTG诱导目的基因的转录。本实验结果表明,IPTG浓度为0.025 mmol·L−1时重组β-AGA酶的水解活力最高,随着诱导剂浓度进一步升高,酶活逐渐降低。这可能是由于IPTG具有一定的细胞毒性,浓度过高时会抑制菌体生长[35],且经高浓度IPTG诱导后重组β-AGA酶的前体蛋白合成速率会大大提高[30,36],容易造成细胞质中无活性包涵体的产生,从而导致正常折叠的蛋白质分泌量显著降低。同时,诱导剂的添加时间也会对重组β-AGA酶的表达量和酶活造成影响,先培养一段时间再添加诱导剂,有利于使菌体增殖到适当的浓度,从而更好地抵御IPTG的毒性,最大限度地诱导细胞产酶。另外,随着发酵时间的延长,菌体在利用营养物质生长繁殖的同时又将更多的代谢产物排出至培养基中,当代谢产物积累到一定量时,会反向阻遏菌体的生产及产酶,因此需要对发酵时间进行严格把控。发酵体系的pH也是影响微生物生长的一个重要因素,适中的pH有利于重组菌的生长、代谢和重组β-AGA酶在发酵期间活力的稳定。除pH外,发酵温度也是调控重组β-AGA酶胞外分泌表达的一项重要指标。本研究结果表明,重组β-AGA酶在25 ℃下发酵活性最高,可能是由于该温度更接近野生型海洋来源β-AGA酶所处的环境温度,此温度下该酶本身稳定性较高,且相对较低的温度也更适合宿主大肠杆菌的生长繁殖和分泌产酶。相反,当发酵温度上升至37 ℃时,宿主菌和重组酶的稳定性均有所下降,在拷贝过程中质粒可能存在较严重的丢失现象,因此胞外酶活显著降低。当诱导温度太低 (20 ℃)时,宿主菌的增殖速率和重组酶的转运速率均会降低,对重组β-AGA酶的胞外生产起到一定的限制作用。碳源和氮源的种类和浓度也是影响微生物生长代谢的重要因素,浓度过低时不足以提供充足的营养,过高则会反向抑制微生物的生长和代谢[37],不同微生物对碳、氮源物质的选择偏好性存在明显差异。
综合上述各因素对重组β-AGA酶进行发酵条件优化后的最终条件为:选取37 ℃下培养5 h的种子液用于发酵产酶,以TB作为发酵出发培养基,设置初始pH为7.0,以6 g·L−1的果糖代替其中的碳源,以30 g·L−1的酵母提取液Ⅱ代替其中的氮源,于25 ℃培养2 h后加入终浓度为0.025 mmol·L−1的IPTG,继续诱导48 h,此条件下发酵所得酶活达到16.72 U·mL−1,约为初始酶活的6倍,在目前已报道的β-琼胶酶中处于领先水平[38-40]。优良的水解活力和分泌稳定性使得重组β-AGA酶具有广阔的应用空间,有望成为定向水解琼脂糖,高效生产特异性新琼寡糖的新型酶解工具。
4. 结论
本研究利用基因挖掘技术从海洋微生物基因库中筛选得到一个具有优良催化潜力的β-琼胶酶基因,采用分子生物学手段实现了其在大肠杆菌系统中的异源表达,并通过发酵调控策略的优化使其酶活提高了约5倍。该酶具有优良的水解活力、稳定的分泌水平和较高的产物特异性,应用前景广阔。研究结果可为开发新型琼胶水解酶和高效制备不同聚合度的新琼寡糖提供参考。
-
图 1 不同光照时长下大口黑鲈精巢组织切片观察
注:a. 6L18D组第60天时的精巢;a1—a3. 6L18D组第0、第30、第60天时的精巢;b. 12L12D组第60天时的精巢;b1—b3. 12L12D组第0、第30、第60天时的精巢;c. 18L6D组第60天时的精巢;c1—c3. 18L6D组第0、第30、第60天时的精巢;Oo3. III期精细胞;Oo4. IV期精细胞; Oo5. V期精细胞;SP. 精子;Ss. 精小囊;L. 小叶腔。
Figure 1. Observation of spermial tissue section of M. salmoides during different photoperiods
Note: a. Testis on 60th day in Group 6L18D; a1−a3. Testis on 0th, 30th, 60th day in Group 6L18D; b. Testis on 60th day in Group 12L12D;b1−b3. Testis on 0th, 30th, 60th day in Group 12L12D; c. Testis on 60th day in Group 18L6D; c1−c3. Testis on 0th, 30th, 60th day in Group 18L6D;Oo3. Stage III spermatocytes; Oo4. Stage IV spermatocytes; Oo5. Stage V spermatocytes; SP. Sperm; Ss. Seminal vesicle; L. Lobular cavity.
图 2 不同光照时长下大口黑鲈卵巢组织切片观察
注: a. 6L18D组第60天的卵巢;a1—a3. 6L18D组第0、第30、第60天时的卵巢;b. 12L12D组第60天时的卵巢;b1—b3. 12L12D组第0、第30、第60天时的卵巢;c. 18L6D组第60天时的卵巢;c1—c3. 18L6D组第0、第30、第60天时的卵巢;Oo2. II期卵母细胞;Oo3. III期卵母细胞; Oo4. IV期卵母细胞; Oo5. V期卵母细胞;yg. 卵黄颗粒;ob. 油球;n. 细胞核。
Figure 2. Observation of ovarian tissue section of M. salmoides during different photoperiods
Note: a. Ovaries on 60th day in Group 6L18D; a1−a3. Ovaries on 0th, 30th, 60th day in Group 6L18D; b. Ovaries on 60th day in Group 12L12D; b1−b3. Ovaries on 0th, 30th, 60th day in Group 12L12D; c. Ovaries on 60th day in Group 18L6D; c1−c3. Ovaries on 0th, 30th, 60th day in Group 18L6D;Oo2. Stage II oocyte; Oo3. Stage III oocyte; Oo4. Stage IV oocyte; Oo5. Stage V oocytes; yg. Yolk granules; ob. Oil globules; n. Nucleus.
图 3 不同光照时长对大口黑鲈血清中睾酮和雌二醇的影响
注:a. 雄鱼T浓度;b. 雌鱼E2浓度;同组中不同字母间存在显著性差异 (p<0.05)。
Figure 3. Effects of different photoperiods on serum testosterone and estradiol of M. salmoides
Note: a. Testosterone content of male fish; b. Estradiol content of female fish; different letters at the top of the bar represent significant differences (p<0.05).
图 4 不同光照时长对大口黑鲈性腺发育相关基因表达的影响
注:a. 雄鱼dmrt1基因表达;b. 雄鱼gsdf基因表达;c. 雌鱼foxl2基因表达;d. 雌鱼cyp19a1a基因表达;柱形图上方字母不同表示差异显著(p<0.05)。
Figure 4. Effects of different photoperiods on gene expression related to gonad development in M. salmoides
Note: a. dmrt1 gene expression in male fish; b. gsdf gene expression in male fish; c. foxl2 gene expression in female fish; d. cyp19a1a gene expression in female fish; different letters at the top of the bar represent significant differences (p<0.05).
表 1 qRT-PCR引物序列
Table 1 qRT-PCR primers sequences
基因
Gene正向引物
Forward primer反向引物
Reverse primerfoxl2 5'-CACTCGACCCAGCCTGTGAG-3' 5'-GCCGTCTCCTCCGAACAAGG-3' cyp19a1a 5'-GTGAGGCAGTGTGTGCTGGA-3' 5'-CAGCCGCAGCTCCACATCT-3' dmrt1 5'-GCTCCCGCTGTAGGAACCAC-3' 5'-CCTGAGCCTGCTGCCTTCTC-3' gsdf 5'-AGCAATGCCCACAGAGCCAT-3' 5'-GGCTGGCAGGATGGATCACC-3' β-actin 5'-AAAGGGAAATCGTGCGTGAC-3' 5'-AAGGAAGGCTGGAAGAGGG-3' 表 2 不同光周期下大口黑鲈体质量、HSI、VSI和GSI的差异
Table 2 Variations in body mass, HSI, VSI and GSI of M. salmoides during different photoperiods
时间
t/d组别
Group成活率
Survival rate/%饵料系数
FCR体质量
Body mass/g肝体比
HSI/%脏体比
VSI/%性腺指数 GSI/% 雄鱼
Male雌鱼
Female0 6L18D 100 1.0±0.08 394.85±14.69 1.37±0.10 7.70±0.92 0.94±0.08 1.45±0.14 12L12D 100 1.0±0.08 399.62±14.34 1.40±0.21 7.70±1.04 0.92±0.05 1.45±0.14 18L6D 100 1.0±0.08 399.14±10.38 1.40±0.20 7.45±0.55 0.96±0.03 1.32±0.35 30 6L18D 100 1.1±0.12a 405.28±19.54a 1.55±0.04a 8.40±0.27a 1.04±0.07a 2.60±0.34a 12L12D 100 0.95±0.07b 432.12±13.42b 1.65±0.04b 8.94±0.25b 1.22±0.13ab 4.35±0.41b 18L6D 100 0.93±0.05b 459.62±17.92c 1.67±0.09b 9.35±0.58c 1.32±0.15b 5.26±0.12c 60 6L18D 96 1.08±0.03a 488.00±14.00a 1.60±0.16a 8.48±0.56a 1.140.11a 4.10±0.27a 12L12D 100 1.0±0.05a 499.66±29.25a 1.68±0.19ab 9.30±1.08ab 1.31±0.05ab 5.60±0.19ab 18L6D 100 1.0±0.05a 511.5±8.98a 1.74±0.08b 9.72±1.04b 1.38±0.14b 5.97±0.14b 注:同列中不同字母间存在显著性差异 (p<0.05)。 Note: Values with different letters within the same column represent significant differences (p<0.05). -
[1] 崔庆奎, 沈志刚, 田宇, 等. 大口黑鲈的卵巢发育周年变化及反季节繁殖研究[J]. 水生生物学报, 2021, 45(1): 76-88. [2] HUSSEIN G H G, CHEN M, QI P P, et al. Aquaculture industry development, annual price analysis and out-of-season spawning in largemouth bass Micropterus salmoides[J]. Aquaculture, 2020, 519: 734901.
[3] 温海深, 林浩然. 环境因子对硬骨鱼类性腺发育成熟及其排卵和产卵的调控[J]. 应用生态学报, 2001, 12(1): 151-155. [4] 刘春晓. 光周期对斑马鱼繁殖力、血浆激素含量及卵巢发育相关基因表达的影响[D]. 上海: 上海海洋大学, 2019: 1-2. [5] 陈戟, 胡炜, 朱作言. 鱼类生殖发育调控研究进展[J]. 科学通报, 2013, 58(2): 103-114. [6] ZOHAR Y, MUNOZ-CUETO J A, ELIZUR A, et al. Neuroendocrinology of reproduction in teleost fish[J]. Gen Comp Endocr, 2010, 165(3): 438-455.
[7] SHEWMON L N, GODWIN J R, MURASHIGE R S, et al. Environmental manipulation of growth and sexual maturation in yellow perch, Perca flavescens[J]. J World Aquac Soc, 2007, 38(3): 383-394. doi: 10.1111/j.1749-7345.2007.00110.x
[8] ACHARJEE A, CHAUBE R, JOY K P. Effects of altered photoperiod and temperature on expression levels of gonadotrophin subunit mRNAs in the female stinging catfish Heteropneustes fossilis[J]. J Fish Biol, 2017, 90(6): 2289-2311.
[9] ZHU D M, YANG K , GUL Y, et al. Effect of photoperiod on growth and gonadal development of juvenile topmouth gudgeon Pseudorasbora parva[J]. Environ Biol Fish, 2014, 97(2): 147-156.
[10] 张明辉. 温度对池蝶蚌性腺发育和内源性性激素的影响及性别相关因子Dmrt1、Foxl2的调控作用[D]. 南昌: 南昌大学, 2023: 9-11. [11] WANG D S, ZHOU L Y, KOBAYASHI T, et al. Doublesex- and Mab-3-related transcription factor-1 repression of aromatase transcription, a possible mechanism favoring the male pathway in tilapia[J]. Endocrinology, 2010, 151(3): 1331-1340.
[12] 樊森. Gsdf在尼罗罗非鱼精子发生和育性中的作用研究[D]. 湛江: 广东海洋大学, 2023: 8-10. [13] KOBAYASHI Y, HORIGUCHI R, NOZU R, et al. Expression and localization of forkhead transcriptional factor 2 (Foxl2) in the gonads of protogynous wrasse, Halichoeres trimaculatus[J]. Biol Sex Differ, 2010, 1(1): 3.
[14] FANG D A, YANG X J, FENG X, et al. FoxL2 combined with Cyp19a1a regulate the spawning upstream migration in Coilia nasus [J]. Gene, 2019, 710: 307-315.
[15] D'COTTA H, FOSTIER A, GUIGUEN Y, et al. Aromatase plays a key role during normal and temperature-induced sex differentiation of tilapia Oreochromis niloticus[J]. Mol Reprod Dev, 2001, 59(3): 265-276.
[16] 郑尧, 王在照, 陈家长. 调控鱼类性腺分化基因的研究进展[J]. 水生生物学报, 2015, 39(4): 798-810. [17] 黄太寿, 白俊杰, 李胜杰, 等. 我国加州鲈养殖现状和绿色发展问题探讨[J]. 中国水产, 2017(12): 44-47. [18] 农业农村部渔业渔政管理局, 全国水产技术推广总站, 中国水产学会. 2024中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2024: 32-33. [19] 杨刚, 李胜杰, 连总强, 等. 温度调控对大口黑鲈性腺发育的影响[J]. 淡水渔业, 2023, 53(4): 26-33. [20] MATTHEWS D. M, STOUT R B. Out-of-season spawning method for Florida largemouth bass to produce advanced-sized fingerlings by early spring[J]. N Am J Aquac, 2013, 75(4): 524-531.
[21] 周家辉, 杜金星, 姜鹏, 等. 17α-甲基睾酮对大口黑鲈生长及性腺发育的影响[J]. 中国水产科学, 2021, 28(9): 1109-1117. doi: 10.12264/JFSC2021-0053 [22] 朱成德. 大银鱼的肝指数及其与卵巢发育的关系[J]. 淡水渔业, 1984(5): 19-21. [23] 侯亚义, 韩晓冬. 光周期影响雌性虹鳟鱼的免疫活动(英文)[J]. 动物学研究, 1999(6): 401-405. [24] 陈慧, 林国文, 刘招坤, 等. 网箱养殖大黄鱼生长特性的研究[J]. 海洋科学, 2010, 34(11): 1-5. [25] FENG C, SARIGAIQIQIGE, LIU W Y, et al. Effect of dark environment on intestinal flora and expression of genes related to liver metabolism in zebrafish (Danio rerio)[J]. Comp Biochem Physiol C: Toxicol Pharmacol, 2021, 249: 109100.
[26] 华元渝, 李廷友, 邹宏海. 养殖型暗纹东方鲀生长与发育特性[J]. 水产学报, 2004, 28(1): 8-14. [27] 张士璀, 孙旭彤, 李红岩. 卵黄蛋白原研究及其进展[J]. 海洋科学, 2002(7): 32-35. [28] 成永旭, 李少菁, 王桂忠, 等. 锯缘青蟹卵黄发生期卵巢和肝胰腺脂类的变化[J]. 海洋学报(中文版), 2001(3): 66-77. [29] 陈文银. 光周期对黄金鲈鱼生长的影响[J]. 上海水产大学学报, 2006(3): 3300-3307. [30] WANG X B, LI Y Y, HU J B, et al. Effects of different photoperiods on growth and ovarian development and maturation of silver pomfret Pampus argenteus[J]. J Fish Biol, 2023, 103(1): 59-72. doi: 10.1111/jfb.15413
[31] BROMAGE N, PORTER M, RANDALL C. The environmental regulation of maturation in farmed finfish with special reference to the role of photoperiod and melatonin[J]. Aquaculture, 2001, 197(1): 63-98.
[32] THORPE J. Maturation responses of salmonids to changing developmental opportunities[J]. Mar Ecol Prog, 2007, 335: 285-288.
[33] GUO T, WANG Y F, JIANG P, et al. miR-155 regulates photoperiod induced gonadal development in Atlantic salmon (Salmo salar) by targeting brain-derived neurotrophic factor[J]. Fishes, 2023, 8(7): 345. doi: 10.3390/fishes8070345
[34] BAEK H J, KIM H B, NOH S, et al. Changes in plasma steroid hormone levels and gonad development by the control of photoperiods and water temperatures on timing of sexual maturity of rockfish (Sebastes schlegeli)[J]. Fish Aquat Sci, 2004, 7(1): 16-22.
[35] LEE C H, PARK Y J, LEE Y D. Effects of photoperiod manipulation on gonadal activity of the damselfish, Chromis notata[J]. Dev Reprod, 2017, 21(2): 223-228. doi: 10.12717/DR.2017.21.2.223
[36] BROWN M L, KASIGA T, SPENGLER D E, et al. Reproductive cycle of northern largemouth bass Micropterus salmoides salmoides[J]. J Exp Zool A, 2019, 331(10): 540-551. doi: 10.1002/jez.2323
[37] NISHIMURA T, TANAKA M. The mechanism of germline sex determination in vertebrates[J]. Biol Reprod, 2016, 95(1): 1-6. doi: 10.1095/biolreprod.116.142240
[38] JIANG D N, YANG H H, LI M H, et al. gsdf is a downstream gene of dmrt1 that functions in the male sex determination pathway of the Nile tilapia[J]. Mol Reprod Dev, 2016, 83(6): 497-508. doi: 10.1002/mrd.22642
[39] YARMOHAMMADI M, POURKAZEMI M, KAZEMI R, et al. Sex steroid level and sexual dimorphism expression of genes in gonads of the great sturgeon Huso huso Linneaus, 1758 during maturity developmental stages[J]. Aquac Res, 2017, 48(4): 1413-1429. doi: 10.1111/are.12977
[40] 祝斐, 张志勇, 徐献明, 等. 光周期对松江鲈鱼生长和性腺发育影响的初步研究[J]. 海洋科学, 2016, 40(9): 32-38. doi: 10.11759//hykx20150606002 [41] QIANG J, HE J, ZHU J H, et al. Optimal combination of temperature and photoperiod for sex steroid hormone secretion and egg development of Oreochromis niloticus as determined by response surface methodology[J]. J Therm Biol, 2021, 97: 102889. doi: 10.1016/j.jtherbio.2021.102889
[42] SAMPATH KUMAR R, IJIRI S, TRANT J M. Changes in the expression of genes encoding steroidogenic enzymes in the channel catfish (Ictalurus punctatus) ovary throughout a reproductive cycle[J]. Biol Reprod, 2000, 63(6): 1676-1682. doi: 10.1095/biolreprod63.6.1676