陆地棉异域×精英群体中株型、发育与种子形态的遗传解析

《The Plant Genome》：Genetic dissection of plant architecture, development, and seed morphology in exotic × elite Upland cotton populations

　　

棉花（Gossypium hirsutum L.）育种长期以来主要强调纤维产量与品质，然而株型结构、发育进程及种子形态性状同样在生产力与栽培管理中发挥关键作用。研究人员利用异域×精英陆地棉群体，对其中4个此类性状的遗传基础进行了解析：株高（HT，plant height）、单果枝节位数（N1B，number of nodes with one branch）、开花天数（DTF，days to flowering）以及单粒种子平均表面积（ASAS，average surface area per seed）。群体特异分析与联合分析共鉴定出17个数量性状位点（QTL，quantitative trait loci），其中ASAS检测到的QTL数量最多。异域系T326在改良ASAS方面显示出最高潜力。尽管HT、N1B与DTF在这些群体中表现出的遗传变异有限，基因本体论（GO，Gene Ontology）富集分析仍指出CCR4–NOT复合体可能是连接ASAS与DTF的潜在共享调控枢纽。此外，ASAS与配套论文中报道的纤维强度及上半部平均长度表现出极显著正相关，并存在QTL区间重叠，提示这些基因组区域可能存在紧密连锁或多效性（pleiotropy）效应，且在商业育种上具有潜在益处。总体而言，数据表明，在改良纤维品质性状过程中，这些次级性状并不存在普遍的负向等位基因权衡或连锁累赘（linkage drag），从而凸显了异域种质在陆地棉育种中同步提升纤维品质与种子形态的应用价值。

该论文发表于《The Plant Genome》，聚焦陆地棉（Gossypium hirsutum L.）次级农艺性状的遗传基础解析。长期以来，陆地棉育种主要围绕纤维产量与纤维品质展开，相关改良已取得显著进展，但株型结构、发育节律和种子形态等性状同样深刻影响群体受光、种植密度适应性、机械化收获效率以及最终生产力。已有研究指出，株高（HT）和分枝分布决定冠层构型与生物量分配；开花天数（DTF）关系到熟期适配与逆境规避；种子表面积则可能影响棉纤维起始与种子外部可利用空间。然而，这些性状在陆地棉改良中相对受到忽视，尤其是平均单粒种子表面积（ASAS）的遗传调控研究更为有限。另一方面，异域种质虽被认为蕴藏丰富遗传多样性，但在精英育种材料中的利用程度仍然不足。因此，开展本研究的必要性在于：在异域×精英遗传背景下系统解析上述次级性状的QTL组成、等位变异来源及其与纤维性状的关系，为陆地棉多性状协同改良提供依据。

研究人员选取5份已完成光周期钝感转育的异域材料T326、T281、T257、T063、T1046，与精英亲本“DES56”构建5个BC₁F₂群体和1个F₃/重组互交（IM，intermated）群体，在2012年和2013年对HT、N1B、DTF和ASAS进行表型鉴定，并结合基因分型测序（GBS，genotyping by sequencing）开展群体特异和联合QTL定位。研究表明，4类性状共检测到17个QTL，其中ASAS的遗传信号最为丰富、广义遗传力（H²）最高，且异域系T326贡献了多数有利等位基因；相比之下，HT、N1B和DTF的遗传变异较小。研究还发现ASAS与DTF在功能富集层面共享CCR4–NOT复合体相关调控通路，而ASAS与纤维强度、上半部平均长度存在显著正相关及QTL共定位。总体结论是，异域种质不仅可用于引入有利纤维品质等位变异，也可在不明显增加负向连锁累赘的情况下促进种子形态改良，具有重要育种应用意义。

主要技术方法概括如下：研究材料来源于6个异域×精英陆地棉分离群体，包括5个BC₁F₂群体和1个F₃/IM群体，2012–2013年在佐治亚大学试验田采用非重复扩增区组设计开展田间表型调查；利用GBS获取全基因组标记，经缺失过滤、次要等位基因频率（MAF，minor allele frequency）筛选、BEAGLE缺失填补后构建物理图谱；采用TASSEL v5.2进行含家系协变量的QTL关联定位，并用Bonferroni校正设定显著性阈值；进一步基于G. hirsutum UTXv3.1参考注释提取候选基因，并采用clusterProfiler开展GO与KEGG富集分析。

在研究结果部分，论文首先在“3.1 Parental analysis and ANOVA”中比较了亲本与群体的表型差异。结果显示，精英亲本DES56在2012和2013年HT与ASAS均低于多数异域亲本，DTF在2013年也低于多数异域系。方差分析表明，在亲本层面仅ASAS达到显著差异；而在分离群体中，环境效应对全部性状和全部群体均极显著，说明年份差异对表型形成具有强烈影响。基因型效应方面，ASAS在所有群体中均达到高度显著，而HT、N1B和DTF只在部分群体中显著，提示后三者在本研究群体中的遗传可分解度较低。

在“3.2 Phenotypic distribution”中，研究人员展示了各群体两年的表型分布。DTF在除群体4和5外的大多数群体中平均值较低；HT在群体3和6中相对较低，且年际差异较大，与方差分析所揭示的环境效应相一致。N1B在不同群体间均值差异不大，而ASAS在两年间总体表现较稳定，其中群体1连续两年均值最高，群体4和6次之。该结果与后续T326来源群体在ASAS上的优势相互印证。

在“3.3 Broad sense heritabilities and correlation analysis”中，ASAS表现出最高的广义遗传力，而其余3个性状由遗传因素解释的表型变异比例较小。性状分布上，ASAS和N1B接近正态连续分布，DTF和HT则略偏离正态。相关分析显示4个性状分为两个关联模块：N1B与HT彼此显著正相关，DTF与ASAS彼此显著正相关，而这两个模块之间则总体呈显著负相关。这提示株型相关性状与发育/种子表面性状可能受到不同但彼此拮抗的调控体系影响。

在“3.4 Correlation between the secondary traits and fiber traits”中，研究人员进一步分析了4个次级性状与2个纤维产量性状及6个HVI（High Volume Instrument，高容量纤维测试）品质性状的关系。结果显示，DTF/ASAS与HT/N1B再次形成两个方向相反的相关模块。例如，DTF和ASAS与马克隆值（MIC）及短纤维含量（SFC）显著负相关，而HT和N1B与这两项指标显著正相关。特别是ASAS与纤维强度（STR）和上半部平均长度（UHM）分别呈显著正相关，为后续QTL共定位分析提供了表型层面的支持。

在“3.5 QTL results and parental contributions”中，研究共检出17个QTL，其中14个为加性效应、3个为显性效应，ASAS检测到10个QTL，数量最多。于“3.5.1 Average surface area per seed”中，全部ASAS QTL均表现为加性效应，T326为9个QTL提供有利等位基因，其中2个位点同时也可由T1046提供有利等位基因，仅qASAS_pop5_2_D09由DES56提供增效等位基因。多个群体特异QTL与联合分析QTL区间高度重叠，尤其是A07、A09与D07上的多个位点在不同分析框架或不同年份中被重复检出，说明这些QTL具有较好的稳定性。研究据此认为T326是改良ASAS最有前景的异域供体。

在“3.5.2 Days to flowering”中，共检测到3个DTF QTL。其中qDTF_pop5_1_D12与qDTF_2_D07为超显性（OD，overdominance）效应，另一个qDTF_pop2_2_D05为加性效应。加性QTL中，T257来源等位基因表现出有利作用；而超显性QTL虽显示较强效应，但其中部分仅在单一分析中出现，仍需进一步验证其稳定性与利用价值。

在“3.5.3 Plant height”中，仅鉴定到2个HT加性QTL：qHT_pop2_1_A11和qHT_pop4_1_D10，分别位于A11和D10染色体。前者中T257等位基因可降低株高，后者的增效等位基因则来自T1046。该结果说明异域材料在局部基因组区段上能够贡献有利的株高调控变异，但整体贡献有限。

在“3.5.4 Number of nodes with one branch”中，仅检测到2个N1B QTL，分别为群体3中的qN1B_pop3_1_D12和联合群体中的qN1B_1_A11。前者为加性效应且有利等位基因来自DES56，后者表现为部分显性（PD，partial dominance）负效应。由此可见，在N1B改良上，本研究所用异域材料并未表现出明显优势。

在“3.6 Putative trait-specific candidate genes in hormonal, RNA, and cytoskeletal pathways”中，研究人员对候选基因及通路进行了功能归纳。N1B相关区间未富集到显著通路。HT相关QTL中，候选基因涉及赤霉素（GA，gibberellin）/DELLA、油菜素内酯（BR，brassinosteroid）、细胞分裂素及离子稳态等调控轴，例如A11区间的木葡聚糖内转糖基酶/水解酶（XTH）基因和D10区间的BR信号及细胞分裂素代谢相关基因，提示株高变异与细胞壁重塑、激素稳态和氧化还原-离子互作有关。DTF和ASAS最突出的共同发现是CCR4–NOT复合体及mRNA降解过程显著富集，且涉及重叠候选基因，提示转录后调控可能同时参与开花转变与种子表面积形成。DTF相关候选基因还包括RIO激酶、FLA蛋白以及茉莉酸、水杨酸和细胞分裂素相关因子，说明花期调控整合了分生组织状态、细胞壁信号与激素通路。ASAS相关候选基因则涉及生长素响应、UDP-糖基转移酶（UGT）、细胞色素P450、细胞骨架组分和COP9信号体（CSN），提示种子表面扩展及纤维起始可能同时受激素代谢、糖基化修饰、蛋白稳态和细胞骨架动态调节。研究还利用泛基因组图谱框架考察了11个唯一QTL区域中的存在/缺失变异（PAV），发现结构变异可能也是这些性状形成的重要遗传来源，但由于研究亲本未被包含在所参照的泛基因组中，这部分证据仍属间接支持。

讨论部分首先强调异域种质的贡献差异。4类性状中，ASAS最具改良潜力，尤其T326在4个独立基因组区域上贡献了多个有利等位基因，且A07上的稳定QTL在不同年份和分析策略下均被重复检出，具有较高的分子标记辅助选择价值。相比之下，T257虽对DTF和HT提供部分有利等位基因，但相关位点仍需进一步验证；N1B则未显示明显的异域利用前景。其次，论文分析了与纤维品质和产量性状的共定位关系。ASAS相关QTL与纤维强度、上半部平均长度QTL存在明显重叠，且相同来源等位基因对相关性状均产生正向影响，这支持多效性或紧密连锁的可能性，并具有商业育种上同步改良多性状的意义。再次，就遗传方式而言，ASAS和HT以加性效应为主，DTF和N1B虽发现少量显性/超显性QTL，但在以自交品种为主的美国棉花生产体系中，显性效应利用空间相对有限。最后，与既往研究比较发现，本研究中的多数ASAS、HT和N1B QTL与前人报道并不重叠，提示可能存在新的遗传变异；而qDTF_pop5_1_D12与先前D12染色体上的花期QTL区间重叠，因而更具后续验证和应用价值。

研究结论可概括为：异域陆地棉种质在纤维品质改良之外，也为种子形态，尤其是ASAS的提升提供了有利等位基因资源；本研究未发现这些次级性状对纤维品质改良存在普遍的负向等位基因权衡或明显连锁累赘，说明将异域种质导入精英基因库不仅有助于缓解等位变异狭窄问题，还有望培育兼具优良纤维性状和更佳种子表型的新型陆地棉品种。