全基因组多样性与Pinna rudis渐渗揭示Thau潟湖是崩溃中的Pinna nobilis保护储备库

《Frontiers in Genetics》：Genome-wide diversity and Pinna rudis admixture reveal Thau lagoon as a conservation reservoir for the collapsing fan mussel Pinna nobilis

　　

Pinna nobilis的保护日益依赖于Thau潟湖种群；该种群在人口统计上保持稳定，并在很大程度上避开了摧毁地中海开放海种群的大规模死亡事件（mass mortality events, MMEs）。在此，研究人员评估法国Thau和Diane潟湖的遗传多样性，并刻画其与Pinna rudis的全基因组混合（admixture）。Thau的线粒体多样性，尤其是在蛋白编码基因上，与先前在开放海种群中记录的结果高度相似；历史样本与当代样本的比较显示16S rRNA单倍型组成存在差异，这与近期线粒体多样性侵蚀相一致。扩展的线粒体和核单核苷酸多态性（single nucleotide polymorphism, SNP）数据证实并细化了关于连续历史种群扩张的早期证据。研究人员还报告了P. nobilis中线粒体异质性（heteroplasmy）的首个证据，该证据在一个同时携带P. nobilis和P. rudis线粒体基因组的个体中鉴定到。相反，五个经确认的F₁杂交个体均仅保留P. nobilis线粒体基因组，表明杂交过程中存在线粒体不对称遗传。核基因组分析揭示与P. rudis的低水平但广泛混合，这与有限但持续基因流的历史相一致。候选渐渗（introgression）区域富集了参与先天免疫识别与信号传导、凋亡、细胞黏附、细胞外基质组织和蛋白稳态（proteostasis）的基因，识别出可能与宿主韧性（host resilience）相关并值得功能研究的通路。然而，混合个体中不存在完全疾病抗性表明，这些区域对疾病韧性的任何潜在贡献都可能较小且依赖情境。这些结果表明，尽管P. nobilis近期发生种群崩溃，其仍保留显著进化潜力；Thau潟湖保存的现存遗传多样性和混合可能构成保护的重要资源。未来将基因组变异与个体感染结局联系起来的研究，对于确定已识别的候选变异是否直接促成疾病耐受和抗性将至关重要。

研究背景与问题：Pinna nobilis是地中海特有的大型双壳类，寿命可超过45年、体长可达1.2 m，作为滤食者可改善水质并为相关物种提供栖息地；其幼虫可在水柱中扩散约20天。过去十年，大规模死亡事件（mass mortality events, MMEs）使开放海种群急剧衰退，而部分潟湖和半封闭环境种群得以持续。尽管病原、环境条件、宿主密度及个体/种群易感性差异可能共同导致空间异质性，其机制仍未解决。反复MMEs预计会造成等位基因和单倍型多样性丧失、稀有变异频率随机变化及有效种群大小（effective population size, Ne）下降；但持续存活的种群也可能保留遗传多样性及与存活相关的变异。Thau潟湖是P. nobilis最后的自然储备之一，研究人员因此评估其是否保留历史开放海遗传多样性、与P. rudis的混合及其人口历史，为保护与恢复提供依据。

研究开展与意义：研究人员比较Thau和Diane潟湖及开放海样本，结合线粒体基因和全基因组SNP，评估遗传多样性、群体结构、与P. rudis的全局和局部混合、线粒体异质性、F₁杂交及候选渐渗区功能，并重建人口历史、检测局部适应和近期正选择。结果表明Thau仍保留可观遗传多样性，存在低水平但广泛的P. rudis混合及与免疫、黏附、应激和蛋白稳态相关的候选渐渗区，但未发现完全疾病抗性；Thau可作为保护遗传资源，未来需连接基因组变异与感染表型。

关键技术方法（≤250字）：样本来自法国Thau潟湖（n=301）和Diane潟湖（n=33）2021—2025年334个成年个体非致死外套膜活检，并整合西班牙、意大利和希腊开放海既往样本45个，共379个；采用Illumina双端（2×150 bp）约10×重测序，检测线粒体基因与核SNP；用ADMIXTURE、RFMix和主成分判别分析（discriminant analysis of principal components, DAPC）分析全局/局部混合，SMC++与BEAST2推断人口历史，pcadapt筛选中性SNP，滑窗遗传分化指数（F_ST）、核苷酸多样性（π）和Tajima’s D检测选择，候选渐渗区行基因本体（Gene Ontology, GO）富集。

研究结果：
3.1 DNA sequencing and variant filtering（DNA测序与变异过滤）：通过Illumina重测序和变异检测，核基因组平均reads为62,072,427（SD=16,510,421），线粒体为73,370（SD=39,736），约10×；识别104,900,089个核变异，过滤后保留50,069,828个双等位SNP，稀疏化后为1,078,378个SNP，说明数据量满足全基因组分析。
3.2 Mitochondrial genetic diversity and structure（线粒体遗传多样性与结构）：通过16S rRNA和蛋白编码基因（protein-coding genes, PCG）单倍型分析，单倍型多样性（haplotype diversity, Hd）高（0.932–1.000）、核苷酸多样性（nucleotide diversity, π）低（0.00286–0.00424），PCG中非同义/同义替换关系显示d_S>d_N，提示纯化选择；各地理组无显著地理结构，单倍型网络呈星状；Thau与开放海多样性相当；历史与当代16S rRNA单倍型比较提示线粒体多样性侵蚀；5个F₁杂交个体均携带P. nobilis线粒体基因组；1个个体检出全线粒体基因组异质性，其核混合仅0.23%。
3.3 Pinna nobilis genetic diversity and structure using neutral SNPs（基于中性SNP的P. nobilis遗传多样性与结构）：基于194,964个中性SNP，观测杂合度（observed heterozygosity, H_O）均值为0.154，近交系数（F_IS）低且多为负；Thau内F_ST均值为0.00061且不显著，Thau与开放海及Diane比较显著但分化低（F_ST<0.014），与西班牙较近，开放海种群间F_ST可达0.026，提示Thau内部近似随机交配并与开放海部分隔离。
3.4 Demographic history（人口历史）：Tajima’s D和Fu and Li’s D多为负；16S rRNA错配分布呈单峰，PCG错配分布较不规则；BEAST2估计历史扩张约800 kya（千年前）发生于Diane和开放海、约500 kya发生于Thau，PCG估计约250 kya；SMC++显示有效种群大小（effective population size, Ne）约200,000—20,000年前上升、随后下降、约30,000—3,000年前瓶颈后快速扩张，提示多阶段人口历史。
3.5 Admixture analyses（混合分析）：ADMIXTURE在K=2时交叉验证误差最低，鉴定5个F₁、20个2.95%—13.43%混合个体、14个0.04%—0.6%低混合个体及339个未检出混合个体；混合比例>2.9%者约占7%，其中Thau为7/301（2.3%）且靠近开放海连接处，Diane为15.15%，意大利开放海为50%，西班牙开放海为25%（不含F₁）。RFMix显示全部P. nobilis个体均携带候选混合片段，41%窗口至少在1个个体中检出混合，个体平均609个候选区、约占基因组2%；269个渐渗热点高度共享；共享热点富集细胞连接、蛋白折叠分子伴侣（chaperone）及清道夫/货物受体活性，Thau另富集细胞黏附、Diane富集受体介导功能；共享热点含先天识别、凋亡/自噬、黏附/细胞外基质及应激/蛋白稳态基因。DAPC在K=8显示主要共享簇及开放海与潟湖间的簇频率差异。
3.6 Selection scans（选择扫描）：通过滑窗F_ST、π和Tajima’s D检测Thau局部适应，共732个窗口显示信号，保留111个由≥3个连续异常窗口构成的区域、含255个编码序列（coding sequences, CDS）；GO未通过多重检验，但候选基因涉及应激/DNA修复（HSP12A/B）、免疫识别（TLR4、DDX58-like、IRF1-like、NACHT结构域及凋亡基因）、溶质运输/代谢及脂质/膜/细胞外基质/信号传导；近期正选择检出LH3-like和EP4-like基因；少数局部适应区与渐渗热点重叠，涉及膜运输、蛋白周转和染色体组织，提示调控性而非直接应激反应途径。

讨论部分总结：保护形势紧迫，H. pinnae仍可检出且其他病原可能参与；残存种群多限于潟湖等环境，连通性受限，保护 increasingly 依赖主动干预。Thau线粒体多样性接近开放海但16S rRNA提示近期侵蚀；核SNP多样性高、内部无显著结构，与开放海低但显著分化、与西班牙较近，提示其为部分隔离的遗传单元。人口历史呈多阶段扩张—瓶颈—恢复；局部适应信号分散于应激、免疫、代谢、运输及膜/细胞外基质基因，并可能受人口历史影响。线粒体异质性极罕见，F₁仅保留P. nobilis线粒体，提示不对称遗传；ADMIXTURE与RFMix结果差异支持古老渐渗而非近期杂交，热点共享且富集免疫、黏附和蛋白稳态相关基因，但混合个体无完全抗病，贡献可能有限且依赖情境。未来需连接基因组变异与感染表型，并评估辅助基因流。

研究结论部分翻译：这些结果表明，尽管P. nobilis近期发生种群崩溃，其仍保留显著进化潜力；Thau潟湖保存的现存遗传多样性和混合可能构成保护的重要资源。未来将基因组变异与个体感染结局联系起来的研究，对于确定已识别的候选变异是否直接促成疾病耐受和抗性将至关重要。