论文标题翻译 Genome sequence of the ornamental plant Digitalis purpureareveals the molecular basis of flower color and morphology variation 中文标题：观赏植物毛地黄（Digitalis purpurea）基因组序列揭示花色与形态变异的分子基础

《BMC Genomics》：Genome sequence of the ornamental plant Digitalis purpurea reveals the molecular basis of flower color and morphology variation

　　

毛地黄（Digitalis purpurea，foxglove）是一种广泛分布的观赏植物。在此，研究人员展示了基于长读长测序的紫红色花毛地黄植株的基因组序列及相应的基因模型预测。N50为4.3 Mbp的高组装连续性及约96%的完整BUSCO基因表明基因组完整性良好。该基因组资源为深入探究毛地黄的花色素沉积奠定了基础。研究人员鉴定了花青素生物合成的结构基因及相应的转录调控因子。通过比较紫红色花与白花植株，发现白花植株的花青素合酶（ANS）基因存在大片段插入，这可能导致该基因功能丧失，从而解释花青素色素沉积的缺失。此外，研究人员发现DpTFL1/CEN基因的大片段插入可能是导致大型顶生花发育的原因。

论文解读

研究背景与立题依据

毛地黄（Digitalis purpurea）是玄参科（Plantaginaceae）二年生植物，兼具药用价值（含强心苷类成分如地高辛）与观赏价值，其花部表型（花色、花序形态）存在显著自然变异——既有典型的紫红色管状花，也存在罕见的白花个体及顶生花突变体。花青素是决定花瓣紫红色的关键色素，其生物合成受结构基因（如CHS、ANS）与转录因子（MYB-bHLH-WD40复合物）共同调控；而花序形态则由TFL1/CEN类基因（磷脂酰乙醇胺结合蛋白家族）调控顶端分生组织命运。尽管前人已解析拟南芥（Arabidopsis thaliana）、金鱼草（Antirrhinum）等模式植物的花青素与花序调控机制，但毛地黄作为非模式植物，其基因组资源匮乏，花色与形态变异的分子基础尚未阐明。为此，研究人员开展本研究，旨在构建毛地黄参考基因组，解析花色丢失与顶生花发育的遗传机制，相关成果发表于《BMC Genomics》。

关键技术方法概述

研究选取1株紫红色花毛地黄（DR1）进行纳米孔（Nanopore）长读长测序，采用NextDenovo组装、NextPolish抛光，结合BUSCO评估完整性；通过GeMoMa基于同源物种（如芝麻、丹参）注释与RNA-seq数据预测基因模型，利用KIPEs、MYB_annotator、bHLH_annotator分别注释花青素通路酶与转录因子；通过LAST、JCVI进行共线性分析与Ks分布拟合以推断全基因组加倍（WGD）事件；利用bcftools进行变异检测与杂合度估算；通过IGV可视化比对红花与白花、顶生花与非顶生花个体的基因组差异，筛选候选变异；最终通过PCR对89株群体植株进行基因型分型验证。

研究结果

基因组序列组装与结构注释

研究人员通过NextDenovo组装获得毛地黄基因组（DR1_v1），大小为940 Mbp，contig N50达4.3 Mbp，包含35,916个蛋白编码基因，BUSCO完整性评估显示96.3%的完整基因（其中55.8%为重复拷贝）。Ks分析显示两个主要峰（Ks=0.21与0.34），表明近期存在全基因组加倍事件；杂合度估算为0.7%，RepeatMasker与EDTA注释显示重复序列占比73.8%（其中LTR反转录转座子占62.4%）。

毛地黄花青素生物合成基因

研究人员通过KIPEs注释鉴定了花青素生物合成全路径结构基因（CHS、CHI、F3H、F3’H、DFR、ANS、arGST、3GT），并通过系统发育分析与表达谱（kallisto定量）验证其功能活性——多数基因在花瓣（尤其紫色斑点区域）高表达。转录调控方面，鉴定到MYB75/PAP1 lineage（DP112203、DP103545、DP109418）、MYB123/TT2 lineage（DP102989）及bHLH42/TT8（DP105472）、WD40/TTG1（DP111282）同源基因，构成经典MYB-bHLH-WD40（MBW）调控复合物。

花色丢失的分子机制

通过比对2株红花（DR1、DR2）与2株白花（DW1、DW2）植株的基因组，研究人员发现白花植株的ANS基因（花青素合酶，催化无色花色素转化为有色花青素的关键酶）第4外显子存在~1.6 kb插入片段。该插入片段经hmmscan比对Pfam数据库，鉴定为LTR反转录转座子（LTR_Finder验证两端LTR序列，Kimura双参数模型估算插入时间约为0.21 Ks对应的进化时期）。插入导致ANS基因阅读框移位，推测使其功能丧失，无法合成花青素，从而导致白花表型。群体PCR验证显示，所有白花植株均携带该插入，而红花植株无此变异。

顶生花发育的遗传基础

针对罕见的顶生花突变体（花序顶端发育为大型辐射对称花，而非正常的管状花），研究人员比对1株顶生花与3株非顶生花个体的基因组，发现顶生花植株的DpTFL1/CEN基因（拟南芥TFL1同源基因，抑制顶端花发育）存在~106 bp缺失突变。该缺失位于基因编码区，导致蛋白C端结构域截短。已知金鱼草CEN、拟南芥TFL1功能丧失会导致顶生花，故推测该缺失破坏了DpTFL1/CEN的蛋白互作或信号转导功能，使顶端分生组织提前启动花发育程序。

讨论与结论

讨论部分指出，毛地黄基因组中高比例的重复序列（73.8%）与近期全基因组加倍事件（Ks=0.21）可能为其表型多样性提供了遗传基础；ANS基因的转座子插入是白花表型的因果变异，且为可逆突变（区别于某些物种的不可逆基因丢失）；DpTFL1/CEN的缺失突变则揭示了花序形态调控的保守性（与拟南芥、金鱼草机制一致）。

研究结论表明：毛地黄参考基因组的构建填补了该物种的基因组资源空白，为药用成分（如强心苷）与观赏性状的遗传改良提供了平台；花色丢失由ANS基因转座子插入导致功能丧失解释，顶生花发育由DpTFL1/CEN基因缺失突变解释，这些发现为植物花部表型多样性的分子机制研究提供了新案例。