《Plant Biotechnology Journal》：A Natural LTR Retrotransposon Insertion in the Promoter of GhNAC140-Dt Boosts Cotton Lint Yield
编辑推荐：
本研究揭示了棉花驯化中一个关键产量性状的分子调控新范式。作者发现，一个自然发生的LTR-Copia型反转座子插入到次级细胞壁NAC转录因子GhNAC140-Dt的启动子区域，其LTR序列中的60 bp cis-调控模块能特异性招募转录因子GhMYB46，从而显著增强下游基因的表达。GhNAC140-Dt的过表达可激活GhCESA4和GhCOBL9等纤维素合成相关基因，促进纤维素沉积，最终提高衣分。群体遗传学分析表明，该插入事件在半野生陆地棉品系中大量出现，并在现代栽培品种中固定率>80%，是棉花驯化过程中人工选择提高产量的直接分子证据。这项工作加深了对“转座子驱动产量演化”的理解，为作物分子设计育种提供了新思路。
文章内容归纳总结：
棉花是全球重要的经济作物，是天然纺织纤维的主要来源，也提供植物油和蛋白质。提高衣分是棉花育种的关键目标之一。转座子（Transposable Elements, TEs）是作物进化和驯化的基本驱动力，但其介导基因活化的具体机制仍不清楚。本文深入剖析了一个天然发生的LTR（Long Terminal Repeat）反转座子插入如何成为提升棉花产量的关键分子开关。
1. GhNAC140-Dt在纤维发育和产量中的重要作用
研究团队在前期工作中发现，NAC转录因子GhNAC140在棉花纤维发育中特异性高表达。进一步分析表明，在陆地棉标准系TM-1中，GhNAC140的表达在开花后15-25天的纤维中达到峰值，这恰好是纤维次生壁（Secondary Cell Wall, SCW）加厚的关键时期，暗示其在SCW形成中扮演要角。值得注意的是，在四倍体陆地棉的At和Dt亚基因组中，GhNAC140的两个同源基因表达差异显著，其中GhNAC140-Dt是优势表达基因。
对TM-1中GhNAC140-At和GhNAC140-Dt启动子的克隆测序揭示，在GhNAC140-Dt翻译起始密码子上游1634 bp处，存在一个5108 bp的插入序列。在线程序RepeatMasker分析确认，该插入序列是一个Ty1/Copia型反转座子。研究人员开发了特异性分子标记，在一个包含225份陆地棉种质的自然群体中进行检测，发现181份材料存在该反转座子插入，44份没有。因此，GhNAC140-Dt存在两个单倍型：带有插入的Hap1和不带插入的Hap2。
利用该群体在15个不同环境下的表型数据，研究人员比较了两种单倍型与棉花产量及纤维品质性状的关联。结果显示，携带反转座子插入的Hap1单倍型，在多个环境中与显著提高的衣分和衣指相关联，而与纤维长度、强度和马克隆值等品质性状的相关性较弱或不显著。这表明GhNAC140-Dt是影响棉花衣分这一关键产量性状的重要基因。
2. 天然LTR插入提升启动子活性
为了探究反转座子插入对基因表达的影响，研究人员构建了包含/不包含该插入的GhNAC140-Dt启动子片段，以及GhNAC140-At的启动子片段，分别驱动GUS（β-glucuronidase）报告基因。通过在20 DPA棉花纤维中进行瞬时转化和GUS活性检测，他们发现含有完整反转座子插入的长启动子片段驱动的GUS表达活性显著高于不含插入的短片段，而两个短片段之间活性无显著差异。这直接证明，反转座子插入增强了GhNAC140-Dt启动子的活性。
对该Copia型反转座子的结构分析显示，其两端各有一个459 bp的长末端重复序列LTR1和LTR2，两者序列高度相似。GUS实验证实，单独的LTR1或LTR2序列均表现出很强的启动子增强活性。进一步的实验表明，完整的反转座子与单独的LTR2在增强启动子活性方面效果相当，这说明LTR2是赋予转录激活功能的关键元件。通过对LTR2进行精细分段，研究人员最终将其核心功能区域定位到一个150 bp的片段LTR2-2上。
3. LTR2-2特异性招募GhMYB46以提升下游基因表达
接下来，研究旨在阐明LTR2-2增强转录的分子机制。由于该150 bp片段在酵母单杂交系统中表现出强烈的自主激活活性，阻碍了互作蛋白的筛选，研究人员转而通过文献调研，收集了十余个在过去十年中被报道与棉花纤维次生壁发育相关的转录因子作为候选。
利用双荧光素酶报告系统进行筛选，发现GhMYBL1、GhMYB7、GhFSN1、GhMYB46、GhFSN5、GhKNL1和GhMYB52-Like这七个转录因子能与LTR2-2序列发生互作。然而，当用不含LTR插入的原始GhNAC140-Dt启动子序列进行检测时，GhMYB46无法与之结合，而其他六个因子仍可结合。这表明，GhMYB46是被LTR2-2特异性招募的转录因子。
进一步的GUS活性对比实验显示，在七个候选转录因子中，GhMYB46蛋白对含有LTR序列的启动子的转录激活增强作用最为显著。通过将150 bp的LTR2-2进一步细分为三个60 bp片段，并结合酵母单杂交、电泳迁移率变动分析（EMSA）等实验，最终确认GhMYB46特异性结合在LTR2-2的第一个60 bp片段上。这些结果共同证明，LTR2-2通过其内部一个60 bp的cis-调控模块特异性招募转录激活因子GhMYB46，从而驱动下游基因的高水平表达。
4. 过表达GhNAC140-Dt增加皮棉产量
亚细胞定位实验证实GhNAC140-Dt是一个定位于细胞核的蛋白质。为了在体内验证GhNAC140-Dt的功能，研究人员构建了该基因的过表达和RNA干扰转基因棉花株系。田间表型调查显示，与野生型相比，GhNAC140-Dt过表达株系的衣分和衣指显著增加，而RNAi株系则明显降低。种子指数在转基因系和野生型间无显著差异，表明该基因主要调控纤维产量而非种子发育。
鉴于NAC转录因子在次生壁发育中的已知功能，研究人员检测了纤维细胞壁厚度。透射电镜分析显示，过表达株系在20 DPA和成熟期的纤维次生壁均显著增厚，而RNAi株系则变薄。相应地，纤维素含量测定也呈现出一致的变化趋势：过表达使纤维素沉积增加，干扰则减少。
然而，在纤维品质方面，GhNAC140-Dt的过表达导致了意料之外的结果：纤维长度变短，强度下降。扫描电镜观察揭示了原因：过表达株系的成熟纤维捻曲度显著降低，且纤维素微纤维排列更为分散和无序；而RNAi株系则表现出更高的捻曲度和更致密有序的微纤维排列。这表明，虽然GhNAC140-Dt的过表达通过增厚细胞壁提高了产量，但却因破坏了细胞壁的有序结构而损害了纤维强度，揭示了产量与品质之间存在权衡关系。
5. GhNAC140-Dt通过直接调控GhCESA4-At/Dt和GhCOBL9-At表达来调控细胞壁组织与生物合成
为阐明GhNAC140-Dt调控纤维发育的下游分子网络，研究人员对转基因株系和野生型的20 DPA纤维进行了转录组测序。基因本体（Gene Ontology, GO）富集分析显示，在过表达株系中共同上调的差异表达基因，以及在RNAi株系中共同下调的基因，均显著富集于“细胞壁组织或生物合成”等相关通路。
同时，研究人员利用DNA亲和纯化测序（DAP-seq）技术在全基因组范围内鉴定了GhNAC140-Dt蛋白的直接结合靶点。GO分析再次确认，这些靶基因显著富集于“植物类型次生壁生物合成”等类别。整合转录组和DAP-seq数据，研究人员锁定了三个共同基因：GH_D07G2262、GH_A07G2317和 GH_A08G0578。它们分别编码纤维素合酶4的Dt和At亚基因组同源基因GhCESA4-Dt、GhCESA4-At，以及一个已知的COBRA-like蛋白编码基因GhCOBL9-At。表达分析证实这三个基因在GhNAC140-Dt过表达株系中上调，在RNAi株系中下调。
通过双荧光素酶报告实验、EMSA和酵母单杂交实验，研究人员最终证实，GhNAC140-Dt蛋白能够直接结合到GhCESA4-Dt、GhCESA4-At和GhCOBL9-At基因的启动子特定区域，并激活它们的转录。已知GhCesA4和GhCOBL9是纤维强度的正调控因子，这解释了为何GhNAC140-Dt能通过激活这些下游基因来促进纤维素合成和细胞壁加厚。
6. GhNAC140-Dt启动子中反转座子插入的进化分析与育种利用
研究人员利用分子标记，在更广泛的棉花种质资源中追踪了这一反转座子插入事件的进化起源与分布。结果显示，在四倍体棉花的可能供体二倍体种（G. arboreum、G. herbaceum、G. raimondii）以及三个野生四倍体棉种（G. tomentosum、G. mustelinum、G. darwinii）中，均只存在不含插入的Hap2单倍型。
在七个陆地棉半野生种族中，最古老的G. hirsutum race yucatanense和G. hirsutum race palmeri同样只携带Hap2。然而，在更接近栽培种的G. hirsutum race latifolium和G. hirsutum race punctatum中，插入频率分别上升到45%和35%。在614份测试的陆地棉栽培种质中，Hap1的频率高达71.01%，显著高于半野生种族。而在269份海岛棉种质中，Hap1的频率仅为1.86%。这些数据表明，该反转座子插入起源于陆地棉种族向栽培种驯化的过渡时期，并在后续的驯化过程中被强烈地正向选择。
选择性清除分析为人工选择提供了直接的群体遗传学证据。对半野生材料、地方品种和现代品种的基因组数据分析显示，在包含该LTR插入位点的染色体Dt-12区域，半野生组与现代品种组、地方品种组之间具有高的遗传分化指数，且现代品种组在该区域的核苷酸多样性显著降低。这符合一个经历了选择性清除的基因组区域特征。
对中国四个主要棉区305个陆地棉商业品种的分析显示，Hap1优良等位基因在各大棉区的频率均超过80%。对中国不同年代育成品种的分析也显示，该优良等位基因的比例在不同育种时期均保持在高位。这表明，在追求高衣分的育种目标驱动下，这一由转座子插入创造的优异等位基因在中国棉花育种史上被持续地选择并保留了下来。
7. 讨论
本研究系统解析了一个Copia型反转座子插入通过重塑基因表达调控网络，从而驱动棉花产量性状进化的完整故事。其核心机制在于：LTR序列中的一个60 bp cis-调控模块，通过招募宿主自身的转录因子GhMYB46，将GhNAC140-Dt置于一个更强的转录调控之下。上调的GhNAC140-Dt进而激活包括GhCESA4和GhCOBL9在内的纤维素合成途径基因，促进纤维次生壁加厚，最终提高衣分。这一机制首次在作物中明确了由转座子衍生的转录开关的具体分子细节。
这项工作不仅加深了我们对“转座子驱动作物驯化”这一理论的理解，揭示了基因组“暗物质”如何通过提供新的调控元件来塑造重要农艺性状，也为棉花乃至其他作物的分子设计育种提供了宝贵的遗传资源和新的策略。例如，该LTR插入及其关联的分子标记可直接用于高衣分品种的辅助选择。此外，研究揭示的产量与品质的权衡关系也提示，未来通过CRISPR等基因编辑技术对包括GhNAC140-Dt在内的基因网络进行精准、协同的调控，可能是实现棉花产量与品质协同提升的可行途径。
