《Biology》：Comparative Plastid Phylogenomics of 29 Newly Sequenced Bougainvillea Cultivars Resolves the Phylogenetic Placement of Double-Bracted Cultivars
编辑推荐：
叶子花（Bougainvillea，紫茉莉科）品种因其鲜艳的苞片而在世界范围内具有极高的园艺价值；然而，频繁的种间杂交和体细胞突变掩盖了其母系起源，尤其是双苞类型。本研究对29个品种（16个双苞、13个单苞）的完整质体基因组进行了测序和组装，并与15个已发表的
叶子花（Bougainvillea，紫茉莉科）品种因其鲜艳的苞片而在世界范围内具有极高的园艺价值；然而，频繁的种间杂交和体细胞突变掩盖了其母系起源，尤其是双苞类型。本研究对29个品种（16个双苞、13个单苞）的完整质体基因组进行了测序和组装，并与15个已发表的分类群进行了比较和系统发育基因组学分析。所有质体基因组均呈现保守的四分体结构（153,994–154,542 bp，包含130–131个功能基因），其中B. × buttiana ‘Mahara Variegated’因经证实的移码诱导提前终止密码子而含有假基因化的ycf1，可作为该品种独特的诊断标记。序列差异主要集中于基因间隔区，鉴定出trnH-GUG–psbA和trnS-GCU–trnG-UCC可作为物种水平鉴定的潜在DNA条形码。相反，种内变异几乎不存在；13个双苞材料序列完全一致。系统发育基因组学分析很好地解析了种间主干（B. spectabilis、B. peruviana和B. glabra），但几乎所有双苞品种与单苞B. × buttiana谱系（包括‘Mrs. Butt’）一起落入未解析的多分支中。这证实了共同的母系起源，并表明双苞表型主要源于反复的体细胞芽变而非独立的杂交事件。此外，不同表型间质体基因组的一致性凸显了母系标记在品种界定中的固有分辨率限制。总体而言，本研究丰富了叶子花基因组资源，澄清了栽培类群间此前模糊的母系谱系，并证明品种水平的区分必须依赖核基因组标记而非质体基因组标记。
论文解读：基于29个新测序叶子花品种的叶绿体系统发育基因组学解析双苞品种的系统发育位置
一、研究背景与科学问题
叶子花（Bougainvillea，紫茉莉科）是全球热带和亚热带地区广泛栽培的观赏植物，以其色彩艳丽的苞片、耐旱性和抗污染能力而备受青睐。该属原产于南美洲，包含18个种，其中B. spectabilis、B. glabra和B. peruviana三个种贡献了绝大多数园艺品种，现代栽培品种主要源于这三个种之间的反复杂交。根据杂交背景，栽培品种通常归入三个杂交群：B. × buttiana（源于B. glabra与粉苞B. peruviana的杂交，以‘Mrs. Butt’为代表）、B. × specto-peruviana（B. spectabilis与B. peruviana的杂交）以及B. × specto-glabra（B. spectabilis与B. glabra的杂交）。在此基础上，通过自然和人工杂交以及自发和诱导突变，已培育出超过500个品种。
按苞片形态，品种分为单苞类型（一轮苞片）和双苞或多苞类型（多轮重叠苞片）。双苞品种最初仅4个——‘Cherry Blossom’、‘Los Banos Beauty’、‘Mahara’和‘Rose Ville Delight’，此后通过进一步突变增至20余个。然而，由于频繁的种间杂交和体细胞突变，双苞品种的母系起源长期存在争议。历史记录相互矛盾：例如‘Cherry Blossom’被认为是‘Los Banos Beauty’的芽变并与B. × buttiana相关，而‘Los Banos Beauty’据称源自B. glabra。此前的系统发育研究要么完全未包含双苞品种（如Bautista等的研究），要么仅取样3个双苞品种（如Wu等对华南品种的研究），远不足以代表该类群的表型多样性。因此，双苞形式究竟源于单一突变事件还是多次独立起源，以及具体由哪些单苞形式衍生而来，这些关键进化问题至今悬而未决。
鉴于双苞品种几乎完全依靠营养繁殖维持，其母系遗传的质体基因组为追踪母系祖先提供了稳定记录。种间杂交和母系替换涉及祖先物种间长期进化分化积累的大量固定序列差异，而体细胞突变在园艺时间尺度上仅产生极少量散在变异。因此，全质体组序列能够有效区分深层母系谱系与品种内体细胞漂变。基于此，研究人员提出两个可检验的假说：（1）双苞品种是源于不同杂交谱系中的多次独立事件，还是共享单一母系背景；（2）完整质体组能否提供足够的系统发育分辨率来澄清这些母系杂交谱系。
二、研究设计与主要技术方法
为检验上述假说，研究人员选取了16个双苞品种和13个单苞品种，覆盖已知双苞形式几乎全部表型范围以及三个主要杂交群的代表性单苞品种，另从GenBank获取15个已发表的叶子花及相关类群质体基因组作为参考。采用改良CTAB法提取总基因组DNA，在DNBSEQ-T7平台进行双端测序；使用NOVOPlasty v4.3.1进行质体基因组组装，GeSeq进行注释，tRNAscan-SE v2.0.7校验tRNA基因。通过Bowtie 2将reads比对回组装结果以验证组装质量。系统发育分析采用最大似然法（IQ-TREE v3.1.1）和贝叶斯推断（MrBayes v3.2），基于完整质体组、trimAl修剪后质体组和串联蛋白编码序列（CDS）三种矩阵。利用DnaSP计算核苷酸多样性（Pi），mVISTA进行序列保守性比较，IRscope比较IR/SC边界。对‘Mahara Variegated’的ycf1假基因化突变，通过高深度read映射、PCR扩增、分子克隆和Sanger测序进行多重验证。
三、主要研究结果
3.1 质体基因组结构特征
29个新组装的质体基因组均呈现典型的四分体结构，包含大单拷贝区（LSC）、小单拷贝区（SSC）和一对反向重复区（IRa和IRb）。基因组长度为153,994–154,542 bp，其中16个品种共享最长154,542 bp。GC含量在全基因组水平（36.5–36.6%）及各区域均变化极小。每个质体基因组携带130–131个功能基因，包括85–86个蛋白编码基因（PCG）、37个tRNA和8个rRNA。IR区包含18个双拷贝基因。14个基因含单一内含子，3个基因（pafI、rps12和clpP1）含两个内含子。所有质体基因组中rps19横跨LSC/IR连接处，ycf1和ndhF横跨IR/SSC连接处。多个品种序列完全一致：‘Flame’与‘Ruby’、‘Chili Orange’与‘Chili Yellow’、五个品种组（‘Elizabeth Angus’、‘Mrs. Eva White’、‘Mrs. Eva’、‘Barbara Karst’和‘Pedro’）以及13个品种的B. × buttiana组。值得注意的是，‘Mahara Variegated’因移码突变引入提前终止密码子而导致ycf1假基因化，高深度read映射（突变位点深度2045×）显示99.85%的reads支持该突变，并经克隆和Sanger测序确认这是真实的、固定的同质突变而非测序伪迹或异质性。
3.2 系统发育关系
基于44个质体组的系统发育分析（ML和BI）显示，叶子花属内分为两个高支持度的主要分支。分支I包含B. pachyphylla、B. spectabilis和B. peruviana，其中B. pachyphylla位于基部，B. spectabilis及其品种（‘Flame’和‘Ruby’）与野生B. peruviana及四个‘Chili’品种形成姐妹亚支。分支II包含其余野生和栽培类群，B. spinosa与所有其他成员姐妹，随后分化为野生种谱系和栽培种谱系。栽培种谱系进一步分为B. glabra分支（含‘Pedro’、‘Barbara Karst’、‘Mrs. Eva’、‘Mrs. Eva White’和‘Elizabeth Angus’）和B. × buttiana分支。关键发现是，所有16个双苞品种均落入B. × buttiana分支，与单苞品种‘Audrey Grey’、‘Mrs. Butt’等聚在一起。然而，B. glabra和B. × buttiana分支内部节点塌缩为多分支，反映序列差异极小。质体单倍型网络分析显示29个品种仅分为5个单倍型，16个双苞品种中15个共享主要单倍型Hap_1，‘Mahara’仅差一个突变步（Hap_2），直观证实双苞品种共享几乎不变的B. × buttiana母系背景。
3.3 序列相似性与结构变异
同种或同杂交群内所有品种SNP差异为零，表明严格的母系谱系保守性。InDel则捕获了密切关联品种间的微小突变：18个B. × buttiana品种间InDel数为1–14个，B. peruviana内最多达51个。B. × buttiana与B. glabra仅相差7个固定SNP和30–38个InDel，显示密切亲缘关系；而两者与B. peruviana和B. spectabilis的差异超过910个SNP和2300个InDel。B. peruviana与B. spectabilis为近缘姐妹类群，相差379个SNP和1100–1150个InDel。序列变异主要集中于LSC和SSC区域的基因间隔区（IGS）。全基因组核苷酸多样性扫描鉴定出十个高变区，其中trnH-GUG–psbA（Pi = 0.02293）和trnS-GCU–trnG-UCC（Pi = 0.02288）为最突出的两个位点，可作为物种水平鉴定的DNA条形码。然而，仅在B. × buttiana分支内计算时，全基因组Pi值降至接近零（0–8.4 × 10?4），表明这些种间热点区域在双苞品种间缺乏多态性，质体标记无法区分该复合群内的大多数品种。
3.4 IR/SC边界比较
17个质体基因组的四个四分体连接处均存在不同程度的IR扩张和收缩。LSC/IRb连接处rps19通常跨越边界（165 bp进入LSC，114 bp进入IRb），但不同类群分裂位置有所变化。IRb/SSC连接处ndhF长度为2223–2259 bp，其中‘Mahara Variegated’和‘Mahara’的ndhF最长（分别有49 bp和47 bp进入IRb）。SSC/IRa连接处位于ycf1基因内部，ycf1延伸进入IRa 1332–1371 bp，导致IRb区域存在相应长度的部分重复假基因。‘Mahara Variegated’表现出独特模式：跨越SSC/IRa连接的完整ycf1拷贝被假基因化，而IRb衍生的重复片段仅为ycf1的部分拷贝。IRa/LSC连接处位于rpl2–trnH间隔区，trnH与边界的距离在不同类群间变化于10–26 bp。
四、讨论与结论
讨论部分要点：质体基因组的高度保守性反映了长期栽培和选择下基本架构和碱基组成的强约束。多个品种组序列完全一致反映了共享母系谱系，也凸显了质体基因组在品种水平的分辨率极限。在木本园艺植物中，不同品种常源于近期体细胞芽变或共享同一母系谱系但花粉供体不同的杂交事件，由于质体基因组进化速率慢且为单亲母系遗传，在如此短期的表型多样化过程中保持不变。‘Mahara Variegated’的ycf1假基因化是固定的谱系特异性突变，可作为该品种的明确分子标记；但其他花叶品种（如‘Mrs. Eusenia Raja Singh’和‘Lizard Double Orange’）保留完整ycf1编码序列，说明ycf1假基因化并非叶子花叶片花叶现象的普遍因果驱动因素。系统发育分析明确将‘Los Banos Beauty’的母系谱系置于B. × buttiana而非此前认为的B. glabra，修正了历史分类学分歧。然而，仅基于母系标记的分类学修订需谨慎解读，叶绿体捕获、不完全谱系分选和细胞核-细胞质不一致等进化过程可能导致质体基因树偏离真实物种系统发育。未来需整合核基因组系统发育学（如Angiosperms353 Hyb-Seq或全基因组重测序）、细胞遗传学和倍性评估以及整合分类学方法。
结论翻译：总之，本研究建立了全面的质体系统发育基因组学框架，阐明了栽培叶子花的母系起源和进化历史，特别聚焦于双苞形式。比较基因组分析表明，双苞表型主要源自B. × buttiana母系谱系通过反复体细胞芽变而非多次独立杂交事件产生。虽然高序列保守性凸显了母系质体标记在精细品种界定中的固有分辨率限制，但特定的基因组特征——尤其是‘Mahara Variegated’中移码诱导的ycf1假基因以及鉴定出的高变基因间隔区——为种质筛选提供了有前景的诊断标记。总体而言，这些发现解决了栽培双苞叶子花母系谱系中长期存在的模糊性，丰富了宝贵的基因组资源，并强调了未来在该观赏属的品种区分、种质保护和分子育种中整合核基因组方法的必要性。
