大珠母贝（Pinctada maxima）肠道微生物组–宿主互作的转录组景观揭示与驯化相关的改变

《Animals》：Transcriptomic Landscape of Gut Microbiota–Host Interactions Reveals Domestication-Related Changes in the Pearl Oyster Pinctada maxima

　　

大珠母贝（Pinctada maxima）的驯化后代在应对近岸复杂多变的环境条件时，表现出多种生理和微生物层面的适应性调整。为支持大珠母贝种质资源恢复，并探究这些表型反应背后的分子机制，研究人员对野生亲本代与驯化后代的肠道转录组和微生物组进行了比较分析。每个样本平均产生43,600,525条高质量读段（clean reads），比对至大珠母贝参考基因组，比对率介于63.91%至79.48%之间，共鉴定出3007个差异表达基因（differentially expressed genes, DEGs）。基因本体论（Gene Ontology, GO）分析显示，差异表达基因富集于有机酸代谢；京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）分析显示，差异表达基因富集于糖鞘脂生物合成及细胞色素P450介导的外源物质代谢途径。微生物谱分析揭示了门水平和属水平的显著组成变化，F1代α多样性增加；优势菌门向螺旋体门（Spirochaetota）和拟杆菌门（Bacteroidota）转变，功能预测提示代谢和免疫能力增强。定量实时PCR（qRT-PCR）验证了免疫相关基因（PmHR96h、PmIAP）的上调以及钙信号相关基因（PmCaM、PmHSP90）的下调，与RNA-seq数据一致。总体而言，这些协调的转录组和微生物组改变反映了宿主–微生物组对近岸环境条件的多方面适应性反应。研究结果为选择性育种和健康监测提供了有价值的分子标记和微生物指标，为在不断变化的环境条件下提高大珠母贝养殖的抗逆性提供了科学依据。

一、研究背景与科学问题

大珠母贝（Pinctada maxima）是分布于中印度-太平洋海域的重要经济贝类，以其大型、光泽优异的海洋珍珠而闻名。然而，过度捕捞、气候变化和栖息地退化已导致该物种在浅于15米水域的种群近乎消失，尤其在南海北部已十分罕见，严重制约了以野生群体为基础的珍珠养殖业的可持续发展。尽管科研人员围绕种质资源分析、人工育苗和育种开展了大量工作，但大珠母贝养殖过程中尤其是壳高3–6厘米阶段的幼苗高死亡率问题仍然突出，已成为全球性关注焦点。已有研究表明，处于死亡高峰期的幼贝经历氧化应激、免疫基因激活和低能量状态恶化等多重生理压力。驯化被认为是实现遗传改良的重要途径，前期研究发现经过连续两代驯化后，养殖幼苗存活率可提高10%以上，并鉴定出多个与钙信号、免疫防御和代谢相关的选择信号基因。然而，关于驯化过程中肠道微生物组与宿主转录组的整合研究仍然有限。肠道在消化营养吸收、病原抵抗和免疫调节中发挥关键作用，多组学整合分析可为理解驯化适应机制提供重要视角。基于此，研究人员提出假设：驯化群体中的高死亡率与快速环境变化引发的强烈自然选择有关，研究目标是鉴定驯化群体的环境适应性标记，为未来驯化策略提供科学依据。

二、研究内容与主要结论

研究人员以2014年采集自海南三亚、养殖于广东湛江北部湾的大珠母贝群体为材料，于2018年随机选取20只健康成熟亲本（F0）通过群体繁殖产生第一代驯化后代（F1），并在相同海区、相同生长阶段（壳长13.43 ± 1.12 cm）对两个世代进行标准化采样。通过整合肠道转录组测序和16S rRNA基因扩增子测序，系统比较了野生亲本与驯化后代的基因表达和肠道微生物组成差异。研究发现，驯化后代肠道中免疫相关基因显著上调，代谢和神经信号通路相关基因表达发生改变，肠道微生物群落结构发生重组，α多样性增加，优势菌门向螺旋体门和拟杆菌门转变，功能预测提示代谢和免疫能力增强。这些协调的宿主–微生物组变化反映了大珠母贝对近岸养殖条件的多方面适应性反应。研究鉴定出的候选基因和微生物标志物为未来选择性育种和健康监测提供了科学依据。该论文发表在《Animals》期刊。

三、主要关键技术方法

研究人员采用了以下关键技术方法：其一，利用Illumina NovaSeq X Plus平台对肠道组织进行转录组测序（RNA-Seq），经质量控制后使用HISAT2软件比对至大珠母贝参考基因组，通过DESeq2进行差异表达基因鉴定；其二，针对肠道内容物提取总DNA，扩增16S rRNA基因V3–V4高变区，在Illumina PE300/PE250平台测序，使用QIIME2中的DADA2插件进行降噪生成扩增子序列变异（ASVs），基于Silva 138数据库进行物种注释，并利用PICRUSt2进行微生物功能预测；其三，采用定量实时PCR（qRT-PCR）对11个差异表达基因进行验证，以β-actin为内参基因，使用2^?ΔΔCt法计算相对表达量。样本来源为广东湛江北部湾养殖海区，F0为野生亲本，F1为驯化后代，每组3个生物学重复。

四、研究结果

1. 1.
   F0与F1转录组分析

六个转录组样本共产生264,094,684条原始读段，质控后保留261,603,152条高质量读段，平均每样本43,600,525条，平均GC含量39.37%，Q20为98.71%，Q30为95.86%。比对率介于63.91%至79.48%之间。主成分分析显示F0与F1组间明显分离，组内生物学重复相关性高。以F0为对照共鉴定出3007个差异表达基因，其中1664个上调、1343个下调。GO富集分析显示差异表达基因显著富集于有机酸代谢过程、羧酸代谢过程和含氧酸代谢过程等生物学过程条目。KEGG富集分析显示差异表达基因显著富集于糖鞘脂生物合成（乳糖和新乳糖系列）以及细胞色素P450介导的外源物质代谢途径，此外还富集于NF-κB信号通路、炎症性肠病和沙门氏菌感染等免疫和信号转导相关通路。

1. 2.
   驯化后代肠道微生物群落组成和多样性差异

16S rRNA高通量测序共获得601,880条高质量读段，鉴定出2293个ASVs。F1个体的ACE、Chao1、Shannon和Sobs多样性指数高于F0，但差异未达统计学显著水平。属水平韦恩图显示两组共有158个ASVs，F0特有124个，F1特有255个。ASV水平主坐标分析显示F0与F1样本明显分离，PC1和PC2分别解释50.71%和25.7%的方差。系统发育树分析进一步支持两组样本分别聚为两个独立分支。在门水平上，F0肠道微生物以变形菌门、厚壁菌门和蓝藻门为主，F1以变形菌门、螺旋体门、拟杆菌门和蓝藻门为主，变形菌门在两组中均为最优势门。在属水平上，F0以支原体属和假单胞菌属为主，F1以螺旋体属、支原体属和假单胞菌属为主。

1. 3.
   肠道微生物功能预测

PICRUSt2功能预测显示F1与F0在代谢、免疫应答、凋亡和细胞信号传导等功能基因丰度上存在差异，与转录组分析结果一致，提示可能存在基因–微生物功能协同机制。F1肠道微生物富集的功能与翻译、核糖体结构和生物发生相关，代谢和环境信息处理相关通路也显著富集。

1. 4.
   RNA-Seq数据的qRT-PCR验证

对11个差异表达基因进行qRT-PCR验证，这些基因分为四类：免疫应答（PmHR96h、PmIAP）、营养代谢（PmSLC46A3、PmPEX7、PmATPase-1、PmATPase-2、PmGGCX）、神经信号（PmPRRT1、PmCAP）和钙离子调节（PmCaM、PmHSP90）。结果显示PmHR96h和PmIAP在F1中显著上调，PmSLC46A3显著上调，PmCAP、PmCaM和PmHSP90显著下调，PmPRRT1显著上调，与RNA-seq结果一致，证实了转录组数据的可靠性。

五、讨论与结论

讨论部分指出，驯化后代的转录组改变涉及免疫、代谢和神经信号通路，反映了对近岸养殖条件的生理适应。神经信号方面，HSP90、钙调蛋白（PmCaM）等基因表达改变提示环境变化可能影响细胞内钙离子稳态和下游钙依赖性代谢。维生素吸收方面，SLC46A3介导支链氨基酸转运，可能影响维生素和代谢物的生物利用度；PEX7参与脂质代谢；GGCX编码维生素K依赖性蛋白的γ-羧化酶。免疫适应方面，IAP家族蛋白通过BIR结构域抑制caspase活性，HR96h编码HSP90家族成员gp96，这些免疫相关基因的上调表明免疫通路被激活以应对近岸环境胁迫。肠道微生物方面，驯化改变了微生物群落结构，拟杆菌门丰度增加可能增强复杂碳水化合物的降解能力，支持宿主能量代谢。研究结论为：近岸驯化重塑了大珠母贝后代的肠道转录组和肠道微生物组，转录组分析揭示3007个差异表达基因，免疫相关基因显著上调，代谢和神经信号通路表达改变，表明对近岸胁迫因子的生理适应；16S rRNA测序揭示微生物群落组成变化和F1个体α多样性增加，功能预测提示代谢和免疫能力增强。这些协调的宿主–微生物组变化反映了多方面的适应性反应，为珍珠贝养殖中不断变化的环境条件下的选择性育种和健康管理提供了有价值的分子标记。