本发明公开了一种预测微藻重碳酸盐利用能力的方法，属于应对气候变化和海洋生物工程领域。不仅能够预测任意重碳酸盐浓度下微藻的重碳酸盐利用能力，而且还能够寻找微藻生长的最适重碳酸盐浓度。测定结果可以量化，具有可比性，为岩溶湖泊微藻碳汇的精确估算提供了科技支撑。
背景技术：
：微藻(microalgae)包括所有生活在水中营浮游生活方式的微小植物，通常就指浮游藻类。微藻结构简单，其生理过程也相对简单，有些种类是科学研究的模式植物，如：莱茵衣藻、小球藻，很多种类还可以人工培养，这为我们的研究提供了便利。微藻对水体无机碳的利用有两种方式，(1)利用大气中的二氧化碳。co2作为线性非极性分子，呈电中性，它可以自由扩散进入细胞双层脂膜，进入细胞中的co2为微藻细胞的光合作用所利用；(2)利用溶液中碳酸氢根离子。碳酸氢根离子既可以直接转运也可间接转运到细胞中为微藻细胞所利用。碳酸氢根离子的直接转运指的是经过细胞质膜表面载体蛋白或阴离子交换蛋白，直接把碳酸氢根离子转运到细胞内，在胞内经碳酸酐酶转化为co2或直接以碳酸氢根离子的形式由叶绿体膜蛋白主动运输到叶绿体内，经碳酸酐酶转化成co2供核酮糖-1,5二磷酸羧化/加氧酶(rubisco)固定；碳酸氢根离子的间接转运是指依赖于胞外碳酸酐酶的碳酸氢根离子的间接转运。水体中存在4种无机碳形式，它们分别是co2,hco3-,h2co3和co32-，这四种形式存在着如下平衡：因此，无论微藻利用二氧化碳还是利用碳酸氢根离子，它们都可能有两个来源，一个是来源于空气中的无机碳，另一个来源于溶液中固有的碳酸氢根离子。无论藻类采用二氧化碳利用途径还是采用碳酸氢根离子利用途径，只要来源空气的无机碳被同化，我们称为藻类直接碳汇，而无论藻类采用二氧化碳利用途径还是采用碳酸氢根离子利用途径，来源水体固有的无机碳被同化称为间接碳汇。直接碳汇直接移去大气中的二氧化碳，而间接碳汇通过移去水体固有的无机碳改变水体无机碳的平衡，间接移去大气中的二氧化碳。以前人们测定或估算的微藻碳汇都为直接碳汇，而间接碳汇为人们所忽视，更谈不上定量了。自然界中碳元素有两种稳定同位素：12c和13c，它们的天然平均丰度分别为98.89％和1.11％。稳定碳同位素组成通常用δ13c(‰)表示，自然界中δ13c的变化为-90‰～+20‰。稳定碳同位素的强烈分馏特征是识别微藻无机碳来源的基础。质量平衡原理以及同位素混合模型和化学计量学方法，是定量识别微藻无机碳来源的基础。喀斯特地区的岩石类型以碳酸岩盐为主，主要成分为caco3和mgco3等。基岩碳酸岩盐在岩溶作用下释放出大量的hco3-和ca2+,进入水体，致使岩溶湖泊水体呈现出高ph和高重碳酸盐的环境。低浓度的重碳酸盐可以促进微藻的生长和重碳酸盐的利用，高浓度的重碳酸盐可以抑制微藻的生长和降低重碳酸盐的利用份额。本发明通过构建微藻蛋白质增加倍数与重碳酸盐浓度、微藻利用重碳酸盐的份额与重碳酸盐浓度的方程，进一步耦合出重碳酸盐利用能力与重碳酸盐浓度的关系方程，。将待测的培养液重碳酸盐浓度带入上述方程中，即可预测出微藻在任意重碳酸盐浓度下的碳酸盐利用能力，为岩溶湖泊微藻的生产力和碳汇的精确估算提供了科技支撑。技术实现要素：本发明要解决的技术问题是，提供一种预测微藻重碳酸盐利用能力的方法，不仅能够预测任意重碳酸盐浓度下微藻的重碳酸盐利用能力，而且还能够寻找微藻生长的最适重碳酸盐浓度。克服了现有技术不能预测和评估微藻间接碳汇能力的缺陷。本发明采取以下技术方案：它包括以下步骤：步骤一，选择两种δ13c值差值大于8‰的碳酸氢钠作为同位素标记1和同位素标记2，设置不同浓度，分别添加到培养液中来培养待测微藻；同位素标记1的碳酸氢钠的δ13c值为δc1，同位素标记2的碳酸氢钠的δ13c值为δc2；步骤二，在培养的起始期测定培养时间为0时藻体蛋白质生物量和各培养液的重碳酸盐浓度c；步骤三，将待测微藻在被考察的培养条件下培养，培养4天后，收获藻体，分别测定添加不同浓度两种同位素标记的碳酸氢钠的培养液培养的相对应的各培养条件下的、被考察微藻的稳定碳同位素组成δ13c的值δt1、δt2；以及藻体微藻蛋白质生物量，计算处理后的微藻蛋白质生物量相对于接种时的增殖倍数，也即藻体蛋白质增加倍数p；步骤四，通过方程计算出微藻不同重碳酸盐浓度下利用重碳酸盐的份额fb；步骤五，以高斯方程拟合微藻蛋白质增加倍数p与重碳酸盐浓度c的方程，获得方程中的各参数值；以高斯方程拟合微藻蛋白质增加倍数p与重碳酸盐浓度c的方程为：其中，参数m指高斯曲线的峰值，c0为其对应的横坐标，n控制着"钟"的宽度；步骤六，以高斯方程拟合微藻利用重碳酸盐的份额fb与重碳酸盐浓度c的方程，获得方程中的各参数值；以高斯方程拟合微藻利用重碳酸盐的份额fb与重碳酸盐浓度c的方程为：其中，参数q指高斯曲线的峰值，c0为其对应的横坐标，r控制着"钟"的宽度；步骤七，进一步获取重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程，重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程为：步骤八，将待测的培养液重碳酸盐浓度带入上述重碳酸盐利用能力与重碳酸盐浓度的关系方程中，即可计算出微藻在待测重碳酸盐浓度下的碳酸盐利用能力。本发明的基本原理为：低浓度的重碳酸盐可以促进微藻的生长和重碳酸盐的利用，高浓度的重碳酸盐可以抑制微藻的生长和降低重碳酸盐的利用份额。微藻利用的无机碳源为空气的无机碳和培养液中添加的无机碳。因此，可以利用两端元的同位素混合模型获取微藻利用来自于空气的无机碳源和利用添加的无机碳源的份额。两端元的同位素混合模型可以表示为：δti＝δai-fbiδai+fbiδbi(i＝1，2，3，------)(1)这里δti为微藻的δ13c值，δai为假定为微藻完全利用空气的无机碳源时藻体的δ13c值，δbi为假定为微藻完全利用添加的无机碳源时藻体的δ13c值，fbi为该考察微藻利用添加的无机碳源所占的份额。很显然，只知道δti很难求出fbi，因此，本发明采用具有较大差异的δ13c值碳酸氢钠分别同时培养微藻，以稳定碳同位素双标记来识别微藻利用添加的无机碳源的份额。对于同位素标记1(i＝1)来说，方程(1)表示如下式：δt1＝δa1-fb1δa1+fb1δb1(2)这里δt1为用第一种已知δ13c值的碳酸氢钠培养的微藻藻体的δ13c值，δa1为假定为微藻完全利用空气的无机碳源时藻体的δ13c值，δb1为假定为微藻完全利用添加的无机碳源时藻体的δ13c值，fb1为该考察微藻利用添加的无机碳为碳源所占的份额。对于同位素标记2(i＝2)来说，方程(1)表示如下式：δt2＝δa2-fb2δa2+fb2δb2(3)这里δt2为用第一种已知δ13c值的碳酸氢钠培养的微藻藻体的δ13c值，δa2为假定为微藻完全利用空气的无机碳源时藻体的δ13c值，δb2为假定为微藻完全利用添加的无机碳源时藻体的δ13c值，fb2为该考察微藻利用添加的无机碳所占的份额。(2)和(3)两个方程中δa1＝δa2，fb＝fbi＝fb1＝fb2，联立求解(4)式中δb1-δb2则可以换算成同位素标记1的碳酸氢钠的δ13c值δc1与同位素标记2的碳酸氢钠的δ13c值δc2的差，则：因此，可以通过测定同位素标记1的碳酸氢钠的δ13c值δc1与同位素标记2的碳酸氢钠的δ13c值δc2，同时测定用对应的标记的碳酸氢钠培养的微藻δ13c值，即测定出δt1和δt2值，依(5)式计算出该考察微藻利用添加的无机碳所占的份额。微藻蛋白质增加倍数p与重碳酸盐浓度c的净光合速率与重碳酸盐浓度的关系可以用高斯方程表示。以高斯方程拟合微藻蛋白质增加倍数p与重碳酸盐浓度c的方程为：其中，参数m指高斯曲线的峰值，c0为其对应的横坐标，n控制着"钟"的宽度。微藻利用重碳酸盐的份额fb与重碳酸盐浓度c的方程可以用高斯方程表示。以高斯方程拟合微藻利用重碳酸盐的份额fb与重碳酸盐浓度c的方程为：其中，参数q指高斯曲线的峰值，c0为其对应的横坐标，r控制着"钟"的宽度。微藻重碳酸盐利用能力则是微藻利用重碳酸盐的份额fb与微藻蛋白质增加倍数p的乘积。因此重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程为：将待测的培养液重碳酸盐浓度带入上述重碳酸盐利用能力与重碳酸盐浓度的关系方程中，即可计算出微藻在待测重碳酸盐浓度下的碳酸盐利用能力。本发明的优点如下：1)本发明能定量预测能够预测任意重碳酸盐浓度下微藻的重碳酸盐利用能力。2)本发明能够便捷、精确定量微藻生长的最适重碳酸盐浓度。3)本发明能够定量研究重碳酸盐对微藻生长发育和碳汇作用的影响。4)本发明克服了现有技术不能预测和评估微藻间接碳汇能力的缺陷，为岩溶湖泊微藻的生产力和碳汇的精确估算提供了科技支撑。具体实施方式本发明的实施例：它包括以下步骤，步骤一，选择两种δ13c值差值大于8‰的碳酸氢钠作为同位素标记1和同位素标记2，设置不同浓度，分别添加到培养液中来培养待测微藻；同位素标记1的碳酸氢钠的δ13c值为δc1，同位素标记2的碳酸氢钠的δ13c值为δc2；步骤二，在培养的起始期测定培养时间为0时藻体蛋白质生物量和各培养液的重碳酸盐浓度c；步骤三，将待测微藻在被考察的培养条件下培养，培养4天后，收获藻体，分别测定添加不同浓度两种同位素标记的碳酸氢钠的培养液培养的相对应的各培养条件下的、被考察微藻的稳定碳同位素组成δ13c的值δt1、δt2；以及藻体微藻蛋白质生物量，计算处理后的微藻蛋白质生物量相对于接种时的增殖倍数，也即藻体蛋白质增加倍数p；步骤四，通过方程计算出微藻不同重碳酸盐浓度下利用重碳酸盐的份额fb；步骤五，以高斯方程拟合微藻蛋白质增加倍数p与重碳酸盐浓度c的方程，获得方程中的各参数值；以高斯方程拟合微藻蛋白质增加倍数p与重碳酸盐浓度c的方程为：其中，参数m指高斯曲线的峰值，c0为其对应的横坐标，n控制着"钟"的宽度；步骤六，以高斯方程拟合微藻利用重碳酸盐的份额fb与重碳酸盐浓度c的方程，获得方程中的各参数值；以高斯方程拟合微藻利用重碳酸盐的份额fb与重碳酸盐浓度c的方程为：其中，参数q指高斯曲线的峰值，c0为其对应的横坐标，r控制着"钟"的宽度；步骤七，进一步获取重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程，重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程为：步骤八，将待测的培养液重碳酸盐浓度带入上述重碳酸盐利用能力与重碳酸盐浓度的关系方程中，即可计算出微藻在待测重碳酸盐浓度下的碳酸盐利用能力。实施例：培养材料为：莱茵衣藻和蛋白核小球藻。基本培养液采用se培养基，基本培养条件为：光周期l/d：12h/12h；温度25℃；光照强度为100μmol·m-2·s-1，ph值8.0(用盐酸和氢氧化钠调节)。添加0、0.5、2.0、4.0、8.0、16.0mmol/l碳酸氢钠到培养液中，添加的碳酸氢钠的δ13c分别为-28.4‰(pdb)(δc1)和-17.4‰(pdb)(δc2)。在培养的起始期测定培养时间为0时藻体蛋白质生物量和各培养液的重碳酸盐浓度c(表1)；收获培养4天后藻体，分别测定藻体蛋白质生物量和藻体的δ13c值(表2)。用本发明方法，得出微藻利用添加的无机碳源的份额fb(表3)，构建出莱茵衣藻和蛋白质核小球藻蛋白质增加倍数与重碳酸盐浓度(p-c)、利用重碳酸盐的份额与重碳酸盐浓度(fb-c)的方程，获得的方程的参数(如表4)，进一步获取重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程，将重碳酸浓度为2.75、3、4.45、5、8、8.30、10、15.70mmol/l带入到重碳酸盐利用能力buc与重碳酸盐浓度c的关系方程中，获得莱茵衣藻和蛋白质核小球藻在重碳酸浓度为2.75、3、4.45、5、8、8.30、10、15.70mmol/l时重碳酸盐利用能力(buc)的预测结果如表5。表1不同浓度重碳酸盐处理下微藻蛋白质增加倍数(p)[hco3-]a[hco3-]b莱茵衣藻蛋白核小球藻mmol/lmmol/lpp0.000.653.51±0.274.24±0.220.501.554.00±0.304.31±0.222.002.754.11±0.184.51±0.234.004.454.26±0.205.28±0.178.008.304.97±0.234.67±0.2116.0015.704.56±0.294.64±0.23a向培养液中添加的碳酸氢钠浓度b培养液中起始期实测的碳酸氢根离子(重碳酸盐)浓度p为处理后的微藻生物量相对于接种时的增殖倍数，即藻体蛋白质增加倍数表2不同浓度重碳酸盐处理下的微藻碳同位素组成b培养液中实测的起始期碳酸氢根离子(重碳酸盐)浓度δt1:添加δ13c为–17.4‰的nahco3的培养液所培养出的微藻δ13c值δt2:添加δ13c为–28.4‰的nahco3的培养液所培养出的微藻δ13c值表3不同浓度重碳酸盐处理下微藻利用重碳酸盐的份额(fb)b培养液中起始期实测的碳酸氢根离子(重碳酸盐)浓度fb为微藻利用重碳酸盐的份额表4莱茵衣藻和蛋白质核小球藻蛋白质增加倍数与重碳酸盐浓度(p-c)、利用重碳酸盐的份额与重碳酸盐浓度(fb-c)的方程参数微藻种类方程类型参数m(q)参数n(r)参数c0莱茵衣藻p-c4.988412.315110.5390莱茵衣藻fb-c0.37254.803211.1324蛋白核小球藻p-c4.991915.69879.2193蛋白核小球藻fb-c0.50954.251911.0310表5莱茵衣藻和蛋白质核小球藻在重碳酸浓度为2.75、3、4.45、5、8、8.30、10、15.70mmol/l时重碳酸盐利用能力(buc)预测结果和实测结果本发明的实施效果如下：从表4中可以看出，莱茵衣藻和蛋白质核小球藻的生长最适重碳酸盐浓度分别为10.5390mmol/l和9.2193mmol/l，利用重碳酸盐份额最大的重碳酸盐浓度分别为11.1324mmol/l和11.0310mmol/l，数值均很接近，表明本发明具有较好的可信性。从表5中可以看出，从重碳酸盐2.75mmol/l到15.7mmol/l范围均具有很好的预测能力，预测误差为0.0051到0.0823buc，表明本发明可以很好地预测岩溶湖泊中微藻的重碳酸盐利用能力，为岩溶湖泊微藻生产力和碳汇的精确估算提供了科技支撑。当前第1页12