期刊：Genomics
影响因子：6.2
文章题目：DenovoassemblytranscriptomeanalysisrevealsthepreliminarymolecularmechanismofpigmentationinjuvenilesofthehardclamMercenariamercenaria
发表期刊：Genomics
技术手段：RNA-Seq+qPCR
派森诺生物与中国科学院海洋研究所携手合作，于近期在《Genomics》上发表了无参转录组分析揭示幼体硬壳蛤（Mercenariamercenaria）色素沉着的初步分子机制相关研究成果。
/研究背景/
颜色是个体重要的表型特征，并且与伪装、免疫、体温调节等生物学功能相关。软体动物具有丰富的外壳颜色，但是，人们对贝壳色素沉着的了解仍然很有限。已有学者利用杂交进行的遗传分析表明，贝壳的颜色可以由基因调控。研究表明，钙信号转导、黑色素合成、卟啉和叶绿素代谢等相关基因参与生物外壳的色素沉着与调控。这些研究大多数以双壳类动物的外被膜作为研究材料，但是，双壳类动物幼体时期是开始色素沉着的重要时期却往往被忽视。
硬壳蛤（Mercenariamercenaria）是一种有着重要商业价值的双壳类动物，其养殖品种具有丰富的壳色变异。本研究根据色素沉着情况将幼体硬壳蛤分为三组，采用RNA-Seq技术进行分析，从中鉴定出参与色素沉着的候选基因和相关代谢通路。
/研究方法/
样本：无色素沉着幼体硬壳蛤（NP）、白色幼体硬壳蛤（W）、红色幼体硬壳蛤（OP）
实验方法：转录组测序
测序平台：IlluminaNovaSeq
技术路线图见下：
/研究结果/
1.转录组测序结果，注释以及差异表达基因（DEGs）分析
借助二代转录组测序技术，对NP、W、OP三组幼体硬壳蛤，每组三个重复，一共9个样品进行转录组测序。通过转录本拼接共得到个unigenes，N50为bp，对unigene功能注释鉴定出很多与信号转导相关的基因。
为了评估在不同分组的基因表达量变化，对三组幼体硬壳蛤的转录表达情况进行两两比较。在OPvsNP比较组中检测到个上调基因，个下调基因。WvsNP比较组中有个上调基因和个下调基因。OPvsW比较组中发现个上调基因和71个下调基因。
为了更好地了解DEGs的功能，作者进行了KEGG富集分析。结果显示，OPvsNP比较组富集了条途径，其中有51条途径显示显著富集（表1），富集最显著的前20个代谢通路如图1A所示。WvsNP和OPvsW比较组分别富集了和条途径，且分别有45和15条途径显示显著富集（表2，表3），富集最显著的前20个代谢通路如图1B和图1C所示。DEGs主要在黑色素瘤、卟啉和叶绿素代谢以及ABC转运通路中富集。
表1OPvsNP显著富集途径条目
图1AOPvsNPKEGG富集分析
表2WvsNP显著富集途径条目
图1BWvsNPKEGG富集分析
表3OPvsW显著富集途径条目
图1COPvsWKEGG富集分析
2.与色素沉着相关的候选基因鉴定
作者研究发现，红色硬壳蛤的黑色素和胡萝卜素浓度明显高于白色硬壳蛤。色素前体分子可以通过ABC转运蛋白运输到色素细胞，并且四吡咯是软体动物中另一种主要的色素。因此，作者挑选了黑色素瘤/黑色素生成、类胡萝卜素相关代谢通路、ABC转运体以及卟啉和叶绿素代谢的DEGs进行聚类热图分析。聚类结果如下：
A.在黑色素瘤/黑色素生成代谢通路途径中，W组首先与OP组聚集，然后与NP组聚集。在该通路中发现了小眼畸形相关转录因子（MITF）、酪氨酸酶（TYR）和钙离子结合相关基因。TYR是合成黑色素所必需的，MITF对TYR具有调控作用。TYR在OP和W组中高表达，MITF在W组中表达相对下调（图2a）。
B.在ABC转运体代谢通路中，W组先与OP组聚类，然后与NP组聚类。ABC转运体可能通过转运色素前体来介导色素沉着。相比于W和NP组，ABCA、ABCC和ABCG在OP组中高表达（图2b）。
C.类胡萝卜的摄取、运输和分解在类胡萝卜色素沉着中起着至关重要的作用。在类胡萝卜素相关代谢通路中清道夫受体、低密度脂蛋白受体相关蛋白和载脂蛋白等的Unigenes差异表达（图2c）。
D.与W组相比，OP组中卟啉和叶绿素代谢显著富集，表明该途径在硬壳蛤中红色素的形成中发挥了作用（图2d）。
图2显示