大蒜全长转录组鉴定AsMYB4作为黄酮代谢与雄性不育的候选调控因子

《Plant Physiology and Biochemistry》：Full-length transcriptome sequencing in garlic identifies AsMYB4 as a candidate regulator of flavonoid metabolism and male sterility

大蒜（Allium sativum L.）是一种全球性栽培作物，但其雄性不育特性严重限制了育种改良。本研究利用PacBio单分子实时（SMRT）测序技术，对大蒜六个器官（根、茎、叶、花薹、花和鳞茎）进行了全长转录组测序，共获得13,260条高质量一致性亚型（isoforms），平均长度2,079.15 bp，其中96.75%的序列获得了功能注释。研究人员鉴定出399个转录因子（TFs）、299个长链非编码RNA（lncRNAs）、1,121个含简单重复序列（SSRs）的转录本以及92个可变剪接（AS）事件。在此基础上，研究人员克隆并鉴定了AsMYB4基因，其编码一个含有两个保守SANT结构域的R2R3-MYB转录因子，定位于细胞核。系统发育分析表明AsMYB4与其他物种的MYB4同源蛋白高度保守。在拟南芥中异源过表达AsMYB4显著抑制了黄酮生物合成基因的表达，导致种子种皮颜色变浅（淡黄色）。更重要的是，转基因植株表现出花丝伸长受损、花粉活力降低（约18.15%的花粉异常）以及种子败育率升高（32.05%），而雌性生殖功能正常。进一步分析发现，AsMYB4过表达在雄蕊中抑制了苯丙烷-黄酮代谢途径基因的表达并降低了总黄酮含量，同时上调了吲哚-3-乙酸（IAA）水平，并下调了生长素转运、生物合成及信号转导相关基因的表达。这些结果表明，AsMYB4作为一个保守的黄酮生物合成抑制因子，通过扰乱苯丙烷-黄酮代谢和生长素稳态进而影响雄性生殖发育。本研究为大蒜提供了全面的转录组资源，并为理解AsMYB4在调控大蒜育性中的功能奠定了基础。

论文解读文章

研究背景与科学问题

大蒜（Allium sativum L.）是全球重要的经济作物，然而绝大多数栽培品种表现出严重的配子体雄性不育，几乎不能产生可育种子，因此长期依赖无性繁殖（鳞茎）。这种繁殖方式导致病毒积累、遗传多样性降低，严重制约了种质创新和新品种选育。尽管已有形态学和细胞学观察，但控制大蒜雄性生殖发育的分子机制仍知之甚少。转录因子，特别是MYB家族成员，在植物次生代谢（如苯丙烷和黄酮合成途径）、胁迫应答及发育调控中发挥核心作用，越来越多的证据也表明MYB基因参与雄配子体发育。然而，大蒜基因组极其庞大（约16 Gb）、高重复和高杂合，使得高质量的基因组组装和注释极为困难。因此，利用全长转录组测序技术绕过基因组限制，系统鉴定调控基因，对于揭示大蒜雄性不育的分子机制具有重要意义。

研究内容与主要结论

为突破大蒜基因组复杂性的瓶颈，研究人员采用PacBio单分子实时（SMRT）测序技术，从大蒜‘Y’品种的六个器官（根、茎、叶、花薹、花和鳞茎）中提取RNA并混合建库，构建了一个全面的全长转录组数据集。基于此资源，研究人员鉴定并克隆了先前通过比较转录组学预测的候选转录因子AsMYB4。通过异源转化拟南芥进行功能验证，研究人员证实AsMYB4是一个保守的R2R3-MYB转录抑制因子，不仅能够抑制黄酮生物合成，还通过扰乱苯丙烷-黄酮代谢和生长素（IAA）稳态，导致花丝伸长缺陷、花粉活力下降和雄性育性降低。这一发现将MYB4的功能从次生代谢拓展到生殖发育调控领域。

关键技术方法

研究人员采用了以下主要技术方法：

1. 1.
   PacBio SMRT全长转录组测序：从大蒜‘Y’品种的六个器官（根、茎、叶、花薹、花和鳞茎）中提取RNA并等量混合建库，在PacBio Sequel II平台上测序。
2. 2.
   生物信息学分析：使用SMRT Link处理原始数据获得环状一致性序列（CCSs）和全长非串联（FLNC）序列，经聚类和校正获得高质量一致性亚型；通过BLASTx比对多个数据库（NR、KOG、Swiss-Prot、KEGG）进行功能注释；使用ANGEL预测开放阅读框（ORFs）；利用PlantTFDB鉴定转录因子家族；使用CNCI和CPC算法预测长链非编码RNA（lncRNAs）；利用MISA检测简单重复序列（SSRs）；使用Cogent和SUPPA鉴定可变剪接（AS）事件。
3. 3.
   系统发育与亚细胞定位：使用ClustalW和MEGA7构建邻接法系统发育树；通过农杆菌介导在本氏烟草叶片中瞬时表达AsMYB4-GFP融合蛋白，结合核标记（AtH2B-mCherry）和膜标记（AtCBL-RFP）进行共聚焦显微镜观察。
4. 4.
   转基因功能验证：采用花浸法将35S::AsMYB4过表达载体转化拟南芥Col-0，筛选获得纯合T₃代株系进行表型分析。
5. 5.
   生理生化测定：采用NaNO₂-AlCl₃-NaOH比色法测定总黄酮含量；使用竞争性酶联免疫吸附测定（ELISA）试剂盒检测吲哚-3-乙酸（IAA）含量；采用亚历山大染色法评估花粉活力。

研究结果

3.1 大蒜全长转录组的构建与特征分析

通过对六个器官混合RNA进行PacBio SMRT测序，研究人员获得了16,419条全长非串联（FLNC）序列，经聚类和校正后得到13,260条高质量一致性亚型，平均长度为2,079.15 bp，N50为2,266 bp，GC含量为41.19%。其中12,830条（96.75%）在至少一个数据库中获得了功能注释，69.66%在所有四个数据库中被一致注释。物种分布分析显示，大多数转录本与天门冬属（Asparagus officinalis）蛋白相似度最高，符合单子叶植物的进化关系。功能分类表明，转录本主要富集于碳水化合物运输与代谢、信号转导以及翻译后修饰等通路。

3.2 调控元件与转录特征的鉴定

研究人员在大蒜转录组中鉴定出399个转录因子（TFs），分属47个家族，其中C3H（34个）、bHLH（29个）和MYB相关（27个）最为丰富。此外，鉴定出299个长链非编码RNA（lncRNAs）。简单重复序列（SSRs）分析发现三核苷酸重复最为丰富（67.71%），其中AAG/CTT基序频率最高（19.4%）。共检测到92个可变剪接（AS）事件，内含子保留（RI）为主要类型（53.26%），且多数基因（53.85%）有两个亚型。进一步分析发现，AS事件和SSRs在调控亚型（TFs和lncRNAs）中的比例显著高于其他功能注释亚型，暗示它们在转录调控中可能发挥更重要的作用。此外，还筛选出27个推定的反义lncRNA-mRNA配对。

3.3 AsMYB4的序列与进化特征

AsMYB4的全长编码区（CDS）为630 bp，编码210个氨基酸，预测分子量为23.82 kDa，等电点为6.32。系统发育分析显示AsMYB4与其他植物MYB4蛋白紧密聚类，且多数物种仅有一个MYB4同源物，表明该蛋白在进化中相对保守。多序列比对发现其N端含有两个高度保守的MYB DNA结合结构域（R2和R3），C端较为多变。亚细胞定位实验证明AsMYB4-GFP荧光信号仅在细胞核中检测到，符合其作为转录因子的预期定位。

3.4 AsMYB4过表达抑制拟南芥黄酮基因表达与种子色素沉积

在AsMYB4过表达的拟南芥株系中，实时定量PCR（RT-qPCR）分析显示黄酮生物合成关键基因（CHS、CHI、F3H、FLS1、DFR、ANS）的转录水平均显著低于野生型。成熟种子表型分析显示，过表达株系的种皮呈淡黄色，而野生型为棕色。用野生型花粉授粉后的过表达母本植株产生的种子仍保持淡黄色表型，表明该表型由母本孢子体组织决定，与胚胎基因型无关。这些结果证实AsMYB4异源过表达抑制了黄酮生物合成基因表达并减少了原花青素的积累。

3.5 AsMYB4过表达抑制拟南芥雄蕊发育

三个独立转基因株系的表型验证排除了插入位置效应。过表达株系的果荚中含有大量皱缩败育的种子（32.05%），显著高于野生型（2.41%）。花形态学观察显示转基因植株花丝伸长受损，导致花药与柱头接触减少。亚历山大染色结果表明，转基因植株中约18.15%的花粉异常。人工授粉实验证实，用野生型花粉给转基因植株授粉后可以正常结籽，表明雌性生殖功能正常，雄性缺陷是导致育性降低的主要原因。在雄蕊中，AsMYB4过表达导致苯丙烷途径基因（PAL1、PAL2、C4H、4CL1）和黄酮合成基因（CHI、F3H）下调，总黄酮含量显著降低。有趣的是，下游分支点基因出现差异化调控：黄酮醇合酶基因FLS1显著下调，而花青素合成基因DFR和ANS显著上调，提示代谢流可能发生了重新分配。同时，IAA含量显著升高，而生长素转运基因（PIN2）、生物合成基因（YUCCA6、TAA1）、信号组分（ARF6、ARF8、IAA4）以及IAA-酰胺合成酶基因（GH3.3、GH3.5、GH3.6）均显著下调，表明生长素稳态受到严重扰动。

讨论与结论

讨论部分总结：本研究利用PacBio SMRT测序成功构建了大蒜全长转录组，获得的转录本平均长度（2,079.15 bp）显著优于此前基于二代测序的研究，为基因结构注释和调控元件鉴定提供了更准确的资源。数据集揭示了多种调控特征（SSRs、AS、lncRNAs、TFs），其中AS和SSRs在调控亚型中的富集提示它们可能在转录调控中扮演重要角色。针对AsMYB4的功能研究表明，其在拟南芥中异源表达后，通过抑制黄酮合成基因（特别是CHS、CHI、F3H、FLS1）导致总黄酮含量下降，并引起种皮原花青素减少（淡黄色表型），这与已知的AtMYB4功能一致，表明MYB4对黄酮合成的抑制作用在进化上是保守的。在雄蕊发育方面，AsMYB4过表达导致花丝缩短和花粉活力降低，这与其对苯丙烷-黄酮途径的抑制有关，尤其是FLS1的下调可能影响了黄酮醇（如槲皮素、山奈酚）的供应，而黄酮醇是花粉活力和花粉管生长所必需的。同时，雄蕊中IAA水平的升高并非源于合成增加，而是可能与生长素转运（PIN2下调）和IAA共轭失活（GH3基因下调）受阻有关，黄酮类化合物本身也能调节PIN依赖的生长素运输。此外，特定黄酮组分的变化可能影响IAA氧化酶的活性，进一步加剧IAA积累。这些代谢和激素信号的共同紊乱最终导致了雄性育性缺陷。

研究结论翻译：本研究为大蒜提供了一个全面的全长转录组资源，并扩展了关于该物种可变剪接事件、简单重复序列、转录因子和长链非编码RNA的信息。利用这一资源，研究人员鉴定并功能表征了AsMYB4，证明其在异源表达时对黄酮生物合成具有保守的抑制作用，并进一步揭示了其通过扰乱苯丙烷-黄酮代谢和生长素稳态来调节雄性生殖发育的功能。总之，这些发现不仅丰富了大蒜现有的分子资源，还将对AsMYB4的功能认知从次生代谢拓展到了生殖发育调控领域。不过，该功能证据来源于拟南芥中的异源过表达，AsMYB4在其原生种大蒜中是否发挥类似的调控作用仍有待证实。未来需要在大蒜中进行稳定的遗传转化以及染色质免疫沉淀等相关研究来确定其天然功能并识别其直接的调控靶标。