科研 | ISME:转录组研究揭示了海洋硅藻温度适应的分子机制

编译:卡德加,编辑:小菌菌、江舜尧。

原创微文,欢迎转发转载。

导读

硅藻是海洋初级生产和碳循环的重要参与者,然而其对于温度变化的调节适应机制仍知之甚少。在本研究中,研究者分别对具有耐冷和耐热特性的两种海洋硅藻进行比较转录组分析,结果显示不同特性的海洋硅藻在基准基因的表达上存在一定的进化趋势。研究者认为这种转录层面的变化趋势,是海洋硅藻用于适应不同生长温度的关键机制。在不同的温度环境下,一些特定的如涉及蛋白处理机制、细胞膜结构、重要的短期生理机制等功能通路将被强化或削弱,这也意味着这些通路对于海洋硅藻的温度适应性起到关键作用。研究结果还显示涉及乙酰辅酶A,包括脂质和支链氨基酸代谢等关联通路,同样受到转录调节表现出进化趋势,这些通路也被海洋硅藻用于实现不同生长温度下光合作用中光捕获和代谢变化的平衡。

论文ID

原名:Molecularmechanisms of temperature acclimation and adaptation in marine diatoms
译名:海洋硅藻温度适应的分子机制
期刊:The ISME Journal
IF: 9.493
发表时间:2019年5月24日(在线发布)
通讯作者:Yue Liang
通讯作者单位:中加拿大达尔豪斯大学海洋学系

实验设计

1. 实验对象及环境

实验对象为两种角毛藻属(chaetoceros sp.),分别为CCMP158和CCMP160。CCMP158分离于加拿大圣玛格丽特湾,自1982年起便在14℃的培养环境下生长;而CCMP160分离于法国马赛湾,同样自1982年后生长于24℃的培养环境。在实验开始前,CCMP158和CCMP160分别在16℃和25℃的环境下生活1-2月。随后以6℃-30℃设置6个梯度区间,每组样设置3个重复,观察CCMP158和CCMP160九个世代后的生长情况(CCMP158的观察周期为9-56个世代,CCMP160的观察周期为16-42个世代),生长速率主要通过叶绿素a荧光实验的结果来判断,随后这两种硅藻分别进行细胞组分和转录组分析。当达到最高生长速率时,定义当时的温度为最适生长温度Tmax生长速率为最高值的50%和86%则分别定义为次适较高温度Tsup和次适较低温度Tsub(图1)。

图1. 角毛藻属CCMP158和CCMP160在不同温度处理下的稳定生长速率

2. 生化指标测定和代谢组学研究

随后每个样取50mL进行生化指标的测定,通过比色法和荧光测定仪等仪器鉴定了叶绿素a、蛋白和RNA的含量。测序平台为Illumina HiSeq 2000,文库片段读长为100 bp的双末端序列,数据过滤采用Trimmomatic。原始数据首先通过Bowtie2比对到PhiX基因组以去除污染序列,随后那些未比对上的reads在进行过滤后,用于后续Trinity软件的de novo组装。随后再次利用Bowtie2,将每个样过滤后的reads比对到Trinity从头拼接组装好的转录组。转录本丰度通过eXpress 1.5.1版本进行估算,随后将那些mapping reads数大于0的基因(3个重复中至少有2个满足此条件)收集,利用DESeq2获得差异表达基因(FDR调整后的p-value < 0.01)。

得到差异基因后,根据它们的表达特点分成9种可能的表达模式(表1)。定义那些仅在冷适应硅藻中且处于次适较高温度,或仅在热适应硅藻中且处于次适较低温度状态下差异表达的基因为潜在强化基因(invested),同理那些具有相反表达模式的差异基因为潜在弱化基因(divested)。通过使用BLAST+、GhostKOALA和TransDecoder等工具,以及现有的KEGG和Swiss-prot数据库,对这些差异基因进行KEGG通路和Brite聚类。这些聚类结果主要涉及代谢、遗传信息处理、环境信息处理和细胞过程等生物学功能。任何一种硅藻中注释结果数量少于10个KO ID的通路将被剔除,并通过超几何分布确定显著富集的KEGG通路(p < 0.05)。

结果

1. 不同生长温度对硅藻生理、生化组成和细胞尺寸的影响

具有冷环境适应的硅藻种(CCMP158,Tmax =16 °C)相较热环境适应硅藻(CCMP160,Tmax = 25°C),具有更慢的生长速率和更低的最适生长温度(Tmax)(图1)。这两类硅藻种在最适生长温度下,表现出光合生理和生化组成的差异,并且这些差异在次适生长温度的条件下(略高于或低于最适温度)同样存在。如NPQ的含量在冷适应硅藻中的含量(0.22-0.3)要显著高于热适应硅藻(0.1-0.16),尽管这种物质在两种硅藻中并未显著表现出含量随生长温度变化的生长趋势。此外,叶绿素a物质的含量在CCMP158中表现出随温度递增的趋势,然而其在CCMP160中并未有此趋势。在最适温度(Tmax)和次适较低温度(Tsub)下,热适应硅藻叶绿素a的相对含量均要高于冷适应硅藻,而在次适较高温度(Tsup)下,该物质含量则差异不大。CCMP160的RNA含量在温度从次适较低温度转向次适较高温度时,发生了显著增长;除此之外,两种硅藻的RNA和蛋白质含量均未有显著变化。在最适温度和次适较低温度下,CCMP160的RNA和蛋白质相对含量均比CCMP158要高。在各自的最适温度下,冷适应硅藻的细胞体积约为热适应硅藻的2.5倍。随着温度的不断降低,CCMP158的细胞体积表现出递增的趋势,然而CCMP160却与CCMP158表现出相反的趋势,尽管这种趋势不具有显著性(p > 0.05)。

2. 相关基因鉴定和注释
通过de novo组装,共在CCMP158和CCMP160中得到42,048和27,580个基因。其中CCMP158中有24.1%的基因成功被KEGG数据库注释,这个比例在CCMP160中为25.8%。两种硅藻共有相同的3,579个KO ID和4,913组OrthoMCL分析结果,其中2,032个KO ID和3,583组OrthoMCL结果在模式海藻物种假微型海链藻(Thalassiosirapseudonana)和三角褐指藻(Phaeodactylumtricornutum)中也有发现。利用GhostKOALA的分类结果显示,这两个de novo组装转录组的结果并未有其它不相关序列的掺入,如来源于细菌、古细菌和真菌的序列。
3. 不同温度影响下的差异基因表达
对于冷适应硅藻CCMP158,在次适较高温度(Tsup)和次适较低温度(Tsub)下,分别共有7%和3%的基因发生差异表达,而这个比例在热适应硅藻CCMP160中,则分别为3%和7%(附表1)。每个硅藻种处于不同温度的所有差异表达基因被分类为9种不同的模式(表1)。在冷适应硅藻CCMP158中,67.4%的差异表达基因仅在次适较高温度表达,而在热适应硅藻CCMP160中,则约有67.8%的差异基因仅在次适较低温度表达。研究者据此推测这些差异表达基因存在被强化或削弱的可能性。
表1. 次适较高温度Tsup和次适较低温度Tsub下表达差异的九种可能模式的描述和解释

4. 与温度适应相关代谢通路的鉴定

基于KO的注释结果,在冷适应硅藻从最适温度至次适较低温度的变化中,糖酵解/葡糖异生和柠檬酸循环代谢通路均发生上调表达,然而伴侣和折叠催化分子,以及内质网中涉及蛋白处理等通路则发生下调表达(基于KO counts和p-value的结果,表3)。涉及细胞周期、DNA修复和重组蛋白的通路发生了上调表达,而光合作用、光合作用蛋白、卟啉和叶绿素代谢、脂肪酸生化合成、脂质生化蛋白、蛋白酶体、核糖体、RNA运输和翻译因子等代谢通路则发生了下调表达.
而在热适应硅藻从最适温度向次适较低温度的变化中,剪接体和核糖体生物合成通路发生上调表达,而涉及脂肪酸降解、缬氨酸、亮氨酸和异亮氨酸降解、伴侣和折叠催化分子、蛋白酶体和肽酶等通路则为下调表达(表3)。当温度变化从最适温度转为次适较高温度后,涉及泛素系统和内质网蛋白处理的相关通路发生上调表达,而赖氨酸生物合成下调表达。在冷适应硅藻中,除了伴侣和折叠催化分子,其余在Tsup温度中下调表达的通路均被鉴定为潜在强化通路。除此之外,涉及不饱和脂肪酸和蛋白酶体生物合成的通路[PATH:ko03050],同样被认为是冷适应硅藻的潜在强化通路,而细胞周期通路则认为是其潜在削弱通路(表3)。研究者认为对于热适应硅藻,所有次适较低温度中下调表达的通路均有被强化的趋势,而剪接体和核糖体生物合成通路则被弱化。
表3. 冷适应硅藻(CCMP158)以及热适应硅藻(CCMP160)显著富集的上调(向上的箭头)和下调(向下的箭头)途径
 

讨论

尽管地球上各式各样的生命整体能适应很大区间的温度变化,绝大多数物种自身的适存温度区间却是很狭窄的。一般而言,最适生长温度较高的物种往往有更强的生长速率,过高或过低的生长温度会导致能量代谢、氧化还原反应不平衡等一系列现象,最终引发自身的补偿应答和压力相应机制。研究者通过以两个亲缘关系很接近且最适生长温度存在差异的角毛藻属为研究对象,开展比较转录组分析,以进一步确定与温度适应相关的具体分子机制,以及硅藻生长速率为何随着温度下降而降低的分子基础。

1. 次适较高或较低温度的转录组数据存在强化或削弱机制
通过转录组注释,两个角毛藻属绝大多数的预测蛋白存在同缘关系,且与模式硅藻种假微型海链藻和三角褐指藻一致(附图1)。在冷适应硅藻CCMP158中,绝大多数差异表达基因出现在次适较高和最适温度的比较中,而对于热适应硅藻CCMP160,大多数差异表达基因则在次适较低和最适温度的比较中出现(表1和附表1)。研究者推测这些转录级别的差异表现,反映出在长时间处于不同模式下,基因开始表现出具有与进化相关的表达趋势。研究者把这些次适高/低温度下的基因表达模式归为9类,并分别赋予它们生物学解释(表1)。其中两种表达模式的变化仅受到温度的影响,而与整体的代谢速率无关,暗示温度是这两种表达模式变化的决定因子,并进而导致酶活性的改变。另外两种模式与生长速率分别表现出正相关和负相关的趋势,并与细胞周期及胁迫相应相关的基因存在关联。另有四类基因表达模式表示仅在一种温度变化模式下出现差异表达(如次适较高或较低温度)。这类在一种温度模式稳定表达,而在另一种温度模式下产生差异的基因,被研究者定义为具有适应性响应和强化/削弱的潜力。举例来说,当一种冷环境适应的物种处于其最适温度时,其基因表达将显著增加,而继续降低温度或升高温度都会让原先最适环境下的基因表达量降低。研究者通过卡方检测,证实每种硅藻不同温度下的差异表达基因均相互独立,如果不同的温度适应性改变了基因表达,则需求不同生长温度的不同物种也将具有不同的差异基因表达。在8组差异基因表达模式中,其中的6组在两种硅藻间存在显著差异,而差异最大的部分就在强化/削弱分类中(即仅在一种温度模式中存在差异基因)(表1),这也表明这些表达模式与温度适应存在强关联性。因此,研究者认为海洋硅藻之所以能适应温度变化,部分原因是通过基因表达的组成型差异,借此进一步增加或减少额外的转录组响应。这个猜想也与之前发现在冷环境适应硅藻中,在次适较高温度而非次适较低温度中发现了更多的差异基因的现象相一致。
2. 被增强或削弱的通路信息中包含已知与温度适应相关的生理机制

之前有研究显示,当处于低温环境时,蛋白质翻译尤其是核糖体的起始翻译速率均会被抑制,而本研究也表明生长温度能显著影响两种硅藻间涉及翻译和蛋白质合成相关基因的差异表达。在冷适应硅藻中,涉及氨酰tRNA合成、核糖体蛋白、RNA转运蛋白和翻译因子通路的基因均发生上调表达,作者据此认为这可能是为了补偿在更冷环境下翻译速率降低的现象,因此硅藻中的一些基准基因表达可能也得到了进一步增强。此外,涉及tRNA合成的酶基因表达同样得到了强化(表3),这也与氨酰tRNA上调表达的结果相一致。而类似的比较则在热适应硅藻中得到了截然不同的结果,如核糖体生物合成和剪切体相关基因的表达量显著降低。热适应硅藻中核糖体合成功能的弱化,不仅能减少材料和能量消耗,也能缓解多肽处理和内质网的压力,伴侣蛋白和折叠催化分子的表达量提高同样能起到降低内质网压力的作用(表3)。这些与蛋白代谢相关的基准基因的表达差异,对于今后的蛋白组学研究能提供很好的理论支持。

研究者同样针对氮磷比与温度的相关性提出了一些假设。之前的研究显示低温能有效促进浮游生物中磷富集rRNA与氮富集蛋白的比例,进而提高整个海平面的氮磷比。然而,研究者的研究结果显示冷适应硅藻似乎比热适应硅藻有更低的RNA/蛋白比率(表2),作者认为这与热适应硅藻有更高的生长速率有关,从而导致更高的RNA相对含量,而RNA/蛋白的比率似乎并不能直接导致海平面的温度变化。

低温环境会降低蛋白质转换的速率并减少膜的流动性,这一结论通过光合作用蛋白和不饱和脂肪酸生物合成的代偿性上调表达获得证实(表3)。具有寒冷适应性的浮游生物,相较那些具有热适性的物种,似乎有自身独到的策略以应对更慢的蛋白转化速率,如生成更多的光合作用蛋白和更稳定的光合色素等。在本文中,研究者便在冷适应硅藻CCMP158中发现了光合作用PSII类蛋白的增强表达(附表4和附图3),以及更多的叶绿素a和更高的NPQ表达水平,这也表明CCMP158能在处于低温时开启自身的补偿机制,以应对获取能量的减少。此外,不饱和脂肪酸生物合成基因的上调表达能调整脂肪酸在脂质双层中的饱和状态,从而起到动态调节膜流动性的目的。

3. 乙酰辅酶A代谢调控对于硅藻温度适应至关重要

在光合生物中,温度变化能改变代谢平衡,从而影响那些核心代谢物的比例,如NADPH:NADP+和乙酰辅酶A:辅酶A等。研究者在两种硅藻中均发现了脂肪酸代谢的上调表达,以及CCMP160中支链氨基酸(BCAA、缬氨酸、亮氨酸和异亮氨酸)的表达量提高,这表明乙酰辅酶A代谢在硅藻温度适应中起到关键作用。脂质能有效缓解不平衡NADPH:NADP+比例所带来的氧化还原压力,而作为脂肪酸代谢和支链氨基酸代谢通路中的核心物质,乙酰辅酶A的不足会降低细胞分裂的速度,而过高程度的细胞溶质乙酰辅酶A则会导致溶质蛋白的过乙酰化,从而抑制酶活性。作者据此认为脂肪酸和支链氨基酸代谢通路在硅藻环境适应中扮演关键角色。处于低温环境时,脂肪酸的生物合成速率提高,这有助于消耗随着温度降低而增多的乙酰辅酶A和NADPH。而当处于高温环境时,这样的情况又有所不同,光合作用下的碳固定可能不能满足随着代谢活动增强,从而导致的对乙酰辅酶A的进一步需求。从另一个角度来看,增强的组蛋白乙酰化过程将导致细胞内辅酶A的数量不断增加,高浓度时可能导致细胞程序性死亡,而乙酰辅酶A的合成能平衡其浓度,使那些热适应物种保持更高的生长速率。

4. 低温下最大生长速率降低的分子机理

整体来看,具有热适应特性的物种往往具有更高的生长速率,然而生长速率和最适温度之间的分子机理长期以来仍处于未知状态。在冷适应硅藻中,研究者发现细胞周期通路[PATH:ko04141]相关的基因发生上调表达,且绝大多数相关基因与有丝分裂和S期相关(附表3和附图2)。在结果中找到了DNA复制许可因子MCM3(K02541),这个因子是微染色体平衡复合体(MCM complex)的一部分,是在S期DNA复制和起始调控的必需因子。在CCMP158中,这些涉及细胞周期相关的基准基因的低表达,可能在有丝分裂之前延长G2期的持续时间,从而导致CCMP158相对CCMP160更大的细胞容积,这个现象在酵母中也有发现。

为了进一步验证结果,研究者通过现有的海洋真核微生物转录组数据集(MMETSP),模拟分析了这两种海洋硅藻在其最适温度下的差异表达基因。分析结果显示CCMP158中蛋白质合成和光合作用相关基因的表达水平均要比CCMP160更高,这个结果也与之前的分析结果一致(附表5),并且尽管蛋白酶体通路在两种硅藻中均得到了强化表达,CCMP158的表达水平要比CCMP160更高。尽管脂肪酸合成和支链氨基酸的差异表达量在CCMP160中并不显著,但两种硅藻中这些通路的KO ID数量整体是一致的(附表5)。这些结果进一步验证了研究者的假设,即对温度刺激的短期反应具有物种特异性,而温度的变化能有效改变关键基准基因和通路的表达水平,以维持代谢平衡。这些关键通路主要包括氨基酸、碳水化合物、能量和脂类代谢,蛋白质折叠和降解、细胞生长和死亡,以及转录和翻译等(表3)。

评论

这篇文章以两种具有不同温度适应性的角毛藻属为研究对象,以温度为梯度,分别选取最适生长温度、次适较高温度和次适较低温度这三个节点,开展比较转录组研究。除了研究对象有一定的特殊性外,实验设计相对并不复杂,分析方法也比较好理解。但有别于一些常规转录组分析流程,本文的讨论部分大大加深了文章的深度,并就乙酰辅酶A等基础物质随温度的变化提出自己的假设,显示出研究者极其出色的分子生物学功底,也使其从一众转录组文章中脱颖而出。此外,精准的采样和实验设计也使得后续的差异基因鉴定取得了很好的效果。这篇文章的不足之处在于还缺少下游实验的验证,尤其是多组学数据的交叉验证,仅从单一转录组和差异基因得到的一些结论,说服力似乎稍显不足。且在文中多次提到基因表达进化的观点,但似乎并没有更细一步关于进化的阐述

(0)

相关推荐