Marine Biology

Effects of salinity on the growth and transcriptome of euryhaline Synechococcus sp. K1*

  • LIAO Ying , 1, 2 ,
  • XIA Xiaomin , 1, 2
Expand
  • 1. Key Laboratory of Tropical Marine Bio-resources and Ecology (South China Sea Institute of Oceanology, Chinese Academy of Sciences), Guangzhou 510301, China
  • 2. University of Chinese Academy of Sciences, Beijing 100049, China
XIA Xiaomin.

Received date: 2021-04-06

  Revised date: 2021-05-01

  Online published: 2021-05-03

Supported by

Nature Science Foundation of China(41906131)

Key Deployment Program of Center for Ocean Mega-Science, Chinese Academy of Sciences(COMS2020Q09)

Guangdong Basic and Applied Basic Research Foundation(2019A1515011340)

Abstract

Synechococcus, which is one of the representative groups of picocyanobacteria, is widely distributed in global oceans and estuaries. It has high phenotypic (pigment) and genetic diversity. According to the level of ability to deal with variation in salinity, Synechococcus can be divided into the strictly marine type and euryhaline type. In this study, we compared the growth of Synechococcus sp. K1 (isolated from the Pearl River estuary) and Synechococcus sp. YX02-1 (isolated from the South China Sea oligotrophic water) under a series of salinity gradient conditions. We also conducted phylogenetic analysis of K1 and YX02-1 based on rpoC1 gene sequences. We found that euryhaline Synechococcus sp. K1 could grow in all salinity levels, while strictly marine Synechococcus sp. YX02-1 was unable to survive below 13‰, which is consistent with the distribution characteristics in the estuary. Transcriptome analysis showed that the expression of genes (ggpS, SPS, stpA) that synthesize osmotic pressure regulator molecules in euryhaline Synechococcus were significantly down-regulated, while the gene expression of the membrane channel protein glzT was significantly up-regulated under low-salinity condition. This suggests that the low-salt tolerance mechanism of euryhaline Synechococcus mainly includes reducing the synthesis of small molecules related to intracellular osmotic pressure and increasing membrane channel proteins to improve the efflux of small molecules. In addition, salinity affected the photosynthesis and metabolism levels of euryhaline Synechococcus. The expression of photosynthesis energy-saving genes (ATPF0B, ATPF0A, ATPF1D, ATPF1A) and pigment protein genes (cpcA, cpcB) were significantly down-regulated under low-salinity condition, and the expression of inorganic nitrogen utilization-related genes were up-regulated. Under low-salinity condition, the demand for small molecules related to osmotic pressure is reduced, which can make more carbon sources to support growth. Meanwhile, the absorption of inorganic nitrogen is enhanced. These two reasons may be responsible for the higher growth of Synechococcus under low-salinity condition.

Cite this article

LIAO Ying , XIA Xiaomin . Effects of salinity on the growth and transcriptome of euryhaline Synechococcus sp. K1*[J]. Journal of Tropical Oceanography, 2022 , 41(2) : 159 -169 . DOI: 10.11978/2021040

海洋聚球藻(cluster 5 Synechococcus)是一类光合自养的微微型原核生物, 直径介于0.6~1.7μm之间, 广泛分布于河口和海洋区域(Moore et al, 2002; 马英 等, 2004)。聚球藻光合膜外排列着许多大分子聚集体——藻胆体(phycobilisomes, PBSs), 其功能是捕获并将已捕获的光能传递给类囊体上的叶绿体a进行光合作用(Yamanaka t al, 1981; 路荣昭 等, 1991)。根据藻胆体组成结构的不同, 聚球藻可分为3种色素类型(Callieri et al, 2002; Six et al, 2007)。Ⅰ类(type1)聚球藻只含藻蓝蛋白(phycocyanobilin, PCB), 吸收红光呈现绿色, 主要分布于河口区域; Ⅱ类(type2)聚球藻同时含有藻红蛋白(phycoerythrobilin, PEB)和PCB, 分布在相对干净的近岸水域; Ⅲ类(type3)聚球藻同时含有藻尿胆素(phycourobilin, PUB)、PEB和PCB, 对蓝光吸收较强, 广泛分布于外海(Xia et al, 2017, 2018)。按照PUB/PEB的比例, Ⅲ类聚球藻又可分为3a型(低PUB)、3b型(中PUB)、3c型(高PUB)、3d型(可变PUB)和3f型(超高PUB)(Everroad et al, 2012; Xia et al, 2017)。另一方面, 根据标记基因的序列差异, 如rpoC1 (Fuller et al, 2003; Mühling et al, 2005)和内源转录间隔区(internally transcribed spacer, ITS)(Rocap et al, 2002; Ahlgren et al, 2006), 海洋聚球藻被分为3个主要进化亚簇(subcluster): subcluster 5.1、subcluster 5.2和subcluster 5.3。其中, subcluster 5.1和subcluster 5.3主要为严格海洋型聚球藻, 广泛分布于外海, 而subcluster 5.2为广盐型聚球藻, 是河口的重要组成部分。近期研究中subcluster 5.1被进一步分为20个遗传型, 主要由Ⅱ类和Ⅲ类聚球藻组成(Xia et al, 2019)。
聚球藻具有丰富的物种多样性, 其分布受多种环境因素影响。前期研究中发现聚球藻的色素类型分布与盐度梯度紧密相关, PEB型聚球藻丰度在高盐区域较高, 而在低盐区域则显著降低(Wang et al, 2011)。此外, Xia等(2017)在香港水域的调查中发现, 随着盐度增加, 聚球藻的优势种逐渐由PCB型聚球藻(Ⅰ型)向PEB型聚球藻(Ⅱ和Ⅲ型)转变(Xia et al, 2017)。而遗传类群中, subcluster 5.2主要分布在低盐海域, 而subcluster 5.1-clade Ⅱ和subcluster 5.1 -clade Ⅲ广泛分布于高盐海域。随着生物技术的发展, 聚球藻基因组分析技术已比较成熟, Palenik等(2003)已破译聚球藻WH8102的全基因组。目前, NCBI数据库中已经有30多个相对完整的聚球藻基因组数据, 但聚球藻中部分基因功能还未知, 而转录组测序技术即RNA-seq (RNA Sequencing), 可反映物种的表达信息, 得到大量功能基因表达数据, 对揭示基因表达与生物学过程的关系具有重要意义。因此, 通过转录组的比较分析, 有助于揭示河口区域聚球藻的盐度适应机制。
本文选取广盐型聚球藻K1和严格海洋型聚球藻YX02-1, 对两者在不同盐度下的生长进行对比, 并进一步通过转录组分析广盐型聚球藻K1的耐低盐机制, 旨在加深人们对聚球藻盐度适应机理的认识, 为厘清聚球藻分布和变化趋势提供新的科学依据和预测方向。

1 材料与方法

1.1 聚球藻的分离纯化

将采集的1mL海水样品转至3mL无硅f/2培养基(表1)(Guillard et al, 1962), 用灭菌海水稀释5倍, 同时加入放线菌酮(20mg·L-1)以抑制真核藻类生长, 并置于25℃、光照:黑暗=12:12、35μmoL·m-2·s-1光照强度下培养数天。当培养基出现淡红色或者淡绿色时, 采取稀释法和平板划线法分离纯化聚球藻K1。稀释法: 取5个12mL细胞培养管, 高温高压灭菌后做标记, 每个培养管中加入9mL用灭菌人工海水配置的无硅f/2培养基, 并加入1mL藻液于1号培养管, 轻摇混匀后, 汲取1mL混合液至2号培养管, 按顺序依次稀释至5号培养管, 最后放置于光照培养箱(Kenton, China)中培养。平板划线法: 在500mL液体无硅培养基中加入5g琼脂, 高温高压灭菌后制成固体培养平板备用, 利用接菌环或涂布棒将藻液轻轻涂抹于平板上, 置于光照培养箱下培养数天并观察生长状况, 以ropC1标记基因利用PCR和流式细胞仪分析为依据验证聚球藻是否纯化成功。
表1 无硅f/2培养基配方

Tab. 1 The medium formula of f/2 without silicon

成分 质量浓度/(g·L-1)
常量元素 NaNO3 7.5×10 -2
NaH2PO4H2O 5.0×10 -3
微量金属 FeCl3·6H2O
Na2EDTA2·H2O
CuSO4·5H2O
Na2MoO4·2H2O
ZnSO4·7H2O
CoCl2·6H2O
MnCl4·4H2O
3.15×10 -3
4.36×10 -3
9.8×10 -3
6.3×10 -3
2.2×10 -2
1.0×10 -2
1.8×10 -1
维生素 Thiamine HCl (vit. B1)
Biotin (vit. H)
Cyanocobalamin (vit. B12)
2.0×10 -4
1.0×10 -3
1.0×10 -3
聚球藻YX02-1分离自南海, 由厦门大学郑强博士提供。

1.2 试验设计及生理测定

试验选用聚球藻YX02-1和K1, 通过调节人工海水中NaCl和MgCl2浓度设置5个盐度: 10‰、13‰、18‰、25‰和33‰, 并设3个重复组, 于40mL细胞培养瓶中连续培养, 培养条件: 25℃、光照:黑暗=12:12、35μmoL·m-2·s-1。培养过程中每24h取1mL藻液, 利用可见分光光度计(A360, China)进行OD440测量, 并绘制生长曲线。

1.3 全吸收光谱和流式细胞分析

采用可见分光光度计测量聚球藻YX02-1和K1的全吸收光谱(400~700nm), 并通过流式细胞仪(Becton-Dickinson FACSCalibur cytometer, Becton, Dickinson and Company, USA) 488nm和635nm的双激光束测定两种藻株的色素类型。

1.4 DNA提取、PCR及测序

采用DNA提取试剂盒(Invitrogen, Thermo Fisher Scientific, Carlsbad, CA, USA)提取聚球藻DNA。将滤膜用无菌剪刀剪成2mm×2mm碎片, 按照试剂盒说明书操作。DNA提取后储存于TE (Tris-EDTA)缓冲液中进行洗脱, 并于-20℃冰箱中保存, 待进一步作特定功能基因rpoC1扩增实验。第一轮PCR反应采用25μL反应体积: 模板1μL、引物(rpoC1-N5和rpoC1-C) 1μL、超纯水9.5μL和taq DNA聚合酶12.5μL。首先将反应体系处于95℃下进行10min高温变性, 然后进行30个循环: 95℃ 45s, 51℃ 45s, 72℃ 45s, 最后, 在Taq酶(72℃)作用下, 以dNTP为原料, 合成与模板互补的DNA链。第二轮采用50μL反应体积: 模板3μL、超纯水18μL、taq DNA聚合酶25μL和引物[rpoC1-39F (5′-GGNATTNGTNTGYGAGCGYTG-3′)和rpoC1- 462R (5′-CGYAGRCFCTGRTCAGCTT-3′)]各2μL (Xia et al, 2015, 2018)。扩增子通过Beckman AWPure Beads 纯化后, 采用Agilent 2100 Bioannlyzer 定量并于Illumina Hiseq 2500进行高通量测序。

1.5 进化树构建

从NCBI数据库下载除K1外的45株聚球藻的rpoC1基因序列。利用MEGA7.0软件进行同源性分析(Kumar et al, 2016), 通过 ClustalW对核苷酸序列进行一致性比对, 并以最大似然数(maximum likelihood)构建系统发育进化树, 以cluster 3 (聚球藻进化簇3)序列为外群, 测试方法Bootstrap, 重复数100。

1.6 样品收集和转录组分析

将培养在33‰和10‰并处于对数生长期的聚球藻K1, 通过0.2μm PC膜(polycarbonate membrane, 聚碳酸酯膜)过滤收集后放于1.5mL离心管中, 加入200μL RNA Later液氮速冻, 于-80℃超低温冰箱中储存等待送样测序。提取工作由广州美格公司承接, 提取后对样品质量进行检测并选取合格样品进行文库构建和初步定量, 并于Illumina Hiseq 2500平台进行文库测序, 获得带接头的原始测序序列, 使用megahit和Prodigal进行基因组装和预测处理原始测序数据(raw reads)后, 利用bowtie计算每个基因丰度, 最后得到高质量数据(clean reads)进行基因表达水平分析, 转录基因丰度升高表明基因表达水平也在提高, 通过序列数即可估计相应基因的表达水平。另外, 将获得的高质量数据进行基因注释并上传至京都基因与基因组数据库(Kyoto encyclopedia of genes and genomes, KEGG)和蛋白相邻类的聚簇数据库(cluster of orthologous groups of proteins, COG)进行序列比对和富集分析。

1.7 差异基因分析

通过R语言DESeq2进行基因差异性分析, 得到高盐和低盐样本中的差异表达基因(differential expression genes, DEGs), 并取p-value<0.05且|log2Foldchange|>1的基因作为差异基因集, 基因上升或下降2倍以上的进一步使用R (version 3.4.0)聚类软件包pheatmap作热图分析, 且所有TPM (Transcript per million, 每千个碱基的转录每百万映射读取的转录本)值进行平方根处理(square root transformed, sqrt), 基因名称、功能等通过KEGG注释。

1.8 数据分析

生长速度通过R (version 3.6.3)和Excel 2016进行统计学分析, 以平均值及标准差形式展示(mean±sd); 通过以下公式计算生长率(μ, 单位: d-1):
$μ[In(N_{t})-In(N_{0})]/(t-t_{0})$
式中: NtN0分别表示tt0时聚球藻在440 nm波长下的吸光值。本研究所获得的序列数据已上传NCBI数据库, 登录号为: PRJNA723138和MW971921, 其他聚球藻序列下载自NCBI数据库。

2 结果

2.1 聚球藻的生理特征

聚球藻K1和YX02-1在不同盐度下的生长曲线如图1所示。K1在10‰、13‰、18‰、25‰和33‰盐度下均生长良好, 13‰时具有最高生长率(对数期), 为0.189, 33‰时生长率最低, 为0.145 (图1b)。而YX02-1在盐度33‰时生长率最高, 为0.2, 10‰和13‰盐度中无法生存, 18‰盐度下的聚球藻生长状态较差, 生长率为0.1 (图1a、图2)。结果表明K1盐度适应性强, 但更偏好中低盐环境, 而YX02-1则在≤13‰盐度下无法生存。
图1 聚球藻YX02-1 (a)和K1 (b)在不同盐度(10‰、13‰、18‰、25‰、33‰)下的生长曲线

Fig. 1 The growth curve of Synechococcus sp. of YX02-1 (a) and K1 (b) at different salinity (10‰, 13‰, 18‰, 25‰, 33‰)

图2 聚球藻YX02-1在不同盐度(10‰、13‰、18‰、25‰、33‰)下培养至第8天时的生长状况

Fig. 2 The growth of Synechococcus sp. YX02-1 cultivated for eight days under different salinity (10‰, 13‰, 18‰, 25‰, 33‰)

2.2 聚球藻鉴定和进化树构建

聚球藻YX02-1 (图3a)和K1 (图3b)的全吸收光谱显示, YX02-1在490nm和550nm处存在两个吸收峰, 分别是PUB和PEB, 而K1在630nm处存在一个吸收峰, 为PCB (Six et al, 2007)。同时, 流式细胞分析显示YX02-1具有很强的藻红蛋白(PEB)荧光信号(图3c), 而K1具有较强的藻蓝蛋白(PCB)荧光信号(图3d)。该结果表明, YX02-1为PEB型聚球藻, K1为PCB型聚球藻。
图3 聚球藻YX02-1 (a)、K1 (b)的全吸收光谱图和YX02-1 (c)、K1 (d)的流式细胞荧光信号图

PEB: 藻红蛋白; PCB: 藻蓝蛋白

Fig. 3 The total absorption spectra of Synechococcus sp. YX02-1 (a) and K1 (b) and the flow cytometric fluorescence signals of Synechococcus sp. YX02-1 (c) and K1 (d)

基于极大似然估计(maximum likelihood)方法对rpoC1功能基因进行比对和分析, 构建包括YX02-1、K1在内的46株聚球藻的系统发育进化树, 进一步分析明确聚球藻的分类地位(图4)。结果显示, 聚球藻被分为3大类: subcluster 5.1、subcluster 5.2和subcluster 5.3, 其中, YX02-1与WH8103亲缘性为100%, 被共同划分为subcluster 5.1-clade Ⅲ, K1与LTW-RED、WH8007共同组成subcluster 5.2-B。结合图3聚球藻色素类型分析, 可鉴定YX02-1为PEB严格海洋型聚球藻, K1为PCB广盐型聚球藻。
图4 聚球藻系统发育进化树

主要分为3个分支: subcluster 5.1、subcluster 5.2和subcluster 5.3。聚球藻K1由红色圆圈标记, 聚球藻YX02-1由蓝色圆圈标记

Fig. 4 The phylogenetic tree of Synechococcus being mainly divided into three groups: subcluster 5.1, subcluster 5.2 and subcluster 5.3. K1 is marked by red dot, while YX02-1 is marked by blue dot

2.3 低盐条件下聚球藻K1在转录组水平的基因表达

聚球藻K1在低盐或高盐培养条件下的差异表达基因分析显示: 以高盐条件为对照组, 低盐条件下显著上调基因共227个, 下调基因70个, 并且差异性基因主要参与碳水化合物转运和代谢、能量产生和转化、信号转导等。在低盐条件下, 与渗透压有关的基因, 包括ggpS (葡萄糖基甘油磷酸合成酶)、SPS (蔗糖磷酸合成酶)、stpA (葡萄糖基甘油3-磷酸酶)等基因丰度发生显著下降(图5b、图6a)。并且聚球藻K1的色素蛋白cpcA和cpcB基因, 以及光合作用中ATPF0BATPF0AATPF1DATPF1A等相关产能基因丰度, 均在低盐条件下显著下降, 说明盐度会影响聚球藻的光合作用(图5b、图6b)。此外, 代谢通路包括碳循环、硫代谢、氮代谢等也均受到盐度影响。其中, cynT基因作为氮代谢的重要参与者, 在低盐培养条件下表达量增加。基因coxA、coxB和coxC是细胞色素亚基单位, 参与氧化磷酸化, 具有将氢离子从基质泵入细胞间隙的能力, 并产生ATP供转运者泵出渗透性成分或铁, 以维持低盐度条件下的渗透压平衡, 这些基因均在低盐条件下明显上调。而热休克蛋白HSP20基因表达量上升, 这类蛋白是非生物胁迫下细胞耐受机制的重要组成部分(图5a)。此外, FTRC基因明显下调, 该基因为铁氧还蛋白硫氧还蛋白还原酶, 其催化生成的还原型烟酰胺腺嘌呤二核苷酸磷酸(nicotinamide adenine dinucleotide phosphate, NADPH)可用于卡尔文循环中的固碳以及硫同化、碳同化等(图5b)。上述基因表达量的丰度变化是聚球藻K1适应低盐环境的重要内在机制。
图5 热图显示低盐条件下聚球藻K1显著上调(a)或显著下调(b)的基因及丰度

图中H1、H2、H3表示高盐组; L1、L2、L3表示低盐组。Sqrt(TPM)表示经过开平方转换的TPM值。图中左侧为基因名称@该基因在KEGG通路中的K号; 右侧字母为基因功能代码, J: 翻译、核糖体结构和生物发生; K: 转录; L: 复制, 重组和修复; D: 细胞周期, 细胞分裂, 染色体分配; V: 防御机制; T: 信号转导机制; M: 细胞壁/膜/包膜生物发生; N: 细胞运动; U: 细胞内运输, 分泌和囊泡转运; O: 翻译后修饰, 蛋白质周转, 分子伴侣; C: 能量产生和转化; G: 碳水化合物转运和代谢; E: 氨基酸转运和代谢; F: 核苷酸转运和代谢; H: 辅酶转运和代谢; I: 脂质转运和代谢; P: 无机离子转运和代谢; Q: 次级代谢产物生物合成, 转运和分解代谢; S: 未知功能

Fig. 5 The heatmap showing that the gene abundance (TPM) of Synechococcus sp. K1 was significantly up-regulated (a) or down-regulated (b) at low salinity

图6 聚球藻K1在高盐或低盐培养条件下渗透压有关的基因SPSggpS、stpA、通道蛋白基因(a)和色素蛋白基因cpcA、cpcB (b)的基因丰度(TPM)比较

**表示差异极显著(p<0.01)

Fig. 6 Comparison of the abundance (TPM) of genes SPS, ggpS, stpA, channel protein (a) and pigment proteins genes cpcA, cpcB (b) under high and low salinity

2.4 盐度对聚球藻K1主要代谢过程的影响

聚球藻K1在低盐和高盐培养条件下均保持着完整的代谢通路过程, 主要包括碳代谢、硫代谢和氮代谢。在碳代谢过程中, 参与葡萄糖降解的有关基因在低盐条件下显著上调。其中, glk基因(glucokinase, 葡萄糖激酶)上调了约1.4倍; G6PD (glucose-6-phosphate-1-dehydrogenase, 6-磷酸葡萄糖1-脱氢酶)上调了约2倍, 该基因是糖酵解和其他葡萄糖分解的重要参与者; PGLS (6-phospho- gluconolactonase, 6-磷酸葡萄糖酸内酯酶)和PGD (6-phosphogluconolactonase dehydrogenase, 6-磷酸葡萄糖酸脱氢酶)分别上调了1.7和1.9倍(图7a)。在碳代谢的其他过程中, 参与合成乙酰辅酶A的PDHA (pyruvate dehydrogenase, 丙酮酸脱氢酶)、PK (pyruvate kinase, 丙酮酸激酶)、PGK (phospho- glycerate kinase, 磷酸甘油酸激酶)等基因表达量基本降低了1~4倍, 其中, PGK下调最多, 约3.5倍(图7b)。硫代谢中, 参与异化硫酸盐还原过程中的基因表达量均下降(该过程为图7b中的硫酸盐还原为硫化物)。而在氮代谢过程中, 硝酸盐还原为铵盐的过程增强, 包括NRT (nitrate/nitrite transporter, 硝酸盐/亚硝酸盐转运子)、narB (nitrate reductase B, 硝酸盐同化还原酶基因B)、narA (nitrate reductase A, 硝酸盐同化还原酶基因A)在内的基因丰度均上调了1倍左右(图7c)。同时, CO2转化为氨基甲酸盐的能力加强, 基因cynT上调了约2.8倍, cynS上调了约1.9倍, 但谷氨酸盐代谢能力减弱。该结果说明盐度对硫代谢过程的消极影响最大, 其碳代谢部分通路在低盐下减弱, 并且还原硫酸盐的基因也受到限制, 但硝酸盐还原为铵的能力增强, 说明K1在河口下的生长过程需要大量的铵盐, 这也进一步解释了分布于河口的聚球藻对氮源的需求。
图7 聚球藻K1碳代谢(a)、硫代谢(b)、氮代谢(c)途径对盐度变化的响应

红色圆圈代表基因上调, 蓝色圆圈代表基因下调。圆圈大小表示调整幅度

Fig. 7 The responses of carbon metabolism (a), sulfur metabolism (b) and nitrogen metabolism (c) of Synechococcus sp. K1 to salinity changes. The red dot represents the genes of up-regulation, and the blue dot represents the genes of down-regulation

2.5 聚球藻K1的其他调控机制

通过对比低盐和高盐培养条件下聚球藻毒素-抗毒素系统(toxin-antitoxin system,TA-system)的相关基因(docmvpA、ftsZ、EARS)丰度, 发现两个试验条件下的基因表达量几乎没有变化(图8)。该结果说明TA-system作为一种重要的防御机制, 能够降低生长、抑制生长或者杀死细胞亚群来抵抗环境压力, 是广盐型聚球藻生活在河口的重要组成, 但并不是其适应盐度的调控机制。
图8 聚球藻K1在低盐或高盐培养条件下TA-system相关基因(docmvpA、ftsZ、EARS)丰度(TPM)的比较

Fig. 8 Comparison of genes (doc, mvpA, ftsZ, EARS) abundance (TPM) about TA-system under low and high salinity

3 讨论

前期针对河口区域聚球藻分布的研究显示, 聚球藻的色素类型组成及群落结构组成在河口区域具有明显的空间差异(Xia et al, 2015, 2017, 2018)。例如在珠江口, 随着盐度增加, 聚球藻的主要类群从淡水类群向subcluster 5.2, 再向subcluster 5.1-clade Ⅲ过渡。研究指出, 河口聚球藻群落组成的空间变化往往是由盐度、营养盐、温度、浊度等环境要素决定(Celepli et al, 2017; Xia et al, 2017)。本研究通过对低盐区域和高盐区域分离得到的聚球藻进行不同盐度条件下的生长对比, 发现K1在盐度10‰~33‰范围内均生长良好, 而YX02-1在盐度≤13‰下无法生长。这一结果进一步验证了盐度是影响聚球藻在河口分布格局的重要因素。此外, 我们发现虽然K1在所有盐度(10‰、13‰、18‰、25‰、33‰)下都能生长, 但是33‰盐度下聚球藻K1的生长率显著低于13‰盐度下的生长率。该结果表明K1更偏好低盐度, 解释了为何广盐型聚球藻的丰度往往在河口中低盐度区域较高(Xia et al, 2017)。进化分析显示, 遗传分类上K1属于subcluster 5.2, 而YX02-1属于subcluster 5.1 (图4), 与Herdman等(2001)的研究结果一致。Herdman等发现subcluster 5.1类群的聚球藻除了clade Ⅷ, 大部分都是严格海洋种, 而subcluster 5.2类群的聚球藻都是能耐受低盐的广盐类群(Herdman et al, 2001)。
通过转录组研究, 我们发现广盐型聚球藻耐受低盐主要是通过减少细胞内渗透压相关小分子合成, 同时增加膜通道蛋白提高小分子物质外排来达到细胞内外渗透压平衡, 这一结果与Xia等(2020)人的研究结论一致, 他们针对PEB型广盐型聚球藻的耐低盐机理开展研究, 发现该类型的聚球藻也是通过减少渗透压调节分子合成、减少小分子向胞内转运, 同时提高小分子外排来实现内外渗透平衡。研究表明, 蓝藻主要利用蔗糖、葡萄糖基甘油和甜菜碱作为渗透压调节分子(Celepli et al, 2017)。我们对生长于低盐和高盐下的聚球藻K1进行转录组比较分析, 发现在低盐条件下, 与糖酵解合成相关的基因ggpS和stpA基因表达量显著下降。ggpS和stpA作为糖酵解催化酶, 能催化ADP-葡萄糖转化为葡萄糖基甘油(GG), 表明葡萄糖基甘油是PCB型广盐型聚球藻适应渗透压变化的关键物质。这与严格海洋型聚球藻及PEB型广盐型聚球藻不同, 这两类聚球藻主要是利用甜菜碱作为渗透压调节分子(Xia et al, 2020)。在波罗的海的研究也发现, 聚球藻甜菜碱合成基因bsmB主要在高盐度站位发现, 而蔗糖合成和葡萄糖基甘油合成的基因SPSggpS在所有盐度都占有较高比例(Celepli et al, 2017)。这些结果表明使用不同的渗透压调节分子有可能是不同聚球藻类群耐受低盐能力差异的一个重要因素(Celepli et al, 2017)。另一方面, 转录组分析结果中还发现一种具有甘氨酸拉链(>fig|1129.219. peg.877, Glycine zipper, GXXXGX XXG)的通道蛋白(channel protein)基因, 其丰度在低盐条件下显著上升(图6a)。通过TMpred网站对该蛋白序列进行预测, 结果显示该蛋白具有跨膜双螺旋结构(图9)。根据最近的研究报道, 该蛋白是一类双向跨膜蛋白, 被命名为glzT, 是渗透压调节机制Env/OmpR系统的重要组成部分(Xia et al, 2020)。该类跨膜蛋白基因glzT, 在PEB型广盐型聚球藻(Xia et al, 2020)和PCB型广盐型聚球藻(本研究)都有发现, 并且在低盐条件下表达量显著上调, 表明该蛋白普遍存在于广盐型聚球藻中, 与低盐耐受密切相关, 是广盐型聚球藻在河口低盐环境中重要的生存机制。
图9 TMpred预测的基因glzT跨膜结构图

图中实线表示膜内向膜外预测; 虚线表示膜外向膜内预测。左和右侧小方框表示膜外区域

Fig. 9 The transmembrane structure diagram of glzT predicted by TMpred

最后, 本研究表明降低盐度可对聚球藻多个代谢通路产生影响。在低盐条件下, 广盐型聚球藻的光合作用可能受到盐度影响, 其编码色素蛋白的基因cpcA、cpcB以及相关产能基因ATPF0BATPF0AATPF1DATPF1A丰度均下调, 这与Deshnium等(1995)报道一致, 该研究将聚球藻从无NaCl培养基转移到含0.4mol·L-1 NaCl (23.4‰)的培养基后, 两天内光合作用减弱、细胞生长缓慢, 说明盐度变化会影响聚球藻的光合作用(Deshnium et al, 1995)。Ludwig等(2012)报道发现ATPF0/F1在高盐下基因丰度上升, 该结果同样说明盐度会影响聚球藻光合作用(Ludwig et al, 2012)。盐度降低对聚球藻的碳代谢、硫代谢、氮代谢也可发生不同程度影响。其中, 硫代谢所受影响最大, 相关基因丰度下调了1~2倍左右, 但氮代谢的相关基因普遍上调, 尤其是硝酸盐还原为铵的能力增强, 说明聚球藻K1在低盐条件下需要较多无机氮, 与低盐下较高的生长率相一致。
Marsan等(2017)研究发现, 在硝酸盐和磷酸盐限制下, 聚球藻CB0101的TA-system相关基因表达量在转录组水平发生明显上升, 这说明TA-system是营养盐限制条件下的重要保护机制, 目前已经在广盐型聚球藻包括CB0101、WH5701、CB0205中均发现了该系统(Marsan et al, 2017)。但通过对比低盐和高盐下聚球藻K1的TA-system相关基因丰度, 发现盐度变化对其丰度影响较小, 该结果说明TA-system是河口多变环境下的重要调节机制, 但可能并不是耐低盐机制的关键调节系统。

4 结论

本研究旨在阐明广盐型聚球藻K1适应低盐的内在分子机制, 通过转录组分析发现, 广盐型聚球藻渗透压调节机制更完善, 其在低盐条件下能通过降低渗透剂(葡萄糖基甘油)的生物合成、蛋白通道转运离子、Env/OmpR系统调节机制等, 改变细胞内渗透压以维持细胞正常生长。同时, 盐度还会影响聚球藻的光合系统及代谢水平。本文深入阐释了广盐型聚球藻的耐低盐机制, 但缺少与严格海洋型聚球藻的比较基因组分析, 后续研究需进一步结合多种研究方法, 从多方面更深入的探讨聚球藻对盐度的适应机制。

*感谢厦门大学郑强博士提供的实验藻株, 感谢编辑老师和审稿专家提供的宝贵意见和帮助。

[1]
路荣昭, 王淑芝, 关志英, 等, 1991. 聚球藻(Synechococcus leopoliensis 625)藻胆体-类囊体膜光谱特性和光能传递的研究[J]. 水生生物学报, 15(4): 368-371.

LU RONGZHAO, WANG SHUZHI, GUAN ZHIYING, et al, 1991. Studies on the spectroscopic properties and energy transfer of phycobilisome-thylakoid of Synechococcus leopoliensis 625[J]. Acta Hydrobiologica Sinica, 15(4): 368-371. (in Chinese with English abstract)

[2]
马英, 焦念志, 2004. 聚球藻(Synechococcus)分子生态学研究进展[J]. 自然科学进展, 14(9): 967-972.

MA YING, JIAO NIANZHI, 2004. Progress in molecular ecology of Synechococcus[J]. Progress in Natural Science, 14(9): 967-972. (in Chinese with English abstract)

[3]
AHLGREN N A, ROCAP G, 2006. Culture isolation and culture-independent clone libraries reveal new marine Synechococcus ecotypes with distinctive light and N physiologies[J]. Applied and Environmental Microbiology, 72(11): 7193-7204.

DOI

[4]
CALLIERI C, STOCKNER J G, 2002. Freshwater autotrophic picoplankton: a review[J]. Journal of Limnology, 61(1): 1-14.

[5]
CELEPLI N, SUNDH J, EKMAN M, et al, 2017. Meta-omic analyses of Baltic Sea cyanobacteria: diversity, community structure and salt acclimation[J]. Environmental Microbiology, 19(2): 673-686.

DOI

[6]
DESHNIUM P, LOS D A, HAYASHI H, et al, 1995. Transformation of Synechococcus with a gene for choline oxidase enhances tolerance to salt stress[J]. Plant Molecular Biology, 29(5): 897-907.

DOI

[7]
EVERROAD R C, WOOD A M, 2012. Phycoerythrin evolution and diversification of spectral phenotype in marine Synechococcus and related picocyanobacteria[J]. Molecular Phylogenetics and Evolution, 64(3): 381-392.

DOI

[8]
FULLER N J, MARIE D, PARTENSKY F, et al, 2003. Clade-specific 16S ribosomal DNA oligonucleotides reveal the predominance of a single marine Synechococcus clade throughout a stratified water column in the Red Sea[J]. Applied and Environmental Microbiology, 69(5): 2430-2443.

DOI

[9]
GUILLARD R R L, RYTHER J H, 1962. Studies of marine planktonic diatoms: Ⅰ. Cyclotella nana hustedt, and Detonula confervacea (cleve) gran[J]. Canadian Journal of Microbiology, 8(2): 229-239.

DOI

[10]
HERDMAN M, CASTENHOLZ RW, WATERNURY J B, et al, 2001. Form-genus ⅩⅢ. Synechococcus[M]//BOONE D R, CASTENHOLZ R W. Bergey’s Manual of Systematic Bacteriology. New York, NY: Springer-Verlag, 508-512.

[11]
KUMAR S, STECHER G, TAMURA K, 2016. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution, 33(7): 1870-1874.

DOI

[12]
LUDWIG M, BRYANT D A, 2012. Synechococcus sp. strain PCC 7002 transcriptome: acclimation to temperature, salinity, oxidative stress, and mixotrophic growth conditions[J]. Frontiers in Microbiology, 3: 354, doi: 10.3389/fmicb.2012. 00354.

DOI

[13]
MARSAN D, PLACE A, FUCICH D, et al, 2017. Toxin-Antitoxin systems in estuarine Synechococcus strain CB0101 and their transcriptomic responses to environmental stressors[J]. Frontiers in Microbiology, 8: 1213, doi: 10.3389/fmicb.2017.01213.

DOI

[14]
MOORE L R, POST A F, ROCAP G, et al, 2002. Utilization of different nitrogen sources by the marine cyanobacteria Prochlorococcus and Synechococcus[J]. Limnology and Oceanography, 47(4): 989-996.

DOI

[15]
MÜHLING M, FULLER N J, MILLARD A, et al, 2005. Genetic diversity of marine Synechococcus and co-occurring cyanophage communities: evidence for viral control of phytoplankton[J]. Environmental Microbiology, 7(4): 499-508.

DOI

[16]
PALENIK B, BRAHAMSHA B, LARIMER F W, et al, 2003. The genome of a motile marine Synechococcus[J]. Nature, 424(6952): 1037-1042.

DOI

[17]
ROCAP G, DISTEL D L, WATERBURY J B, et al, 2002. Resolution of Prochlorococcus and Synechococcus ecotypes by using 16S-23S ribosomal DNA internal transcribed spacer sequences[J]. Applied and Environmental Microbiology, 68(3): 1180-1191.

DOI

[18]
SIX C, THOMAS J C, GARCZAREK L, et al, 2007. Diversity and evolution of phycobilisomes in marine Synechococcus spp.: a comparative genomics study[J]. Genome Biology, 8(12): R259.

DOI

[19]
WANG KUI, WOMMACK K E, CHEN FENG, 2011. Abundance and distribution of Synechococcus spp. and cyanophages in the Chesapeake Bay[J]. Applied and Environmental Microbiology, 77(21): 7459-7468.

DOI

[20]
XIA XIAOMIN, VIDYARATHNA N K, PALENIK B, et al, 2015. Comparison of the seasonal variations of Synechococcus assemblage structures in estuarine waters and coastal waters of Hong Kong[J]. Applied and Environmental Microbiology, 81(21): 7644-7655.

DOI

[21]
XIA XIAOMIN, GUO WANG, TAN SHANGJIN, et al, 2017. Synechococcus assemblages across the salinity gradient in a salt wedge estuary[J]. Frontiers in Microbiology, 8: 1254, doi: 10.3389/fmicb.2017.01254.

DOI

[22]
XIA XIAOMIN, LIU HONGBIN, CHOI D, et al, 2018. Variation of Synechococcus pigment genetic diversity along two turbidity gradients in the China Seas[J]. Microbial Ecology, 75(1): 10-21.

DOI

[23]
XIA XIAOMIN, CHEUNG S, ENDO H, et al, 2019. Latitudinal and vertical variation of Synechococcus assemblage composition along 170° W transect from the South Pacific to the Arctic Ocean[J]. Microbial Ecology, 77(2): 333-342.

DOI

[24]
XIA XIAOMIN, LEE P, CHEUNG S, et al, 2020. Discovery of euryhaline phycoerythrobilin-containing Synechococcus and its mechanisms for adaptation to estuarine environments[J]. mSystems, 5(6): e00842-20.

[25]
YAMANAKA G, GLAZER A N, 1981. Dynamic aspects of phycobilisome structure: modulation of phycocyanin content of Synechococcus phycobilisomes[J]. Archives of Microbiology, 130(1): 23-30.

DOI

Outlines

/