超级基因与基因组嵌合驱动加拉帕戈斯树雀的表型演化

《Nature Communications》：A supergene and genomic mosaicism are associated with phenotypic evolution in Galápagos tree finches

　　

适应辐射（adaptive radiations）为研究表型新颖性和物种多样性的起源提供了强大模型。达尔文雀（Darwin's finches）是适应辐射的经典例子，其中喙形态（beak morphology）是与适合度相关的关键性状。在此，研究人员聚焦于Camarhynchus树雀，通过分析全基因组并结合形态学数据，重新审视了在遗传证据出现之前提出的关于其起源的长期问题。Camarhynchus在过去10万年间发生多样化，尽管遗传分化程度较浅，却产生了显著令人瞩目的形态差异。研究人员鉴定出多个高度分化的位点，其中包括G03——一个影响喙和体大小的超级基因（supergene）。研究人员揭示了一个显著的基因组嵌合（genomic mosaicism）案例，其中两个同域物种之间的常染色体分化主要限于一个接近固定的单一基因座——G03——使其从一个表型最适状态转移到另一个。在更广泛的系统发育框架中，研究人员还揭示了G03单倍型的意外多样性，并表明小与大G03等位基因在辐射早期即发生分化，提示这是达尔文雀物种间实现生态扩张的关键遗传创新。研究结果进一步将Camarhynchus确立为伴随基因流的物种形成（speciation with gene flow）的模型范例，强调了杂交在产生生物多样性中的广泛作用。

研究人员在《自然·通讯》（Nature Communications）上发表了关于加拉帕戈斯树雀（Camarhynchus）的基因组研究，结合全基因组测序与博物馆标本形态数据，揭示了超级基因和基因组嵌合现象如何驱动表型演化。适应辐射（adaptive radiation，指一个祖先物种快速分化为多种生态类型）是研究新表型起源和物种多样性的理想体系。达尔文雀是适应辐射的经典代表，喙形态（beak morphology）是与适合度直接相关的关键性状。以往研究多聚焦于Geospiza地雀，而树雀的早期演化推测长期停留在形态学层面，例如David Lack提出的谱系关系和扩散路径缺乏遗传证据检验。为此，研究人员对树雀进行了迄今最全面的基因组分析。

研究人员分析了415只达尔文雀的全基因组数据，其中包括77只Camarhynchus树雀个体，覆盖所有五个物种和大部分亚种；并加入来自小安的列斯群岛的外群Loxigilla noctis（n=3）和L. barbadensis（n=5）。通过比对到Camarhynchus parvulus参考基因组，获得1120万个双等位基因SNP。主要结论为：树雀在过去约10万年内发生多样化，形态差异显著但遗传分化极浅；在多个高度分化位点中，G03超级基因（supergene，指一个区域内多个连锁基因共同控制复杂表型）是影响喙与体大小的核心；西部psittacula（C. p. affinis）是一个镶嵌谱系，其常染色体背景几乎完全来自parvulus，但G03位点固定了大型（L）单倍型；G03的L与S单倍型在辐射早期（约0.9–1.1百万年前）即已分化。这些发现将Camarhynchus确立为基因流（gene flow，指种群间遗传物质交流）伴随物种形成的典范，表明杂交在产生生物多样性中的重要作用。

关键技术方法包括：全基因组重测序（平均深度15.8×）和Sentieon流程变异检测；利用tsinfer构建祖先重组图（ancestral recombination graph, ARG），并计算累积分布统计（genealogical nearest neighbor, GNN）；基于10 kb窗口计算F_ST、d_xy和π评估分化与多样性；采用SplitsTree构建系统发育网络、IQ-TREE构建最大似然树；通过Dsuite进行ABBA-BABA检验和距离分数（distance fraction, df）扫描检测局部基因流；使用GEMMA进行混合映射；使用正常混合模型（normal mixture model, NMM）对348个博物馆标本的形态数据进行无先验分组。

结果部分如下：

**尖喙树雀的进化历史——woodpecker（C. pallidus）和mangrove（C. heliobates）雀**：系统发育网络显示pallidus三个亚种与heliobates相互嵌套，分支模式对应从东到西的地理分布；GNN和ABBA-BABA分析检测到西部pallidus与heliobates之间的杂交个体与近期基因流。

**钝喙树雀的进化历史——小、中、大树雀**：基因组扫描表明钝喙树雀间的遗传分化主要局限于1A、5、6和Z染色体上的少数位点，其中G03（位于1A）的ZF_ST偏离基因组平均值约25个标准差。西部psittacula和parvulus之间的常染色体F_ST仅1.7%，但它们在G03上几乎固定不同等位基因。

**西部psittacula是一个镶嵌谱系**：G03区系统树和TWISST分析显示，西部psittacula在G03区域与北部psittacula聚类，而在去除G03后的常染色体树上则与parvulus聚类；df分析支持西部与北部psittacula在G03共享大量等位基因，其L单倍型序列与尖喙雀更为相似，提示该单倍型可能来源于尖喙雀。

**混合映射揭示树雀遗传分化位点**：用GEMMA对parvulus、psittacula（排除西部psittacula）和尖喙雀进行两两混合映射，G03是每个比较中最显著的关联位点，同时鉴定出15个新数量性状位点（quantitative trait locus, QTL），命名为C01–C15。

**达尔文雀中G03单倍型广泛变异**：G03单倍型网络揭示多种L和S单倍型；L与S主支的分化时间估计为0.9–1.1百万年，先于多数达尔文雀物种分化。Cocos雀在G03局部系统树上与小地雀和小树雀更近，但在全基因组树上位于较早分支。

**树雀在达尔文雀系统发育历史中的位置**：全基因组最大似然系统树支持北部psittacula最早分化，西部和中部psittacula与parvulus、pauper形成单系；Z染色体树则显示pauper与西部psittacula为姐妹群，表明基因流和不完全谱系分选影响了系统发育推断。

讨论部分：研究人员认为树雀间的遗传分化水平通常处于亚种之间（F_ST约2%–15%），这使传统分类学界定变得困难。例如，全基因组树提示北部psittacula可能应被视为独立物种C. habeli，而C. pallidus San Cristóbal也可能达到种级分化。G03超级基因在树雀和地雀中都控制喙与体大小，是跨生态位重复利用的关键遗传元件。与地雀研究不同，西部psittacula的例子表明，两个同域物种可以在仅一个常染色体位点（G03）上高度分化而其余基因组趋于同质化，这归因于反复的基因流与选择共同作用。研究结论指出：Camarhynchus辐射提供了伴随基因流的物种形成的典型实例，杂交和分化不是对立过程，它们在物种形成早期可相互协同产生演化新颖性，从而加深了对生物多样性生成机制的理解。