生物标志物梯烷脂指示胶州湾厌氧氨氧化活动❋

史明炜， 肖晓彤❋❋， 刘 珂， 李 莉， 王玉珏， 赵美训

(1.中国海洋大学深海圈层与地球系统前沿科学中心和海洋化学理论与工程技术教育部重点实验室， 山东 青岛 266100;2.青岛海洋科学与技术试点国家实验室 海洋生态与环境科学功能实验室， 山东 青岛 266237;3.华东师范大学河口海岸学国家重点实验室，上海 200241)

沿海区域是氮元素海洋地球化学循环的重要场所，同时也被认为是厌氧氨氧化作用的重要地区[13-14]。尤其是在较易受人类活动影响的海湾区域，沿海人口增加及土地利用扩张导致的营养盐荷载输入增加为厌氧氨氧化活动提供了充足的反应底物；同时，营养盐增加会引发赤潮和绿潮等生态灾害，可能导致水体低氧区域增多，为厌氧氨氧化反应提供重要的环境条件。目前厌氧氨氧化作用主要的分析方法有生物标志物法、同位素标记法和分子信息学。在边缘海区域，前人主要使用氮同位素标记法和分子信息学方法对厌氧氨氧化活动进行研究。如Song等[15]通过同位素标记法发现长江口有13%～50%的固定氮流失来源于厌氧氨氧化反应，Hou 等[16]使用生物信息学方法发现随盐度的增加，长江口沉积物中厌氧氨氧化种群结构由CandidatusBrocadia为主转变为CandidatusScalindua为主。以上两种方法主要用于研究正在进行的厌氧氨氧化活动，而使用生物标志物法可以区分厌氧氨氧化菌的来源并获得较长时间尺度的厌氧氨氧化活动记录。

梯烷脂是厌氧氨氧化菌产生的特征生物标志物，可以指示海洋水体及沉积物中的厌氧氨氧化活动。梯烷脂主要存在于厌氧氨氧化菌中的厌氧氨氧化体(Anammoxosome)上，核心结构为3～5个链接起来的环丁烷基团(见图1)。这种膜脂结构特殊，为厌氧氨氧化菌特有，且在自然条件下较为稳定，易于辨识，指示性专一。通过对海洋沉积梯烷脂进行分析，能有效地指示该区域厌氧氨氧化活动的发生和相关菌种的来源，并进行氮循环历史记录的回溯[17-22]。国内对梯烷脂的研究主要集中在河口的典型缺氧区域， Hu等[21]通过对长江口水体颗粒物的分析，指出梯烷脂主要来自于水体原位产生；Zhao等[17]对长江口沉积柱的研究发现梯烷脂能较好地在百年尺度上重建长江口-东海区域的缺氧及富营养化历史。但目前鲜有在与人类活动密切相关的海湾的梯烷脂研究。

图1 核心梯烷脂结构示意图Fig.1 The structure of ladderane fatty acids

胶州湾位于中国山东半岛南部黄海之滨，形成于全新世时期黄海对山东半岛的海侵运动，是典型的半封闭-封闭型海湾。其海域面积367.4 km2，但湾口入黄海最窄处仅有3.3 km。胶州湾除了与黄海进行物质交换，还受到陆源输入和人类活动排放的强烈影响。海湾西部有大沽河、洋河，东部有墨水河、白沙河、李村河等河流汇入，河流入湾总径流量约9×108m3[23]。胶州湾周围设置了数十个污水排放口，主要集中在东岸。海湾内水文状况主要受潮流和余流的影响，加之河流径流的输入，导致海湾水文状况较复杂，整体易受人类活动排放的影响。胶州湾被中国北部沿海的大型城市青岛所环绕，其海域状况与青岛市海洋经济发展水平、污水排放及城市化程度密切相关。已有的研究表明，胶州湾水质状况从20世纪末呈明显的恶化趋势，但在21世纪初随着海洋国家自然保护区的划定有所改善。营养盐浓度在2002年左右呈中度营养和富营养化水平，随后海水富营养化有所缓解[24]。溶解氧浓度从1998年的9 mg/L迅速下降至2002年的5 mg/L，随着随后十年中水质的改善，溶解氧浓度得以上升[25]。本研究使用生物标志物法，运用梯烷脂从有机碳循环的角度研究胶州湾的厌氧氨氧化活动。本文主要通过生物标志物和总有机质指标(TOC, TOC/TN，δ13C)来研究沉积物中陆海源有机碳的贡献，通过梯烷脂的分布来讨论厌氧氨氧化活动的来源，并探讨厌氧氨氧化活动对海洋环境因素的响应与指示作用。

1 研究区域及方法

1.1 样品采集

(蓝色圆点为采样站位。Blue points represent sampling sites)图2 研究区域及采样站位Fig.2 Map of study area and sampling sites

1.2 沉积物粒度及总有机碳分析

有机质会附着在沉积颗粒表面，使沉积颗粒的粒径增加，为保证绝大多数有机质能被去除，又不至于高温破坏颗粒结构，将冻干沉积物样品置于马弗炉350 ℃高温去除有机碳(12 h)。降至室温后加入分散剂(六偏磷酸钠溶液，1 mol/L)，随后使用粒度自动测定仪(Mastersizer 3000)进行测定，测定结果取平均粒度进行分析。

称取1 g左右冻干沉积物并研磨均匀，使用4 mol/L的盐酸除去沉积物样品中的无机碳组分后将样品在55 ℃的烘箱中烘干，平衡24 h后, 使用CHN元素分析仪(Thermo Flash EA2000)测定TOC和TN。得出原始数据后进行校正以避免酸洗过程中的损失对TOC和TN的影响。校正方法如下：

酸洗前称量冻干沉积物总质量m1，酸洗后称量冻干沉积物总质量m2，酸洗后实测样品TOC含量为a。求得校正后的样品TOC含量为 TOC=a×m2/m1；TN校正方法同TOC。

TOC和TN的标准品为Atropina Standard(TOC=70.56%)和Low Organic Content Soil Standard OAS(TOC=(1.61±0.09)%)，TOC与TN标准偏差使用标准Low Organic Content Soil Standard OAS测得，分别为0.02% (n=6)及0.002% (n=6)。

上述预处理沉积物使用稳定同位素比质谱仪(Elemental Analyzer-Isotope Ratio Mass Spectrometer, EA-IRMS)在连续流模式下进行有机碳的δ13C测定，使用稳定同位素比质谱仪测定δ13C的标准品为USGS-40(δ13C =-26.39‰), IAEA-600(δ13C =-27.77‰)和IAEA-CH3(δ13C =-24.72‰),标准偏差使用标准USGS-40(δ13C =-26.39‰)测得为0.6‰(n=6)。

标准偏差公式如下：

式中：Xi为标准品实测值；X为测定均值；n=6。

1.3 生物标志物分析

本研究计算陆源有机质贡献指标所使用的类脂生物标志物主要包括正构烷烃、甾醇和支链甘油二烷基甘油四醚(br-GDGTs)，使用萃取分离法进行富集。将冻干并研磨均匀的沉积物样品加入内标(C19-醇, C24-氘烷和C46-GDGT)，使用二氯甲烷∶甲醇(3∶1)溶剂进行超声萃取。萃取后混合液使用氮吹富集，并用6%氢氧化钾-甲醇溶液水解。水解12 h后使用正己烷萃取，先使用正己烷洗脱正构烷烃组分，再使用二氯甲烷∶甲醇(95∶5)洗脱甾醇和GDGTs组分。正构烷烃和甾醇组分使用气相色谱仪(FID-GC，Agilent 7890N，色谱柱为HP-1毛细管柱)进行测定，GDGTs组分使用高效液相色谱-质谱联用仪(HPLC-MS，Agilent1200/Waters-Quattro Ultima Pt mass，色谱柱为Prevail Cyano 正相色谱柱)进行测定。

本研究使用梯烷脂作为厌氧氨氧化的标志物，使用萃取分离法进行富集。冻干并研磨均匀的沉积物样品加入二氯甲烷∶甲醇(2∶1)溶剂进行萃取，然后使用活化铜粉除去残余的硫，氮吹富集后用6%氢氧化钾-甲醇溶液水解。水解12 h后使用正己烷萃取4次，合并萃取液获得中性组分。在剩余的混合液中加入6 mol/L盐酸调节pH，继续用正己烷萃取，得到的萃取液经无水硫酸钠脱水后加入三氟化硼-甲醇溶液甲基化处理。甲基化后的脂肪酸甲酯(FAMEs)经过正己烷的萃取后加入内标2-羟基-4-甲氧基二苯丙酮，使用高效液相色谱-质谱联用仪(HPLC-MS，Agilent1200/Waters-Quattro Ultima Pt mass，色谱柱为Zorbax Eclipse XDB-C8反相键合色谱柱)进行测定。

1.4 TMBR、BIT和NL5的计算

相对陆源指标TMBR[28]和BIT[29]的计算公式分别如下：

式中：∑(C27+ C29+ C31n-alkanes)为陆地高等植物生物标志物长直链烷烃含量；∑(B + D)为海洋浮游植物生物标志物；B代表菜子甾醇含量；D代表甲藻甾醇含量。

式中：Ia + IIa + IIIa代表不同的支链甘油二烷基甘油四醚：br-GDGT-I、 br-GDGT-II、br-GDGT-III含量之和，IV代表海洋中泉古菌醇含量。

梯烷脂的NL5指标与厌氧氨氧化菌所生活的温度相关[30]，公式如下：

式中：ladderane I代表甲酯化处理后的梯烷脂C18-[5]-ladderane FAME(I)含量；ladderane III代表甲酯化处理后的梯烷脂C20-[5]-ladderane FAME(III) 含量。

2 结果

2.1 胶州湾水体及沉积物的环境参数

胶州湾夏季底层水体温度差异较小，为26.0～27.5 ℃，平均26.7 ℃，胶州湾东北角及西部海域水温较高(见图3(a))，是由于这些区域水深较浅，底层水易受太阳辐射影响。底层水盐度为26.9～31.6，平均为30.8，胶州湾东北角及洋河河口附近盐度较低，其余站位盐度均为31.0附近(见图3(b))。胶州湾底层水三种无机氮盐(亚硝氮、硝氮、氨氮)空间分布上较为类似，主要的高值区出现在湾东北部及西部河口区域(见图3(c)，3(d)，3(e))。底层水体溶解氧含量为5.2～8.3 mg/L，平均6.3 mg/L，低值区同样集中于胶州湾东北角及西部近岸河口区域(见图3(f))。

((a)胶州湾底层水体温度(℃)； (b)胶州湾底层水体盐度； (c)胶州湾底层水体含量(μmol/L)； (d)胶州湾底层水体含量(μmol/L)； (e)胶州湾底层水体含量(μmol/L)；(f)胶州湾底层水体溶解氧含量(mg/L)。 (a)Bottom sea water temperature (℃) of Jiaozhou Bay; (b)Bottom sea water salinity (psu) of Jiaozhou Bay; (c)Bottom content (μmol/L) of Jiaozhou Bay; (f)Bottom DO content (mg/L) of Jiaozhou Bay.)图3 胶州湾海底水体参数Fig.3 Bottom sea water parameters in Jiaozhou Bay

2.2 总有机碳参数及粒度

胶州湾表层沉积物中的TOC含量为0.06%～0.94%，平均0.48%，高值区集中在胶州湾东部近岸处(见图4(a))。表层沉积物总有机质δ13C值为-22.4‰～-19.7‰，δ13C值普遍偏高且极差较小，但是空间分布上呈较明显的西高东低的趋势，最低值出现在胶州湾东北角及湾南部出口处(见图4(b))。表层沉积物的碳氮比TOC/TN为3.7～12.3，平均8.4，碳氮比普遍偏低，空间上呈西低东高的趋势(见图4(c))。

表层沉积物的平均粒径为24～187 μm，大部分底泥为粉砂质(Silt)，而粒径较高的砂质(Sand)区集中在三个区域，分别是海湾东北角，西部大沽河/洋河河口，及胶州湾南部湾口处(见图4(d))。

((a) 沉积物TOC含量(%)；(b) 沉积物TOC的δ13C值(‰)；(c) 沉积物碳氮比(mol/mol)；(d) 沉积物平均粒径(μm)。 (a) TOC of sediments (%); (b) TOC-δ13C value of sediments (‰); (c) TOC/TN ratio of sediments (mol/mol); (d) mean grain size of sediments (μm).)图4 胶州湾沉积物总有机质参数及平均粒径Fig.4 Bulk organic matter and mean grain sizes of sediment samples

2.3 生物标志物和比值指标

胶州湾表层沉积物生物标志物的空间分布上，指示陆源有机质贡献的指标TMBR和BIT都在胶州湾东北部及西部比值较高(见图5(a)，5(b))，TMBR比值指标为0.23～0.69，平均为0.46，BIT比值指标为0.10～0.29，平均0.17。

((a) 沉积物有机质TBMR指标，数值由长链烷烃及甾醇含量比值得出，代表陆源有机质的贡献比例；(b) 沉积物有机质BIT指标，数值由不同GDGT及泉古菌醇含量比值得出，代表陆源有机质的贡献比例。 (a) Distribution of TMBR index, which is calculated by contents of long chain n-alkanes, brassicasterol and dinosterol; (b) Distribution of BIT index, which is calculated by contents of br-GDGTs and crenachaeol.)图5 胶州湾沉积有机质陆源贡献指标Fig.5 Distribution of terrestrial OC contribution indices

胶州湾所有站位的表层沉积物中都检出了梯烷脂，含量为4～574 μg/g OC，平均含量为163 μg/g OC。在胶州湾西部含量较低，高值站位主要分布胶州湾东部120.3°E附近及湾口处(见图6(a))。本研究聚焦四种梯烷核心脂，C18-[5]-ladderane FAME(I)，C18-[3]-ladderane FAME(II)，C20-[5]-ladderane FAME(III)和C20-[3]-ladderane FAME(IV)。胶州湾表层沉积物中C20-[5]-ladderane FAME(III)和C20-[3]-ladderane FAME(IV)为主要组分，分别占同站位总梯烷脂含量的23%～66%(平均38%)和24%～77%(平均56%)，且不同组分梯烷脂的空间分布和总梯烷脂空间分布状况基本一致。

(根据山东省生态环境厅发布的省控重点污染源排放监控数据库，红点代表胶州湾周围污水主要输入；蓝点代表胶州湾沿岸河流主要输入；橙点代表河流-污水混合输入。 (a) ∑ ladderane/TOC含量区间分布图(μg/g OC)；(b) ladderane I/TOC含量分布图(μg/g OC)；(c) ladderane II/TOC 含量分布图(μg/g OC)；(d) ladderane III/TOC含量分布图(μg/g OC)；(e) ladderane IV/TOC含量分布图(μg/g OC)。 Red dots represent main sewage input, according to the emission monitoring database of key pollution sources released by Shandong Provincial Department of Ecology and Environment. Blue dots represent main river water input of Jiaozhou Bay. Orange dots represent mixing inputs of sewage and river water. Red dots represent main sewage input around Jiaozhou Bay. (a) Distribution of summed TOC-normalized ladderane contents (μg/g OC); (b) Distribution of TOC-normalized ladderane I contents (μg/g OC); (c) Distribution of ladderane II contents (μg/g OC); (d) Distribution of TOC-normalized ladderane III contents (μg/g OC); (e) Distribution of TOC-normalized ladderane IV contents (μg/g OC).)图6 胶州湾沉积物梯烷脂含量分布图Fig.6 The spatial distribution of different ladderane FAMEs contentsin sediments

3 讨论

3.1 胶州湾有机质来源与贡献

利用梯烷脂研究海湾厌氧氨氧化活动的关键是厘清梯烷脂是由海洋产生还是来自于陆源输入，因此首先对胶州湾沉积物中有机质的来源进行判断与讨论，并后续建立有机碳与梯烷脂的相关关系，进一步对梯烷脂的来源及响应进行分析。

总有机质参数碳氮比值TOC/TN及稳定碳同位素δ13C能从总有机碳角度对有机质的来源进行分析。胶州湾沉积物有机质TOC/TN值与δ13C值存在良好的相关关系(r2= 0.82)(见图7(b))，表明在胶州湾沉积有机质中有机质来源的影响因素相似，可以通过两者来有效判断有机质来源。碳氮比高值区(TOC/TN>9，mol/mol)和δ13C低值区(δ13C <-21‰)明显集中于海湾东部沿岸区域，表明与西北部区域相比，东部沿岸区域沉积物中陆源有机质的贡献更多。一般海源浮游植物碳氮比约为5～7(mol/mol)，接近Redfield比值6.6(mol/mol)，而陆地高等植被的碳氮比高达15～1 000(mol/mol)[31]。海源浮游植物δ13C值一般在-23‰～-12‰[32]，陆地高等植物的δ13C值则更低，为-33‰～-23‰[32]。胶州湾沉积物中碳氮比值平均为8.4(mol/mol)， δ13C值平均为-21.4‰。总的来说，胶州湾沉积物有机质TOC/TN值较低、δ13C值较高，表明在胶州湾沉积有机质中海源初级生产力可能占主导，原因可能为近年来胶州湾附近各河流中下游设置较多的拦水坝，使得陆源物质如高等植物碎屑较难直接输入到海湾中。TOC-TN的拟合直线显示两者相关性良好，斜率为8.38(高于Redfield比值)，且回归线落在TN轴的正截距上(见图7(a))，表现出沉积物对无机氮的吸附作用[33]，可能与胶州湾陆源输入较多的无机氮盐相关。胶州湾西北大沽河口区域δ13C值较高(见图4(b))，可能受沿海河口区域C4植物(δ13C值-16‰～-9‰)的影响[34]，另外，土壤来源有机质及农肥也可能是导致海湾地区δ13C值偏高的原因[35-36]。综合而言，总有机质参数表明胶州湾沉积有机质以海源为主导。

((a) 沉积物TOC和TN线性关系图；(b) 沉积物总有机质的δ13C和TOC/TN线性关系图。 (a) Plot of sediment TOC and TN; (b) Plot of sediment TOC-δ13C and TOC/TN.)图7 胶州湾沉积物总有机质参数线性关系图Fig.7 Plot of sediment OC parameters

考虑到胶州湾总有机质参数难以精确指示不同来源有机质贡献程度，因此使用生物标志物指标对不同来源有机质贡献进行进一步解析。为了避免沉积物本身性质(如粒度等)对有机质保存的影响，可使用相对陆源指标TMBR和BIT衡量湾底沉积物中海/陆源有机质的贡献比例。整体来说，TMBR(平均值0.46)和BIT值(平均值0.17)都表明胶州湾沉积物中海源有机质占主要地位，其贡献约为54%(TMBR得出)和83%(BIT得出)。TMBR指标和BIT指标所指示陆源相对贡献的空间分布类似，陆源贡献较高的区域位于海湾东北部和西部，主要受河流输入的影响。湾东北部区域陆源贡献指标尤其高，虽然此处附近的河流流量较小，仅占胶州湾整体河流入海通量的不到20%[23]，但陆源贡献指标明显高于西部河口区域(见图5(a)，5(b))。原因可能是海湾东北部区域较小且较狭窄，不同于西部河口直接面向开阔海域，海湾东北部多条河流输入的陆源有机质能较多的聚集于此区域。在本文研究区域，BIT值普遍低于TMBR值，这可能是由于TMBR所使用的生标(长链烷烃)与BIT使用的生标(br-GDGTs)在传输机制及存在方式上有所不同。BIT指标使用的br-GDGTs主要由土壤细菌产生， 只能指示河流输入土壤对海湾沉积有机质的贡献；同时br-GDGTs在河口区域的降解也可能是BIT值偏低的原因[37]；而TMBR指标使用的长链烷烃存在于高等植物碎屑中，除了随河流输入海湾，还可以通过大气沉降输送入海[38]。

3.2 胶州湾梯烷脂的成分与来源

胶州湾表层沉积物中普遍检出梯烷脂。由于胶州湾沉积物平均粒径的空间分布存在较大差异，为了部分排除粒径对有机质保存产生的影响，因此在讨论梯烷脂的空间分布时使用总有机碳均一化后的梯烷脂含量(ladderane/TOC)来分析。前人研究指出沉积物中的梯烷脂主要来自于活体的厌氧氨氧化细菌，也可能来自于已经凋亡或失活的厌氧氨氧化细菌[39]。与其他研究方式不同，沉积物中的梯烷脂不仅能代表正在进行的厌氧氨氧化活动，还能记录曾经进行的厌氧氨氧化活动。通过追溯沉积物中梯烷脂的来源并分析其成分，可以推断梯烷脂是由外部(包括河流和污水口)输入海湾还是由海洋原位的厌氧氨氧化活动产生。因此，运用梯烷脂能更好地对胶州湾脱氮的潜力及贡献进行评估。

胶州湾沉积梯烷脂空间分布不均匀，海湾西部梯烷脂含量普遍偏低，高值区主要位于胶州湾东北角，北岸的红岛附近，以及团岛和前湾港附近的站位(见图6(a))。对胶州湾沉积梯烷脂与陆源有机质贡献指标进行相关性比对，结果表明梯烷脂与陆源贡献指标均无显著线性关系(见图8)。这说明沉积物中的梯烷脂不单受陆源或海源主导，需要考虑不同方式的混合对梯烷脂分布的影响。

((a)沉积物BIT指标与总梯烷脂含量线性关系；(b) 沉积物TMBR指标与总梯烷脂含量线性关系。 (a) Plot of BIT index and TOC-normalized ladderane contents; (b) Plot of TMBR index and TOC-normalized ladderane contents.)图8 胶州湾沉积物有机质陆源贡献指标与总梯烷脂含量线性关系图Fig.8 Plot of terrestrial OC contribution indices versus summed TOC-normalized ladderane contents

研究表明，在不同温度下，厌氧氨氧化菌会通过改变合成的梯烷脂链长来保持细胞膜的流动性[30]，一般来说厌氧氨氧化菌所处环境温度越高，膜脂上所合成的链长较长的ladderane III含量占比会增加，链长较短的ladderane I则减少，而基于ladderane III和ladderane I比值的温度指标NL5值越大。根据NL5-T经验公式的函数性质，NL5值反演梯烷脂合成温度的最适区间为12～20 ℃[19]，反演温度的NL5最大值为0.9，对应的温度约为25 ℃。胶州湾所有站位表层沉积物中梯烷脂NL5值均大于0.9，无法进行具体温度反演计算，可能是较高的合成温度使得厌氧氨氧化菌膜脂上所合成的ladderane III含量远高于ladderane I。根据胶州湾沉积物NL5的数值，定性上可以判断梯烷脂合成的环境温度高于25 ℃，与胶州湾夏季平均水温26.7 ℃较为相符。同样，祝孟玲等[54]和Zhao等[17]分别对珠江口和东海陆架沉积物中的梯烷脂来源进行了分析，并通过NL5值反演温度推断沉积物中梯烷脂可能主要来自于上层水体颗粒物沉降。另外，虽然海洋沉积物中存在厌氧氨氧化活动[27,39,49]，但胶州湾沉积物夏季温度约为21 ℃[51]，明显低于夏季水体温度，与梯烷脂较高的NL5值不符，且夏季表层沉积物间隙水并未明显缺氧(平均4.5 mg/L)[27]，因此胶州湾沉积物中的梯烷脂可能主要来自于夏季水体颗粒物的沉降。但由于NL5-T经验公式的温度最佳适用范围有限(12～20 ℃)，并不能据此完全排除胶州湾沉积物原位产生梯烷脂的可能性。

厌氧氨氧化菌在海洋水体中主要附着于颗粒物上进行厌氧氨氧化活动[21]。厌氧氨氧化菌主要有五种不同属，分别为Candidatusbrocadia、Candidatusscalindua、Candidatuskuenenia、Candidatusjettenia、Candidatusanammoxoglobus，不同属厌氧氨氧化菌种有不同的生理性质及形态特征，会产生不同种类的梯烷脂，因此可以分析各成分梯烷脂的含量及比例来探究胶州湾沉积物中梯烷脂的不同来源。近海及河流环境中主要发现的菌群为Candidatusscalindua，Candidatusbrocadia和Candidatuskuenenia[26-27,40]。其中Scalindua为高盐种，最常见于大洋环境；而Brocadia及Kuenenia为低盐种，主要存在于污水处理反应器及河流/河口区域[41]。培养试验的结果显示，菌种Scalindua的梯烷膜脂中ladderane IV为单独主导，ladderane III/ladderane IV约为0.6；菌种Brocadia和Kuenenia中ladderane III和IV共同主导，ladderane III/ladderane IV均为1.0左右[42]。胶州湾中各站位梯烷脂含量差异较大且空间分布不均匀，因此将梯烷脂高值站位(>200 μg/g OC)与低值站位的ladderane III/ladderane IV进行双样本t检验，以分辨不同站位中梯烷脂的来源。在排除异常值后，高值站位与低值站位ladderane III/ladderane IV存在显著性差异(p=0.001)。具体来说，梯烷脂高值站位中的ladderane III/ladderane IV较高(见图9(a))，平均达到了0.9(见图9(c))，比例更接近淡水菌种Brocadia 和 Kuenenia；而低值站位中ladderane III/ladderane IV平均值为0.7，排除异常值后平均值为0.6(见图9(c))，比例更接近高盐种厌氧氨氧化菌Scalindua。需要注意的是，梯烷脂高值站位也能检出Scalindua的特有膜脂ladderane II[42]，因此这些站位可能同样存在Scalindua菌种。以上结果表明，在梯烷脂低值站位，高盐种厌氧氨氧化菌Scalindua可能是梯烷脂的主要贡献者，即梯烷脂主要由海水原位产生；而在梯烷脂高值站位，除了Scalindua的贡献之外，还可能有较多来自于淡水排放菌种Brocadia和Kuenenia的贡献，即梯烷脂主要来自于河流/污水输送的颗粒物。空间分布上，高值站位也集中位于胶州湾沿岸的污水处理厂附近(见图6(a))。而湾西部的河流(大沽河、洋河)虽然能贡献胶州湾80%以上的河流入海流量及70%以上的泥沙入海通量[23]，但相比于人口密集的东岸，西部河流流域受人类排放的影响较少[48]，且沉积物平均粒径较高不利于生物标志物的保存，因此梯烷脂在湾西部无明显的高值站位。未来在典型区域的端元(如污水端，河流淡水端及现场生产悬浮颗粒物)以及不同季节对梯烷脂进行进一步的结构及分布研究能更好的校验梯烷脂对海湾区域厌氧氨氧化活动的指示意义。

((a) 梯烷脂低值站位中(<200 μg/g OC)的不同组分含量占比； (b) 梯烷脂高值站位中(> 200 μg/g OC)的不同组分含量占比； (c) 梯烷脂低值及高值站位ladderane III与ladderane IV比值。(a) Compositions of ladderane FAMEs in sites with lower ladderane contents(< 200 μg/g OC); (b) Compositions of ladderane FAMEs in sites with higher ladderane contents(> 200 μg/g OC) ; (c). ladderane III/ladderane IV ratio in sites with lower ladderane contents and sites with higher ladderane contents.)图9 不同站位梯烷脂不同组分含量占比箱式图Fig.9 Box plot of ladderane FAMEs in different sites

3.3 胶州湾沉积物梯烷脂的环境指示意义

目前，厌氧氨氧化活动已经在不同环境中被广泛发现，包括大洋、近海、淡水、河口和污水等。根据青岛市海洋环境公报，青岛市夏季(8月)海水水质状况为一年中最佳，主要表现在无机氮盐平均浓度较低，富营养化海域面积较小且主要集中在东部沿岸，与本次采样数据相吻合。若使用分子信息学方法研究沉积物的厌氧氨氧化活动，由于采样时间和季节的不同，沉积物中存在凋亡或非活性态的厌氧氨氧化菌，可能因此引入误差，错估了研究海域中厌氧氨氧化对氮脱除的贡献。胶州湾沉积物沉积速率为0.28～0.64 cm/a[45]，因此本研究采集的表层沉积物(0～2 cm)反映的是近几年的沉积状况。梯烷脂作为百年尺度上稳定的生物标志物，在确定其来源后，能够反映研究海域的厌氧氨氧化活动状况(包括已经进行的和正在进行的；河流/人为排放的和海水中原位合成的)，从而更好地对氮循环及脱氮过程进行评价，并综合指示海洋中的厌氧氨氧化状况。

表1 胶州湾海底环境各参数皮尔森相关关系abTable 1 Pearson correlations of different parameters in Jiaozhou Bay bottom environment

((a)主成分分析载荷图；(b)主成分分析得分图，其中红色圆点代表梯烷脂高值站位，黑色方点代表梯烷脂低值站位。 (a) Loadings plot from principle component analysis; (b) Scores plot from principle component analysis, red dots represent higher ladderane contents sites and black squares represent lower ladderane contents sites.)图10 胶州湾沉积物梯烷脂及其他参数的主成分载荷图及得分图Fig.10 Loadings plot and scores plot from principle component analysis (PCA) containing ladderanes and other parameters

4 结论

(1)胶州湾表层沉积物有机质中海源贡献普遍大于陆源贡献，部分区域的高陆源有机质贡献可能受河流输入的影响。

(2)沉积物中梯烷脂含量的空间分布及其各组分的比例表明，胶州湾区域存在广泛的厌氧氨氧化活动，且主要夏季水体颗粒物的沉降而非在沉积物中产生。大部分站位沉积梯烷脂主要来自于海水水体中高盐种Candidatusscalindua，在梯烷脂高值站位也有菌种Candidatusbrocadia和Candidatuskuenenia的贡献，主要来源于河流/污水输入的颗粒物。

致谢：本文相关样品的采集由马海力、朋鹏和丁杨提供协助，样品总有机质参数的测定由张海龙和金贵娥给予帮助，作者对以上人员表示诚挚谢意。