微囊藻CO2-浓缩机制基因型的动态变化及其对CO2的竞争效应*

杜雨欣,郑保海,李佳欣,李玉鑫,黄立成,施军琼,吴忠兴

(1:西南大学三峡库区生态环境教育部重点实验室,重庆市三峡库区植物生态与资源重点实验室,重庆 400715) (2:昆明市滇池高原湖泊研究院,昆明 650228)

自工业革命开始,化石燃料的使用快速增加,同时伴随着乱砍滥伐、过度开垦等环境问题的发生,使得大气中的CO2浓度以前所未有的速度持续增加,预计将在21世纪末超过0.08%(V/V)[1]。CO2作为光合作用的主要原料,在藻类的繁殖和生长中发挥着重要作用[2]。而大气中CO2浓度的升高会影响海洋和湖泊等水体表面的CO2气体交换,使得水体环境发生变化,影响水体中浮游植物细胞的新陈代谢过程,进而对水生生态系统和生物地球化学循环产生深远影响[3]。因而,关于浮游植物应对大气CO2浓度升高的生理和生态响应一直是热点问题[4-6],但是目前的相关研究较多集中于海洋浮游植物,而关于大气CO2浓度升高对淡水藻类生长的具体作用研究甚少[7]。

因此,本研究在云南滇池对微囊藻水华过程进行跟踪采样,分析微囊藻水华前、中、后和末期CCM基因型的动态变化特征,并在室内模拟低、中和高3种CO2浓度条件,分析不同CCM基因型的微囊藻对外界环境中CO2浓度变化的适应策略,旨在探究自然水体条件下CCM微进化与微囊藻种群演替的关系,进而为揭示环境CO2升高对微囊藻水华演替和发展以及水华形成机制提供重要的理论和实验依据。

1 材料与方法

1.1 研究区域概况与样品采集

滇池(24°40′～25°02′N,102°36′～102°40′E) 位于云贵高原,属亚热带高原季风气候。年平均气温15℃。平均水域面积310 km2,平均水深4.7 m,流域面积2920 km2。滇池被海埂大堤划分为草海和外海两部分。其中外海占湖区总面积的96.7%。本研究分别在外海的4个方位设置了4个采样点(图1)。分别是东大河(DDH)、观音山(GYS)、洛龙河(LLH)和生态所(STS)。于2021年5-11月的每月采集水样1～2次。使用水样采集器采集表层水样(距水面0.5 m),重复采样3次,混合后立即放入冰盒中于当天带回实验室做进一步分析。

图1 滇池采样点分布示意Fig.1 Distributions of sampling sites at Lake Dianchi

1.2 水体理化指标的测定

1.3 基因组DNA的提取

取250 mL混匀后的水样经0.22 μm滤膜抽滤,含藻滤膜于-80℃冰箱中冻存。滤膜剪碎后采用改良CTAB法提取DNA[11]。

1.4 扩增引物及荧光定量PCR标准曲线的构建

本研究所选用的引物如表1所示。其中引物对16S-F/16S-R和引物对sbtA-F2/sbtA-R2由Giovanni等设计[12-13],分别用于特异性扩增微囊藻的16S rDNA基因和sbtA基因。引物对bicA-F/bicA/R通过软件DNAMAN设计,用于特异性扩增微囊藻bicA基因。所有引物由生工生物技术有限公司合成。3个基因片段的PCR扩增产物经琼脂糖凝胶电泳后,对目的条带进行胶回收。纯化后的产物通过TA克隆与PMD-19T载体连接,并转入感受态E.coli细胞DH5α 中。挑取阳性克隆过夜摇菌,采用质粒小体试剂盒(天根,北京)提取质粒。将质粒线性化片段溶于ddH2O制备成标准品。采用微量分光光度计NanoDrop 1000(Thermo,美国)测定标准品中的质粒浓度。标准品基因拷贝数(copies/mL)=质粒浓度×10-9×6.02×1023/[质粒长度(bp)+目的片段(bp)×660]。将制备的标准品10倍梯度稀释构建标准曲线。3条标准曲线扩增效率分别为85.71%、94.61%和89.34%。

表1 本研究所用引物序列Tab.1 Primer sequences used in this study

1.5 荧光定量PCR

定量PCR反应体系为20 μL,其中含有TB GreenPremix Ex TaqII 10 μL、ROX Reference Dye II 0.4 μL,上下游引物各0.8 μL、DNA模板2 μL,超纯水6 μL。反应程序如下:95℃预变性2 min,95℃变性10 s,52℃退火30 s,72℃ 15 s,40个循环。每个样品重复测定3次。将得到的Ct值代入标准曲线计算样品中16S rDNA、bicA和sbtA基因的拷贝数。以微囊藻16S rDNA基因的拷贝数表征水体微囊藻总量,以bicA和sbtA基因拷贝数与16S rDNA拷贝数的比值表征两种基因型微囊藻的相对丰度。

1.6 藻种培养与室内实验设计

实验所用微囊藻藻株FACHB-908(sbtA型)来自中国科学院淡水藻种库,A5(bicA型)分离自三峡水库,纯化后的微囊藻藻种于MA培养基中进行扩大培养。培养温度(25±1)℃,光照强度30 μmol protons /(m2·s),光照周期12 h∶12 h(光照∶黑暗)。

实验将等生物量的两种CCM基因型的微囊藻分别在浓度为0.02%、0.04%和0.08%的CO2条件下进行混合培养。培养体积为400 mL,初始藻细胞接入量为OD680=0.07,每个处理设置3个重复。不同浓度的CO2通过空气及纯氮气以一定比例的流量混合实现。所有气体进入之前均经过0.22 μm无菌滤膜过滤,经复合气体分析仪检测校正CO2浓度后用于实验。第0、3、7、11和15天对各处理进行了取样,样品的DNA提取及定量PCR方法同前。

1.7 数据处理及分析

采用Excel和SPSS 24(IBM Inc., Chicago, IL, 美国)进行数据的统计分析,P<0.05被认为具有统计学意义,并使用Origin 2021(Origin Lab Inc.,Northampton, MA, 美国)进行绘图。不同基因型微囊藻和环境因子Mantel 检验分析在R 4.2.2中使用linkET包和dplyr包进行分析和绘制[14],其中Bray-curtis距离用于物种分析,欧几里得距离(Euclidean distance)用于环境因子分析。

2 结果

2.1 滇池各样点叶绿素a浓度的动态变化

为确定滇池水华的发生状况,测定了4个样点5-11月份的叶绿素a浓度。结果显示,滇池的叶绿素a浓度在时间上变化较大,且不同样点的变化不一致(图2)。根据叶绿素a浓度的变化,发现各样点均有水华的发生,其发生的时间和程度在不同样点间有着一定的差异,但最大值都超过了150 μg/L(图2)。根据Wu[15]的研究,将叶绿素浓度高于150 μg/L的月份划定为水华中期,8-9月为水华后期,10-11月为水华末期。

图2 不同时期滇池各采样点的叶绿素a浓度Fig.2 Chl.a concentration at different sampling points in Lake Dianchi at different periods

2.2 滇池不同样点无机碳浓度的动态变化

图3 水华过程中滇池各样点的无机碳浓度变化(a:东大河(DDH); b:观音山(GYS); c:洛龙河(LLH); d:生态所(STS))Fig.3 Changes in inorganic carbon concentration at different locations in Lake Dianchi during bloom

2.3 滇池不同基因型微囊藻在水华过程中的竞争

对滇池中不同基因型微囊藻鉴定的结果显示。在滇池的各样点中,不同基因型微囊藻的相对丰度在空间上具有一定的差异,并且在整个水华过程中sbtA基因型微囊藻和bicA基因型微囊藻的变化趋势相反。相同的是,sbtA基因型微囊藻在各样点中始终占据绝对优势,相对丰度远远高于bicA基因型微囊藻(图4)。在整个水华过程中,除STS外,滇池的bicA基因型微囊藻的相对丰度表现出了先降低后升高的趋势,DDH与LLH的结果相一致,最小值均出现在水华中期,且显著低于其它3个时期,水华后期显著低于末期(P<0.05,图4a和c);GYS的最小值出现在水华后期,与水华中期无显著差异,但显著低于水华前和末期(P<0.05,图4b);STS的变化趋势与其他样点不同,bicA基因型微囊藻的相对丰度逐渐降低,水华后、末期具有最小值,且显著低于前、中期(P<0.05,图4d)。对于sbtA基因型微囊藻,其相对丰度表现出的趋势则恰恰相反,DDH和STS表现一致,在水华中期达到最大并显著高于另外3个时期(P<0.05,图4a和d);GYS则表现为水华中、后期显著高于水华末期(P<0.05,图4b);LLH同样也是在水华中期达到最大值,且显著高于水华后期(P<0.05,图4c)。

2.4 理化因子对水华过程中不同基因型微囊藻的影响

为探究水体理化因子与bicA基因型和sbtA基因型微囊藻之间的关系,经Mantel 检验分析发现,在不同的水华时期,与两种基因型微囊藻显著相关的理化因子并不一致,且多表现出与sbtA基因型微囊藻显著相关(图5)。结果显示,在水华前期,bicA基因型微囊藻与Chl.a和TP的相关性系数较高,分别为0.59和-0.41,但并无显著性差异(图5a),sbtA基因型微囊藻则与Chl.a呈显著正相关(r=0.57,P<0.05,图5a);水华中期的bicA基因型微囊藻仅与Alk相关性较高但不显著(r=0.36,图5b),sbtA基因型微囊藻则与Alk呈显著正相关(P<0.05,图5b);在水华后期,bicA基因型微囊藻与TN显著正相关(P<0.05,图5c),而无理化因子与sbtA基因型微囊藻显著相关(图5c);水华末期与sbtA基因型微囊藻显著相关的理化因子有pH和CO2(aq)(P<0.05,图5d),并无与bicA基因型微囊藻显著相关的理化因子。对整个水华期两种基因型微囊藻和理化因子的Mantel 检验分析结果显示sbtA基因型微囊藻与pH和CO2(aq)显著正相关,bicA基因型微囊藻与Chl.a显著正相关(P<0.05,图5e);结果还表明CT与sbtA基因型微囊藻负相关,与bicA基因型微囊藻正相关。

图5 在不同水华期理化因子与两种基因型微囊藻的相关性分析(a、b、c、d、e分别代表水华前、中、后、末期、整个水华时期)Fig.5 Correlation analysis between physicochemical factors and two genotypes of Microcystis at different bloom stages

2.5 不同基因型微囊藻对不同浓度CO2的响应

为了进一步探究不同基因型微囊藻对CO2的响应,使用针对微囊藻的特异性引物对不同组合微囊藻的DNA进行荧光定量PCR。结果显示,在不同浓度CO2处理下bicA和sbtA两种基因型微囊藻的基因拷贝数在培养过程中均出现了明显分异(图6)。在低浓度CO2处理组中,sbtA型的微囊藻基因拷贝数从共培养开始至结束逐渐升高且始终显著高于bicA基因型微囊藻(P<0.001),具有绝对优势,而bicA基因型微囊藻一直保持比较低的数量(图6);在中浓度CO2处理组中,sbtA型的微囊藻随着共培养时间的增加,其基因拷贝数也增加,同样显著高于bicA基因型微囊藻(P<0.001)(图6);在高浓度CO2处理组中,两种基因型微囊藻的基因拷贝数均随着共培养时间的增加而增加,并且从第3天开始发生明显的分异,第3天以前是sbtA型微囊藻占多数,第3天以后则是bicA基因型微囊藻占绝对优势,其基因拷贝数要显著高于sbtA型微囊藻(P<0.05)。因此,随着CO2浓度的降低,处于竞争中的不同基因型微囊藻会发生从bicA基因型占优势到sbtA基因型占优势的转变。

图6 bicA 和sbtA 基因型微囊藻在不同浓度CO2处理下的竞争响应(a、b、c分别为两种基因型藻株在0.08%、0.04%、0.02% CO2浓度下的竞争;*表示0.01

3 讨论

CO2作为蓝藻光合作用的重要原料,对蓝藻尤其是具有不同CCM机制的微囊藻水华的形成具有重要的影响。滇池作为云南省内最大的淡水湖频繁暴发蓝藻水华[16-17],已经引起了广泛的关注。

3.1 滇池不同水华期CCM基因型微囊藻的动态变化

不同水华期滇池CCM微囊藻特异性引物荧光定量PCR结果显示(图4),滇池中所有采样点的sbtA基因型相对丰度均高于bicA基因型微囊藻,占绝对优势。这与曾佳颖等[20]对滇池不同CCM基因型微囊藻的研究结果一致。产生这种现象的原因可能与水体中低CO2浓度有关,已有研究发现低浓度CO2条件下更有利于sbtA基因型微囊藻获得竞争优势[21],促使bicA基因型向sbtA基因型微囊藻的转化[22]。

3.2 环境因子对不同CCM基因型微囊藻的影响

滇池水华前中后末期理化因素对不同CCM基因型微囊藻的相关性分析结果显示(图5),在野外条件下微囊藻的生长受到除无机碳外的其他理化因素的影响。在水华早期,Chl.a与bicA基因型正相关,与sbtA基因型显著正相关。由于Chl.a浓度能反映浮游植物的生物量[27],说明滇池中的浮游植物以微囊藻为主,尤其是sbtA基因型微囊藻(图5a),这与2.3节的结果中滇池sbtA基因型微囊藻在微囊藻中占绝对优势一致。

sbtA基因型微囊藻在水华中期与Alk显著正相关,bicA基因型微囊藻在水华中后期也与Alk相关程度很高,并且与TN显著正相关(图5b和5c)。偏碱性的环境和充足的氮源更有利于微囊藻的生长[19,28]。浮游植物在水华中期的大量生长导致了CO2浓度的大幅下降,CO2和氮磷营养盐等营养物质的同化增加了水华发展过程中的碱度和pH值[29-30]。pH值和碱度的增加使无机碳成分向碳酸氢盐转移[31],由于sbtA基因型微囊藻对碳酸氢盐具有高亲和力,因此在水华中后期sbtA基因型微囊藻的相对丰度高于水华早期和末期。

pH能够影响微囊藻的生长[25]和水体中的无机碳[26],曾佳颖等[20]对滇池的研究发现,pH与sbtA基因型的相对优势度呈正比,而与bicA基因型呈反比。而本研究发现这种现象出现在滇池水华末期,并且pH与sbtA基因型微囊藻呈正比(图5d)。在滇池水华末期水体中的pH均值(8.6)低于水华中后期(8.7和8.8),CO2浓度的增加引起了水体pH的下降[32],使水华末期sbtA基因型微囊藻丰度降低(图4)。同样的现象也出现在对滇微囊藻和理化因子的整体Mantel 检验分析中。同时在对滇池整体的分析中还发现CT的浓度和sbtA基因型微囊藻负相关,与bicA基因型微囊藻正相关(图5e)。这和先前研究发现sbtA基因型微囊藻在低CT条件下会更好的生长,但在较高的CT环境中则生长较差;反之,bicA基因型微囊藻在低CT环境下生长较差,但在高CT条件下表现良好的实验结果一致[20]。同时在研究中还发现水体中CO2浓度的升高同样也有利于sbtA基因型微囊藻的生长,可能的原因是CO2作为光合作用的重要原料,其浓度的升高对微囊藻的生长和繁殖具有促进作用。通过对两种基因型微囊藻和环境变量之间的相关性分析结果可知,不同CCM基因型微囊藻对不同无机碳的利用相对复杂,基本和先前的研究结果一致[12,21],但也出现像bicA基因型在生态所的4个水华时期比较特殊的情况,本研究认为这可能和环境因子、物种之间的交互作用有关。

3.3 不同CCM基因型微囊藻对不同CO2浓度响应

4 结论

云南滇池不同采样点之间空间差异较小。不同水华时期,sbtA基因型微囊藻始终占据绝对优势,随滇池水体无机碳浓度的变化促使微囊藻基因型呈现bicA→sbtA→bicA转变的趋势,使得具有CCM机制的微囊藻能够在生态系统中处于竞争优势。环境变量中,Chl.a浓度对水华早期的微囊藻有较强的相关性,总碱度和TN在水华中后期对两种基因型影响较大,而pH和CO2浓度对水华后期sbtA基因型产生了显著影响。室内模拟实验进一步验证了高浓度CO2适合bicA生长,而sbtA基因型对低浓度CO2具有更强的适应性。因此,CO2浓度对两种CCM机制微囊藻的演替起着重要的调控作用。